A three-dimensional Magnetic Resonance Spin Tomography in Time-domain protocol for high-resolution multiparametric quantitative magnetic resonance imaging.

Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT) is a multiparametric quantitative MR framework, which allows for simultaneously acquiring quantitative tissue parameters such as T1, T2, and proton density from one single short scan. A typical two-dimensional (2D) MR-STAT acquisition uses a gradient-spoiled, gradient-echo sequence with a slowly varying RF flip-angle train and Cartesian readouts, and the quantitative tissue maps are reconstructed by an iterative, model-based optimization algorithm. In this work, we design a three-dimensional (3D) MR-STAT framework based on previous 2D work, in order to achieve better image signal-to-noise ratio, higher though-plane resolution, and better tissue characterization. Specifically, we design a 7-min, high-resolution 3D MR-STAT sequence, and the corresponding two-step reconstruction algorithm for the large-scale dataset. To reduce the long acquisition time, Cartesian undersampling strategies such as SENSE are adopted in our transient-state quantitative framework. To reduce the computational burden, a data-splitting scheme is designed for decoupling the 3D reconstruction problem into independent 2D reconstructions. The proposed 3D framework is validated by numerical simulations, phantom experiments, and in vivo experiments. High-quality knee quantitative maps with 0.8 × 0.8 × 1.5 mm3 resolution and bilateral lower leg maps with 1.6 mm isotropic resolution can be acquired using the proposed 7-min acquisition sequence and the 3-min-per-slice decoupled reconstruction algorithm. The proposed 3D MR-STAT framework could have wide clinical applications in the future.

Prospective longitudinal cohort study of quantitative muscle magnetic resonance imaging in a healthy control population.

Quantitative muscle magnetic resonance imaging (qMRI) is a valuable methodology for assessing muscular injuries and neuromuscular disorders. Notably, muscle diffusion tensor imaging (DTI) gives insights into muscle microstructural and macrostructural characteristics. However, the long-term reproducibility and robustness of these measurements remain relatively unexplored. The purpose of this prospective longitudinal cohort study was to assess the long-term robustness and range of variation of qMRI parameters, especially DTI metrics, in the lower extremity muscles of healthy controls under real-life conditions. Twelve volunteers (seven females, age 44.1 ± 12.1 years, body mass index 23.3 ± 2.0 kg/m2) underwent five leg muscle MRI sessions every 20 ± 4 weeks over a total period of 1.5 years. A multiecho gradient-echo Dixon-based sequence, a multiecho spin-echo T2-mapping sequence, and a spin-echo echo planar imaging diffusion-weighted sequence were acquired bilaterally with a Philips 3-T Achieva MR System using a 16-channel torso coil. Fifteen leg muscles were segmented in both lower extremities. qMRI parameters, including fat fraction (FF), water T2 relaxation time, and the diffusion metrics fractional anisotropy (FA) and mean diffusivity (MD), were evaluated. Coefficients of variance (wsCV) and intraclass correlation coefficients (ICCs) were calculated to assess the reproducibility of qMRI parameters. The standard error of measurement (SEM) and the minimal detectable change (MDC) were calculated to determine the range of variation. All tests were applied to all muscles and, subsequently, to each muscle separately. wsCV showed good reproducibility (≤ 10%) for all qMRI parameters in all muscles. The ICCs revealed excellent agreement between time points (FF = 0.980, water T2 = 0.941, FA = 0.952, MD = 0.948). Random measurement errors assessed by SEM and the MDC were low (< 12%). In conclusion, in this study, we showed that qMRI parameters in healthy volunteers living normal lives are stable over 18 months, thereby defining a benchmark for the expected range of variation over time.

Quantitative muscle magnetic resonance imaging in limb-girdle muscular dystrophy type R1 (LGMDR1): A prospective longitudinal cohort study.

