Magnetic Resonance Imaging

We explore MRI physics, signal formulation, image reconstruction, and tissue modeling. Our goal is to develop innovative MRI methods for understanding fundamental pathophysiological processes of diseases and providing tools to assess tissue biological structures and activities in the brain, joints, and body.

A full list of our publications in this category can be found here.


Research Highlight

Quantitative MRI for Tissue Characterization

MRI is a versatile and powerful imaging method used in clinical diagnosis and scientific research. It provides detailed soft-tissue images and various contrast options without using radiation. This enables the capture of functional, hemodynamic, metabolic, and high-resolution anatomical images for a comprehensive examination. MRI also allows for the quantitative measurement of tissue properties in living organisms. There is a growing interest in using quantitative MRI for routine clinical assessments as it offers increased sensitivity to different diseases and may facilitate early identification of pathologies. Additionally, quantitative measurements can provide valuable information about tissue composition and microstructure. Our research group focuses on quantitative MRI, particularly in the areas of multi-contrast spin-lattice relaxation time (T1), spin-spin relaxation time (T2), spin-lattice relaxation time in the rotating frame (T1ρ), magnetization transfer imaging, ultra-short echo time imaging, and susceptibility imaging in the brain, body, and musculoskeletal system applications.

T1, T2, and T1ρ Relaxation Mapping

The T2 and T1ρ relaxation parameters are extensively studied for assessing tissue degeneration in knee and brain diseases. T2 relaxation time is the most commonly used method and has the most extensive body of literature. It can detect subtle changes in the macromolecular and water content as well as the ultrastructure of tissue extracellular matrix, which are associated with early tissue degeneration. These changes occur before the onset of visible tissue damage. Cartilage T1ρ can indicate changes in the cartilage extracellular matrix with high sensitivity to the loss of important macromolecules, making it a potential tool for early detection of Osteoarthritis. T1 time is also actively researched and has been found to correlate with mechanical property changes of cartilage, different from T2 and T1ρ. It also shows high sensitivity to brain tissue degeneration.

We have developed imaging techniques to thoroughly assess tissue microstructure and composition using advanced quantitative MR relaxometry. Our MR imaging protocol involves performing bi-component T2 mapping of cartilage and brain using two steady-state imaging sequences. We acquire Spoiled gradient-echo (SPGR) scans over multiple flip angles to infer T1 information, and Balanced steady-state free precession (bSSFP) scans over multiple flip angles to infer combined T2/T1 information. These SPGR and bSSFP images are input into a bi-component tissue model to decouple two water pools. When evaluating cartilage matrix, bi-component T2 maps can be used to assess proteoglycan and collagen content. In brain tissue assessment, this technique can evaluate myelin integrity through myelin water characterization.

Bi-component analysis of water fraction (FPG) and T2 relaxation time (T2PG) of water bound to the proteoglycan component is sensitive to the proteoglycan content of cartilage. FPG shows a positive correlation with proteoglycan content, while T2PG shows a positive correlation with proteoglycan denaturation. Trypsin is an enzyme to specifically degrades proteoglycan in the cartilage specimen. The second harmonic generation image shows clumped collagen fibers (d) in the degraded cartilage due to proteoglycan loss caused by trypsin.
FPG and T2PG are sensitive biomarkers to detect early cartilage degeneration. Arrows indicate a cartilage lesion in the trochlea of the cartilage of a 52-year-old male knee Osteoarthritis patient. The FPG has higher diagnostic values than all other T2 parameters for differentiating normal versus abnormal cartilage.
Magnetization Transfer Mapping

