Measurement Guidance in Diffusion Models: Insight from Medical Image Synthesis.

In the field of healthcare, the acquisition of sample is usually restricted by multiple considerations, including cost, labor- intensive annotation, privacy concerns, and radiation hazards, therefore, synthesizing images-of-interest is an important tool to data augmentation. Diffusion models have recently attained state-of-the-art results in various synthesis tasks, and embedding energy functions has been proved that can effectively guide the pre-trained model to synthesize target samples. However, we notice that current method development and validation are still limited to improving indicators, such as Fréchet Inception Distance score (FID) and Inception Score (IS), and have not provided deeper investigations on downstream tasks, like disease grading and diagnosis. Moreover, existing classifier guidance which can be regarded as a special case of energy function can only has a singular effect on altering the distribution of the synthetic dataset. This may contribute to in-distribution synthetic sample that has limited help to downstream model optimization. All these limitations remind that we still have a long way to go to achieve controllable generation. In this work, we first conducted an analysis on previous guidance as well as its contributions on further applications from the perspective of data distribution. To synthesize samples which can help downstream applications, we then introduce uncertainty guidance in each sampling step and design an uncertainty-guided diffusion models. Extensive experiments on four medical datasets, with ten classic networks trained on the augmented sample sets provided a comprehensive evaluation on the practical contributions of our methodology. Furthermore, we provide a theoretical guarantee for general gradient guidance in diffusion models, which would benefit future research on investigating other forms of measurement guidance for specific generative tasks. Codes and models are available at: https://github.com/yangqy1110/MGDM.

Exercise-induced calf muscle hyperemia quantified with dynamic blood oxygen level-dependent (BOLD) imaging.

Muscle hyperemia in exercise is usually the combined result of increased cardiac output and local muscle vasodilation, with the latter reflecting muscle's capacity for increased blood perfusion to support exercise. In this study, we aim to quantify muscle's vasodilation capability with dynamic BOLD imaging. A deoxyhemoglobin-kinetics model is proposed to analyze dynamic BOLD signals acquired during exercise recovery, deriving a hyperemia index (HI) for a muscle group of interest. We demonstrated the method's validity with calf muscles of healthy subjects who performed plantar flexion for muscle stimulation. In a test with exercise load incrementally increasing from 0 to 16 lbs., gastrocnemius HI showed considerable variance among the 4 subjects, but with a consistent trend, i.e. low at light load (e.g. 0-6 lbs) and linearly increasing at heavy load. The high variability among different subjects was confirmed with the other 10 subjects who exercised with a same moderate load of 8 lbs., with coefficient of variance among subjects' medial gastrocnemius 87.8%, lateral gastrocnemius 111.8% and soleus 132.3%. These findings align with the fact that intensive exercise induces high muscle hyperemia, but a comparison among different subjects is hard to make, presumably due to the subjects' different rate of oxygen utilization. For the same 10 subjects who exercised with load of 8 lbs., we also performed dynamic contrast enhanced (DCE) MRI to measure muscle perfusion (F). With a moderate correlation of 0.654, HI and F displayed three distinctive responses of calf muscles: soleus of all the subjects were in the cluster of low F and low HI, and gastrocnemius of most subjects had high F and either low or high HI. This finding suggests that parameter F encapsulates blood flow through vessels of all sizes, but BOLD-derived HI focuses on capillary flow and therefore is a more specific indicator of muscle vasodilation. In conclusion, the proposed hyperemia index has the potential of quantitatively assessing muscle vasodilation induced with exercise.

Utilization of CRISPR-Cas genome editing technology in filamentous fungi: function and advancement potentiality.

Filamentous fungi play a crucial role in environmental pollution control, protein secretion, and the production of active secondary metabolites. The evolution of gene editing technology has significantly improved the study of filamentous fungi, which in the past was laborious and time-consuming. But recently, CRISPR-Cas systems, which utilize small guide RNA (sgRNA) to mediate clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated proteins (Cas), have demonstrated considerable promise in research and application for filamentous fungi. The principle, function, and classification of CRISPR-Cas, along with its application strategies and research progress in filamentous fungi, will all be covered in the review. Additionally, we will go over general matters to take into account when editing a genome with the CRISPR-Cas system, including the creation of vectors, different transformation methodologies, multiple editing approaches, CRISPR-mediated transcriptional activation (CRISPRa) or interference (CRISPRi), base editors (BEs), and Prime editors (PEs).

