Muscle fiber strain rates in the lower leg during ankle dorsi-/plantarflexion exercise.

Static quantitative magnetic resonance imaging (MRI) provides readouts of structural changes in diseased muscle, but current approaches lack the ability to fully explain the loss of contractile function. Muscle contractile function can be assessed using various techniques including phase-contrast MRI (PC-MRI), where strain rates are quantified. However, current two-dimensional implementations are limited in capturing the complex motion of contracting muscle in the context of its three-dimensional (3D) fiber architecture. The MR acquisitions (chemical shift-encoded water-fat separation scan, spin echo-echoplanar imaging with diffusion weighting, and two time-resolved 3D PC-MRI) wereperformed at 3 T. PC-MRI acquisitions and performed with and without load at 7.5% of the maximum voluntary dorsiflexion contraction force. Acquisitions (3 T, chemical shift-encoded water-fat separation scan, spin echo-echo planar imaging with diffusion weighting, and two time-resolved 3D PC-MRI) were performed with and without load at 7.5% of the maximum voluntary dorsiflexion contraction force. Strain rates and diffusion tensors were calculated and combined to obtain strain rates along and perpendicular to the muscle fibers in seven lower leg muscles during the dynamic dorsi-/plantarflexion movement cycle. To evaluate strain rates along the proximodistal muscle axis, muscles were divided into five equal segments. t-tests were used to test if cyclic strain rate patterns (amplitude > 0) were present along and perpendicular to the muscle fibers. The effects of proximal-distal location and load were evaluated using repeated measures ANOVAs. Cyclic temporal strain rate patterns along and perpendicular to the fiber were found in all muscles involved in dorsi-/plantarflexion movement (p < 0.0017). Strain rates along and perpendicular to the fiber were heterogeneously distributed over the length of most muscles (p < 0.003). Additional loading reduced strain rates of the extensor digitorum longus and gastrocnemius lateralis muscle (p < 0.001). In conclusion, the lower leg muscles involved in cyclic dorsi-/plantarflexion exercise showed cyclic fiber strain rate patterns with amplitudes that varied between muscles and between the proximodistal segments within the majority of muscles.

Image registration and mutual thresholding enable low interimage variability across dynamic MRI measurements of supraclavicular brown adipose tissue during mild cold exposure.

PURPOSE: Activated brown adipose tissue (BAT) enhances lipid catabolism and improves cardiometabolic health. Quantitative MRI of the fat fraction (FF) of supraclavicular BAT (scBAT) is a promising noninvasive measure to assess BAT activity but suffers from high scan variability. We aimed to test the effects of coregistration and mutual thresholding on the scan variability in a fast (1 min) time-resolution MRI protocol for assessing scBAT FF changes during cold exposure. METHODS: Ten volunteers (age 24.8 ± 3.0 years; body mass index 21.2 ± 2.1 kg/m2 ) were scanned during thermoneutrality (32°C; 10 min) and mild cold exposure (18°C; 60 min) using a 12-point gradient-echo sequence (70 consecutive scans with breath-holds, 1.03 min per dynamic). Dynamics were coregistered to the first thermoneutral scan, which enabled drawing of single regions of interest in the scBAT depot. Voxel-wise FF changes were calculated at each time point and averaged across regions of interest. We applied mutual FF thresholding, in which voxels were included if their FF was greater than 30% FF in the reference scan and the registered dynamic. The efficacy of the coregistration was determined by using a moving average and comparing the mean squared error of residuals between registered and nonregistered data. Registered scBAT ΔFF was compared with single-scan thresholding using the moving average method. RESULTS: Registered scBAT ΔFF had lower mean square error values than nonregistered data (0.07 ± 0.05% vs. 0.16 ± 0.14%; p < 0.05), and mutual thresholding reduced the scBAT ΔFF variability by 30%. CONCLUSION: We demonstrate that coregistration and mutual thresholding improve stability of the data 2-fold, enabling assessment of small changes in FF following cold exposure.