Cross-relaxation imaging is a promising technique that can investigate the tissue extracellular matrix by utilizing the magnetization transfer (MT) effect between the water protons and the macromolecular protons of the tissue. This imaging method can offer detailed measurements of the MT parameters at the voxel level, including the fraction of the macromolecular bound protons (f) and the T2 relaxation time of macromolecular bound protons. Our studies have demonstrated the sensitivity of MT parameters to changes in the cartilage extracellular matrix through both ex-vivo and in-vivo research. Additionally, we have successfully imaged the whole brain’s multi-component T2 and MT parameters using an efficient imaging acquisition protocol. The MT parameters obtained through cross-relaxation imaging may provide valuable insights into the collagen fiber network of cartilage, the myelin content of brain tissue, and various macromolecule content in other parts of the body. By combining advanced relaxation imaging and MT imaging, we can develop sensitive and specific imaging biomarkers to comprehensively assess the contents of macromolecules in disease formation and progression.

Ongoing development in the lab aims to acquire multiple contrast MT and relaxation parameters simultaneously. This is achieved through an integrated acquisition framework that compensates for system imperfections and motion, while maintaining a super high spatial resolution, low signal-to-noise ratio, and short scan time.

The myelin water fraction (MWF) map obtained from mcDESPOT and mcRISE (proposed) and f map obtained from mcRISE for one slice of healthy volunteer is shown in a. The corresponding histograms in b show the overestimated MWF in mcDESPOT compared to mcRISE.
(a) Macromolecular fraction f, (b) fundamental rate constant k and (c) T1 maps obtained from MT phantom. (d) Plotting f as a function of Agar percentage shows a linear relationship, as expected. (e) The fundamental rate constant is independent of Agar percentage, corroborating previous studies. (f) Comparing T1 estimated using previous vFA and our proposed BTS shows that vFA consistently underestimates due to overlooking MT effects, and the degree of bias increases with macromolecular fraction. Representative fits of (g) 2%, (h) 4%, (i) 8% Agar and (j) egg white taken from pixels indicated by circles in (a). The difference between BL (blue) and BTS (orange) curves increases with increasing f.
Ultra-short Echo Time Imaging

Musculoskeletal tissues like tendons, ligaments, and meniscus contain a lot of organized collagen fibers, which make them appear dark in regular MRI scans because their signal decays very quickly. Recently, ultrashort echo time (UTE) techniques have been developed to capture this rapid signal decay in these tissues. By using multiple very short echoes (as short as 0.008 msec), these techniques can calculate the ultra-short echo time T2* (UTE-T2*) relaxation time of tendons, ligaments, meniscus, and cortical bone. We have created new multi-component UTE-T2* acquisition and image processing methods to improve the imaging ability of UTE for better understanding tissue properties despite low signal and various system confounding factors.

Sagittal images through the patellar tendon in a 27-year-old healthy male volunteer at all 16 echoes acquired using the 3D-Cones UTE-T2* mapping sequence. There is little MRI signal in the patellar tendon on images with TEs of 4.3 msec and longer. b,c: Signal intensity curves for a homogenous region of interest placed in a sagittal image through the central patellar tendon in a 27-year-old healthy volunteer and a 21-year-old patient with patellar tendinopathy, respectively. Note that there is a visibly improved curve fit of the signal intensity values when using a bicomponent exponential signal model.
UTE-T2* parameter maps in a 27-year-old healthy male volunteer and 32-year-old high-level recreational male runner and volleyball player with grade 2 patellar tendinopathy (ie, the increased signal intensity between 25% and 50% of the axial cross-sectional tendon width). Note the large focal area of increased T2F and decreased FF in the proximal patellar tendon in the patient with patellar tendinopathy with no visible change in T2S (arrows).
Rapid MRI using Compressed Sensing, Parallel Imaging and Non-Cartesian Acquisition

We are researching different techniques to speed up MRI scans, such as compressed sensing, parallel imaging, and a combination of both. Additionally, we are developing new MRI sequences to improve non-Cartesian k-space sampling, aiming for efficient scanning, fewer imaging artifacts, three-dimensional volumetric imaging, and better reconstruction for brain, knee, and body imaging.

Rapid 3D Knee MRI using Compressed Sensing