Corrigendum: Identification of ACKR4 as an immune checkpoint in pulmonary arterial hypertension.

[This corrects the article DOI: 10.3389/fimmu.2023.1153573.].

A motion-corrected deep-learning reconstruction framework for accelerating whole-heart magnetic resonance imaging in patients with congenital heart disease.

BACKGROUND: Cardiovascular magnetic resonance (CMR) is an important imaging modality for the assessment and management of adult patients with congenital heart disease (CHD). However, conventional techniques for three-dimensional (3D) whole-heart acquisition involve long and unpredictable scan times and methods that accelerate scans via k-space undersampling often rely on long iterative reconstructions. Deep-learning-based reconstruction methods have recently attracted much interest due to their capacity to provide fast reconstructions while often outperforming existing state-of-the-art methods. In this study, we sought to adapt and validate a non-rigid motion-corrected model-based deep learning (MoCo-MoDL) reconstruction framework for 3D whole-heart MRI in a CHD patient cohort. METHODS: The previously proposed deep-learning reconstruction framework MoCo-MoDL, which incorporates a non-rigid motion-estimation network and a denoising regularization network within an unrolled iterative reconstruction, was trained in an end-to-end manner using 39 CHD patient datasets. Once trained, the framework was evaluated in eight CHD patient datasets acquired with seven-fold prospective undersampling. Reconstruction quality was compared with the state-of-the-art non-rigid motion-corrected patch-based low-rank reconstruction method (NR-PROST) and against reference images (acquired with three-or-four-fold undersampling and reconstructed with NR-PROST). RESULTS: Seven-fold undersampled scan times were 2.1 ± 0.3 minutes and reconstruction times were ∼30 seconds, approximately 240 times faster than an NR-PROST reconstruction. Image quality comparable to the reference images was achieved using the proposed MoCo-MoDL framework, with no statistically significant differences found in any of the assessed quantitative or qualitative image quality measures. Additionally, expert image quality scores indicated the MoCo-MoDL reconstructions were consistently of a higher quality than the NR-PROST reconstructions of the same data, with the differences in 12 of the 22 scores measured for individual vascular structures found to be statistically significant. CONCLUSION: The MoCo-MoDL framework was applied to an adult CHD patient cohort, achieving good quality 3D whole-heart images from ∼2-minute scans with reconstruction times of ∼30 seconds.

Recurrent neural network-based simultaneous cardiac T1, T2, and T1ρ mapping.

The purpose of the current study was to explore the feasibility of training a deep neural network to accelerate the process of generating T1, T2, and T1ρ maps for a recently proposed free-breathing cardiac multiparametric mapping technique, where a recurrent neural network (RNN) was utilized to exploit the temporal correlation among the multicontrast images. The RNN-based model was developed for rapid and accurate T1, T2, and T1ρ estimation. Bloch simulation was performed to simulate a dataset of more than 10 million signals and time correspondences with different noise levels for network training. The proposed RNN-based method was compared with a dictionary-matching method and a conventional mapping method to evaluate the model's effectiveness in phantom and in vivo studies at 3 T, respectively. In phantom studies, the RNN-based method and the dictionary-matching method achieved similar accuracy and precision in T1, T2, and T1ρ estimations. In in vivo studies, the estimated T1, T2, and T1ρ values obtained by the two methods achieved similar accuracy and precision for 10 healthy volunteers (T1: 1228.70 ± 53.80 vs. 1228.34 ± 52.91 ms, p > 0.1; T2: 40.70 ± 2.89 vs. 41.19 ± 2.91 ms, p > 0.1; T1ρ: 45.09 ± 4.47 vs. 45.23 ± 4.65 ms, p > 0.1). The RNN-based method can generate cardiac multiparameter quantitative maps simultaneously in just 2 s, achieving 60-fold acceleration compared with the dictionary-matching method. The RNN-accelerated method offers an almost instantaneous approach for reconstructing accurate T1, T2, and T1ρ maps, being much more efficient than the dictionary-matching method for the free-breathing multiparametric cardiac mapping technique, which may pave the way for inline mapping in clinical applications.

Stop moving: MR motion correction as an opportunity for artificial intelligence.