Diagnosing myasthenia gravis using orthoptic measurements: assessing extraocular muscle fatiguability.

INTRODUCTION: Diagnosing ocular myasthenia gravis (MG) can be challenging because serum antibodies are often not detected. We aimed to explore whether determining extraocular muscle (EOM) weakness using orthoptic measures, including an adapted Hess chart examination, can aid in diagnosing MG. METHODS: We conducted a prospective study among patients with acetylcholine receptor antibody positive MG (20 recently diagnosed, 19 chronic) and 14 seronegative MG patients. We compared orthoptic measures to 19 healthy and 18 disease controls with Graves orbitopathy, chronic progressive external ophthalmoplegia or oculopharyngeal muscular dystrophy. Maximal eye duction angles were measured using a synoptophore. Gaze deviations between eyes were measured using standard Hess chart examination with addition of 1 min persistent gaze to assess MG-associated fatiguability. Receiver operating characteristics curve analysis was performed. RESULTS: For duction angles, the area under the curve (AUC) was 0.73 comparing MG to healthy, and 0.69 comparing to patient controls. For the outer field of the Hess chart, the AUC was 0.89 comparing to healthy and 0.54 to patient controls. For drift, the AUC was 0.93 comparing to healthy and 0.93 to patient controls. The sensitivity and specificity of the presence of drift was 81% and 100%. DISCUSSION: Orthoptic measurements can be used to diagnose MG by quantifying EOM weakness and fatiguability. Drift during persistent gaze on a Hess chart is specific for MG and could be used for diagnostic purposes. The Hess chart examination is widely available, inexpensive and fast. Moreover, orthoptic measurements may be a clinically relevant outcome measure for clinical trials.

Effects of muscle damage on 31 P-MRS indices of energetic status and sarcolemma integrity in young mdx mice.

The article of Lopez et al describes the use of a multi-parametric MR approach to study muscle T2 relaxation times and 31 P-MRS indices of energetics and sarcolemma integrity in a mouse model of DMD, the mdx mouse. Muscular dystrophies have a multi-factorial disease cascade, and there are several MR methods used to assess these. Aspects that reflect disease progression are outlined on the left, while features that are more related to disease activity are outlined on the right.

Improved detection limits of J-coupled neurometabolites in the human brain at 7 T with a J-refocused sLASER sequence.

In a standard spin echo, the time evolution due to homonuclear couplings is not reversed, leading to echo time (TE)-dependent modulation of the signal amplitude and signal loss in the case of overlapping multiplet resonances. This has an adverse effect on quantification of several important metabolites such as glutamate and glutamine. Here, we propose a J-refocused variant of the sLASER sequence (J-sLASER) to improve quantification of J-coupled metabolites at ultrahigh field (UHF). The use of the sLASER sequence is particularly advantageous at UHF as it minimizes chemical shift displacement error and results in relatively homogenous refocusing. We simulated the MRS signal from brain metabolites over a broad range of TE values with sLASER and J-sLASER, and showed that the signal of J-coupled metabolites was increased with J-sLASER with TE values up to ~80 ms. We further simulated "brain-like" spectra with both sequences at the shortest TE available on our scanner. We showed that, despite the slightly longer TE, the J-sLASER sequence results in significantly lower Cramer-Rao lower bounds (CRLBs) for J-coupled metabolites compared with those obtained with sLASER. Following phantom validation, we acquired spectra from two brain regions in 10 healthy volunteers (age 38 ± 15 years) using both sequences. We showed that using J-sLASER results in a decrease of CRLBs for J-coupled metabolites. In particular, we measured a robust ~38% decrease in the mean CRLB (glutamine) in parietal white matter and posterior cingulate cortex (PCC). We further showed, in 10 additional healthy volunteers (age 34 ± 15 years), that metabolite quantification following two separate acquisitions with J-sLASER in the PCC was repeatable. The improvement in quantification of glutamine may in turn improve the independent quantification of glutamate, the main excitatory neurotransmitter in the brain, and will simultaneously help to track possible modulations of glutamine, which is a key player in the glutamatergic cycle in astrocytes.

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.

