Physics-guided self-supervised learning for retrospective T1 and T2 mapping from conventional weighted brain MRI: Technical developments and initial validation in glioblastoma.

PURPOSE: To develop a self-supervised learning method to retrospectively estimate T1 and T2 values from clinical weighted MRI. METHODS: A self-supervised learning approach was constructed to estimate T1, T2, and proton density maps from conventional T1- and T2-weighted images. MR physics models were employed to regenerate the weighted images from the network outputs, and the network was optimized based on loss calculated between the synthesized and input weighted images, alongside additional constraints based on prior information. The method was evaluated on healthy volunteer data, with conventional mapping as references. The reproducibility was examined on two 3.0T scanners. Performance in tumor characterization was inspected by applying the method to a public glioblastoma dataset. RESULTS: For T1 and T2 estimation from three weighted images (T1 MPRAGE, T1 gradient echo sequences, and T2 turbo spin echo), the deep learning method achieved global voxel-wise error ≤9% in brain parenchyma and regional error ≤12.2% in six types of brain tissues. The regional measurements obtained from two scanners showed mean differences ≤2.4% and correlation coefficients >0.98, demonstrating excellent reproducibility. In the 50 glioblastoma patients, the retrospective quantification results were in line with literature reports from prospective methods, and the T2 values were found to be higher in tumor regions, with sensitivity of 0.90 and specificity of 0.92 in a voxel-wise classification task between normal and abnormal regions. CONCLUSION: The self-supervised learning method is promising for retrospective T1 and T2 quantification from clinical MR images, with the potential to improve the availability of quantitative MRI and facilitate brain tumor characterization.

Retrospective motion correction for cardiac multi-parametric mapping with dictionary matching-based image synthesis and a low-rank constraint.

PURPOSE: To develop a model-based motion correction (MoCo) method that does not need an analytical signal model to improve the quality of cardiac multi-parametric mapping. METHODS: The proposed method constructs a hybrid loss that includes a dictionary-matching loss and a signal low-rankness loss, where the former registers the multi-contrast original images to a set of motion-free synthetic images and the latter forces the deformed images to be spatiotemporally coherent. We compared the proposed method with non-MoCo, a pairwise registration method (Pairwise-MI), and a groupwise registration method (pTVreg) via a free-breathing Multimapping dataset of 15 healthy subjects, both quantitatively and qualitatively. RESULTS: The proposed method achieved the lowest contour tracking errors (epicardium: 2.00 ± 0.39 mm vs 4.93 ± 2.29 mm, 3.50 ± 1.26 mm, and 2.61 ± 1.00 mm, and endocardium: 1.84 ± 0.34 mm vs 4.93 ± 2.40 mm, 3.43 ± 1.27 mm, and 2.55 ± 1.09 mm for the proposed method, non-MoCo, Pairwise-MI, and pTVreg, respectively; all p < 0.01) and the lowest dictionary matching errors among all methods. The proposed method also achieved the highest scores on the visual quality of mapping (T1: 4.74 ± 0.33 vs 2.91 ± 0.82, 3.58 ± 0.87, and 3.97 ± 1.05, and T2: 4.48 ± 0.56 vs 2.59 ± 0.81, 3.56 ± 0.93, and 4.14 ± 0.80 for the proposed method, non-MoCo, Pairwise-MI, and pTVreg, respectively; all p < 0.01). Finally, the proposed method had similar T1 and T2 mean values and SDs relative to the breath-hold reference in nearly all myocardial segments, whereas all other methods led to significantly different T1 and T2 measures and increases of SDs in multiple segments. CONCLUSION: The proposed method significantly improves the motion correction accuracy and mapping quality compared with non-MoCo and alternative image-based methods.

3D B1+ corrected simultaneous myocardial T1 and T1ρ mapping with subject-specific respiratory motion correction and water-fat separation.