Limb-girdle muscular dystrophy (LGMD) type R1 (LGMDR1) is the most common subtype of LGMD in Europe. Prospective longitudinal data, including clinical assessments and new biomarkers such as quantitative magnetic resonance imaging (qMRI), are needed to evaluate the natural course of the disease and therapeutic options. We evaluated eight thigh and seven leg muscles of 13 LGMDR1 patients (seven females, mean age 36.7 years, body mass index 23.9 kg/m2) and 13 healthy age- and gender-matched controls in a prospective longitudinal design over 1 year. Clinical assessment included testing for muscle strength with quick motor function measure (QMFM), gait analysis and patient questionnaires (neuromuscular symptom score, activity limitation [ACTIVLIM]). MRI scans were performed on a 3-T MRI scanner, including a Dixon-based sequence, T2 mapping and diffusion tensor imaging. The qMRI values of fat fraction (FF), water T2 relaxation time (T2), fractional anisotropy, mean diffusivity, axial diffusivity and radial diffusivity were analysed. Within the clinical outcome measures, significant deterioration between baseline and follow-up was found for ACTIVLIM (p = 0.029), QMFM (p = 0.012). Analysis of qMRI parameters of the patient group revealed differences between time points for both FF and T2 when analysing all muscles (FF: p < 0.001; T2: p = 0.016). The highest increase of fat replacement was found in muscles with an FF of between 10% and 50% at baseline. T2 in muscles with low-fat replacement increased significantly. No significant differences were found for the diffusion metrics. Significant correlations between qMRI metrics and clinical assessments were found at baseline and follow-up, while only T2 changes in thigh muscles correlated with changes in ACTIVLIM over time (ρ = -0.621, p < 0.05). Clinical assessments can show deterioration of the general condition of LGMDR1 patients. qMRI measures can give additional information about underlying pathophysiology. Further research is needed to establish qMRI outcome measures for clinical trials.

31 P MR Spectroscopy in the Pancreas: Repeatability, Comparison With Liver, and Pilot Pancreatic Cancer Data.

BACKGROUND: Non-invasive evaluation of phosphomonoesters (PMEs) and phosphodiesters (PDEs) by 31-phosphorus MR spectroscopy (31 P MRS) may have potential for early therapy (non-)response assessment in cancer. However, 31 P MRS has not yet been applied to investigate the human pancreas in vivo. PURPOSE: To assess the technical feasibility and repeatability of 31 P MR spectroscopic imaging (MRSI) of the pancreas, compare 31 P metabolite levels between pancreas and liver, and determine the feasibility of 31 P MRSI in pancreatic cancer. STUDY TYPE: Prospective cohort study. POPULATION: 10 healthy subjects (age 34 ± 12 years, four females) and one patient (73-year-old female) with pancreatic ductal adenocarcinoma. FIELD STRENGTH/SEQUENCE: 7-T, 31 P FID-MRSI, 1 H gradient-echo MRI. ASSESSMENT: 31 P FID-MRSI of the abdomen (including the pancreas and liver) was performed with a nominal voxel size of 20 mm (isotropic). For repeatability measurements, healthy subjects were scanned twice on the same day. The patient was only scanned once. Test-retest 31 P MRSI data of pancreas and liver voxels (segmented on 1 H MRI) of healthy subjects were quantified by fitting in the time domain and signal amplitudes were normalized to γ-adenosine triphosphate. In addition, the PME/PDE ratio was calculated. Metabolite levels were averaged over all voxels within the pancreas, right liver lobe and left liver lobe, respectively. STATISTICAL TESTS: Repeatability of test-retest data from healthy pancreas was assessed by paired t-tests, Bland-Altman analyses, and calculation of the intrasubject coefficients of variation (CoVs). Significant differences between healthy pancreas and right and left liver lobes were assessed with a two-way analysis of variance (ANOVA) for repeated measures. A P-value <0.05 was considered statistically significant. RESULTS: The intrasubject CoVs for PME, PDE, and PME/PDE in healthy pancreas were below 20%. Furthermore, PME and PME/PDE were significantly higher in pancreas compared to liver. In the patient with pancreatic cancer, qualitatively, elevated relative PME signals were observed in comparison with healthy pancreas. DATA CONCLUSION: In vivo 31 P MRSI of the human healthy pancreas and in pancreatic cancer may be feasible at 7 T. EVIDENCE LEVEL: 3 TECHNICAL EFFICACY: Stage 2.

Pre-Participation Screenings Frequently Miss Occult Cardiovascular Conditions in Apparently Healthy Male Middle-Aged First-Time Marathon Runners.