Subject motion is a long-standing problem of magnetic resonance imaging (MRI), which can seriously deteriorate the image quality. Various prospective and retrospective methods have been proposed for MRI motion correction, among which deep learning approaches have achieved state-of-the-art motion correction performance. This survey paper aims to provide a comprehensive review of deep learning-based MRI motion correction methods. Neural networks used for motion artifacts reduction and motion estimation in the image domain or frequency domain are detailed. Furthermore, besides motion-corrected MRI reconstruction, how estimated motion is applied in other downstream tasks is briefly introduced, aiming to strengthen the interaction between different research areas. Finally, we identify current limitations and point out future directions of deep learning-based MRI motion correction.

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.

Accurate exclusion of kidney regions affected by susceptibility artifact in blood oxygenation level-dependent (BOLD) images using deep-learning-based segmentation.

Susceptibility artifact (SA) is common in renal blood oxygenation level-dependent (BOLD) images, and including the SA-affected region could induce much error in renal oxygenation quantification. In this paper, we propose to exclude kidney regions affected by SA in gradient echo images with different echo times (TE), based on a deep-learning segmentation approach. For kidney segmentation, a ResUNet was trained with 4000 CT images and then tuned with 60 BOLD images. Verified by a Monte Carlo simulation, the presence of SA leads to a bilinear pattern for the segmented area of kidney as function of TE, and the segmented kidney in the image of turning point's TE would exclude SA-affected regions. To evaluate the accuracy of excluding SA-affected regions, we compared the SA-free segmentations by the proposed method against manual segmentation by an experienced user for BOLD images of 35 subjects, and found DICE of 93.9% ± 3.4%. For 10 kidneys with severe SA, the DICE was 94.5% ± 1.7%, for 14 with moderate SA, 92.8% ± 4.7%, and for 46 with mild or no SA, 94.3% ± 3.8%. For the three sub-groups of kidneys, correction of SA led to a decrease of R2* of 8.5 ± 2.8, 4.7 ± 1.8, and 1.6 ± 0.9 s-1, respectively. In conclusion, the proposed method is capable of segmenting kidneys in BOLD images and at the same time excluding SA-affected region in a fully automatic way, therefore can potentially improve both speed and accuracy of the quantification procedure of renal BOLD data.

Identification of ACKR4 as an immune checkpoint in pulmonary arterial hypertension.

Objective: Inflammation is recognized as a contributor in the development of pulmonary arterial hypertension (PAH), and the recruitment and functional capacity of immune cells are well-orchestrated by chemokines and their receptors. This study is aimed at identification of critical chemokines in the progression of PAH via transcriptomic analysis. Methods: Differentially expressed genes (DEGs) from lungs of PAH patients were achieved compared to controls based on Gene Expression Omnibus (GEO) database. Gene set enrichment analysis (GSEA) was applied for functional annotation and pathway enrichement. The abundance of immune cells was estimated by the xCell algorithm. Weighted correlation network analysis (WGCNA) was used to construct a gene expression network, based on which a diagnostic model was generated to determine its accuracy to distinguish PAH from control subjects. Target genes were then validated in lung of hypoxia-induce pulmonary hypertension (PH) mouse model. Results: ACKR4 (atypical chemokine receptor 4) was downregulated in PAH lung tissues in multiple datasets. PAH relevant biological functions and pathways were enriched in patients with low-ACKR4 level according to GSEA enrichment analysis. Immuno-infiltration analysis revealed a negative correlation of activated dendritic cells, Th1 and macrophage infiltration with ACKR4 expression. Three gene modules were associated with PAH via WGCNA analysis, and a model for PAH diagnosis was generated using CXCL12, COL18A1 and TSHZ2, all of which correlated with ACKR4. The ACKR4 expression was also downregulated in lung tissues of our experimental PH mice compared to that of controls. Conclusions: The reduction of ACKR4 in lung tissues of human PAH based on transcriptomic data is consistent with the alteration observed in our rodent PH. The correlation with immune cell infiltration and functional annotation indicated that ACKR4 might serve as a protective immune checkpoint for PAH.

Fast Multi-Contrast MRI Acquisition by Optimal Sampling of Information Complementary to Pre-Acquired MRI Contrast.