PURPOSE: To develop a 3D free-breathing cardiac multi-parametric mapping framework that is robust to confounders of respiratory motion, fat, and B1+ inhomogeneities and validate it for joint myocardial T1 and T1ρ mapping at 3T. METHODS: An electrocardiogram-triggered sequence with dual-echo Dixon readout was developed, where nine cardiac cycles were repeatedly acquired with inversion recovery and T1ρ preparation pulses for T1 and T1ρ sensitization. A subject-specific respiratory motion model relating the 1D diaphragmatic navigator to the respiration-induced 3D translational motion of the heart was constructed followed by respiratory motion binning and intra-bin 3D translational and inter-bin non-rigid motion correction. Spin history B1+ inhomogeneities were corrected with optimized dual flip angle strategy. After water-fat separation, the water images were matched to the simulated dictionary for T1 and T1ρ quantification. Phantoms and 10 heathy subjects were imaged to validate the proposed technique. RESULTS: The proposed technique achieved strong correlation (T1: R2 = 0.99; T1ρ: R2 = 0.98) with the reference measurements in phantoms. 3D cardiac T1 and T1ρ maps with spatial resolution of 2 × 2 × 4 mm were obtained with scan time of 5.4 ± 0.5 min, demonstrating comparable T1 (1236 ± 59 ms) and T1ρ (50.2 ± 2.4 ms) measurements to 2D separate breath-hold mapping techniques. The estimated B1+ maps showed spatial variations across the left ventricle with the septal and inferior regions being 10%-25% lower than the anterior and septal regions. CONCLUSION: The proposed technique achieved efficient 3D joint myocardial T1 and T1ρ mapping at 3T with respiratory motion correction, spin history B1+ correction and water-fat separation.

A Novel Workflow for Electrophysiology Studies in Patients with Brugada Syndrome.

Brugada Syndrome (BrS) is a primary electrical epicardial disease characterized by ST-segment elevation followed by a negative T-wave in the right precordial leads on the surface electrocardiogram (ECG), also known as the 'type 1' ECG pattern. The risk stratification of asymptomatic individuals with spontaneous type 1 ECG pattern remains challenging. Clinical and electrocardiographic prognostic markers are known. As none of these predictors alone is highly reliable in terms of arrhythmic prognosis, several multi-factor risk scores have been proposed for this purpose. This article presents a new workflow for processing endocardial signals acquired with high-density RV electro-anatomical mapping (HDEAM) from BrS patients. The workflow, which relies solely on Matlab software, calculates various electrical parameters and creates multi-parametric maps of the right ventricle. The workflow, but it has already been employed in several research studies involving patients carried out by our group, showing its potential positive impact in clinical studies. Here, we will provide a technical description of its functionalities, along with the results obtained on a BrS patient who underwent an endocardial HDEAM.

Free-breathing simultaneous native myocardial T1, T2 and T1ρ mapping with Cartesian acquisition and dictionary matching.

BACKGROUND: T1, T2 and T1ρ are well-recognized parameters for quantitative cardiac MRI. Simultaneous estimation of these parameters allows for comprehensive myocardial tissue characterization, such as myocardial fibrosis and edema. However, conventional techniques either quantify the parameters individually with separate breath-hold acquisitions, which may result in unregistered parameter maps, or estimate multiple parameters in a prolonged breath-hold acquisition, which may be intolerable to patients. We propose a free-breathing multi-parametric mapping (FB-MultiMap) technique that provides co-registered myocardial T1, T2 and T1ρ maps in a single efficient acquisition. METHODS: The proposed FB-MultiMap performs electrocardiogram-triggered single-shot Cartesian acquisition over 16 consecutive cardiac cycles, where inversion, T2 and T1ρ preparations are introduced for varying contrasts. A diaphragmatic navigator was used for prospective through-plane motion correction and the in-plane motion was corrected retrospectively with a group-wise image registration method. Quantitative mapping was conducted through dictionary matching of the motion corrected images, where the subject-specific dictionary was created using Bloch simulations for a range of T1, T2 and T1ρ values, as well as B1 factors to account for B1 inhomogeneities. The FB-MultiMap was optimized and validated in numerical simulations, phantom experiments, and in vivo imaging of 15 healthy subjects and six patients with suspected cardiac diseases. RESULTS: The phantom T1, T2 and T1ρ values estimated with FB-MultiMap agreed well with reference measurements with no dependency on heart rate. In healthy subjects, FB-MultiMap T1 was higher than MOLLI T1 mapping (1218 ± 50 ms vs. 1166 ± 38 ms, p < 0.001). The myocardial T2 and T1ρ estimated with FB-MultiMap were lower compared to the mapping with T2- or T1ρ-prepared 2D balanced steady-state free precession (T2: 41.2 ± 2.8 ms vs. 42.5 ± 3.1 ms, p = 0.06; T1ρ: 45.3 ± 4.4 ms vs. 50.2 ± 4.0, p < 0.001). The pathological changes in myocardial parameters measured with FB-MultiMap were consistent with conventional techniques in all patients. CONCLUSION: The proposed free-breathing multi-parametric mapping technique provides co-registered myocardial T1, T2 and T1ρ maps in 16 heartbeats, achieving similar mapping quality to conventional breath-hold mapping methods.