Microvascular response to exercise varies along the length of the tibialis anterior muscle.

Microvascular function is an important component in the physiology of muscle. One of the major parameters, blood perfusion, can be measured noninvasively and quantitatively by arterial spin labeling (ASL) MRI. Most studies using ASL in muscle have only reported data from a single slice, thereby assuming that muscle perfusion is homogeneous within muscle, whereas recent literature has reported proximodistal differences in oxidative capacity and perfusion. Here, we acquired pulsed ASL data in 12 healthy volunteers after dorsiflexion exercise in two slices separated distally by 7 cm. We combined this with a Look-Locker scheme to acquire images at multiple postlabeling delays (PLDs) and with a multiecho readout to measure T2 *. This enabled the simultaneous evaluation of quantitative muscle blood flow (MBF), arterial transit time (ATT), and T2 * relaxation time in the tibialis anterior muscle during recovery. Using repeated measures analyses of variance we tested the effect of time, slice location, and their interaction on MBF, ATT, and T2 *. Our results showed a significant difference as a function of time postexercise for all three parameters (MBF: F = 34.0, p < .0001; T2 *: F = 73.7, p < .0001; ATT: F = 13.6, p < .001) and no average differences between slices over the total time postexercise were observed. The interaction effect between time postexercise and slice location was significant for MBF and T2 * (F = 5.5, p = 0.02, F = 6.1, p = 0.02, respectively), but not for ATT (F = 2.2, p = .16). The proximal slice showed a higher MBF and a lower ATT than the distal slice during the first 2 min of recovery, and T2 * showed a delayed response in the distal slice. These results imply a higher perfusion and faster microvascular response to exercise in the proximal slice, in line with previous literature. Moreover, the differences in ATT indicate that it is difficult to correctly determine perfusion based on a single PLD as is commonly performed in the muscle literature.

Baseline fat fraction is a strong predictor of disease progression in Becker muscular dystrophy.

In Becker muscular dystrophy (BMD), muscle weakness progresses relatively slowly, with a highly variable rate among patients. This complicates clinical trials, as clinically relevant changes are difficult to capture within the typical duration of a trial. Therefore, predictors for disease progression are needed. We assessed if temporal increase of fat fraction (FF) in BMD follows a sigmoidal trajectory and whether fat fraction at baseline (FFbase) could therefore predict FF increase after 2 years (ΔFF). Thereafter, for two different MR-based parameters, we tested the additional predictive value to FFbase. We used 3-T Dixon data from the upper and lower leg, and multiecho spin-echo MRI and 7-T 31 P MRS datasets from the lower leg, acquired in 24 BMD patients (age: 41.4 [SD 12.8] years). We assessed the pattern of increase in FF using mixed-effects modelling. Subsequently, we tested if indicators of muscle damage like standard deviation in water T2 (stdT2 ) and the phosphodiester (PDE) over ATP ratio at baseline had additional value to FFbase for predicting ∆FF. The association between FFbase and ΔFF was described by the derivative of a sigmoid function and resulted in a peak ΔFF around 0.45 FFbase (fourth-order polynomial term: t = 3.7, p < .001). StdT2 and PDE/ATP were not significantly associated with ∆FF if FFbase was included in the model. The relationship between FFbase and ∆FF suggests a sigmoidal trajectory of the increase in FF over time in BMD, similar to that described for Duchenne muscular dystrophy. Our results can be used to identify muscles (or patients) that are in the fast progressing stage of the disease, thereby facilitating the conduct of clinical trials.

Characterization of patients with Becker muscular dystrophy by histology, magnetic resonance imaging, function, and strength assessments.