Recent studies on multi-contrast MRI reconstruction have demonstrated the potential of further accelerating MRI acquisition by exploiting correlation between contrasts. Most of the state-of-the-art approaches have achieved improvement through the development of network architectures for fixed under-sampling patterns, without considering inter-contrast correlation in the under-sampling pattern design. On the other hand, sampling pattern learning methods have shown better reconstruction performance than those with fixed under-sampling patterns. However, most under-sampling pattern learning algorithms are designed for single contrast MRI without exploiting complementary information between contrasts. To this end, we propose a framework to optimize the under-sampling pattern of a target MRI contrast which complements the acquired fully-sampled reference contrast. Specifically, a novel image synthesis network is introduced to extract the redundant information contained in the reference contrast, which is exploited in the subsequent joint pattern optimization and reconstruction network. We have demonstrated superior performance of our learned under-sampling patterns on both public and in-house datasets, compared to the commonly used under-sampling patterns and state-of-the-art methods that jointly optimize the reconstruction network and the under-sampling patterns, up to 8-fold under-sampling factor.

Free-breathing simultaneous water-fat separation and T1 mapping of the whole liver (SWALI) with isotropic resolution using 3D golden-angle radial trajectory.

Background: Conventional liver T1 mapping techniques are typically performed under breath-holding conditions; they have limited slice coverage and often rely on multiple acquisitions. Furthermore, liver fat affects the accuracy of T1 quantification. Therefore, we aim to propose a free-breathing technique for simultaneous water-fat separation and T1 mapping of the whole liver (SWALI) in a single scan. Methods: The proposed SWALI sequence included an inversion recovery (IR) preparation pulse followed by a series of multiecho three-dimensional (3D) golden-angle radial acquisitions. For each echo time (TE), a series of images containing a mix of water and fat were reconstructed using a sliding window method. For each inversion time (TI), water and fat were separated, and then water and fat T1 estimation was conducted. The fat fraction (FF) was calculated based on the last TI image. The FF and water T1 quantification accuracy were compared with the gold standard sequences in the phantom. The in vivo feasibility was tested in 9 healthy volunteers, 2 patients with fatty liver, and 3 patients with hepatocellular carcinoma (HCC). The reproducibility was evaluated in the patients with fatty liver and in the healthy volunteers. Results: The mean FF and the mean water T1 values obtained by the SWALI sequence showed good agreements with chemical shift-encoded magnetic resonance imaging (CSE-MRI; r=0.998; P<0.001) and fat-suppressed (FS) IR-spin echo (SE; r=0.997; P<0.001) in the phantom. For the patients with fatty liver and the healthy volunteers, the SWALI sequence showed no significant difference with CSE-MRI in FF quantification (P=0.53). In T1 quantification, comparable T1 values were obtained with the SWALI sequence and modified Look-Locker inversion recovery (MOLLI; P=0.10) in healthy volunteers, while the water T1 estimated by the SWALI sequence was significantly lower than the water-fat compound T1 estimated by MOLLI (P<0.001) in patients with fatty liver. In the reproducibility study, the intraclass correlation coefficients (ICCs) for the estimated FF and water T1 were 0.997 and 0.943, respectively. Water T1 of the patients with HCC calculated using the SWALI sequence showed a significant reduction after the contrast administration (P<0.001). Conclusions: Free-breathing water-fat separation and T1 mapping of the whole liver with 2.5 mm isotropic spatial resolution were achieved simultaneously using the SWALI sequence in a 5-min scan.

Simultaneous T1, T2, and T2* Mapping of Carotid Plaque: The SIMPLE* Technique.