INTRODUCTION: The optimal pre-participation screening strategy to identify athletes at risk for exercise-induced cardiovascular events is unknown. We therefore aimed to compare the American College of Sports Medicine (ACSM) and European Society of Cardiology (ESC) pre-participation screening strategies against extensive cardiovascular evaluations in identifying high-risk individuals among 35-50-year-old apparently healthy men. METHODS: We applied ACSM and ESC pre-participation screenings to 25 men participating in a study on first-time marathon running. We compared screening outcomes against medical history, physical examination, electrocardiography, blood tests, echocardiography, cardiopulmonary exercise testing, and magnetic resonance imaging. RESULTS: ACSM screening classified all participants as "medical clearance not necessary." ESC screening classified two participants as "high-risk." Extensive cardiovascular evaluations revealed ≥1 minor abnormality and/or cardiovascular condition in 17 participants, including three subjects with mitral regurgitation and one with a small atrial septal defect. Eleven participants had dyslipidaemia, six had hypertension, and two had premature atherosclerosis. Ultimately, three (12%) subjects had a serious cardiovascular condition warranting sports restrictions: aortic aneurysm, hypertrophic cardiomyopathy (HCM), and myocardial fibrosis post-myocarditis. Of these three participants, only one had been identified as "high-risk" by the ESC screening (for dyslipidaemia, not HCM) and none by the ACSM screening. CONCLUSION: Numerous occult cardiovascular conditions are missed when applying current ACSM/ESC screening strategies to apparently healthy middle-aged men engaging in their first high-intensity endurance sports event.

Intravoxel incoherent motion (IVIM)-derived perfusion fraction mapping for the visual evaluation of MR-guided high intensity focused ultrasound (MR-HIFU) ablation of uterine fibroids.

BACKGROUND: A method for periprocedural contrast agent-free visualization of uterine fibroid perfusion could potentially shorten magnetic resonance-guided high intensity focused ultrasound (MR-HIFU) treatment times and improve outcomes. Our goal was to test feasibility of perfusion fraction mapping by intravoxel incoherent motion (IVIM) modeling using diffusion-weighted MRI as method for visual evaluation of MR-HIFU treatment progression. METHODS: Conventional and T2-corrected IVIM-derived perfusion fraction maps were retrospectively calculated by applying two fitting methods to diffusion-weighted MRI data (b = 0, 50, 100, 200, 400, 600 and 800 s/mm2 at 1.5 T) from forty-four premenopausal women who underwent MR-HIFU ablation treatment of uterine fibroids. Contrast in perfusion fraction maps between areas with low perfusion fraction and surrounding tissue in the target uterine fibroid immediately following MR-HIFU treatment was evaluated. Additionally, the Dice similarity coefficient (DSC) was calculated between delineated areas with low IVIM-derived perfusion fraction and hypoperfusion based on CE-T1w. RESULTS: Average perfusion fraction ranged between 0.068 and 0.083 in areas with low perfusion fraction based on visual assessment, and between 0.256 and 0.335 in surrounding tissues (all p < 0.001). DSCs ranged from 0.714 to 0.734 between areas with low perfusion fraction and the CE-T1w derived non-perfused areas, with excellent intraobserver reliability of the delineated areas (ICC 0.97). CONCLUSION: The MR-HIFU treatment effect in uterine fibroids can be visualized using IVIM perfusion fraction mapping, in moderate concordance with contrast enhanced MRI. IVIM perfusion fraction mapping has therefore the potential to serve as a contrast agent-free imaging method to visualize the MR-HIFU treatment progression in uterine fibroids.

Diffusion tensor imaging and quantitative T2 mapping to monitor muscle recovery following hamstring injury.

MRI examinations are accurate for diagnosing sports-related acute hamstring injuries. However, sensitive imaging methods for assessing recovery of these injuries are lacking. Diffusion tensor imaging (DTI) and quantitative T2 (qT2) mapping have both shown promise for assessing recovery of muscle micro trauma and exercise effects. The purpose of this study was to explore the potential of DTI and qT2 mapping for monitoring the muscle recovery processes after acute hamstring injury. In this prospective study, athletes with an acute hamstring injury underwent a 3-T MRI examination of the injured and contralateral hamstrings including DTI and qT2 measurements at three time points: (1) within 1 week after sustaining the injury, (2) 2 weeks after time point 1, and (3) return to play (RTP). A linear mixed model was used for time-effect analysis and paired t-tests for the detection of differences between injured and uninjured muscles. Forty-one athletes (age 27.8 ± 7 years; two females and 39 males) were included. Mean RTP time was 50 (range 12-169) days. A significant time effect was found for mean diffusivity, radial diffusivity, and the second and third eigenvalues (p ≤ 0.001) in the injured muscles. Fractional anisotropy (p = 0.40), first eigenvalue (p = 0.02), and qT2 (p = 0.61) showed no significant time effect. All DTI indices, except for fractional anisotropy, were significantly elevated compared with control muscles right after the injury (p < 0.001). Values normalized during the recovery period, with no significant differences between control and injured muscles at RTP (p values ranged from 0.08 to 0.51). Mean qT2 relaxation times in injured muscles were not significantly elevated compared with control muscles at any time point (p > 0.04). In conclusion, DTI can be used to monitor recovery after an acute hamstring injury. Future work should explore the potential of DTI indices to predict RTP and recovery times in athletes after an acute strain injury.