INTRODUCTION/AIMS: Becker muscular dystrophy (BMD) is characterized by variable disease severity and progression, prompting the identification of biomarkers for clinical trials. We used data from an ongoing phase II study to provide a comprehensive characterization of a cohort of patients with BMD, and to assess correlations between histological and magnetic resonance imaging (MRI) markers with muscle function and strength. METHODS: Eligible patients were ambulatory males with BMD, aged 18 to 65 years (200 to 450 meters on 6-minute walk test). The following data were obtained: function test results, strength, fat-fraction quantification using chemical shift-encoded MRI (whole thigh and quadriceps), and fibrosis and muscle fiber area (MFA) of the brachial biceps. RESULTS: Of 70 patients screened, 51 entered the study. There was substantial heterogeneity between patients in muscle morphology (histology and MRI), with high fat replacement. Total fibrosis correlated significantly and mostly moderately with all functional endpoints, including both upper arm strength assessments (left and right elbow flexion rho -.574 and -.588, respectively [both P < .0001]), as did MRI fat fraction (whole thigh and quadriceps), for example, with four-stair-climb velocity -.554 and -.550, respectively (both P < .0001). Total fibrosis correlated significantly and moderately with both MRI fat fraction assessments (.500 [P = .0003] and .423 [.0024], respectively). DISCUSSION: In this BMD cohort, micro- and macroscopic morphological muscle parameters correlated moderately with each other and with functional parameters, potentially supporting the use of MRI fat fraction and histology as surrogate outcome measures in patients with BMD, although additional research is required to validate this.

Association of shivering threshold time with body composition and brown adipose tissue in young adults.

PURPOSE: Brown adipose tissue (BAT) increases metabolic heat production in response to cold exposure. Body size and composition are involved in the human cold response, yet the influence of BAT herein have not fully been explored. Here, we aimed to study the association of the cold-induced shivering threshold time with body composition, BAT, the perception of shivering and skin temperature in young adults. METHODS: 110 young healthy adults (81 females; age = 21.7 ± 2.1 years, BMI = 24.2 ± 4.3 kg/m2) underwent 2 h of individualized cooling, followed by the quantification of BAT using a18F-fluorodeoxyglucose ([18F]FDG) positron emission tomography-computed tomography (PET-CT) scan. Body mass index (BMI), lean mass, fat mass and body surface area (BSA) were also measured. Shivering threshold time was defined as the time until shivering occurred using an individualized cooling protocol. RESULTS: The shivering threshold time was on average 116.1 min for males and 125.8 min for females, and was positively associated to BMI (ß = 3.106; R2 = 0.141; p = 0.001), lean mass (ß = 2.295; R2 = 0.128; p = 0.001) and fat mass (ß = 1.492; R2 = 0.121; p = 0.001) in females, but not in males (all p ≥ 0.409). The shivering threshold time was positively associated with BSA in males (p = 0.047) and females (p = 0.001), but it was not associated with BAT volume or [18F]FDG uptake nor with the perception of shivering and skin temperature perception in both sexes. CONCLUSION: The shivering threshold time is positively associated with whole-body adiposity and lean mass in females, but not in males. The shivering threshold time was positively associated with BSA, but no association was observed with BAT nor with the perception of shivering or skin temperature. Future research should consider the influence of body composition when applying cooling protocols among individuals with different phenotypical features.

The Infrared Thermography Toolbox: An Open-access Semi-automated Segmentation Tool for Extracting Skin Temperatures in the Thoracic Region including Supraclavicular Brown Adipose Tissue.