Commercially available on most MRI vendor platforms, three-dimensional fast spin-echo (3D-FSE-CUBE) sequences can acquire thin continuous slices through joints. These slices can be reformatted in any orientation, eliminating the need to repeat sequences with identical tissue contrast in multiple planes. The use of 3D-FSE-CUBE sequences in clinical practice could significantly decrease MRI examination times, improving patient comfort and increasing the clinical efficiency of the MRI scanner.

Currently, 3D-FSE-CUBE sequences have long scan times to achieve high isotropic resolution. Compressed sensing (CS) is a method that could reduce the scan time by acquiring less image data through k-space undersampling. We conducted studies to explore the possibility of using CS to speed up 3D-FSE-CUBE imaging of the knee and to determine the best imaging parameters of CS for improved image quality in assessing various joint structures.

(a) Estimated low-frequency noise standard deviation maps for two healthy volunteers (#1 a central sagittal slice through the middle of the knee joint in a 30-year-old male and #2 a sagittal slice through the lateral femoral condyle in a 29-year-old female) for CUBE and CUBE-CS and corresponding (b) Color-coded noise amplification NA factor maps superimposed on top of the CUBE source image. (c) NA factor histograms show different noise amplification distributions for these two subjects, but both histograms have a mean value close to 1, indicating no noise amplification.
(a,b) CUBE and CUBE-CS images in a 26-year-old male show a similar appearance of an anterior cruciate ligament tear (arrows). (c,d) CUBE and CUBE-CS images in a 46-year-old female show a similar appearance of a posterior horn medial meniscus tear (arrows).
Novel Non-Cartesian MRI

Our group is working closely with Dr. Li Feng from New York University to develop, optimize, and evaluate non-Cartesian imaging in various clinical applications. We are particularly interested in the GRASP imaging technique, originally developed by Dr. Feng and his colleagues at New York University. We are exploring the integration of GRASP and its variants with our AI-based reconstruction techniques, including SANTIS, MANTIS, and RELAX, to achieve faster image generation, improved quality, accurate quantification, and enhanced spatial-temporal resolution while ensuring algorithm convergence.

The GRASP project represents a decade of innovation by our team consisting of imaging scientists, clinicians, and our industry partners. The GRASP paper was announced as the third most-cited MRM paper at the 2017 ISMRM annual meeting, and the XD-GRASP paper was announced as the top most-cited MRM paper at the 2019 ISMRM annual meeting. With the rise of Artificial Intelligence in recent years, we are now aiming to integrate GRASP MRI with deep learning approaches to enable further improvement in reconstruction quality and speed, as well as new uses of this imaging framework. The initial feasibility of deep-learning-enabled golden-angle radial MRI has been demonstrated by us with a technique called SANTIS (see Deep Learning for Rapid MRI), and we are in the process of developing new quantitative imaging methods based on a combination of GRASP with deep learning.


by Dr. Li Feng
IR-prepared stack-of-stars acquisitions. The imaging sequence was developed based on a stack-of-stars 3D GRE sequence (RAVE). A, An adiabatic non-selective 180° IR pulse is periodically played-out to achieve magnetization preparation. After each IR pulse, a series of radial stacks rotated by a pre-defined rotation scheme is acquired until the magnetization reaches a steady state. B, After synchronizing all the acquired repetitions, a composite IR-prepared dynamic image series can be generated where N consecutive golden-angle rotations form k-space at each time point to ensure uniform coverage. C, The IR-prepared stack-of-stars sequence can also be performed for multi-echo acquisitions, where the rotating angles for different echoes are the same, and the user can select the number of echoes.
Comparison of brain T1 maps obtained from MP2RAGE and MP-GRASP (proposed ) in one volunteer. The T1 maps are visually comparable except for the CSF and the skull region. The linear regression shows that mean T1 values across all the subjects exhibit a good correlation (R2 = 0.955). The Bland-Altman plot suggests that MP2RAGE yielded lower T1 values than MP-GRASP.