In vivo phosphorus magnetic resonance spectroscopic imaging of the whole human liver at 7 T using a phosphorus whole-body transmit coil and 16-channel receive array: Repeatability and effects of principal component analysis-based denoising.

Quantitative three-dimensional (3D) imaging of phosphorus (31 P) metabolites is potentially a promising technique with which to assess the progression of liver disease and monitor therapy response. However, 31 P magnetic resonance spectroscopy has a low sensitivity and commonly used 31 P surface coils do not provide full coverage of the liver. This study aimed to overcome these limitations by using a 31 P whole-body transmit coil in combination with a 16-channel 31 P receive array at 7 T. Using this setup, we determined the repeatability of whole-liver 31 P magnetic resonance spectroscopic imaging (31 P MRSI) in healthy subjects and assessed the effects of principal component analysis (PCA)-based denoising on the repeatability parameters. In addition, spatial variations of 31 P metabolites within the liver were analyzed. 3D 31 P MRSI data of the liver were acquired with a nominal voxel size of 20 mm isotropic in 10 healthy volunteers twice on the same day. Data were reconstructed without denoising, and with PCA-based denoising before or after channel combination. From the test-retest data, repeatability parameters for metabolite level quantification were determined for 12 31 P metabolite signals. On average, 31 P MR spectra from 100 ± 25 voxels in the liver were analyzed. Only voxels with contamination from skeletal muscle or the gall bladder were excluded and no voxels were discarded based on (low) signal-to-noise ratio (SNR). Repeatability for most quantified 31 P metabolite levels in the liver was good to excellent, with an intrasubject variability below 10%. PCA-based denoising increased the SNR ~ 3-fold, but did not improve the repeatability for mean liver 31 P metabolite quantification with the fitting constraints used. Significant spatial heterogeneity of various 31 P metabolite levels within the liver was observed, with marked differences for the phosphomonoester and phosphodiester metabolites between the left and right lobe. In conclusion, using a 31 P whole-body transmit coil in combination with a 16-channel 31 P receive array at 7 T allowed 31 P MRSI acquisitions with full liver coverage and good to excellent repeatability.

Diffusion-weighted MRI with deep learning for visualizing treatment results of MR-guided HIFU ablation of uterine fibroids.

OBJECTIVES: No method is available to determine the non-perfused volume (NPV) repeatedly during magnetic resonance-guided high-intensity focused ultrasound (MR-HIFU) ablations of uterine fibroids, as repeated acquisition of contrast-enhanced T1-weighted (CE-T1w) scans is inhibited by safety concerns. The objective of this study was to develop and test a deep learning-based method for translation of diffusion-weighted imaging (DWI) into synthetic CE-T1w scans, for monitoring MR-HIFU treatment progression. METHODS: The algorithm was retrospectively trained and validated on data from 33 and 20 patients respectively who underwent an MR-HIFU treatment of uterine fibroids between June 2017 and January 2019. Postablation synthetic CE-T1w images were generated by a deep learning network trained on paired DWI and reference CE-T1w scans acquired during the treatment procedure. Quantitative analysis included calculation of the Dice coefficient of NPVs delineated on synthetic and reference CE-T1w scans. Four MR-HIFU radiologists assessed the outcome of MR-HIFU treatments and NPV ratio based on the synthetic and reference CE-T1w scans. RESULTS: Dice coefficient of NPVs was 71% (± 22%). The mean difference in NPV ratio was 1.4% (± 22%) and not statistically significant (p = 0.79). Absolute agreement of the radiologists on technical treatment success on synthetic and reference CE-T1w scans was 83%. NPV ratio estimations on synthetic and reference CE-T1w scans were not significantly different (p = 0.27). CONCLUSIONS: Deep learning-based synthetic CE-T1w scans derived from intraprocedural DWI allow gadolinium-free visualization of the predicted NPV, and can potentially be used for repeated gadolinium-free monitoring of treatment progression during MR-HIFU therapy for uterine fibroids. KEY POINTS: • Synthetic CE-T1w scans can be derived from diffusion-weighted imaging using deep learning. • Synthetic CE-T1w scans may be used for visualization of the NPV without using a contrast agent directly after MR-HIFU ablations of uterine fibroids.

Effect of two eccentric hamstring exercises on muscle architectural characteristics assessed with diffusion tensor MRI.