Infrared thermography (IRT) is widely used to assess skin temperature in response to physiological changes. Yet, it remains challenging to standardize skin temperature measurements over repeated datasets. We developed an open-access semi-automated segmentation tool (the IRT-toolbox) for measuring skin temperatures in the thoracic area to estimate supraclavicular brown adipose tissue (scBAT) activity, and compared it to manual segmentations. The IRT-toolbox, designed in Python, consisted of image pre-alignment and non-rigid image registration. The toolbox was tested using datasets of 10 individuals (BMI = 22.1 ± 2.1 kg/m2, age = 22.0 ± 3.7 years) who underwent two cooling procedures, yielding four images per individual. Regions of interest (ROIs) were delineated by two raters in the scBAT and deltoid areas on baseline images. The toolbox enabled direct transfer of baseline ROIs to the registered follow-up images. For comparison, both raters also manually drew ROIs in all follow-up images. Spatial ROI overlap between methods and raters was determined using the Dice coefficient. Mean bias and 95% limits of agreement in mean skin temperature between methods and raters were assessed using Bland-Altman analyses. ROI delineation time was four times faster with the IRT-toolbox (01:04 min) than with manual delineations (04:12 min). In both anatomical areas, there was a large variability in ROI placement between methods. Yet, relatively small skin temperature differences were found between methods (scBAT: 0.10 °C, 95%LoA[-0.13 to 0.33 °C] and deltoid: 0.05 °C, 95%LoA[-0.46 to 0.55 °C]). The variability in skin temperature between raters was comparable between methods. The IRT-toolbox enables faster ROI delineations, while maintaining inter-user reliability compared to manual delineations. (Trial registration number (ClinicalTrials.gov): NCT04406922, [May 29, 2020]).

Compartmental diffusion and microstructural properties of human brain gray and white matter studied with double diffusion encoding magnetic resonance spectroscopy of metabolites and water.

Double diffusion encoding (DDE) of the water signal offers a unique ability to separate the effect of microscopic anisotropic diffusion in structural units of tissue from the overall macroscopic orientational distribution of cells. However, the specificity in detected microscopic anisotropy is limited as the signal is averaged over different cell types and across tissue compartments. Performing side-by-side water and metabolite DDE spectroscopic (DDES) experiments provides complementary measures from which intracellular and extracellular microscopic fractional anisotropies (µFA) and diffusivities can be estimated. Metabolites are largely confined to the intracellular space and therefore provide a benchmark for intracellular µFA and diffusivities of specific cell types. By contrast, water DDES measurements allow examination of the separate contributions to water µFA and diffusivity from the intra- and extracellular spaces, by using a wide range of b values to gradually eliminate the extracellular contribution. Here, we aimed to estimate tissue and compartment specific human brain microstructure by combining water and metabolites DDES experiments. We performed our DDES measurements in two brain regions that contain widely different amounts of white matter (WM) and gray matter (GM): parietal white matter (PWM) and occipital gray matter (OGM) in a total of 20 healthy volunteers at 7 Tesla. Metabolite DDES measurements were performed at b = 7199 s/mm2, while water DDES measurements were performed with a range of b values from 918 to 7199 s/mm2. The experimental framework we employed here resulted in a set of insights pertaining to the morphology of the intracellular and extracellular spaces in both gray and white matter. Results of the metabolite DDES experiments in both PWM and OGM suggest a highly anisotropic intracellular space within neurons and glia, with the possible exception of gray matter glia. The water µFA obtained from the DDES results at high b values in both regions converged with that of the metabolite DDES, suggesting that the signal from the extracellular space is indeed effectively suppressed at the highest b value. The µFA measured in the OGM significantly decreased at lower b values, suggesting a considerably lower anisotropy of the extracellular space in GM compared to WM. In PWM, the water µFA remained high even at the lowest b value, indicating a high degree of organization in the interstitial space in WM. Tortuosity values in the cytoplasm for water and tNAA, obtained with correlation analysis of microscopic parallel diffusivity with respect to GM/WM tissue fraction in the volume of interest, are remarkably similar for both molecules, while exhibiting a clear difference between gray and white matter, suggesting a more crowded cytoplasm and more complex cytomorphology of neuronal cell bodies and dendrites in GM than those found in long-range axons in WM.

A split-label design for simultaneous measurements of perfusion in distant slices by pulsed arterial spin labeling.