Background Quantitative T1, T2, and T2* measurements of carotid atherosclerotic plaque are important in evaluating plaque vulnerability and monitoring its progression. Purpose To develop a sequence to simultaneously quantify T1, T2, and T2* of carotid plaque. Materials and Methods The simultaneous T1, T2, and T2* mapping of carotid plaque (SIMPLE*) sequence is composed of three modules with different T2 preparation pulses, inversion-recovery pulses, and acquisition schemas. Single-echo data were used for T1 and T2 quantification, while the multiecho (ME) data were used for T2* quantification. The quantitative accuracy of SIMPLE* was tested in a phantom study by comparing its measurements with those of reference standard sequences. In vivo feasibility of the technique was prospectively evaluated between November 2020 and February 2022 in healthy volunteers and participants with carotid atherosclerotic plaque. The Pearson or Spearman correlation test, Student t test, and Wilcoxon rank-sum test were used. Results T1, T2, and T2* estimated with SIMPLE* strongly correlated with inversion-recovery spin-echo (SE) (correlation coefficient [r] = 0.99), ME-SE (r = 0.99), and ME gradient-echo (r = 0.99) sequences in the phantom study. In five healthy volunteers (mean age, 25 years ± 3 [SD]; three women), measurements were similar between SIMPLE* and modified Look-Locker inversion recovery, or MOLLI (1151 msec ± 71 vs 1098 msec ± 64; P = .14), ME turbo SE (31 msec ± 1 vs 31 msec ± 1; P = .32), and ME turbo field echo (24 msec ± 2 vs 25 msec ± 2; P = .19). In 18 participants with carotid plaque (mean age, 65 years ± 9; 16 men), quantitative T1, T2, and T2* of plaque components were consistent with their signal characteristics on multicontrast images. Conclusion A quantitative technique for simultaneous T1, T2, and T2* mapping of carotid plaque with 100-mm3 coverage and 0.8-mm3 resolution was developed using the proposed SIMPLE* sequence and demonstrated high accuracy and in vivo feasibility. © RSNA, 2023 Supplemental material is available for this article.

Free-running 3D whole-heart T1 and T2 mapping and cine MRI using low-rank reconstruction with non-rigid cardiac motion correction.

PURPOSE: To introduce non-rigid cardiac motion correction into a novel free-running framework for the simultaneous acquisition of 3D whole-heart myocardial T 1 $$ {T}_1 $$ and T 2 $$ {T}_2 $$ maps and cine images, enabling a ∼ $$ \sim $$ 3-min scan. METHODS: Data were acquired using a free-running 3D golden-angle radial readout interleaved with inversion recovery and T 2 $$ {T}_2 $$ -preparation pulses. After correction for translational respiratory motion, non-rigid cardiac-motion-corrected reconstruction with dictionary-based low-rank compression and patch-based regularization enabled 3D whole-heart T 1 $$ {T}_1 $$ and T 2 $$ {T}_2 $$ mapping at any given cardiac phase as well as whole-heart cardiac cine imaging. The framework was validated and compared with established methods in 11 healthy subjects. RESULTS: Good quality 3D T 1 $$ {T}_1 $$ and T 2 $$ {T}_2 $$ maps and cine images were reconstructed for all subjects. Septal T 1 $$ {T}_1 $$ values using the proposed approach ( 1200 ± 50 $$ 1200\pm 50 $$ ms) were higher than those from a 2D MOLLI sequence ( 1063 ± 33 $$ 1063\pm 33 $$ ms), which is known to underestimate T 1 $$ {T}_1 $$ , while T 2 $$ {T}_2 $$ values from the proposed approach ( 51 ± 4 $$ 51\pm 4 $$ ms) were in good agreement with those from a 2D GraSE sequence ( 51 ± 2 $$ 51\pm 2 $$ ms). CONCLUSION: The proposed technique provides 3D T 1 $$ {T}_1 $$ and T 2 $$ {T}_2 $$ maps and cine images with isotropic spatial resolution in a single ∼ $$ \sim $$ 3.3-min scan.

Motion detection and correction for carotid MRI using a markerless optical system.

PURPOSE: Motion related artifact is a challenge for MRI, especially when imaging regions like the carotid artery where complex motion (abrupt and bulk motion) may occur. This study aims to develop a non-contact motion detection and correction system for carotid MRI using a markerless optical tracking system. METHODS: The proposed markerless optical tracking system consisted of a cross-line laser, an MRI-compatible camera and plastic holders mounted inside the scanner bore. The neck motion of the subject can be captured by monitoring the change of the projected laser position in real-time. The system was used to correct both abrupt motion and bulk motion for carotid MRI. The abrupt motion (e.g. coughing) was compensated by discarding the corrupted k-space lines and re-estimating the missing lines using SPIRiT algorithm. The bulk motion was corrected by phase adjustment of k-space lines according to the measured 1D-translational bulk motion (along anterior-posterior direction) and optimized in-plane translation parameters. Ten volunteers underwent carotid MRI with real-time neck motion detection and retrospective motion correction. Artery sharpness, vessel wall thickness and overall image quality score were compared between the motion-corrupted image and motion-corrected images of different correction strategies. RESULTS: Both the abrupt motion and the bulk motion during carotid scanning were successfully detected and corrected. The results of ten volunteers demonstrated significant improvement in carotid artery sharpness, vessel wall thickness measurement, and overall image quality score using the proposed markerless optical tracking system and motion correction strategies. CONCLUSION: The proposed markerless structured light based motion detection and correction system can sensitively detect both abrupt and bulk motion during carotid MR scans. By correcting for both abrupt and bulk motion, vessel wall delineation was improved in carotid MR images, which could potentially facilitate carotid plaque identification and atherosclerosis diagnosis in the future.

