Correlated spectroscopic imaging of calf muscle in three spatial dimensions using group sparse reconstruction of undersampled single and multichannel data

Lipid Deposition in Skeletal Muscle Tissues and Its Correlation with Intra-Abdominal Fat: A Pilot Investigation in Type 2 Diabetes Mellitus

Background/Objectives: This study evaluated metabolites and lipid composition in the calf muscles of Type 2 diabetes mellitus (T2DM) patients and age-matched healthy controls using multi-dimensional MR spectroscopic imaging. We also explored the association between muscle metabolites, lipids, and intra-abdominal fat in T2DM. Methods: Participants included 12 T2DM patients (60.3 ± 8.6 years), 9 age-matched healthy controls (AMHC) (60.9 ± 7.8 years), and 10 young healthy controls (YHC) (28.3 ± 1.8 years). We acquired the 2D MR spectra of calf muscles using an enhanced accelerated 5D echo-planar correlated spectroscopic imaging (EP-COSI) technique and abdominal MRI with breath-hold 6-point Dixon sequence. Results: In YHC, choline levels were lower in the gastrocnemius (GAS) and soleus (SOL) muscles but higher in the tibialis anterior (TA) compared to AMHC. YHC also showed a higher unsaturation index (U.I.) of extramyocellular lipids (EMCL) in TA, intramyocellular lipids (IMCL) in GAS, carnosine in SOL, and taurine and creatine in TA. T2DM patients exhibited higher choline in TA and myo-inositol in SOL than AMHC, while triglyceride fat (TGFR2) levels in TA were lower. Correlation analyses indicated associations between IMCL U.I. and various metabolites in muscles with liver, pancreas, and abdominal fat estimates in T2DM. Conclusions: This study highlights distinct muscle metabolite and lipid composition patterns across YHC, AMHC, and T2DM subjects. Associations between IMCL U.I. and abdominal fat depots underscore the interplay between muscle metabolism and adiposity in T2DM. These findings provide new insights into metabolic changes in T2DM and emphasize the utility of advanced MR spectroscopic imaging in characterizing muscle-lipid interactions.

Myosteatosis: diagnostic significance and assessment by imaging approaches

Myosteatosis has emerged as an important concept in muscle health as it is associated with an increased risk of adverse health outcomes, a higher rate of complications, and increased mortality associated with ageing, chronic systemic and neuromuscular diseases, cancer, metabolic syndromes, degenerative events, and trauma. Myosteatosis involves ectopic infiltration of fat into skeletal muscle, and it exhibits a negative correlation with muscle mass, strength, and mobility representing a contributing factor to decreased muscle quality. While myosteatosis serves as an additional biomarker for sarcopenia, cachexia, and metabolic syndromes, it is not synonymous with sarcopenia. Myosteatosis induces proinflammatory changes that contribute to decreased muscle function, compromise mitochondrial function, and increase inflammatory response in muscles. Imaging techniques such as computed tomography (CT), particularly opportunistic abdominal CT scans, and magnetic resonance imaging (MRI) or magnetic resonance spectroscopy (MRS), have been used in both clinical practice and research. And in recent years, ultrasound has emerged as a promising bedside tool for measuring changes in muscle tissue. Various techniques, including CT-based muscle attenuation (MA) and intermuscular adipose tissue (IMAT) quantification, MRI-based proton density fat fraction (PDFF) and T1-T2 mapping, and musculoskeletal ultrasound (MSUS)-based echo intensity (EI) and shear wave elastography (SWE), are accessible in clinical practice and can be used as adjunct biomarkers of myosteatosis to assess various debilitating muscle health conditions. However, a stan¬dard definition of myosteatosis with a thorough understanding of the pathophysiological mechanisms, and a consensus in assessment methods and clinical outcomes has not yet been established. Recent developments in image acquisition and quantification have attempted to develop an appropriate muscle quality index for the assessment of myosteatosis. Additionally, emerging studies on artificial intelligence (AI) may provide further insights into quantification and automated assessment, including MRS analysis. In this review, we discuss the pathophysiological aspects of myosteatosis, all the current imaging techniques and recent advances in imaging assessment as potential biomarkers of myosteatosis, and the most common clinical conditions involved.

Enhanced spectral resolution for correlated spectroscopic imaging using inner-product and covariance transform: a pilot analysis of metabolites and lipids in breast cancer in vivo