Assessment of MRF for simultaneous T1 and T2 quantification and water-fat separation in the liver at 0.55 T.

OBJECTIVE: The goal of this work was to assess the feasibility of performing MRF in the liver on a 0.55 T scanner and to examine the feasibility of water-fat separation using rosette MRF at 0.55 T. MATERIALS AND METHODS: Spiral and rosette MRF sequences were implemented on a commercial 0.55 T scanner. The accuracy of both sequences in T1 and T2 quantification was validated in the ISMRM/NIST system phantom. The efficacy of rosette MRF in water-fat separation was evaluated in simulations and water/oil phantoms. Both spiral and rosette MRF were performed in the liver of healthy subjects. RESULTS: In the ISMRM/NIST phantom, both spiral and rosette MRF achieved good agreement with reference values in T1 and T2 measurements. In addition, rosette MRF enables water-fat separation and can generate water- and fat- specific T1 maps, T2 maps, and proton density images from the same dataset for a spatial resolution of 1.56 × 1.56 × 5mm3 within the acquisition time of 15 s. CONCLUSION: It is feasible to measure T1 and T2 simultaneously in the liver using MRF on a 0.55 T system with lower performance gradients compared to state-of-the-art 1.5 T and 3 T systems within an acquisition time of 15 s. In addition, rosette MRF enables water-fat separation along with T1 and T2 quantification with no time penalty.

Single-shot multi-parametric mapping based on multiple overlapping-echo detachment (MOLED) imaging.

Multi-parametric quantitative magnetic resonance imaging (mqMRI) allows the characterization of multiple tissue properties non-invasively and has shown great potential to enhance the sensitivity of MRI measurements. However, real-time mqMRI during dynamic physiological processes or general motions remains challenging. To overcome this bottleneck, we propose a novel mqMRI technique based on multiple overlapping-echo detachment (MOLED) imaging, termed MQMOLED, to enable mqMRI in a single shot. In the data acquisition of MQMOLED, multiple MR echo signals with different multi-parametric weightings and phase modulations are generated and acquired in the same k-space. The k-space data is Fourier transformed and fed into a well-trained neural network for the reconstruction of multi-parametric maps. We demonstrated the accuracy and repeatability of MQMOLED in simultaneous mapping apparent proton density (APD) and any two parameters among T2, T2*, and apparent diffusion coefficient (ADC) in 130-170 ms. The abundant information delivered by the multiple overlapping-echo signals in MQMOLED makes the technique potentially robust to system imperfections, such as inhomogeneity of static magnetic field or radiofrequency field. Benefitting from the single-shot feature, MQMOLED exhibits a strong motion tolerance to the continuous movements of subjects. For the first time, it captured the synchronous changes of ADC, T2, and T1-weighted APD in contrast-enhanced perfusion imaging on patients with brain tumors, providing additional information about vascular density to the hemodynamic parametric maps. We expect that MQMOLED would promote the development of mqMRI technology and greatly benefit the applications of mqMRI, including therapeutics and analysis of metabolic/functional processes.

Multi-parametric artificial neural network fitting of phase-cycled balanced steady-state free precession data.

PURPOSE: Standard relaxation time quantification using phase-cycled balanced steady-state free precession (bSSFP), eg, motion-insensitive rapid configuration relaxometry (MIRACLE), is subject to a considerable underestimation of tissue T1 and T2 due to asymmetric intra-voxel frequency distributions. In this work, an artificial neural network (ANN) fitting approach is proposed to simultaneously extract accurate reference relaxation times (T1 , T2 ) and robust field map estimates ( B1+ , ΔB0 ) from the bSSFP profile. METHODS: Whole-brain bSSFP data acquired at 3T were used for the training of a feedforward ANN with N = 12, 6, and 4 phase-cycles. The magnitude and phase of the Fourier transformed complex bSSFP frequency response served as input and the multi-parametric reference set [T1 , T2 , B1+ , ∆B0 ] as target. The ANN predicted relaxation times were validated against the target and MIRACLE. RESULTS: The ANN prediction of T1 and T2 for trained and untrained data agreed well with the reference, even for only four acquired phase-cycles. In contrast, relaxometry based on 4-point MIRACLE was prone to severe off-resonance-related artifacts. ANN predicted B1+ and ∆B0 maps showed the expected spatial inhomogeneity patterns in high agreement with the reference measurements for 12-point, 6-point, and 4-point bSSFP phase-cycling schemes. CONCLUSION: ANNs show promise to provide accurate brain tissue T1 and T2 values as well as reliable field map estimates. Moreover, the bSSFP acquisition can be accelerated by reducing the number of phase-cycles while still delivering robust T1 , T2 , B1+ , and ∆B0 estimates.