OBJECTIVES: To evaluate the effect of a Nordic hamstring exercise or Diver hamstring exercise intervention on biceps femoris long head, semitendinosus and semimembranosus muscle's fascicle length and orientation through diffusion tensor imaging (DTI) with magnetic resonance imaging. METHODS: In this three-arm, single-center, randomized controlled trial, injury-free male basketball players were randomly assigned to a Nordic, Diver hamstring exercise intervention or control group. The primary outcome was the DTI-derived fascicle length and orientation of muscles over 12 weeks. RESULTS: Fifty-three participants were included for analysis (mean age 22 ± 7 years). Fascicle length in the semitendinosus over 12 weeks significantly increased in the Nordic-group (mean [M]: 20.8 mm, 95% confidence interval [95% CI]: 7.8 to 33.8) compared with the Control-group (M: 0.9 mm, 95% CI: -7.1 to 8.9), mean between-groups difference: 19.9 mm, 95% CI: 1.9 to 37.9, p = 0.026. Fascicle orientation in the biceps femoris long head over 12 weeks significantly decreased in the Diver-group (mean: -2.6°, 95% CI: -4.1 to -1.0) compared with the Control-group (mean: -0.2°, 95% CI: -1.4 to 1.0), mean between-groups difference: -2.4°, 95% CI: -4.7 to -0.1, p = 0.039. CONCLUSION: The Nordic hamstring exercise intervention did significantly increase the fascicle length of the semitendinosus and the Diver hamstring exercise intervention did significantly change the orientation of fascicles of the biceps femoris long head. As both exercises are complementary to each other, the combination is relevant for preventing hamstring injuries.

Dynamic brain ADC variations over the cardiac cycle and their relation to tissue strain assessed with DENSE at high-field MRI.

PURPOSE: The ADC of brain tissue slightly varies over the cardiac cycle. This variation could reflect physiology, including mixing of the interstitial fluid, relevant for brain waste clearance. However, it is known from cardiac diffusion imaging that tissue deformation by itself affects the magnitude of the MRI signal, leading to artificial ADC variations as well. This study investigates to what extent tissue deformation causes artificial ADC variations in the brain. THEORY AND METHODS: We implemented a high-field MRI sequence with stimulated echo acquisition mode that simultaneously measures brain tissue deformation and ADC. Based on the measured tissue deformation, we simulated the artificial ADC variation by combining established theoretical frameworks and compared the results with the measured ADC variation. We acquired data in 8 healthy volunteers with diffusion weighting b = 300 and b = 1000 s/mm2 . RESULTS: Apparent diffusion coefficient variation was largest in the feet-to-head direction and showed the largest deviation from the mean ADC at peak systole. Artificial ADC variation estimated from tissue deformation was 1.3 ± 0.37·10-5  mm2 /s in the feet-to-head direction for gray matter, and 0.75 ± 0.29·10-5  mm2 /s for white matter. The measured ADC variation in the feet-to-head direction was 5.6·10-5 ± 1.5·10-5  mm2 /s for gray matter and 3.2·10-5 ± 1.0·10-5  mm2 /s for white matter, which was a factor of 3.5 ± 0.82 and 3.4 ± 0.57 larger than the artificial diffusion variations. The measured diffusion variations in the right-to-left/anterior-to-posterior direction were a factor of 1.5 ± 1.0/1.7 ± 1.4 and 2.0 ± 0.91/2.5 ± 0.94 larger than the artificial diffusion variations for gray matter and white matter, respectively. CONCLUSION: Apparent diffusion coefficient variations in the brain likely largely reflect physiology.

Residual quadrupolar couplings observed in 7 Tesla deuterium MR spectra of skeletal muscle.

PURPOSE: Deuterium metabolic imaging could potentially be used to investigate metabolism in skeletal muscle noninvasively. However, skeletal muscle is a tissue with a high degree of spatial organization. In this study, we investigated the effect of incomplete motional averaging on the naturally abundant deuterated water signal in 7 Tesla deuterium spectra of the lower leg muscles and the dependence on the angle between the muscle fibers and the main magnetic field B0 , as determined by DTI. METHODS: Natural abundance deuterium MRSI measurements of the right lower leg muscles were performed at 7 Tesla. Three subjects were scanned in a supine position, with the right leg parallel with the B0 field. One subject was scanned twice; during the second scan, the subject was laying on his right side and the right knee was bent such that the angle between the right lower leg and B0 was approximately 45°. DTI was performed in the same subjects in the same positions at 3 Tesla to determine muscle fiber angles. RESULTS: We observed splittings in the natural abundance deuterated water signal. The size of the splittings varied between different muscles in the lower leg but were mostly similar among subjects for each muscle. The splittings depended on the orientation of the muscle fibers with respect to the main magnetic field B0 . CONCLUSION: Partial molecular alignment in skeletal muscle leads to residual deuteron quadrupolar couplings in deuterated water, the size of which depends on the angle between the muscle fibers and B0 .