Accelerated 3D free-breathing high-resolution myocardial T1ρ mapping at 3 Tesla.

PURPOSE: To develop a fast free-breathing whole-heart high-resolution myocardial T1ρ mapping technique with robust spin-lock preparation that can be performed at 3 Tesla. METHODS: An adiabatically excited continuous-wave spin-lock module, insensitive to field inhomogeneities, was implemented with an electrocardiogram-triggered low-flip angle spoiled gradient echo sequence with variable-density 3D Cartesian undersampling at a 3 Tesla whole-body scanner. A saturation pulse was performed at the beginning of each cardiac cycle to null the magnetization before T1ρ preparation. Multiple T1ρ -weighted images were acquired with T1ρ preparations with different spin-lock times in an interleaved fashion. Respiratory self-gating approach was adopted along with localized autofocus to enable 3D translational motion correction of the data acquired in each heartbeat. After motion correction, multi-contrast locally low-rank reconstruction was performed to reduce undersampling artifacts. The accuracy and feasibility of the 3D T1ρ mapping technique was investigated in phantoms and in vivo in 10 healthy subjects compared with the 2D T1ρ mapping. RESULTS: The 3D T1ρ mapping technique provided similar phantom T1ρ measurements in the range of 25-120 ms to the 2D T1ρ mapping reference over a wide range of simulated heart rates. With the robust adiabatically excited continuous-wave spin-lock preparation, good quality 2D and 3D in vivo T1ρ -weighted images and T1ρ maps were obtained. Myocardial T1ρ values with the 3D T1ρ mapping were slightly longer than 2D breath-hold measurements (septal T1ρ : 52.7 ± 1.4 ms vs. 50.2 ± 1.8 ms, P < 0.01). CONCLUSION: A fast 3D free-breathing whole-heart T1ρ mapping technique was proposed for T1ρ quantification at 3 T with isotropic spatial resolution (2 mm3 ) and short scan time of ∼4.5 min.

Quantification of changes in myocardial T1 * values with exercise cardiac MRI using a free-breathing non-electrocardiograph radial imaging.

PURPOSE: To develop and evaluate a free breathing non-electrocardiograph (ECG) myocardial T1 * mapping sequence using radial imaging to quantify the changes in myocardial T1 * between rest and exercise (T1 *reactivity ) in exercise cardiac MRI (Ex-CMR). METHODS: A free-running T1 * sequence was developed using a saturation pulse followed by three Look-Locker inversion-recovery experiments. Each Look-Locker continuously acquired data as radial trajectory using a low flip-angle spoiled gradient-echo readout. Self-navigation was performed with a temporal resolution of ∼100 ms for retrospectively extracting respiratory motion. The mid-diastole phase for every cardiac cycle was retrospectively detected on the recorded electrocardiogram signal using an empirical model. Multiple measurements were performed to obtain mean value to reduce effects from the free-breathing acquisition. Finally, data acquired at both mid-diastole and end-expiration are picked and reconstructed by a low-rank plus sparsity constraint algorithm. The performance of this sequence was evaluated by simulations, phantoms, and in vivo studies at rest and after physiological exercise. RESULTS: Numerical simulation demonstrated that changes in T1 * are related to the changes in T1 ; however, other factors such as breathing motion could influence T1 * measurements. Phantom T1 * values measured using free-running T1 * mapping sequence had good correlation with spin-echo T1 values and was insensitive to heart rate. In the Ex-CMR study, the measured T1 * reactivity was 10% immediately after exercise and declined over time. CONCLUSION: The free-running T1 * mapping sequence allows free-breathing non-ECG quantification of changes in myocardial T1 * with physiological exercise. Although, absolute myocardial T1 * value is sensitive to various confounders such as B1 and B0 inhomogeneity, quantification of its change may be useful in revealing myocardial tissue properties with exercise.