Acquisition duration of correlated spectroscopy in vivo can be longer due to a large number of t1 increments along the indirect (F1) dimension. Limited number of t1 increments on the other hand leads to poor spectral resolution along F1. Covariance transformation (CT) instead of Fourier transform along t1 is an alternative way of increasing the resolution of the 2D COSY spectrum. Prospectively undersampled five-dimensional echo-planar correlated spectroscopic imaging (EP-COSI) data from ten malignant patients and ten healthy women were acquired and reconstructed using compressed sensing. The COSY spectrum at each voxel location was then generated using FFT, CT and a variant of CT called Inner Product (IP). Metabolite and lipid ratios were computed with respect to water from unsuppressed one-dimensional spectrum. The effects of t1-ridging artifacts commonly seen with FFT were not observed with CT/IP. Statistically significant differences were observed in the fat cross peaks measured with CT/IP/FFT. Spectral resolution was increased ~ 8.5 times (~ 19.53 Hz in FFT, ~ 2.32 Hz in CT/IP) without affecting the spectral width along F1 was possible with CT/IP. CT and IP enabled substantially increased F1 resolution effectively with significant gain in scan time and reliable measure of unsaturation index as a biomarker for malignant breast cancer.

Ensemble Learning for Breast Cancer Lesion Classification: A Pilot Validation Using Correlated Spectroscopic Imaging and Diffusion-Weighted Imaging

The main objective of this work was to evaluate the application of individual and ensemble machine learning models to classify malignant and benign breast masses using features from two-dimensional (2D) correlated spectroscopy spectra extracted from five-dimensional echo-planar correlated spectroscopic imaging (5D EP-COSI) and diffusion-weighted imaging (DWI). Twenty-four different metabolite and lipid ratios with respect to diagonal fat peaks (1.4 ppm, 5.4 ppm) from 2D spectra, and water and fat peaks (4.7 ppm, 1.4 ppm) from one-dimensional non-water-suppressed (NWS) spectra were used as the features. Additionally, water fraction, fat fraction and water-to-fat ratios from NWS spectra and apparent diffusion coefficients (ADC) from DWI were included. The nine most important features were identified using recursive feature elimination, sequential forward selection and correlation analysis. XGBoost (AUC: 93.0%, Accuracy: 85.7%, F1-score: 88.9%, Precision: 88.2%, Sensitivity: 90.4%, Specificity: 84.6%) and GradientBoost (AUC: 94.3%, Accuracy: 89.3%, F1-score: 90.7%, Precision: 87.9%, Sensitivity: 94.2%, Specificity: 83.4%) were the best-performing models. Conventional biomarkers like choline, myo-Inositol, and glycine were statistically significant predictors. Key features contributing to the classification were ADC, 2D diagonal peaks at 0.9 ppm, 2.1 ppm, 3.5 ppm, and 5.4 ppm, cross peaks between 1.4 and 0.9 ppm, 4.3 and 4.1 ppm, 2.3 and 1.6 ppm, and the triglyceryl-fat cross peak. The results highlight the contribution of the 2D spectral peaks to the model, and they demonstrate the potential of 5D EP-COSI for early breast cancer detection.

Correlated MR spectroscopic imaging of breast cancer to investigate metabolites and lipids: acceleration and compressed sensing reconstruction

Objectives

The main objective of this work was to detect novel biomarkers in breast cancer by spreading the MR spectra over two dimensions in multiple spatial locations using an accelerated 5D EP-COSI technology.

Methods

The 5D EP-COSI data were non-uniformly undersampled with an acceleration factor of 8 and reconstructed using group sparsity-based compressed sensing reconstruction. Different metabolite and lipid ratios were then quantified and statistically analyzed for significance. Linear discriminant models based on the quantified metabolite and lipid ratios were generated. Spectroscopic images of the quantified metabolite and lipid ratios were also reconstructed.

Results

The 2D COSY spectra generated using the 5D EP-COSI technique showed differences among healthy, benign, and malignant tissues in terms of their mean values of metabolite and lipid ratios, especially the ratios of potential novel biomarkers based on unsaturated fatty acids, myo-inositol, and glycine. It is further shown the potential of choline and unsaturated lipid ratio maps, generated from the quantified COSY signals across multiple locations in the breast, to serve as complementary markers of malignancy that can be added to the multiparametric MR protocol. Discriminant models using metabolite and lipid ratios were found to be statistically significant for classifying benign and malignant tumor from healthy tissues.