Robustness and stability of volume-based tractography in a multicenter setting.

Muscle diffusion tensor imaging (mDTI)-based tractography is a promising tool with which to detect subclinical changes in muscle injuries and to evaluate pathophysiology in neuromuscular diseases. Classic region of interest (ROI)-based tractography is very time-consuming and requires an examiner with extensive experience. (Semi)automatic approaches such as volume-based tractography (VBT) can diminish this problem but its robustness and stability are unknown. The aim of the current study was to assess the performance of VBT in a multicenter setting and to evaluate semiautomatic segmentation approaches in the analysis of VBT-derived data in terms of the comparability of the outcome measures. Five traveling volunteers underwent 3-T mDTI of seven calf muscles of both legs at six different MR sites. Tract properties and diffusion metrics were calculated using VBT. Within-subject coefficients of variance (wsCVs) and intraclass correlation coefficients (ICCs) were calculated to assess the multicenter reproducibility of tract properties such as tract density (TD), mean tract length, volume and tract propagation angle, and diffusion metrics such as fractional anisotropy, mean diffusivity, axial diffusivity (λ1 ) and radial diffusivity in traveling subjects. Furthermore, 50 individual datasets from five different centers (10 datasets per center) were pooled to assess the feasibility of VBT with manual and semiautomatic segmentation. To assess the differences of tract properties and diffusion metrics between segmentation approaches an ANOVA was performed, and ICC and Bland-Altman plots were analyzed. wsCVs and ICCs showed good reproducibility of the tract properties TD and volume, as well as diffusion metrics. ANOVA showed no significant differences between manual and semiautomatic approaches. ICCs were excellent (≥ 0.992) and Bland-Altman analysis did not reveal any systemic bias between the methods. Tract properties and diffusion metrics derived from VBT showed good comparability among centers. Semiautomatic approaches revealed excellent agreement with gold standard of manual segmentation. These findings suggest that pooling data from different centers to construct a reference database for tractography results is feasible using semiautomatic segmentation approaches.

Multi-parametric quantitative magnetic resonance imaging of the upper arm muscles of patients with spinal muscular atrophy.

Quantitative magnetic resonance imaging (qMRI) is frequently used to map the disease state and disease progression in the lower extremity muscles of patients with spinal muscular atrophy (SMA). This is in stark contrast to the almost complete lack of data on the upper extremity muscles, which are essential for carrying out daily activities. The aim of this study was therefore to assess the disease state in the upper arm muscles of patients with SMA in comparison with healthy controls by quantitative assessment of fat fraction, diffusion indices, and water T2 relaxation times, and to relate these measures to muscle force. We evaluated 13 patients with SMA and 15 healthy controls with a 3-T MRI protocol consisting of DIXON, diffusion tensor imaging, and T2 sequences. qMRI measures were compared between groups and related to muscle force measured with quantitative myometry. Fat fraction was significantly increased in all upper arm muscles of the patients with SMA compared with healthy controls and correlated negatively with muscle force. Additionally, fat fraction was heterogeneously distributed within the triceps brachii (TB) and brachialis muscle, but not in the biceps brachii muscle. Diffusion indices and water T2 relaxation times were similar between patients with SMA and healthy controls, but we did find a slightly reduced mean diffusivity (MD), λ1, and λ3 in the TB of patients with SMA. Furthermore, MD was positively correlated with muscle force in the TB of patients with SMA. The variation in fat fraction further substantiates the selective vulnerability of muscles. The reduced diffusion tensor imaging indices, along with the positive correlation of MD with muscle force, point to myofiber atrophy. Our results show the feasibility of qMRI to map the disease state in the upper arm muscles of patients with SMA. Longitudinal data in a larger cohort are needed to further explore qMRI to map disease progression and to capture the possible effects of therapeutic interventions.

Post-exercise intramuscular O2 supply is tightly coupled with a higher proximal-to-distal ATP synthesis rate in human tibialis anterior.