Accelerating 3D MTC-BOOST in patients with congenital heart disease using a joint multi-scale variational neural network reconstruction.

PURPOSE: Free-breathing Magnetization Transfer Contrast Bright blOOd phase SensiTive (MTC-BOOST) is a prototype balanced-Steady-State Free Precession sequence for 3D whole-heart imaging, that employs the endogenous magnetisation transfer contrast mechanism. This achieves reduction of flow and off-resonance artefacts, that often arise with the clinical T2prepared balanced-Steady-State Free Precession sequence, enabling high quality, contrast-agent free imaging of the thoracic cardiovascular anatomy. Fully-sampled MTC-BOOST acquisition requires long scan times (~10-24 min) and therefore acceleration is needed to permit its clinical incorporation. The aim of this study is to enable and clinically validate the 5-fold accelerated MTC-BOOST acquisition with joint Multi-Scale Variational Neural Network (jMS-VNN) reconstruction. METHODS: Thirty-six patients underwent free-breathing, 3D whole-heart imaging with the MTC-BOOST sequence, which is combined with variable density spiral-like Cartesian sampling and 2D image navigators for translational motion estimation. This sequence acquires two differently weighted bright-blood volumes in an interleaved fashion, which are then joined in a phase sensitive inversion recovery reconstruction to obtain a complementary fully co-registered black-blood volume. Data from eighteen patients were used for training, whereas data from the remaining eighteen patients were used for testing/evaluation. The proposed deep-learning based approach adopts a supervised multi-scale variational neural network for joint reconstruction of the two differently weighted bright-blood volumes acquired with the 5-fold accelerated MTC-BOOST. The two contrast images are stacked as different channels in the network to exploit the shared information. The proposed approach is compared to the fully-sampled MTC-BOOST and 5-fold undersampled MTC-BOOST acquisition with Compressed Sensing (CS) reconstruction in terms of scan/reconstruction time and bright-blood image quality. Comparison against conventional 2-fold undersampled T2-prepared 3D bright-blood whole-heart clinical sequence (T2prep-3DWH) is also included. RESULTS: Acquisition time was 3.0 ±â€¯1.0 min for the 5-fold accelerated MTC-BOOST versus 9.0 ±â€¯1.1 min for the fully-sampled MTC-BOOST and 11.1 ±â€¯2.6 min for the T2prep-3DWH (p < 0.001 and p < 0.001, respectively). Reconstruction time was significantly lower with the jMS-VNN method compared to CS (10 ±â€¯0.5 min vs 20 ±â€¯2 s, p < 0.001). Image quality was higher for the proposed 5-fold undersampled jMS-VNN versus conventional CS, comparable or higher to the corresponding T2prep-3DWH dataset and similar to the fully-sampled MTC-BOOST. CONCLUSION: The proposed 5-fold accelerated jMS-VNN MTC-BOOST framework provides efficient 3D whole-heart bright-blood imaging in fast acquisition and reconstruction time with concomitant reduction of flow and off-resonance artefacts, that are frequently encountered with the clinical sequence. Image quality of the cardiac anatomy and thoracic vasculature is comparable or superior to the clinical scan and 5-fold CS reconstruction in faster reconstruction time, promising potential clinical adoption.

End-to-End Deep Learning of Non-rigid Groupwise Registration and Reconstruction of Dynamic MRI.

Temporal correlation has been exploited for accelerated dynamic MRI reconstruction. Some methods have modeled inter-frame motion into the reconstruction process to produce temporally aligned image series and higher reconstruction quality. However, traditional motion-compensated approaches requiring iterative optimization of registration and reconstruction are time-consuming, while most deep learning-based methods neglect motion in the reconstruction process. We propose an unrolled deep learning framework with each iteration consisting of a groupwise diffeomorphic registration network (GRN) and a motion-augmented reconstruction network. Specifically, the whole dynamic sequence is registered at once to an implicit template which is used to generate a new set of dynamic images to efficiently exploit the full temporal information of the acquired data via the GRN. The generated dynamic sequence is then incorporated into the reconstruction network to augment the reconstruction performance. The registration and reconstruction networks are optimized in an end-to-end fashion for simultaneous motion estimation and reconstruction of dynamic images. The effectiveness of the proposed method is validated in highly accelerated cardiac cine MRI by comparing with other state-of-the-art approaches.