Conclusions

Accelerated 5D EP-COSI technique demonstrates the potential to detect novel biomarkers such as glycine, myo-inositol, and unsaturated fatty acids in addition to commonly reported choline in breast cancer, and facilitates metabolite and lipid ratio maps which have the potential to play a significant role in breast cancer detection.

Advances in knowledge

This study presents the first evaluation of a multidimensional MR spectroscopic imaging technique for the detection of potentially novel biomarkers based on glycine, myo-inositol, and unsaturated fatty acids, in addition to commonly reported choline. Spatial mapping of choline and unsaturated fatty acid ratios with respect to water in malignant and benign breast masses are also shown. These metabolic characteristics may serve as additional biomarkers for improving the diagnostic and therapeutic evaluation of breast cancer.

Accelerated MR spectroscopic imaging-a review of current and emerging techniques

Over more than 30 years in vivo MR spectroscopic imaging (MRSI) has undergone an enormous evolution from theoretical concepts in the early 1980s to the robust imaging technique that it is today. The development of both fast and efficient sampling and reconstruction techniques has played a fundamental role in this process. State-of-the-art MRSI has grown from a slow purely phase-encoded acquisition technique to a method that today combines the benefits of different acceleration techniques. These include shortening of repetition times, spatial-spectral encoding, undersampling of k-space and time domain, and use of spatial-spectral prior knowledge in the reconstruction. In this way in vivo MRSI has considerably advanced in terms of spatial coverage, spatial resolution, acquisition speed, artifact suppression, number of detectable metabolites and quantification precision. Acceleration not only has been the enabling factor in high-resolution whole-brain 1 H-MRSI, but today is also common in non-proton MRSI (31 P, 2 H and 13 C) and applied in many different organs. In this process, MRSI techniques had to constantly adapt, but have also benefitted from the significant increase of magnetic field strength boosting the signal-to-noise ratio along with high gradient fidelity and high-density receive arrays. In combination with recent trends in image reconstruction and much improved computation power, these advances led to a number of novel developments with respect to MRSI acceleration. Today MRSI allows for non-invasive and non-ionizing mapping of the spatial distribution of various metabolites' tissue concentrations in animals or humans, is applied for clinical diagnostics and has been established as an important tool for neuro-scientific and metabolism research. This review highlights the developments of the last five years and puts them into the context of earlier MRSI acceleration techniques. In addition to 1 H-MRSI it also includes other relevant nuclei and is not limited to certain body regions or specific applications.

Deep learning can accelerate and quantify simulated localized correlated spectroscopy

Nuclear magnetic resonance spectroscopy (MRS) allows for the determination of atomic structures and concentrations of different chemicals in a biochemical sample of interest. MRS is used in vivo clinically to aid in the diagnosis of several pathologies that affect metabolic pathways in the body. Typically, this experiment produces a one dimensional (1D) 1H spectrum containing several peaks that are well associated with biochemicals, or metabolites. However, since many of these peaks overlap, distinguishing chemicals with similar atomic structures becomes much more challenging. One technique capable of overcoming this issue is the localized correlated spectroscopy (L-COSY) experiment, which acquires a second spectral dimension and spreads overlapping signal across this second dimension. Unfortunately, the acquisition of a two dimensional (2D) spectroscopy experiment is extremely time consuming. Furthermore, quantitation of a 2D spectrum is more complex. Recently, artificial intelligence has emerged in the field of medicine as a powerful force capable of diagnosing disease, aiding in treatment, and even predicting treatment outcome. In this study, we utilize deep learning to: (1) accelerate the L-COSY experiment and (2) quantify L-COSY spectra. All training and testing samples were produced using simulated metabolite spectra for chemicals found in the human body. We demonstrate that our deep learning model greatly outperforms compressed sensing based reconstruction of L-COSY spectra at higher acceleration factors. Specifically, at four-fold acceleration, our method has less than 5% normalized mean squared error, whereas compressed sensing yields 20% normalized mean squared error. We also show that at low SNR (25% noise compared to maximum signal), our deep learning model has less than 8% normalized mean squared error for quantitation of L-COSY spectra. These pilot simulation results appear promising and may help improve the efficiency and accuracy of L-COSY experiments in the future.

Phased-array combination of 2D MRS for lipid composition quantification in patients with breast cancer