Stimulated echo based mapping (STEM) of T1 , T2 , and apparent diffusion coefficient: validation and protocol optimization.

PURPOSE: To present a stimulated-echo based mapping (STEM) approach for simultaneous T1 , T2 , and ADC mapping. METHODS: Diffusion-weighted stimulated-echo images with various combinations of mixing time (TM), TE, and b-value were acquired to enable simultaneous mapping of T1 , T2 , and ADC. The proposed STEM method was performed by densely sampling the TM-TE-b space in a phantom and in brain and prostate of healthy volunteers. T1 , T2 , and ADC from STEM were compared to reference mapping methods. Additionally, protocol optimization was performed to enable rapid STEM acquisition within 2 min by sparsely sampling the TM-TE-b space. The T1 , T2 , and ADC measurements from rapid acquisitions were compared to the densely sampled STEM for evaluation. Finally, a patient with biopsy-proven high-risk prostate cancer was imaged to demonstrate the ability of STEM to differentiate cancer and healthy tissues. RESULTS: Relative to the reference measurements, densely sampled STEM provided accurate quantitative T1 , T2 , and ADC mapping in phantoms (R2 = 0.999, slope between 0.97-1.03), as well as in brain and prostate. Further, the T1 , T2 , and ADC measurements from the optimized rapid STEM acquisitions agreed closely with densely sampled STEM. Finally, STEM showed decreased T2 and ADC in prostate cancer compared to healthy prostate tissue. CONCLUSION: STEM provides accurate simultaneous mapping of T1 , T2 , and ADC. This method may enable rapid and accurate multi-parametric tissue characterization for clinical and research applications.

Specific white matter tissue microstructure changes associated with obesity.

Obesity-related structural brain alterations point to a consistent reduction in gray matter with increasing body mass index (BMI) but changes in white matter have proven to be more complex and less conclusive. Hence, more recently diffusion tensor imaging (DTI) has been employed to investigate microstructural changes in white matter structure. Altogether, these studies have mostly shown a loss of white matter integrity with obesity-related factors in several brain regions. However, the variety of these obesity-related factors, including inflammation and dyslipidemia, resulted in competing influences on the DTI indices. To increase the specificity of DTI results, we explored specific brain tissue properties by combining DTI with quantitative multi-parameter mapping in lean, overweight and obese young adults. By means of multi-parameter mapping, white matter structures showed differences in MRI parameters consistent with reduced myelin, increased water and altered iron content with increasing BMI in the superior longitudinal fasciculus, anterior thalamic radiation, internal capsule and corpus callosum. BMI-related changes in DTI parameters revealed mainly alterations in mean and axial diffusivity with increasing BMI in the corticospinal tract, anterior thalamic radiation and superior longitudinal fasciculus. These alterations, including mainly fiber tracts linking limbic structures with prefrontal regions, could potentially promote accelerated aging in obese individuals leading to an increased risk for cognitive decline.

Automatic Brain Tissue and Lesion Segmentation and Multi-Parametric Mapping of Contrast-Enhancing Gliomas without the Injection of Contrast Agents: A Preliminary Study.

This study aimed to develop a rapid, 1 mm3 isotropic resolution, whole-brain MRI technique for automatic lesion segmentation and multi-parametric mapping without using contrast by continuously applying balanced steady-state free precession with inversion pulses throughout incomplete inversion recovery in a single 6 min scan. Modified k-means clustering was performed for automatic brain tissue and lesion segmentation using distinct signal evolutions that contained mixed T1/T2/magnetization transfer properties. Multi-compartment modeling was used to derive quantitative multi-parametric maps for tissue characterization. Fourteen patients with contrast-enhancing gliomas were scanned with this sequence prior to the injection of a contrast agent, and their segmented lesions were compared to conventionally defined manual segmentations of T2-hyperintense and contrast-enhancing lesions. Simultaneous T1, T2, and macromolecular proton fraction maps were generated and compared to conventional 2D T1 and T2 mapping and myelination water fraction mapping acquired with MAGiC. The lesion volumes defined with the new method were comparable to the manual segmentations (r = 0.70, p < 0.01; t-test p > 0.05). The T1, T2, and macromolecular proton fraction mapping values of the whole brain were comparable to the reference values and could distinguish different brain tissues and lesion types (p < 0.05), including infiltrating tumor regions within the T2-lesion. Highly efficient, whole-brain, multi-contrast imaging facilitated automatic lesion segmentation and quantitative multi-parametric mapping without contrast, highlighting its potential value in the clinic when gadolinium is contraindicated.