KEY POINTS: The post-exercise recovery of phosphocreatine, a measure of the oxidative capacity of muscles, as assessed by 31 P MR spectroscopy, shows a striking increase from distal to proximal along the human tibialis anterior muscle. To investigate why this muscle exhibits a greater oxidative capacity proximally, we tested whether the spatial variation in phosphocreatine recovery rate is related to oxygen supply, muscle fibre type or type of exercise. We revealed that oxygen supply also increases from distal to proximal along the tibialis anterior, and that it strongly correlated with phosphocreatine recovery. Carnosine level, a surrogate measure for muscle fibre type was not different between proximal and distal, and type of exercise did not affect the gradient in phosphocreatine recovery rate. Taken together, the findings of this study suggest that the post-exercise spatial gradients in oxygen supply and phosphocreatine recovery are driven by a higher intrinsic mitochondrial oxidative capacity proximally. ABSTRACT: Phosphorus magnetic resonance spectroscopy (31 P MRS) of human tibialis anterior (TA) revealed a strong proximo-distal gradient in the post-exercise phosphocreatine (PCr) recovery rate constant (kPCr ), a measure of muscle oxidative capacity. The aim of this study was to investigate whether this kPCr gradient is related to O2 supply, resting phosphorylation potential, muscle fibre type, or type of exercise. Fifteen male volunteers performed continuous isometric ankle dorsiflexion at 30% maximum force until exhaustion. At multiple locations along the TA, we measured the oxidative PCr resynthesis rate (VPCr = kPCr × PCr depletion) by 31 P MRS, the oxyhaemoglobin recovery rate constant (kO2Hb ) by near infrared spectroscopy, and muscle perfusion with MR intravoxel incoherent motion imaging. The kO2Hb , kPCr , VPCr and muscle perfusion depended on measurement location (P < 0.001, P < 0.001, P = 0.032 and P = 0.003, respectively), all being greater proximally. The kO2Hb and muscle perfusion correlated with kPCr (r = 0.956 and r = 0.852, respectively) and VPCr (r = 0.932 and r = 0.985, respectively), the latter reflecting metabolic O2 consumption. Resting phosphorylation potential (PCr/inorganic phosphate) was also higher proximally (P < 0.001). The surrogate for fibre type, carnosine content measured by 1 H MRS, did not differ between distal and proximal TA (P = 0.884). Performing intermittent exercise to avoid exercise ischaemia, still led to larger kPCr proximally than distally (P = 0.013). In conclusion, the spatial kPCr gradient is strongly associated with the spatial variation in O2 supply. It cannot be explained by exercise-induced ischaemia nor by fibre type. Our findings suggest it is driven by a higher proximal intrinsic mitochondrial oxidative capacity, apparently to support contractile performance of the TA.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge.

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

PCA denoising and Wiener deconvolution of 31 P 3D CSI data to enhance effective SNR and improve point spread function.

PURPOSE: This study evaluates the performance of 2 processing methods, that is, principal component analysis-based denoising and Wiener deconvolution, to enhance the quality of phosphorus 3D chemical shift imaging data. METHODS: Principal component analysis-based denoising increases the SNR while maintaining spectral information. Wiener deconvolution reduces the FWHM of the voxel point spread function, which is increased by Hamming filtering or Hamming-weighted acquisition. The proposed methods are evaluated using simulated and in vivo 3D phosphorus chemical shift imaging data by 1) visual inspection of the spatial signal distribution; 2) SNR calculation of the PCr peak; and 3) fitting of metabolite basis functions. RESULTS: With the optimal order of processing steps, we show that the effective SNR of in vivo phosphorus 3D chemical shift imaging data can be increased. In simulations, we show we can preserve phosphorus-containing metabolite peaks that had an SNR < 1 before denoising. Furthermore, using Wiener deconvolution, we were able to reduce the FWHM of the voxel point spread function with only partially reintroducing Gibb-ringing artifacts while maintaining the SNR. After data processing, fitting of the phosphorus-containing metabolite signals improved. CONCLUSION: In this study, we have shown that principal component analysis-based denoising in combination with regularized Wiener deconvolution allows increasing the effective spectral SNR of in vivo phosphorus 3D chemical shift imaging data, with reduction of the FWHM of the voxel point spread function. Processing increased the effective SNR by at least threefold compared to Hamming weighted acquired data and minimized voxel bleeding. With these methods, fitting of metabolite amplitudes became more robust with decreased fitting residuals.

No need to detune transmitters in 32-channel receiver arrays at 7 T.