Lipid composition in breast cancer, a central marker of disease progression, can be non-invasively quantified using 2D MRS method of double quantum filtered correlation spectroscopy (DQF-COSY). The low signal to noise ratio (SNR), arising from signal retention of only 25% and depleted lipids within tumour, demands improvement approaches beyond signal averaging for clinically viable applications. We therefore adapted and examined combination algorithms, designed for 1D MRS, for 2D MRS with both internal and external references. Lipid composition spectra were acquired from 17 breast tumour specimens, 15 healthy female volunteers and 25 patients with breast cancer on a clinical 3 T MRI scanner. Whitened singular value decomposition (WSVD) with internal reference yielded maximal SNR with an improvement of 53.3% (40.3-106.9%) in specimens, 84.4 ± 40.6% in volunteers, 96.9 ± 54.2% in peritumoural adipose tissue and 52.4% (25.1-108.0%) in tumours in vivo. Non-uniformity, as variance of improvement across peaks, was low at 21.1% (13.7-28.1%) in specimens, 5.5% (4.2-7.2%) in volunteers, 6.1% (5.0-9.0%) in peritumoural tissue, and 20.7% (17.4-31.7%) in tumours in vivo. The bias (slope) in improvement ranged from - 1.08 to 0.21%/ppm along the diagonal directions. WSVD is therefore the optimal algorithm for lipid composition spectra with highest SNR uniformly across peaks, reducing acquisition time by up to 70% in patients, enabling clinical applications.

Fast data acquisition techniques in magnetic resonance spectroscopic imaging

Magnetic resonance spectroscopic imaging (MRSI) is an important technique for assessing the spatial variation of metabolites in vivo. The long scan times in MRSI limit clinical applicability due to patient discomfort, increased costs, motion artifacts, and limited protocol flexibility. Faster acquisition strategies can address these limitations and could potentially facilitate increased adoption of MRSI into routine clinical protocols with minimal addition to the current anatomical and functional acquisition protocols in terms of imaging time. Not surprisingly, a lot of effort has been devoted to the development of faster MRSI techniques that aim to capture the same underlying metabolic information (relative metabolite peak areas and spatial distribution) as obtained by conventional MRSI, in greatly reduced time. The gain in imaging time results, in some cases, in a loss of signal-to-noise ratio and/or in spatial and spectral blurring. This review examines the current techniques and advances in fast MRSI in two and three spatial dimensions and their applications. This review categorizes the acceleration techniques according to their strategy for acquisition of the k-space. Techniques such as fast/turbo-spin echo MRSI, echo-planar spectroscopic imaging, and non-Cartesian MRSI effectively cover the full k-space in a more efficient manner per T_R . On the other hand, techniques such as parallel imaging and compressed sensing acquire fewer k-space points and employ advanced reconstruction algorithms to recreate the spatial-spectral information, which maintains statistical fidelity in test conditions (ie no statistically significant differences on voxel-wise comparisions) with the fully sampled data. The advantages and limitations of each state-of-the-art technique are reviewed in detail, concluding with a note on future directions and challenges in the field of fast spectroscopic imaging.

Prior-knowledge Fitting of Accelerated Five-dimensional Echo Planar J-resolved Spectroscopic Imaging: Effect of Nonlinear Reconstruction on Quantitation

¹H Magnetic Resonance Spectroscopic imaging (SI) is a powerful tool capable of investigating metabolism in vivo from mul- tiple regions. However, SI techniques are time consuming, and are therefore difficult to implement clinically. By applying non-uniform sampling (NUS) and compressed sensing (CS) reconstruction, it is possible to accelerate these scans while re- taining key spectral information. One recently developed method that utilizes this type of acceleration is the five-dimensional echo planar J-resolved spectroscopic imaging (5D EP-JRESI) sequence, which is capable of obtaining two-dimensional (2D) spectra from three spatial dimensions. The prior-knowledge fitting (ProFit) algorithm is typically used to quantify 2D spectra in vivo, however the effects of NUS and CS reconstruction on the quantitation results are unknown. This study utilized a simulated brain phantom to investigate the errors introduced through the acceleration methods. Errors (normalized root mean square error >15%) were found between metabolite concentrations after twelve-fold acceleration for several low concentra- tion (<2 mM) metabolites. The Cramér Rao lower bound% (CRLB%) values, which are typically used for quality control, were not reflective of the increased quantitation error arising from acceleration. Finally, occipital white (OWM) and gray (OGM) human brain matter were quantified in vivo using the 5D EP-JRESI sequence with eight-fold acceleration.