PURPOSE: Multislice arterial spin labeling (ASL) MRI acquisitions are currently challenging in skeletal muscle because of long transit times, translating into low-perfusion SNR in distal slices when large spatial coverage is required. However, fiber type and oxidative capacity vary along the length of healthy muscles, calling for multislice acquisitions in clinical studies. We propose a new variant of flow alternating inversion recovery (FAIR) that generates sufficient ASL signal to monitor exercise-induced perfusion changes in muscle in two distant slices. METHODS: Label around and between two 7-cm distant slices was created by applying the presaturation/postsaturation and selective inversion modules selectively to each slice (split-label multislice FAIR). Images were acquired using simultaneous multislice EPI. We validated our approach in the brain to take advantage of the high resting-state perfusion, and applied it in the lower leg muscle during and after exercise, interleaved with a single-slice FAIR as a reference. RESULTS: We show that standard multislice FAIR leads to an underestimation of perfusion, while the proposed split-label multislice approach shows good agreement with separate single-slice FAIR acquisitions in brain, as well as in muscle following exercise. CONCLUSION: Split-label FAIR allows measuring muscle perfusion in two distant slices simultaneously without losing sensitivity in the distal slice.

The feasibility of quantitative MRI of extra-ocular muscles in myasthenia gravis and Graves' orbitopathy.

Although quantitative MRI can be instrumental in the diagnosis and assessment of disease progression in orbital diseases involving the extra-ocular muscles (EOM), acquisition can be challenging as EOM are small and prone to eye-motion artefacts. We explored the feasibility of assessing fat fractions (FF), muscle volumes and water T2 (T2water ) of EOM in healthy controls (HC), myasthenia gravis (MG) and Graves' orbitopathy (GO) patients. FF, EOM volumes and T2water values were determined in 12 HC (aged 22-65 years), 11 MG (aged 28-71 years) and six GO (aged 28-64 years) patients at 7 T using Dixon and multi-echo spin-echo sequences. The EOM were semi-automatically 3D-segmented by two independent observers. MANOVA and t-tests were used to assess differences in FF, T2water and volume of EOM between groups (P < .05). Bland-Altman limits of agreement (LoA) were used to assess the reproducibility of segmentations and Dixon scans. The scans were well tolerated by all subjects. The bias in FF between the repeated Dixon scans was -0.7% (LoA: ±2.1%) for the different observers; the bias in FF was -0.3% (LoA: ±2.8%) and 0.03 cm3 (LoA: ± 0.36 cm3 ) for volume. Mean FF of EOM in MG (14.1% ± 1.6%) was higher than in HC (10.4% ± 2.5%). Mean muscle volume was higher in both GO (1.2 ± 0.4 cm3 ) and MG (0.8 ± 0.2 cm3 ) compared with HC (0.6 ± 0.2 cm3 ). The average T2water for all EOM was 24.6 ± 4.0 ms for HC, 24.0 ± 4.7 ms for MG patients and 27.4 ± 4.2 ms for the GO patient. Quantitative MRI at 7 T is feasible for measuring FF and muscle volumes of EOM in HC, MG and GO patients. The measured T2water was on average comparable with skeletal muscle, although with higher variation between subjects. The increased FF in the EOM in MG patients suggests that EOM involvement in MG is accompanied by fat replacement. The unexpected EOM volume increase in MG may provide novel insights into underlying pathophysiological processes.

Muscle architecture is associated with muscle fat replacement in Duchenne and Becker muscular dystrophies.

INTRODUCTION/AIMS: Duchenne and Becker muscular dystrophies (DMD and BMD, respectively) are characterized by fat replacement of different skeletal muscles in a specific temporal order. Given the structural role of dystrophin in skeletal muscle mechanics, muscle architecture could be important in the progressive pathophysiology of muscle degeneration. Therefore, the aim of this study was to assess the role of muscle architecture in the progression of fat replacement in DMD and BMD. METHODS: We assessed the association between literature-based leg muscle architectural characteristics and muscle fat fraction from 22 DMD and 24 BMD patients. Dixon-based magnetic resonance imaging estimates of fat fractions at baseline and 12 (only DMD) and 24 months were related to fiber length and physiological cross-sectional area (PCSA) using age-controlled linear mixed modeling. RESULTS: DMD and BMD muscles with long fibers and BMD muscles with large PCSAs were associated with increased fat fraction. The effect of fiber length was stronger in muscles with larger PCSA. DISCUSSION: Muscle architecture may explain the pathophysiology of muscle degeneration in dystrophinopathies, in which proximal muscles with a larger mass (fiber length × PCSA) are more susceptible, confirming the clinical observation of a temporal proximal-to-distal progression. These results give more insight into the mechanical role in the pathophysiology of muscular dystrophies. Ultimately, this new information can be used to help support the selection of current and the development of future therapies.