Myocardial tissue imaging with cardiovascular magnetic resonance.

Alteration in myocardial tissue, such as myocardial fibrosis, edema, inflammation, or accumulation with amyloid, lipids, or iron, has an important role in the cardiac remodeling that leads to diastolic and/or systolic dysfunction and the development of chronic heart failure, increasing the risk of adverse cardiovascular events. Thus, the early detection of changes at myocardial tissue level has great diagnostic and prognostic potential. The gold standard technique to assess these myocardial alterations is endomyocardial biopsy. However, this has been limited to a few patients due to the invasive nature, sampling errors, and its inability to assess the entire myocardium. Cardiovascular magnetic resonance (CMR) has emerged as the gold standard imaging not only for assessing cardiac volume, function quantification, and viability but also for noninvasive myocardial tissue characterization over the past decade. Its ability to characterize myocardial tissue composition is unique among noninvasive imaging modalities in cardiovascular disease. Currently, multi-parametric myocardial characterization with T1, T2, and extracellular volume has the potential to identify and track diffuse pathology in various diseases. In this review article, we present the role of established and emerging CMR techniques in myocardial tissue characterization, with an emphasis on T1 and T2 mapping, in clinical practice.

Brain iron content in systemic iron overload: A beta-thalassemia quantitative MRI study.

OBJECTIVE: Multisystem iron poisoning is a major concern for long-term beta-thalassemia management. Quantitative MRI-based techniques routinely show iron overload in heart, liver, endocrine glands and kidneys. However, data on the brain are conflicting and monitoring of brain iron content is still matter of debate. METHODS: This 3T-MRI study applied a well validated high-resolution whole-brain quantitative MRI assessment of iron content on 47 transfusion-dependent (mean-age: 36.9 ±â€¯10.3 years, 63% females), 23 non-transfusion dependent (mean-age: 29.2 ±â€¯11.7 years, 56% females) and 57 healthy controls (mean-age: 33.9 ±â€¯10.8 years, 65% females). Clinical data, Wechsler Adult Intelligence Scale scores and treatment regimens were recorded. Beside whole-brain R2* analyses, regional R2*-values were extracted in putamen, globus pallidum, caudate nucleus, thalamus and red nucleus; hippocampal volumes were also determined. RESULTS: Regional analyses yielded no significant differences between patients and controls, except in those treated with deferiprone that showed lower R2*-values (p<0.05). Whole-brain analyses of R2*-maps revealed strong age-R2* correlations (r2=0.51) in both groups and clusters of significantly increased R2*-values in beta-thalassemia patients in the hippocampal formations and around the Luschka foramina; transfusion treatment was associated with additional R2* increase in dorsal thalami. Hippocampal formation R2*-values did not correlate with hippocampal volume; hippocampal volume did not differ between patients and controls. All regions with increased R2*-values shared a strict anatomical contiguity with choroid plexuses suggesting a blooming effect as the likely cause of R2* increase, in agreement with the available histopathologic literature evidence. CONCLUSION: According to our MRI findings and the available histopathologic literature evidence, concerns about neural tissue iron overload in beta-thalassemia appear to be unjustified.

Clinical applications of multiparametric CMR in left ventricular hypertrophy.

There are a number of diseases which can increase left ventricular myocardial wall thickness through a number of different mechanisms. Multi-parametric mapping techniques are a new addition to the cardiovascular magnetic resonance (CMR) armoury with a number of potential clinical roles. In this review article, we will explore the role of imaging in left ventricular hypertrophy, and particularly developments in CMR. We focus on ability of CMR to characterize myocardial tissue using multiparametric mapping (native T1, T2 and extracellular volume mapping), to bridge from the microscopic histological domain and into the clinical domain of non-invasive imaging.