Ultrahigh field magnetic resonance imaging facilitates high spatiotemporal resolution that benefits from increasing the number of receiver elements. Because high-density receiver arrays have a relatively small element size compared with the transmitter, a side effect is that such setups cause low flux coupling between the transmitter and receiver. Moreover, when transmitters are designed in a multitransmit configuration, their relative size is much smaller than the sample, reducing coupling to the sample and thereby potentially also the coupling to the receivers. Transmitters are traditionally detuned during reception. In this study, we investigate, for a 32-channel receiver head array at 7 T, if transmitter detuning of a quadrature birdcage or of an eight-channel transmit coil can be omitted without substantially sacrificing signal-to-noise ratio (SNR). The transmit elements are operated once with and once without detuning and, in the latter, the received signals are either merged with the array or excluded for image reconstruction. For each of the three measurements, SNR and 1/g-factor maps are investigated. The tuning of the quadrature and eight-channel transmit coils during signal reception introduced a 10.1% and 6.5% penalty in SNR, respectively, relative to the SNR received with detuned transmitters. When also incorporating the signal of the transmit coils, the SNR was regained to 98.5% or 101.4% for the quadrature and eight-channel coil, respectively, relative to the detuned transmitters, while the 1/g-factor maps improved slightly. For the 32-channel receive coil used the SNR penalty can become negligible when omitting detuning of the transmit coils. This not only simplifies transmit coil designs, potentially increasing their efficiency, but also enables the transmitters to be used as receivers in parallel to the receiver array, thus increasing parallel imaging performance.

Evaluation of interrater reliability of different muscle segmentation techniques in diffusion tensor imaging.

INTRODUCTION: Muscle diffusion tensor imaging (mDTI) is a quantitative MRI technique that can provide information about muscular microstructure and integrity. Ultrasound and DTI studies have shown intramuscular differences, and therefore separation of different muscles for analysis is essential. The commonly used methods to assess DTI metrics in muscles are manual segmentation and tract-based analysis. Recently methods such as volume-based tractography have been applied to optimize muscle architecture estimation, but can also be used to assess DTI metrics. PURPOSE: To evaluate diffusion metrics obtained using three different methods-volume-based tractography, manual segmentation-based analysis and tract-based analysis-with respect to their interrater reliability and their ability to detect intramuscular variance. MATERIALS AND METHODS: 30 volunteers underwent an MRI examination in a 3 T scanner using a 16-channel Torso XL coil. Diffusion-weighted images were acquired to obtain DTI metrics. These metrics were evaluated in six thigh muscles using volume-based tractography, manual segmentation and standard tractography. All three methods were performed by two independent raters to assess interrater reliability by ICC analysis and Bland-Altman plots. Ability to assess intramuscular variance was compared using an ANOVA with muscle as a between-subjects factor. RESULTS: Interrater reliability for all methods was found to be excellent. The highest interrater reliability was found for volume-based tractography (ICC ≥ 0.967). Significant differences for the factor muscle in all examined diffusion parameters were shown in muscles using all methods (main effect p < 0.001). CONCLUSIONS: Diffusion data can be assessed by volume tractography, standard tractography and manual segmentation with high interrater reliability. Each method produces different results for the investigated DTI parameters. Volume-based tractography was superior to conventional manual segmentation and tractography regarding interrater reliability and detection of intramuscular variance, while tract-based analysis showed the lowest coefficients of variation.

Quantification of disease progression in spinal muscular atrophy with muscle MRI-a pilot study.

OBJECTIVES: Quantitative MRI (qMRI) of muscles is a promising tool to measure disease progression or to assess therapeutic effects in neuromuscular diseases. Longitudinal imaging studies are needed to show sensitivity of qMRI in detecting disease progression in spinal muscular atrophy (SMA). In this pilot study we therefore studied one-year changes in quantitative MR parameters in relation to clinical scores. METHODS: We repeated quantitative 3 T MR analysis of thigh muscles and clinical testing one year after baseline in 10 treatment-naïve patients with SMA, 5 with Type 2 (21.6 ± 7.0 years) and 5 with Type 3 (33.4 ± 11.9 years). MR protocol consisted of Dixon, T2 mapping and diffusion tensor imaging (DTI). The temporal relation of parameters was examined with a mixed model. RESULTS: We detected a significant increase in fat fraction (baseline, 38.2% SE 0.6; follow-up, 39.5% SE 0.6; +1.3%, p = 0.001) in all muscles. Muscles with moderate to high fat infiltration at baseline show a larger increase over time (+1.6%, p < 0.001). We did not find any changes in DTI parameters except for low fat-infiltration muscles (m. adductor longus and m. biceps femoris (short head)). The T2 of muscles decreased from 28.2 ms to 28.0 ms (p = 0.07). Muscle strength and motor function scores were not significantly different between follow-up and baseline. CONCLUSION: Longitudinal imaging data show slow disease progression in skeletal muscle of the thigh of (young-) adult patients with SMA despite stable strength and motor function scores. Quantitative muscle imaging demonstrates potential as a biomarker for disease activity and monitoring of therapy response.