Recommendations for the analysis of gene expression data to identify intrinsic differences between similar tissues.

Identifying genes involved in functional differences between similar tissues from expression profiles is challenging, because the expected differences in expression levels are small. To exemplify this challenge, we studied the expression profiles of two skeletal muscles, deltoid and biceps, in healthy individuals. We provide a series of guides and recommendations for the analysis of this type of studies. These include how to account for batch effects and inter-individual differences to optimize the detection of gene signatures associated with tissue function. We provide guidance on the selection of optimal settings for constructing gene co-expression networks through parameter sweeps of settings and calculation of the overlap with an established knowledge network. Our main recommendation is to use a combination of the data-driven approaches, such as differential gene expression analysis and gene co-expression network analysis, and hypothesis-driven approaches, such as gene set connectivity analysis. Accordingly, we detected differences in metabolic gene expression between deltoid and biceps that were supported by both data- and hypothesis-driven approaches. Finally, we provide a bioinformatic framework that support the biological interpretation of expression profiles from related tissues from this combination of approaches, which is available at github.com/tabbassidaloii/AnalysisFrameworkSimilarTissues.

Water-fat separation in spiral magnetic resonance fingerprinting for high temporal resolution tissue relaxation time quantification in muscle.

PURPOSE: To minimize the known biases introduced by fat in rapid T1 and T2 quantification in muscle using a single-run magnetic resonance fingerprinting (MRF) water-fat separation sequence. METHODS: The single-run MRF acquisition uses an alternating in-phase/out-of-phase TE pattern to achieve water-fat separation based on a 2-point DIXON method. Conjugate phase reconstruction and fat deblurring were applied to correct for B0 inhomogeneities and chemical shift blurring. Water and fat signals were matched to the on-resonance MRF dictionary. The method was first tested in a multicompartment phantom. To test whether the approach is capable of measuring small in vivo dynamic changes in relaxation times, experiments were run in 9 healthy volunteers; parameter values were compared with and without water-fat separation during muscle recovery after plantar flexion exercise. RESULTS: Phantom results show the robustness of the water-fat resolving MRF approach to undersampling. Parameter maps in volunteers show a significant (P < .01) increase in T1 (105 ± 94 ms) and decrease in T2 (14 ± 6 ms) when using water-fat-separated MRF, suggesting improved parameter quantification by reducing the well-known biases introduced by fat. Exercise results showed smooth T1 and T2 recovery curves. CONCLUSION: Water-fat separation using conjugate phase reconstruction is possible within a single-run MRF scan. This technique can be used to rapidly map relaxation times in studies requiring dynamic scanning, in which the presence of fat is problematic.

T2 relaxation-time mapping in healthy and diseased skeletal muscle using extended phase graph algorithms.

PURPOSE: Multi-echo spin-echo (MSE) transverse relaxometry mapping using multi-component models is used to study disease activity in neuromuscular disease by assessing the T2 of the myocytic component (T2water ). Current extended phase graph algorithms are not optimized for fat fractions above 50% and the effects of inaccuracies in the T2fat calibration remain unexplored. Hence, we aimed to improve the performance of extended phase graph fitting methods over a large range of fat fractions, by including the slice-selection flip angle profile, a through-plane chemical-shift displacement correction, and optimized calibration of T2fat . METHODS: Simulation experiments were used to study the influence of the slice flip-angle profile with chemical-shift and T2fat estimations. Next, in vivo data from four neuromuscular disease cohorts were studied for different T2fat calibration methods and T2water estimations. RESULTS: Excluding slice flip-angle profiles or chemical-shift displacement resulted in a bias in T2water up to 10 ms. Furthermore, a wrongly calibrated T2fat caused a bias of up to 4 ms in T2water . For the in vivo data, one-component calibration led to a lower T2fat compared with a two-component method, and T2water decreased with increasing fat fractions. CONCLUSION: In vivo data showed a decline in T2water for increasing fat fractions, which has important implications for clinical studies, especially in multicenter settings. We recommend using an extended phase graph-based model for fitting T2water from MSE sequences with two-component T2fat calibration. Moreover, we recommend including the slice flip-angle profile in the model with correction for through-plane chemical-shift displacements.