Quantifying fat replacement of muscle by quantitative MRI in muscular dystrophy

The muscular dystrophies are rare orphan diseases, characterized by progressive muscle weakness: the most common and well known is Duchenne muscular dystrophy which affects young boys and progresses quickly during childhood. However, over 70 distinct variants have been identified to date, with different rates of progression, implications for morbidity, mortality, and quality of life. There are presently no curative therapies for these diseases, but a range of potential therapies are presently reaching the stage of multi-centre, multi-national first-in-man clinical trials. There is a need for sensitive, objective end-points to assess the efficacy of the proposed therapies. Present clinical measurements are often too dependent on patient effort or motivation, and lack sensitivity to small changes, or are invasive. Quantitative MRI to measure the fat replacement of skeletal muscle by either chemical shift imaging methods (Dixon or IDEAL) or spectroscopy has been demonstrated to provide such a sensitive, objective end-point in a number of studies. This review considers the importance of the outcome measures, discusses the considerations required to make robust measurements and appropriate quality assurance measures, and draws together the existing literature for cross-sectional and longitudinal cohort studies using these methods in muscular dystrophy.

Pilot Assessment of Brain Metabolism in Perinatally HIV-Infected Youths Using Accelerated 5D Echo Planar J-Resolved Spectroscopic Imaging

Purpose

To measure cerebral metabolite levels in perinatally HIV-infected youths and healthy controls using the accelerated five dimensional (5D) echo planar J-resolved spectroscopic imaging (EP-JRESI) sequence, which is capable of obtaining two dimensional (2D) J-resolved spectra from three spatial dimensions (3D).

Materials and Methods

After acquisition and reconstruction of the 5D EP-JRESI data, T₁-weighted MRIs were used to classify brain regions of interest for HIV patients and healthy controls: right frontal white (FW), medial frontal gray (FG), right basal ganglia (BG), right occipital white (OW), and medial occipital gray (OG). From these locations, respective J-resolved and TE-averaged spectra were extracted and fit using two different quantitation methods. The J-resolved spectra were fit using prior knowledge fitting (ProFit) while the TE-averaged spectra were fit using the advanced method for accurate robust and efficient spectral fitting (AMARES).

Results

Quantitation of the 5D EP-JRESI data using the ProFit algorithm yielded significant metabolic differences in two spatial locations of the perinatally HIV-infected youths compared to controls: elevated NAA/(Cr+Ch) in the FW and elevated Asp/(Cr+Ch) in the BG. Using the TE-averaged data quantified by AMARES, an increase of Glu/(Cr+Ch) was shown in the FW region. A strong negative correlation (r < -0.6) was shown between tCh/(Cr+Ch) quantified using ProFit in the FW and CD4 counts. Also, strong positive correlations (r > 0.6) were shown between Asp/(Cr+Ch) and CD4 counts in the FG and BG.

Conclusion

The complimentary results using ProFit fitting of J-resolved spectra and AMARES fitting of TE-averaged spectra, which are a subset of the 5D EP-JRESI acquisition, demonstrate an abnormal energy metabolism in the brains of perinatally HIV-infected youths. This may be a result of the HIV pathology and long-term combinational anti-retroviral therapy (cART). Further studies of larger perinatally HIV-infected cohorts are necessary to confirm these findings.

Sparse Reconstruction Techniques in Magnetic Resonance Imaging: Methods, Applications, and Challenges to Clinical Adoption

The family of sparse reconstruction techniques, including the recently introduced compressed sensing framework, has been extensively explored to reduce scan times in magnetic resonance imaging (MRI). While there are many different methods that fall under the general umbrella of sparse reconstructions, they all rely on the idea that a priori information about the sparsity of MR images can be used to reconstruct full images from undersampled data. This review describes the basic ideas behind sparse reconstruction techniques, how they could be applied to improve MRI, and the open challenges to their general adoption in a clinical setting. The fundamental principles underlying different classes of sparse reconstructions techniques are examined, and the requirements that each make on the undersampled data outlined. Applications that could potentially benefit from the accelerations that sparse reconstructions could provide are described, and clinical studies using sparse reconstructions reviewed. Lastly, technical and clinical challenges to widespread implementation of sparse reconstruction techniques, including optimization, reconstruction times, artifact appearance, and comparison with current gold standards, are discussed.