Multi-parametric MR in Becker muscular dystrophy patients.

Quantitative MRI and MRS of muscle are increasingly being used to measure individual pathophysiological processes in Becker muscular dystrophy (BMD). In particular, muscle fat fraction was shown to be highly associated with functional tests in BMD. However, the muscle strength per unit of contractile cross-sectional area is lower in patients with BMD compared with healthy controls. This suggests that the quality of the non-fat-replaced (NFR) muscle tissue is lower than in healthy controls. Consequently, a measure that reflects changes in muscle tissue itself is needed. Here, we explore the potential of water T2 relaxation times, diffusion parameters and phosphorus metabolic indices as early disease markers in patients with BMD. For this purpose, we examined these measures in fat-replaced (FR) and NFR lower leg muscles in patients with BMD and compared these values with those in healthy controls. Quantitative proton MRI (three-point Dixon, multi-spin-echo and diffusion-weighted spin-echo echo planar imaging) and 2D chemical shift imaging 31 P MRS data were acquired in 24 patients with BMD (age 18.8-66.2 years) and 13 healthy controls (age 21.3-63.6 years). Muscle fat fractions, phosphorus metabolic indices, and averages and standard deviations (SDs) of the water T2 relaxation times and diffusion tensor imaging (DTI) parameters were assessed in six individual leg muscles. Phosphodiester levels were increased in the NFR and FR tibialis anterior, FR peroneus and FR gastrocnemius lateralis muscles. No clear pattern was visible for the other metabolic indices. Increased T2 SD was found in the majority of FR muscles compared with NFR and healthy control muscles. No differences in average water T2 relaxation times or DTI indices were found between groups. Overall, our results indicate that primarily muscles that are further along in the disease process showed increases in T2 heterogeneity and changes in some metabolic indices. No clear differences were found for the DTI indices between groups.

31 P magnetic resonance spectroscopy in skeletal muscle: Experts' consensus recommendations.

Skeletal muscle phosphorus-31 31 P MRS is the oldest MRS methodology to be applied to in vivo metabolic research. The technical requirements of 31 P MRS in skeletal muscle depend on the research question, and to assess those questions requires understanding both the relevant muscle physiology, and how 31 P MRS methods can probe it. Here we consider basic signal-acquisition parameters related to radio frequency excitation, TR, TE, spectral resolution, shim and localisation. We make specific recommendations for studies of resting and exercising muscle, including magnetisation transfer, and for data processing. We summarise the metabolic information that can be quantitatively assessed with 31 P MRS, either measured directly or derived by calculations that depend on particular metabolic models, and we give advice on potential problems of interpretation. We give expected values and tolerable ranges for some measured quantities, and minimum requirements for reporting acquisition parameters and experimental results in publications. Reliable examination depends on a reproducible setup, standardised preconditioning of the subject, and careful control of potential difficulties, and we summarise some important considerations and potential confounders. Our recommendations include the quantification and standardisation of contraction intensity, and how best to account for heterogeneous muscle recruitment. We highlight some pitfalls in the assessment of mitochondrial function by analysis of phosphocreatine (PCr) recovery kinetics. Finally, we outline how complementary techniques (near-infrared spectroscopy, arterial spin labelling, BOLD and various other MRI and 1 H MRS measurements) can help in the physiological/metabolic interpretation of 31 P MRS studies by providing information about blood flow and oxygen delivery/utilisation. Our recommendations will assist in achieving the fullest possible reliable picture of muscle physiology and pathophysiology.