Normal diastolic and systolic myocardial T1 values at 1.5-T MR imaging: correlations and blood normalization

Estimating synthetic hematocrit and extracellular volume from native blood pool T1 times at 3 Tesla CMR: Derivation of a conversion equation, accuracy and comparison with published formulas

PURPOSE

Calculation of extracellular volume fraction (ECV), a marker of myocardial fibrosis in cardiac magnetic resonance imaging (CMR), requires hematocrit (Hct). We aimed to correlate Hct levels with native blood T1 times, to derive a formula for estimating synthetic Hct (Hctsyn) and synthetic ECV (ECVsyn), to assess accuracy of ECVsyn and to compare our model with published formulas.

METHOD

In this retrospective study, a cohort of 250 CMR scans with T1 mapping (3T, MOLLI 5(3)3, endsystolic aquisition), was divided into a derivation and validation cohort. Native T1 times of the left ventricular blood pool (T1native,midLV) were correlated with Hct levels from blood sampling within 24 h (Hct24h) and a formula for calculation of Hctsyn was derived by linear regression.

RESULTS

In the derivation cohort (n = 167), Hct24h showed a good association with T1native,midLV (r = -0.711, p < 0.001). The resulting regression equation was Hctsyn = 1/T1native,midLV * 1355.52-0.310. In the validation cohort (n = 83), Hctsyn and Hct24h showed good correlation (r = 0.726, p < 0.001), while ECVsyn, and ECV24h demonstrated excellent correlation (r = 0.940, p < 0.001). ECVsyn had a minimal bias of 0.28 % and the misclassification rate (8.8 %) was comparable to the variability introduced by repeated Hct measurements (misclassification in 7.5 %). Applying published formulas in our cohort resulted in incorrect classification in up to 60 %.

CONCLUSION

We provide a formula for estimating Hctsyn from native blood T1 on a 3T scanner. The measurement error of ECVsyn is low and comparable to the error due to retest variability of conventional Hct. Scanner- and sequence-specific formulas should be used.

Conditions for late gadolinium enhancement MRI in myocardial infarction model rats that better reflect microscopic tissue staining

Late gadolinium enhancement (LGE) is a widely used magnetic resonance imaging method for assessing cardiac disease. However, the relationship between different LGE signal thresholds and microscopic tissue staining images is unclear. In this study, we performed cardiovascular MRI on myocardial infarction (MI) model rats and evaluated the relationship between LGE with different signal thresholding methods and tissue staining images. We prepared 16 rats that underwent MRI 14-18 days following a surgery to create an MI model. We captured cine and LGE images of the cardiac short-axis and longitudinal two- and four-chamber views. The mean ± 2SD, ± 3SD, and ± 5SD of the pixel values in the non-infarcted area were defined as the LGE area. We compared areas of Sirius red staining, determined by the color tone, with their respective LGE areas at end-diastole and end-systole. We observed that the LGE area calculated as the mean ± 2SD of the non-infarcted area at end-diastole demonstrated a significant positive correlation with the area of Sirius red staining (Pearson's correlation coefficient in both: 0.81 [p < 0.01]). Therefore, the LGE area calculated as the mean ± 2SD of the non-infarcted area at end-diastole best reflected the MI area in tissue staining.

Sex-specific ventricular morphology, function, and tissue characteristics in arterial hypertension: a magnetic resonance study of the Hamburg city health cohort

OBJECTIVE

To determine the influence of arterial hypertension (AHT), sex, and the interaction between both left- and right ventricular (LV, RV) morphology, function, and tissue characteristics.

METHODS

The Hamburg City Health Study (HCHS) is a population-based, prospective, monocentric study. 1972 individuals without a history of cardiac diseases/ interventions underwent 3 T cardiac MR imaging (CMR). Generalized linear models were conducted, including AHT, sex (and the interaction if significant), age, body mass index, place of birth, diabetes mellitus, smoking, hyperlipoproteinemia, atrial fibrillation, and medication.

RESULTS

Of 1972 subjects, 68% suffered from AHT. 42% with AHT and 49% controls were female. Females overall showed a higher ejection fraction (EF) (LV: regression coefficient +2.4% [95% confidence interval: 1.7; 3.1]), lower volumes and LV mass (-19.8% [-21.3; -18.5]), and prolonged native septal T1 (+22.1 ms [18.3; 25.9])/T2 relaxation times (+1.1 ms [0.9; 1.3]) (all p < 0.001) compared to males. Subjects with AHT showed a higher EF (LV: +1.2% [0.3; 2.0], p = 0.009) and LV mass (+6.6% [4.3; 9.0], p < 0.001) than controls. The interaction between sex and AHT influenced mapping. After excluding segments with LGE, males (-0.7 ms [-1.0; -0.3 | ) and females with AHT (-1.1 ms [-1.6; -0.6]) showed shorter T2 relaxation times than the sex-respective controls (p < 0.001), but the effect was stronger in females.

CONCLUSION

In the HCHS, female and male subjects with AHT likewise showed a higher EF and LV mass than controls, independent of sex. However, differences in tissue characteristics between subjects with AHT and controls appeared to be sex-specific.

CLINICAL RELEVANCE STATEMENT

The interaction between sex and cardiac risk factors is an underestimated factor that should be considered when comparing tissue characteristics between hypertensive subjects and controls, and when establishing cut-off values for normal and pathological relaxation times.

KEY POINTS

There are sex-dependent differences in arterial hypertension, but it is unclear if cardiac MR parameters are sex-specific. Differences in cardiac MR parameters between hypertensive subjects and healthy controls appeared to be sex-specific for tissue characteristics. Sex needs to be considered when comparing tissue characteristics in patients with arterial hypertension to healthy controls.

Deep image prior cine MR fingerprinting with B1 + spin history correction

PURPOSE

To develop a deep image prior (DIP) reconstruction for B1 + -corrected 2D cine MR fingerprinting (MRF).

METHODS

The proposed method combines low-rank (LR) modeling with a DIP to generate cardiac phase-resolved parameter maps without motion correction, employing self-supervised training to enforce consistency with undersampled spiral k-space data. Two implementations were tested: one approach (DIP) for cine T1 , T2 , and M0 mapping, and a second approach (DIP with effective B1 + estimation [DIP-B1]) that also generated an effective B1 + map to correct for errors due to RF transmit inhomogeneities, through-plane motion, and blood flow. Cine MRF data were acquired in 14 healthy subjects and four reconstructions were compared: LR, low-rank motion-corrected (LRMC), DIP, and DIP-B1. Results were compared to diastolic ECG-triggered MRF, MOLLI, and T2 -prep bSSFP. Additionally, bright-blood and dark-blood images calculated from cine MRF maps were used to quantify ventricular function and compared to reference cine measurements.

RESULTS

DIP and DIP-B1 outperformed other cine MRF reconstructions with improved noise suppression and delineation of high-resolution details. Within-segment variability in the myocardium (reported as the coefficient of variation for T1 /T2 ) was lowest for DIP-B1 (2.3/8.3%) followed by DIP (2.7/8.7%), LRMC (3.5/10.5%), and LR (15.3/39.6%). Spatial homogeneity improved with DIP-B1 having the lowest intersegment variability (2.6/4.1%). The mean bias in ejection fraction was -1.1% compared to reference cine scans.

CONCLUSION

A DIP reconstruction for 2D cine MRF enabled cardiac phase-resolved mapping of T1 , T2 , M0 , and the effective B1 + with improved noise suppression and precision compared to LR and LRMC.

STEAM-SASHA: a novel approach for blood- and fat-suppressed native T1 measurement in the right ventricular myocardium

OBJECTIVE

The excellent blood and fat suppression of stimulated echo acquisition mode (STEAM) can be combined with saturation recovery single-shot acquisition (SASHA) in a novel STEAM-SASHA sequence for right ventricular (RV) native T1 mapping.

MATERIALS AND METHODS

STEAM-SASHA splits magnetization preparation over two cardiac cycles, nulling blood signal and allowing fat signal to decay. Breath-hold T1 mapping was performed in a T1 phantom and twice in 10 volunteers using STEAM-SASHA and a modified Look-Locker sequence at peak systole at 3T. T1 was measured in 3 RV regions, the septum and left ventricle (LV).

RESULTS

In phantoms, MOLLI under-estimated while STEAM-SASHA over-estimated T1, on average by 3.0% and 7.0% respectively, although at typical 3T myocardial T1 (T1 > 1200 ms) STEAM-SASHA was more accurate. In volunteers, T1 was higher using STEAM-SASHA than MOLLI in the LV and septum (p = 0.03, p = 0.006, respectively), but lower in RV regions (p > 0.05). Inter-study, inter-observer and intra-observer coefficients of variation in all regions were < 15%. Blood suppression was excellent with STEAM-SASHA and noise floor effects were minimal.

DISCUSSION

STEAM-SASHA provides accurate and reproducible T1 in the RV with excellent blood and fat suppression. STEAM-SASHA has potential to provide new insights into pathological changes in the RV in future studies.

Misclassification of females and males in cardiovascular magnetic resonance parametric mapping: the importance of sex-specific normal ranges for diagnosis of health vs. disease

AIMS

Cardiovascular magnetic resonance parametric mapping enables non-invasive quantitative myocardial tissue characterization. Human myocardium has normal ranges of T1 and T2 values, deviation from which may indicate disease or change in physiology. Normal myocardial T1 and T2 values are affected by biological sex. Consequently, normal ranges created with insufficient numbers of each sex may result in sampling biases, misclassification of healthy values vs. disease, and even misdiagnoses. In this study, we investigated the impact of using male normal ranges for classifying female cases as normal or abnormal (and vice versa).

METHODS AND RESULTS

One hundred and forty-two healthy volunteers (male and female) were scanned on two Siemens 3T MR systems, providing averaged global myocardial T1 and T2 values on a per-subject basis. The Monte Carlo method was used to generate simulated normal ranges from these values to estimate the statistical accuracy of classifying healthy female or male cases correctly as 'normal' when using sex-specific vs. mixed-sex normal ranges. The normal male and female T1- and T2-mapping values were significantly different by sex, after adjusting for age and heart rate.

CONCLUSION

Using 15 healthy volunteers who are not sex specific to establish a normal range resulted in a typical misclassification of up to 36% of healthy females and 37% of healthy males as having abnormal T1 values and up to 16% of healthy females and 12% of healthy males as having abnormal T2 values. This paper highlights the potential adverse impact on diagnostic accuracy that can occur when local normal ranges contain insufficient numbers of both sexes. Sex-specific reference ranges should thus be routinely adopted in clinical practice.

Increased cardiac involvement in Fabry disease using blood-corrected native T1 mapping

Fabry disease (FD) is a rare lysosomal storage disorder resulting in myocardial sphingolipid accumulation which is detectable by cardiovascular magnetic resonance as low native T1. However, myocardial T1 contains signal from intramyocardial blood which affects variability and consequently measurement precision and accuracy. Correction of myocardial T1 by blood T1 increases precision. We therefore deployed a multicenter study of FD patients (n = 218) and healthy controls (n = 117) to investigate if blood-correction of myocardial native T1 increases the number of FD patients with low T1, and thus reclassifies FD patients as having cardiac involvement. Cardiac involvement was defined as a native T1 value 2 standard deviations below site-specific means in healthy controls for both corrected and uncorrected measures. Overall low T1 was 135/218 (62%) uncorrected vs. 145/218 (67%) corrected (p = 0.02). With blood-correction, 13/83 previously normal patients were reclassified to low T1. This reclassification appears clinically relevant as 6/13 (46%) of reclassified had focal late gadolinium enhancement or left ventricular hypertrophy as signs of cardiac involvement. Blood-correction of myocardial native T1 increases the proportion of FD subjects with low myocardial T1, with blood-corrected results tracking other markers of cardiac involvement. Blood-correction may potentially offer earlier detection and therapy initiation, but merits further prospective studies.

Reference values of myocardial native T1 and T2 mapping values in normal Indian population at 1.5 Tesla scanner

T1 and T2 mapping techniques on cardiovascular magnetic resonance (CMR) provide insights into the myocardial tissue characterisation. We sought to establish the normal reference values of native T1 and T2 mapping in Indian population which can be used subsequently in clinical practice for addressing various cardiac pathologies. This prospective study included consecutive healthy volunteers (18-60 years) who underwent CMR on a 1.5 Tesla scanner using standard protocol. T1 mapping sequence was performed using MOLLI sequence with two different flip angles (FA) (35° and 50°). T2 mapping was performed using a hybrid gradient and spin-echo sequence sequence with two different FA (70° and 12°). Images were analysed with ROIs drawn in all the 16 AHA myocardial segments. 50 volunteers (average age-34 years, males-72%) were included. All the scans were normal. The mean T1 value at 35-degree FA was 946.86 + 14.16 ms and at 50-degree FA was 941.60 + 11.89 ms. The mean T2 mapping value at 70-degree FA was 45.67 + 1.39 ms and at 12-degree FA was 45.61 + 1.47 ms. The mapping values were not statistically different between males and females (all p > 0.2). The T1 and T2 mapping values did not show any significant correlation with LVEF, age, BMI or heart rate (all r < 0.33). The T1 mapping values significantly differ at 35- and 50-degree FAs (p = 0.002). The results establish the normal reference T1 and T2 mapping value for Indian population in institutes using the same protocol and parameters at 1.5 Tesla and may guide future research.

Cardiac phase-resolved late gadolinium enhancement imaging

Late gadolinium enhancement (LGE) with cardiac magnetic resonance (CMR) imaging is the clinical reference for assessment of myocardial scar and focal fibrosis. However, current LGE techniques are confined to imaging of a single cardiac phase, which hampers assessment of scar motility and does not allow cross-comparison between multiple phases. In this work, we investigate a three step approach to obtain cardiac phase-resolved LGE images: (1) Acquisition of cardiac phase-resolved imaging data with varying T₁ weighting. (2) Generation of semi-quantitative T1* maps for each cardiac phase. (3) Synthetization of LGE contrast to obtain functional LGE images. The proposed method is evaluated in phantom imaging, six healthy subjects at 3T and 20 patients at 1.5T. Phantom imaging at 3T demonstrates consistent contrast throughout the cardiac cycle with a coefficient of variation of 2.55 ± 0.42%. In-vivo results show reliable LGE contrast with thorough suppression of the myocardial tissue is healthy subjects. The contrast between blood and myocardium showed moderate variation throughout the cardiac cycle in healthy subjects (coefficient of variation 18.2 ± 3.51%). Images were acquired at 40–60 ms and 80 ms temporal resolution, at 3T and 1.5, respectively. Functional LGE images acquired in patients with myocardial scar visualized scar tissue throughout the cardiac cycle, albeit at noticeably lower imaging resolution and noise resilience than the reference technique. The proposed technique bears the promise of integrating the advantages of phase-resolved CMR with LGE imaging, but further improvements in the acquisition quality are warranted for clinical use.

Saturation-pulse prepared heart-rate independent inversion-recovery (SAPPHIRE) biventricular T1 mapping: inter-field strength, head-to-head comparison of diastolic, systolic and dark-blood measurements

Background

To assess the feasibility of biventricular SAPPHIRE T₁ mapping in vivo across field strengths using diastolic, systolic and dark-blood (DB) approaches.

Methods

10 healthy volunteers underwent same-day non-contrast cardiovascular magnetic resonance at 1.5 Tesla (T) and 3 T. Left and right ventricular (LV, RV) T₁ mapping was performed in the basal, mid and apical short axis using 4-variants of SAPPHIRE: diastolic, systolic, 0th and 2nd order motion-sensitized DB and conventional modified Look-Locker inversion recovery (MOLLI).

Results

LV global myocardial T₁ times (1.5 T then 3 T results) were significantly longer by diastolic SAPPHIRE (1283 ± 11|1600 ± 17 ms) than any of the other SAPPHIRE variants: systolic (1239 ± 9|1595 ± 13 ms), 0th order DB (1241 ± 10|1596 ± 12) and 2nd order DB (1251 ± 11|1560 ± 20 ms, all p < 0.05). In the mid septum MOLLI and diastolic SAPPHIRE exhibited significant T₁ signal contamination (longer T₁) at the blood-myocardial interface not seen with the other 3 SAPPHIRE variants (all p < 0.025). Additionally, systolic, 0th order and 2nd order DB SAPPHIRE showed narrower dispersion of myocardial T₁ times across the mid septum when compared to diastolic SAPPHIRE (interquartile ranges respectively: 25 ms, 71 ms, 73 ms vs 143 ms, all p < 0.05). RV T₁ mapping was achievable using systolic, 0th and 2nd order DB SAPPHIRE but not with MOLLI or diastolic SAPPHIRE. All 4 SAPPHIRE variants showed excellent re-read reproducibility (intraclass correlation coefficients 0.953 to 0.996).

Conclusion

These small-scale preliminary healthy volunteer data suggest that DB SAPPHIRE has the potential to reduce partial volume effects at the blood-myocardial interface, and that systolic SAPPHIRE could be a feasible solution for right ventricular T₁ mapping. Further work is needed to understand the robustness of these sequences and their potential clinical utility.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12880-022-00843-0.

Phantom-based correction for standardization of myocardial native T1 and extracellular volume fraction in healthy subjects at 3-Tesla cardiac magnetic resonance imaging

OBJECTIVES

To investigate the effect of the phantom-based correction method for standardizing myocardial native T1 and extracellular volume fraction (ECV) in healthy subjects.

METHODS

Seventy-one healthy asymptomatic adult (≥ 20 years) volunteers of five different age groups (34 men and 37 women, 45.5 ± 15.5 years) were prospectively enrolled in three academic hospitals. Cardiac MRI including Modified Look - Locker Inversion recovery T1 mapping sequence was performed using a 3-Tesla system with a different type of scanner for each hospital. Native T1 and ECV were measured in the short-axis T1 map and analyzed for mean values of the 16 entire segments. The myocardial T1 value of each subject was corrected based on the site-specific equation derived from the T1 Mapping and ECV Standardization phantom. The global native T1 and ECV were compared between institutions before and after phantom-based correction, and the variation in native T1 and ECV among institutions was assessed using a coefficient of variation (CoV).

RESULTS

The global native T1 value significantly differed between the institutions (1198.7 ± 32.1 ms, institution A; 1217.7 ± 39.9 ms, institution B; 1232.7 ± 31.1 ms, institution C; p = 0.002), but the mean ECV did not (26.6-27.5%, p = 0.355). After phantom-based correction, the global native T1 and ECV were 1289.7 ± 32.4 ms and 25.0 ± 2.7%, respectively, and CoV for native T1 between the three institutions decreased from 3.0 to 2.5%. The corrected native T1 value did not significantly differ between institutions (1284.5 ± 31.5 ms, institution A; 1296.5 ± 39.1 ms, institution B; 1291.3 ± 29.3 ms, institution C; p = 0.440), and neither did the ECV (24.4-25.9%, p = 0.078).

CONCLUSIONS

The phantom-based correction method can provide standardized reference T1 values in healthy subjects.

KEY POINTS

• After phantom-based correction, the global native T1 of 16 entire myocardial segments on 3-T cardiac MRI is 1289.4 ± 32.4 ms, and the extracellular volume fraction was 25.0 ± 2.7% for healthy subjects. • After phantom - based correction was applied, the differences in the global native T1 among institutions became insignificant, and the CoV also decreased from 3.0 to 2.5%.

Metabolite-cycled echo-planar spectroscopic imaging of the human heart

Purpose

Spectroscopic imaging could provide insights into regional cardiac triglyceride variations, but is hampered by relatively long scan times. It is proposed to synergistically combine echo‐planar spectroscopic imaging (EPSI) with motion‐adapted gating, weighted acquisition and metabolite cycling to reduce scan times to less than 10 min while preserving spatial‐spectral quality. The method is compared to single‐voxel measurements and to metabolite‐cycled EPSI with conventional acquisition for assessing triglyceride‐to‐water (TG/W) ratios in the human heart.

Methods

Measurements were performed on 10 healthy volunteers using a clinical 1.5T system. EPSI data was acquired both without and with motion‐adapted gating in combination with weighted acquisition to assess TG/W ratios and relative Cramér‐Rao lower bounds (CRLB) of TG. For comparison, single‐voxel (PRESS) spectra were acquired in the interventricular septum.

Results

Bland–Altman analyses did not show a significant bias in TG/W when comparing both metabolite‐cycled EPSI methods to PRESS for any of the cardiac segments. Scan time was 8.05 ± 2.06 min and 17.91 ± 3.93 min for metabolite‐cycled EPSI with and without motion‐adapted gating and weighted acquisition, respectively, while relative CRLB of TG did not differ significantly between the two methods for any of the cardiac segments.

Conclusions

Metabolite‐cycled EPSI with motion‐adapted gating and weighted acquisition allows detecting TG/W ratios in different regions of the in vivo human heart. Scan time is reduced by more than 2‐fold to less than 10 min as compared to conventional acquisition, while keeping the quality of TG fitting constant.

Click here for author‐reader discussions

T2* assessment of the three coronary artery territories of the left ventricular wall by different monoexponential truncation methods

Objectives

This study aimed at evaluating left ventricular myocardial pixel-wise T2* using two truncation methods for different iron deposition T2* ranges and comparison of segmental T2* in different coronary artery territories.

Material and methods

Bright blood multi-gradient echo data of 30 patients were quantified by pixel-wise monoexponential T2* fitting with its R2 and SNR truncation. T2* was analyzed at different iron classifications. At low iron classification, T2* values were also analyzed by coronary artery territories.

Results

The right coronary artery has a significantly higher T2* value than the other coronary artery territories. No significant difference was found in classifying severe iron by the two truncation methods in any myocardial region, whereas in moderate iron, it is only apparent at septal segments. The R2 truncation produces a significantly higher T2* value than the SNR method when low iron is indicated.

Conclusion

Clear T2* differentiation between the three coronary territories by the two truncation methods is demonstrated. The two truncation methods can be used interchangeably in classifying severe and moderate iron deposition at the recommended septal region. However, in patients with low iron indication, different results by the two truncation methods can mislead the investigation of early iron level progression.

Myocardial T1 Values at 1.5 T: Normal Values for General Electric Scanners and Sex-Related Differences

BACKGROUND

No data are available about normal ranges for native T1 in human myocardium using General Electric (GE) scanners.

PURPOSE

To establish normal ranges for myocardial T1 values and evaluate regional variability and the influence of physiological factors.

STUDY TYPE

Prospective.

SUBJECTS

One hundred healthy volunteers with normal electrocardiogram, no cardiovascular/systemic diseases, or risk factors (age range: 20-70 years; 50 females).

FIELD STRENGTH/SEQUENCE

1.5 T/Steady-state free precession cine and a modified Look-Locker inversion recovery sequence in diastole (also in systole for 61 volunteers).

ASSESSMENT

Image analysis was performed by operators with >10 years experience in cardiac MR using commercially available software. T1 values were calculated for 16 myocardial segments, and the global value was the mean. Segments were grouped according to circumferential region (anterior, septal, inferior, and lateral) and to level (basal, medial, apical). Twenty images were analyzed twice by the same operator and by a different operator to assess reproducibility.

STATISTICAL TESTS

Independent-samples t-test or Mann-Whitney test; paired sample t-test or Wilcoxon signed-rank test; one-way repeated measures ANOVA or Friedman tests; Pearson's or Spearman's correlation. Reproducibility evaluated using coefficient of variability (CoV).

RESULTS

Due to artifacts and/or partial-volume effects, 45/1600 (2.8%) segments were excluded. A good intra- and inter-operator reproducibility was detected (CoV < 5%). There were significant differences in segmental T1 values (P < 0.05). A significant circumferential variability was present (P < 0.05): the mean native T1 value over the lateral region was significantly lower than in the other three regions. An increasing gradient from basal to apical slices was detected (P < 0.05). Segmental and global T1 values were not associated with age (range P = 0.052-0.911) but were significantly lower in males than in females (global: 993 ± 32 vs. 1037 ± 27 ms; P < 0.05) and significantly correlated with heart rate (range R for segmental values = 0.247-0.920; P < 0.05). Almost all segmental T1 values were inversely correlated with wall thickness (R from -0.233 to -0.514; P < 0.05). Systolic T1 values were significantly lower than diastolic values in basal anteroseptal segment, in all medial segments except the inferior one, and in all apical segments (P < 0.05).

DATA CONCLUSION

Myocardial T1 values differ among myocardial regions, are influenced by sex, heart rate, and wall thickness and vary according to the cardiac cycle in healthy adults.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

Impact of the Choice of Native T1 in Pixelwise Myocardial Blood Flow Quantification

BACKGROUND

Quantification of myocardial blood flow (MBF) from dynamic contrast-enhanced (DCE) MRI can be performed using a signal intensity model that incorporates T1 values of blood and myocardium.

PURPOSE

To assess the impact of T1 values on pixelwise MBF quantification, specifically to evaluate the influence of 1) study population-averaged vs. subject-specific, 2) diastolic vs. systolic, and 3) regional vs. global myocardial T1 values.

STUDY TYPE

Prospective.

SUBJECTS

Fifteen patients with chronic coronary heart disease.

FIELD STRENGTH/SEQUENCE

3T; modified Look-Locker inversion recovery for T1 mapping and saturation recovery gradient echo for DCE imaging, both acquired in a mid-ventricular short-axis slice in systole and diastole.

ASSESSMENT

MBF was estimated using Fermi modeling and signal intensity nonlinearity correction with different T1 values: study population-averaged blood and myocardial, subject-specific systolic and diastolic, and segmental T1 values. Myocardial segments with perfusion deficits were identified visually from DCE series.

STATISTICAL TESTS

The relationships between MBF parameters derived by different methods were analyzed by Bland-Altman analysis; corresponding mean values were compared by t-test.

RESULTS

Using subject-specific diastolic T1 values, global diastolic MBF was 0.61 ± 0.13 mL/(min·g). It did not differ from global MBF derived from the study population-averaged T1 (P = 0.88), but the standard deviation of differences was large (0.07 mL/(min·g), 11% of mean MBF). Global diastolic and systolic MBF did not differ (P = 0.12), whereas global diastolic MBF using systolic (0.62 ± 0.13 mL/(min·g)) and diastolic T1 values differed (P < 0.05). If regional instead of global T1 values were used, segmental MBF was lower in segments with perfusion deficits (bias = -0.03 mL/(min·g), -7% of mean MBF, P < 0.05) but higher in segments without perfusion deficits (bias = 0.01 mL/(min·g), 1% of mean MBF, P < 0.05).

DATA CONCLUSION

Whereas cardiac phase-specific T1 values have a minor impact on MBF estimates, subject-specific and myocardial segment-specific T1 values substantially affect MBF quantification.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY STAGE: 3.

Reference ranges ("normal values") for cardiovascular magnetic resonance (CMR) in adults and children: 2020 update

Cardiovascular magnetic resonance (CMR) enables assessment and quantification of morphological and functional parameters of the heart, including chamber size and function, diameters of the aorta and pulmonary arteries, flow and myocardial relaxation times. Knowledge of reference ranges ("normal values") for quantitative CMR is crucial to interpretation of results and to distinguish normal from disease. Compared to the previous version of this review published in 2015, we present updated and expanded reference values for morphological and functional CMR parameters of the cardiovascular system based on the peer-reviewed literature and current CMR techniques. Further, databases and references for deep learning methods are included.

Systolic modified Look-Locker inversion recovery myocardial T1 mapping improves the accuracy of T1 and extracellular volume fraction measurements of patients with high heart rate or atrial fibrillation

Image data for T1 mapping are generally acquired during mid-diastole period. However, T1 mapping tends to fail for patients with high heart rate or atrial fibrillation because of short or irregular R-R interval. Focusing on the evidence that the timing of systole is more stable than that of diastole from the R wave, we compared systolic T1 mapping with conventional diastolic T1 mapping for all participants (n = 58) by visual scoring of T1 calculation error. The systolic scores were significantly better than the diastolic scores (p < 0.05). This advantage of the systolic scores was confirmed selectively for patients with atrial fibrillation (p < 0.05, n = 19). The successful number of nonrigid image registration alignment for extracellular volume fraction (ECV) analysis also increased significantly for systolic images compared with diastolic images (p < 0.05). Thus, systolic T1 mapping improves the accuracy of T1 values and ECV analysis.

Automated quantification of myocardial tissue characteristics from native T1 mapping using neural networks with uncertainty-based quality-control

BACKGROUND

Tissue characterisation with cardiovascular magnetic resonance (CMR) parametric mapping has the potential to detect and quantify both focal and diffuse alterations in myocardial structure not assessable by late gadolinium enhancement. Native T1 mapping in particular has shown promise as a useful biomarker to support diagnostic, therapeutic and prognostic decision-making in ischaemic and non-ischaemic cardiomyopathies.

METHODS

Convolutional neural networks (CNNs) with Bayesian inference are a category of artificial neural networks which model the uncertainty of the network output. This study presents an automated framework for tissue characterisation from native shortened modified Look-Locker inversion recovery ShMOLLI T1 mapping at 1.5 T using a Probabilistic Hierarchical Segmentation (PHiSeg) network (PHCUMIS 119-127, 2019). In addition, we use the uncertainty information provided by the PHiSeg network in a novel automated quality control (QC) step to identify uncertain T1 values. The PHiSeg network and QC were validated against manual analysis on a cohort of the UK Biobank containing healthy subjects and chronic cardiomyopathy patients (N=100 for the PHiSeg network and N=700 for the QC). We used the proposed method to obtain reference T1 ranges for the left ventricular (LV) myocardium in healthy subjects as well as common clinical cardiac conditions.

RESULTS

T1 values computed from automatic and manual segmentations were highly correlated (r=0.97). Bland-Altman analysis showed good agreement between the automated and manual measurements. The average Dice metric was 0.84 for the LV myocardium. The sensitivity of detection of erroneous outputs was 91%. Finally, T1 values were automatically derived from 11,882 CMR exams from the UK Biobank. For the healthy cohort, the mean (SD) corrected T1 values were 926.61 (45.26), 934.39 (43.25) and 927.56 (50.36) for global, interventricular septum and free-wall respectively.

CONCLUSIONS

The proposed pipeline allows for automatic analysis of myocardial native T1 mapping and includes a QC process to detect potentially erroneous results. T1 reference values were presented for healthy subjects and common clinical cardiac conditions from the largest cohort to date using T1-mapping images.

Free breathing lung T1 mapping using image registration in patients with idiopathic pulmonary fibrosis

PURPOSE

To assess the use of image registration for correcting respiratory motion in free breathing lung T₁ mapping acquisition in patients with idiopathic pulmonary fibrosis (IPF).

THEORY AND METHODS

The method presented used image registration to synthetic images during postprocessing to remove respiratory motion. Synthetic images were generated from a model of the inversion recovery signal of the acquired images that incorporated a periodic lung motion model. Ten healthy volunteers and 19 patients with IPF underwent 2D Look-Locker T₁ mapping acquisition at 1.5T during inspiratory breath-hold and free breathing. Eight healthy volunteers and seven patients with IPF underwent T₁ mapping acquisition during expiratory breath-hold. Fourteen patients had follow-up scanning at 6 months. Dice similarity coefficient (DSC) was used to evaluate registration efficacy.

RESULTS

Image registration increased image DSC (P < .001) in the free breathing inversion recovery images. Lung T₁ measured during a free breathing acquisition was lower in patients with IPF when compared with healthy controls (inspiration: P = .238; expiration: P = .261; free breathing: P = .021). Measured lung T₁ was higher in expiration breath-hold than inspiration breath-hold in healthy volunteers (P < .001) but not in patients with IPF (P = .645). There were no other significant differences between lung T₁ values within subject groups.

CONCLUSIONS

The registration technique significantly reduced motion in the Look-Locker images acquired during free breathing and may improve the robustness of lung T₁ mapping in patients who struggle to hold their breath. Lung T₁ measured during a free breathing acquisition was significantly lower in patients with IPF when compared with healthy controls.

Navigator-free metabolite-cycled proton spectroscopy of the heart

PURPOSE

Respiratory gating in cardiac water-suppressed (WS) proton spectroscopy leads to long and unpredictable scan times. Metabolite cycling allows to perform frequency and phase correction on the water signal and, hence, offers an approach to navigator-free cardiac spectroscopy with fixed scan time. The objective of the present study was to develop and implement navigator-free metabolite-cycled cardiac proton spectroscopy (MC nonav) and compare it with standard navigator-gated WS (WS nav) and navigator-free WS (WS nonav) measurements for the assessment of triglyceride-to-water ratios (TG/W) and creatine-to-water ratios (CR/W) in the intraventricular septum of the in vivo heart.

METHODS

Navigator-free metabolite-cycled spectroscopy was implemented on a clinical 1.5T system. In vivo measurements were performed on 10 young and 5 older healthy volunteers to assess signal-to-noise ratio efficiency as well as TG/W and CR/W and the relative Cramér-Rao lower bounds for CR. The performance of the metabolite-cycled sequence was verified using simulations.

RESULTS

On average, scan times of MC nonav were 3.4 times shorter compared with WS nav, while no significant bias for TG/W was observed (coefficient of variation = 14.0%). signal-to-noise ratio efficiency of both TG and CR increased for MC nonav compared with WS nav. Relative Cramér-Rao lower bounds of CR decreased for MC nonav. Overall spectral quality was found comparable between MC nonav and WS nav, while it was inferior for WS nonav.

CONCLUSION

Navigator-free metabolite-cycled cardiac proton spectroscopy offers 3.4-fold accelerated assessment of TG/W and CR/W in the heart with preserved spectral quality when compared with navigator-gated WS scans.

Early detection of heart function abnormality by native T1: a comparison of two T1 quantification methods

Objective

To compare the robustness of native T1 mapping using mean and median pixel-wise quantification methods.

Methods

Fifty-seven consecutive patients without overt signs of heart failure were examined in clinical routine for suspicion of cardiomyopathy. MRI included the acquisition of native T1 maps by a motion-corrected modified Look-Locker inversion recovery sequence at 1.5 T. Heart function status according to four established volumetric left ventricular (LV) cardio MRI parameter thresholds was used for retrospective separation into subgroups of normal (n = 26) or abnormal heart function (n = 31). Statistical normality of pixel-wise T1 was tested on each myocardial segment and mean and median segmental T1 values were assessed.

Results

Segments with normally distributed pixel-wise T1 (57/58%) showed no difference between mean and median quantification in either patient group, while differences were highly significant (p < 0.001) for the respective 43/42% non-normally distributed segments. Heart function differentiation between two patient groups was significant in 14 myocardial segments (p < 0.001–0.040) by median quantification compared with six (p < 0.001–0.042) by using the mean. The differences by median quantification were observed between the native T1 values of the three coronary artery territories of normal heart function patients (p = 0.023) and insignificantly in the abnormal patients (p = 0.053).

Conclusion

Median quantification increases the robustness of myocardial native T1 definition, regardless of statistical normality of the data. Compared with the currently prevailing method of mean quantification, differentiation between LV segments and coronary artery territories is better and allows for earlier detection of heart function impairment.

Key Points

• Median pixel-wise quantification of native T1 maps is robust and can be applied regardless of the statistical distribution of data points.

• Median quantification is more sensitive to early heart function abnormality compared with mean quantification.

• The new method yields significant native T1 value differentiation between the three coronary artery territories.

Electronic supplementary material

The online version of this article (10.1007/s00330-019-06364-9) contains supplementary material, which is available to authorized users.

Functional LGE Imaging: Cardiac Phase-Resolved Assessment of Focal Fibrosis

Cardiac Magnetic Resonance Imaging (CMR) is a central tool for diagnosis of various ischemic and non-ischemic cardiomyopathies. CMR protocols commonly comprise assessment of functional properties using cardiac phase-resolved CINE MRI and characterization of myocardial viability using late gadolinium enhancement (LGE) imaging. Conventional LGE imaging requires inversion recovery preparation with a specific inversion time to null the healthy myocardium, which restricts the acquisition to a single cardiac phase. In turn, this necessitates separate scans for cardiac function and viability. In this work, we develop a new method for functional LGE imaging in a single breath-hold using a three-step approach: 1) ECG-triggered multi-contrast data is acquired for each cardiac phase, 2) semi-quantitative relaxation maps are generated, 3) LGE imaging contrast is synthesized based on the semi-quantitative maps. The proposed functional LGE method is evaluated in four healthy subject and 20 patients at 1.5T and 3T. Thorough suppression of the healthy myocardium, as well as 40-80ms temporal resolution are achieved, with no visually apparent temporal blurring at tissue interfaces. Functional LGE in patients with focal scar demonstrates robust hyperenhancement in the scar area throughout all cardiac phases, allowing for visual assessment of scar motility. The proposed technique bears the potential to simplify and speedup common cardiac imaging protocols, while enabling improved data fusion of functional and viability information for improved evaluation of CMR.

Retrospective phase-based gating for cardiac proton spectroscopy with fixed scan time

BACKGROUND

Respiratory motion is a major limiting factor for spectral quality and duration of in vivo proton MR spectroscopy of the heart. Prospective navigator gating is frequently applied to minimize the effects of respiratory motion, but scan durations are subject-dependent and hence difficult to predict.

PURPOSE

To implement cardiac proton MRS with fixed scan time by employing retrospective phase-based gating and to compare the proposed method to conventional navigator-gated MRS.

STUDY TYPE

Prospective.

SUBJECTS

Eighteen healthy volunteers (29.7 ± 7.8 years).

FIELD STRENGTH/SEQUENCE

1.5, navigator-gated (16 averages without, 96 with water suppression [WS]) data acquisition as reference and navigator-free data acquisition with a fixed scan time (48 without WS, 304 with WS), cardiac-triggered point-resolved spectroscopy (PRESS).

ASSESSMENT

Navigator-free data acquisition with retrospective phase-based gating was compared with prospective navigator-gating. Navigator-free acquisition was repeated in 10 subjects to assess reproducibility. Scan time was assessed for prospective and retrospective gating. Retrospective phase-based gating was performed using a threshold based on the standard deviation (SD) of individual water (W) and triglyceride (TG) phases.

STATISTICAL TESTS

T-tests and Bland-Altman analysis.

RESULTS

The duration of the prospective navigator-gated scans ranged from 6:09 minutes to 21:50 minutes (mean 10:05 ± 3:46 min, gating efficiency 40.4 ± 10.5%), while data acquisition for retrospective phase-based gating had a fixed scan time of 11:44 minutes. Retrospective phase-based gating using a threshold of 1 × SD yielded a gating efficiency of 72.7 ± 4.3% and a coefficient of variation (CoV) of triglyceride-to-water ratios of 9.8% compared with the navigated reference. The intrasubject reproducibility of retrospective gating revealed a CoV of 9.5%.

DATA CONCLUSION

Cardiac proton MRS employing retrospective phase-based gating is feasible and provides reproducible assessment of TG/W in a fixed scan time. Since scan time is independent of respiratory motion, retrospective phase-based gating offers an approach to motion compensation with predictable exam time for proton MRS of the heart.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;50:1973-1981.

Physiological and pathological skeletal muscle T1 changes quantified using a fast inversion-recovery radial NMR imaging sequence

We investigated the response of skeletal muscle global T1 under different physiological and pathological conditions using an inversion-recovery radial T1 mapping sequence. Thirty five healthy volunteers, seven patients with Becker muscular dystrophy (BMD) and seven patients with sporadic inclusion body myositis (IBM) were investigated in order to evaluate the effects of gender, age, muscle group, exercise and pathological processes on global T1 values. In addition, the intramuscular fat content was measured using 3-point Dixon and the global T2 and water T2 (T2_H2O) were determined with a multi-spin-echo sequence. In the muscles of healthy volunteers, there was no impact of age on global T1. However, we measured a significant effect of sex and muscle group. After exercise, a significant 7.7% increase of global T1 was measured in the recruited muscles, and global T1 variations were highly correlated to T2_H2O variations (R = 0.91). In pathologies, global T1 values were reduced in fat infiltrated muscles. When fat fraction was taken into account, global T1 values were higher in IBM patients compared to BMD. Global T1 variations are a sensitive indicator of tissue changes in skeletal muscle related to several physiological and pathological events.

Assessment of sympathetic reinnervation after cardiac transplantation using hybrid cardiac PET/MRI: A pilot study

Background

Sympathetic reinnervation after heart transplantation (HTX) is a known phenomenon, which has an impact on patient heart rate variability and exercise capacity. The impact of reinnervation on myocardial structure has not been evaluated yet.

Propose

To evaluate the feasibility of simultaneous imaging of cardiac reinnervation and cardiac structure using a hybrid PET/MRI system.

Study type

Prospective / pilot study.

Subjects

Ten patients, 4–21 years after cardiac transplantation.

Field Strength/Sequence

3 T hybrid PET/MRI system. Cine SSFP, T₁ mapping (modified Look–Locker inversion recovery sequence) pre/postcontrast as well as dynamic [¹¹C]meta‐hydroxyephedrine ([¹¹C]mHED) PET.

Assessment

All MRI and PET parameters were evaluated by experienced readers using dedicated postprocessing software packages for cardiac MRI and PET. For all parameters a 16‐segment model for the left ventricle was applied.

Statistical Tests

Mann–Whitney U‐test; Spearman correlations.

Results

Thirty‐six of 160 myocardial segments showed evidence of reinnervation by PET. On a segment‐based analysis, mean native T₁ relaxation times were nonsignificantly altered in segments with evidence of reinnervation (1305 ± 151 msec vs. 1270 ± 112 msec; P = 0.1), whereas mean extracellular volume (ECV) was significantly higher in segments with evidence of reinnervation (35.8 ± 11% vs. 30.9 ± 7%; P = 0.019). There were no significant differences in wall motion (WM) and wall thickening (WT) between segments with or without reinnervation (mean WM: 7.6 ± 4 mm vs. group B: 9.3 ± 7 mm [P = 0.13]; WT: 79 ± 63% vs. 94 ± 74% [P = 0.27]) under resting conditions.

Data Conclusion

The assessment of cardiac reinnervation using a hybrid PET/MRI system is feasible. Segments with evidence of reinnervation by PET showed nonsignificantly higher T₁ relaxation times and a significantly higher ECV, suggesting a higher percentage of diffuse fibrosis in these segments, without impairment of rest WM and WT.

Level of Evidence: 3

Technical Efficacy: Stage 3

J. Magn. Reson. Imaging 2019;50:1326–1335.

Reference range determination for imaging biomarkers: Myocardial T1

BACKGROUND

Imaging biomarkers, such as the T₁ relaxation time of the myocardium using MRI, can be valuable in cardiac medicine if they are properly validated. Consensus statements recommend that for myocardial T₁ , each investigator should establish a reference range.

PURPOSE

To describe a statistically valid method for determining and reporting the reference range in each center, which simultaneously minimizes the twin risks of undersampling, leading to a uselessly uncertain range, and oversampling, which exposes volunteers to unnecessary scanning and wastes resources.

STUDY TYPE

Cohort.

POPULATION

In all, 278 normal human subjects without cardiac disease from two cardiac MR centers.

FIELD STRENGTH/SEQUENCE

1.5 T and 3 T; Modified Look-Locker Inversion recovery sequence.

ASSESSMENT

The T₁ relaxation time was estimated from multiple samples of tissue magnetization after inversion. A valid method for calculating a reference range was used.

STATISTICAL TESTS

Shapiro-Wilk test for normality; Tukey robust approach for identification of outliers; reference range calculation with confidence intervals.

RESULTS

Reference ranges for measurement of myocardial T₁ were calculated, with confidence intervals, enabling comparison with clinically important differences. At 3 T: 1129 to 1301 msec at site 1 (n = 21) and 1160 to 1309 msec at site 2 (n = 59), and at 1.5 T at site 2: 933 to 1020 msec (male, n = 130) and 965 to 1054 msec (female, n = 68). The 3 T reference range from site 1 was successfully benchmarked against the 3 T reference range at site 2.

DATA CONCLUSION

Myocardial T₁ reference ranges can be properly characterized, enabling clinical comparison to a valid reference range with known confidence intervals, using methodology similar to that described in this report.

LEVEL OF EVIDENCE

3 Technical Efficacy Stage: 3 J. Magn. Reson. Imaging 2019;50:771-778.

Real-time Triggered RAdial Single-Shot Inversion recovery for arrhythmia-insensitive myocardial T1 mapping: motion phantom validation and in vivo comparison

PURPOSE

Cardiac T₁ mapping has become an increasingly important imaging technique, contributing novel diagnostic options. However, currently utilized methods are often associated with accuracy problems because of heart rate variations and cardiac arrhythmia, limiting their value in clinical routine. This study aimed to introduce an improved arrhythmia-related robust T₁ mapping sequence called RT-TRASSI (real-time Triggered RAdial Single-Shot Inversion recovery).

METHODS

All measurements were performed on a 3.0T whole-body imaging system. A real-time feedback algorithm for arrhythmia detection was implemented into the previously described pulse sequence. A programmable motion phantom was constructed and measurements with different simulated arrhythmias arranged. T₁ mapping accuracy and susceptibility to artifacts were analyzed. In addition, in vivo measurements and comparisons with 3 prevailing T₁ mapping sequences (MOLLI, ShMOLLI, and SASHA) were carried out to investigate the occurrence of artifacts.

RESULTS

In the motion phantom measurements, RT-TRASSI showed excellent agreement with predetermined reference T₁ values. Percentage scattering of the T₁ values ranged from -0.6% to +1.9% in sinus rhythm and -1.0% to +3.1% for high-grade arrhythmias. In vivo, RT-TRASSI showed diagnostic image quality with only 6% of the acquired T₁ maps including image artifacts. In contrast, more than 40% of the T₁ maps acquired with MOLLI, ShMOLLI, or SASHA included motion artifacts.

CONCLUSION

Accuracy issues because of heart rate variability and arrhythmia are a prevailing problem in current cardiac T₁ mapping techniques. With RT-TRASSI, artifacts can be minimized because of the short acquisition time and effective real-time feedback, avoiding potential data acquisition during systolic heart phase.

Cardiac- versus diaphragm-based respiratory navigation for proton spectroscopy of the heart

OBJECTIVES

To study inter-individual differences of the relation between diaphragm and heart motion, the objective of the present study was to implement respiratory navigation on the heart and compare it against the established method of navigator gating on the diaphragm for single-voxel cardiac ¹H-MRS.

MATERIALS AND METHODS

¹H-MRS was performed on a 1.5T system in 19 healthy volunteers of mixed age (range 24-75 years). Spectra were recorded in a 6-8 ml voxel in the ventricular septum using a PRESS (point-resolved spectroscopy) sequence and ECG gating. Water-unsuppressed data acquired with pencil beam navigation on the heart were compared to data with navigation on the diaphragm. Water-suppressed data were obtained to assess triglyceride-to-water ratios.

RESULTS

Water phase and amplitude fluctuations for cardiac versus diaphragm navigation did not reveal significant differences. Both navigator positions provided comparable triglyceride-to-water ratios and gating efficiencies (coefficient of variation (CoV) 7.0%). The cardiac navigator showed a good reproducibility (CoV 5.2%).

DISCUSSION

Respiratory navigation on the heart does not convey an advantage over diaphragm-based navigator gating for cardiac ¹H-MRS, but also no disadvantage. Consequently, cardiac and diaphragm respiratory navigation may be used interchangeably.

Cardiac magnetic resonance T1 mapping. Part 1: Aspects of acquisition and evaluation

While an enormous number of studies have documented pathological alterations of the myocardial native longitudinal relaxation time (T1) and the fraction of the extracellular myocardial volume (ECV), it has also become clear that continuously evolving T1 mapping sequence, acquisition and evaluation techniques have a substantial impact on quantitative results, making the translation of reported findings into routine clinical use particularly challenging. To provide a basis for the discussion of pathological myocardial T1 and ECV alterations, the present review aims to summarize the methodological aspects of myocardial T1 mapping along with technical and physiological factors influencing results and normal ranges of myocardial native T1 and ECV reported across studies.

Accurate and robust systolic myocardial T1 mapping using saturation recovery with individualized delay time: comparison with diastolic T1 mapping

T₁ mapping data are generally acquired in patients' diastolic phase, wherein their myocardium is the thinnest in the cardiac cycle. However, the analysis of the thin myocardium may cause errors in image registrations and settings related to the region of interest. In this study, we validated systolic T₁ mapping using the saturation recovery with individualized delay time (SR-IDT) method and compared it with conventional diastolic T₁ mapping. Both diastolic and systolic T₁ mappings were performed in the mid-ventricular plane in 10 healthy volunteers (35 ± 9 years, 9 males) and 29 consecutive patients with cardiac diseases (68 ± 14 years, 19 males). Comparison of the myocardial T₁ value at diastole and systole was performed with both the Pearson correlation coefficient (r) and the Bland-Altman analysis. Additionally, the systolic myocardial T₁ value was compared between the volunteers and patients by using Tukey's test. Pearson correlation analysis demonstrated a strong positive correlation between diastolic and systolic T₁ values (r = 0.88, P < 0.001). The Bland-Altman plot suggested that left ventricular T₁ values in the diastole and systole showed high agreement (mean difference and 95% limits of agreement = 17 ± 104 ms). Further, systolic T₁ values with SR-IDT in patients in the late gadolinium enhancement (LGE) group were significantly higher than those in the control group (1585 ± 118 ms vs 1469 ± 69 ms; P = 0.024). Therefore, the proposed systolic T₁ mapping with the SR-IDT, which was validated with respect to the conventional diastolic method, is a useful clinical tool for the quantitative characterization of the myocardium.

A functional form for a representative individual arterial input function measured from a population using high temporal resolution DCE MRI

PURPOSE

To measure the arterial input function (AIF), an essential component of tracer kinetic analysis, in a population of patients using an optimized dynamic contrast-enhanced (DCE) imaging sequence and to estimate inter- and intrapatient variability. From these data, a representative AIF that may be used for realistic simulation studies can be extracted.

METHODS

Thirty-nine female patients were imaged on multiple visits before and during a course of neoadjuvant chemotherapy for breast cancer. A total of 97 T₁ -weighted DCE studies were analyzed including bookend estimates of T₁ and model-fitting to each individual AIF. Area under the curve and cardiac output were estimated from each first pass peak, and these data were used to assess inter- and intrapatient variability of the AIF.

RESULTS

Interpatient variability exceeded intrapatient variability of the AIF. There was no change in cardiac output as a function of MR visit (mean value 5.6 ± 1.1 L/min) but baseline blood T₁ increased significantly following the start of chemotherapy (which was accompanied by a decrease in hematocrit).

CONCLUSION

The AIF in an individual patient can be measured reproducibly but the variability of AIFs between patients suggests that use of a population AIF will decrease the precision of tracer kinetic analysis performed in cross-patient comparison studies. A representative AIF is presented that is typical of the population but retains the characteristics of an individually measured AIF.

T1 and T2 mapping in the identification of acute myocardial injury in patients with NSTEMI

AIMS

To test T1 and T2 mapping in the assessment of acute myocardial injury in patients with non-ST-segment elevation myocardial infarction (NSTEMI), evaluated before revascularization.

METHODS

Forty-seven patients with acute NSTEMI underwent cardiac magnetic resonance (CMR) at 1.5 T, including T1 and T2 mapping.

RESULTS

Coronary angiography (CA) evidenced an obstructive coronary artery disease (CAD) in 36 patients (80%) and a non-obstructive CAD in 11 patients (20%). Edema was detected in 51.1/65.9% of patients in T1/T2 maps, respectively. This difference was due to artifacts in T1 maps. T1/T2 values were significantly higher in the infarcted myocardium (IM) compared with the remote myocardium (RM) (in T1: 1151.6 ± 53.5 ms vs. 958.2 ± 38.6 ms, respectively; in T2: 69 ± 6 ms vs. 51.9 ± 2.9 ms, respectively; p < 0.0001 for both). We found both an obstructive CAD at CA and myocardial edema at CMR in 53.2% of patients, while 8.5% of patients had a non-obstructive CAD and no edema. However, 25.5% of patients had an obstructive CAD without edema, while 12.8% of patients showed edema despite a non-obstructive CAD. Furthermore, in 6 of the edema-positive patients with multi-vessels obstructive CAD, CMR identified myocardial edema in a vascular territory different from that of the lesion supposed to be the culprit at CA.

CONCLUSIONS

In a non-negligible percentage of NSTEMI patients, T1 and T2 mapping detect myocardial edema without significant stenosis at CA and vice versa. Therefore, CA and CMR edema imaging might provide complementary information in the evaluation of NSTEMI.

Inverse association of MRI-derived native myocardial T1 and perfusion reserve index in women with evidence of ischemia and no obstructive CAD: A pilot study

BACKGROUND

It has recently been shown that magnetic resonance (MR) "native T1" mapping is capable of characterizing abnormal microcirculation in patients with obstructive coronary artery disease (CAD). In studies involving women with signs and symptoms of ischemia and no obstructive CAD (INOCA), however, the potential role of native T1 as an imaging marker and its association with indices of diastolic function or vasodilator-induced myocardial ischemia have not been explored. We investigated whether native T1 in INOCA is associated with reduced myocardial perfusion reserve index (MPRI) or with diastolic dysfunction.

METHODS

Twenty-two female patients with INOCA and twelve female reference controls with matching age and body-mass index were studied. The patients had evidence of vasodilator-induced ischemia without obstructive CAD or any prior infarction. All 34 subjects underwent stress/rest MR including native T1 mapping (MOLLI 5(3)3) at 1.5-Tesla.

RESULTS

Compared with controls, patients had similar morphology/function. As expected, MPRI was significantly reduced in patients compared to controls (1.78 ± 0.39 vs. 2.49 ± 0.41, p < 0.0001). Native T1 was significantly elevated in patients (1040.1 ± 29.3 ms vs. 1003.8 ± 18.5 ms, p < 0.001) and the increased T1 showed a significant inverse correlation with MPRI (r = -0.481, p = 0.004), but was not correlated with reduced diastolic strain rate.

CONCLUSIONS

Symptomatic women with INOCA have elevated native T1 compared to matched reference controls and there is a significant association between elevated native T1 and impaired MPRI, considered a surrogate measure of ischemia severity in this cohort. Future studies in a larger cohort are needed to elucidate the mechanism underlying this inverse relationship.

Evaluation of Modified Look-Locker Inversion Recovery and Arrhythmia-Insensitive Rapid Cardiac T1 Mapping Pulse Sequences in Cardiomyopathy Patients

OBJECTIVE

The aim of this study was to compare the performance of arrhythmia-insensitive rapid (AIR) and modified Look-Locker inversion recovery (MOLLI) T1 mapping in patients with cardiomyopathies.

METHODS

In 58 patients referred for clinical cardiac magnetic resonance imaging at 1.5 T, we compared MOLLI and AIR native and postcontrast T1 measurements. Two readers independently analyzed myocardial and blood T1 values. Agreement between techniques, interreader agreement per technique, and intrascan agreement per technique were evaluated.

RESULTS

The MOLLI and AIR T1 values were strongly correlated (r = 0.98); however, statistically significantly different T1 values were derived (bias 80 milliseconds, pooled data, P < 0.01). Both techniques demonstrated high repeatability (MOLLI, r = 1.00 and coefficient of repeatability [CR] = 72 milliseconds; AIR, r = 0.99 and CR = 184.2 milliseconds) and produced high interreader agreement (MOLLI, r = 1.00 and CR = 51.7 milliseconds; AIR, r = 0.99 and CR = 183.5 milliseconds).

CONCLUSIONS

Arrhythmia-insensitive rapid and MOLLI sequences produced significantly different T1 values in a diverse patient cohort.

A Framework for the Generation of Realistic Synthetic Cardiac Ultrasound and Magnetic Resonance Imaging Sequences From the Same Virtual Patients

The use of synthetic sequences is one of the most promising tools for advanced in silico evaluation of the quantification of cardiac deformation and strain through 3-D ultrasound (US) and magnetic resonance (MR) imaging. In this paper, we propose the first simulation framework which allows the generation of realistic 3-D synthetic cardiac US and MR (both cine and tagging) image sequences from the same virtual patient. A state-of-the-art electromechanical (E/M) model was exploited for simulating groundtruth cardiac motion fields ranging from healthy to various pathological cases, including both ventricular dyssynchrony and myocardial ischemia. The E/M groundtruth along with template MR/US images and physical simulators were combined in a unified framework for generating synthetic data. We efficiently merged several warping strategies to keep the full control of myocardial deformations while preserving realistic image texture. In total, we generated 18 virtual patients, each with synthetic 3-D US, cine MR, and tagged MR sequences. The simulated images were evaluated both qualitatively by showing realistic textures and quantitatively by observing myocardial intensity distributions similar to real data. In particular, the US simulation showed a smoother myocardium/background interface than the state-of-the-art. We also assessed the mechanical properties. The pathological subjects were discriminated from the healthy ones by both global indexes (ejection fraction and the global circumferential strain) and regional strain curves. The synthetic database is comprehensive in terms of both pathology and modality, and has a level of realism sufficient for validation purposes. All the 90 sequences are made publicly available to the research community via an open-access database.

Temporally resolved parametric assessment of Z-magnetization recovery (TOPAZ): Dynamic myocardial T1 mapping using a cine steady-state look-locker approach

PURPOSE

To develop and evaluate a cardiac phase-resolved myocardial T₁ mapping sequence.

METHODS

The proposed method for temporally resolved parametric assessment of Z-magnetization recovery (TOPAZ) is based on contiguous fast low-angle shot imaging readout after magnetization inversion from the pulsed steady state. Thereby, segmented k-space data are acquired over multiple heartbeats, before reaching steady state. This results in sampling of the inversion-recovery curve for each heart phase at multiple points separated by an R-R interval. Joint T₁ and B1+ estimation is performed for reconstruction of cardiac phase-resolved T₁ and B1+ maps. Sequence parameters are optimized using numerical simulations. Phantom and in vivo imaging are performed to compare the proposed sequence to a spin-echo reference and saturation pulse prepared heart rate-independent inversion-recovery (SAPPHIRE) T₁ mapping sequence in terms of accuracy and precision.

RESULTS

In phantom, TOPAZ T₁ values with integrated B1+ correction are in good agreement with spin-echo T₁ values (normalized root mean square error = 4.2%) and consistent across the cardiac cycle (coefficient of variation = 1.4 ± 0.78%) and different heart rates (coefficient of variation = 1.2 ± 1.9%). In vivo imaging shows no significant difference in TOPAZ T₁ times between the cardiac phases (analysis of variance: P = 0.14, coefficient of variation = 3.2 ± 0.8%), but underestimation compared with SAPPHIRE (T₁ time ± precision: 1431 ± 56 ms versus 1569 ± 65 ms). In vivo precision is comparable to SAPPHIRE T₁ mapping until middiastole (P > 0.07), but deteriorates in the later phases.

CONCLUSIONS

The proposed sequence allows cardiac phase-resolved T₁ mapping with integrated B1+ assessment at a temporal resolution of 40 ms. Magn Reson Med 79:2087-2100, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Segment-by-segment assessment of left ventricular myocardial affection in Anderson-Fabry disease by non-enhanced T1-mapping

Background Anderson-Fabry disease (AFD) is an X-linked lysosomal enzyme disorder associated with an intracellular accumulation of sphingolipids, which shorten myocardial T1 relaxation times. Myocardial affection, however, varies between different segments. Purpose To evaluate the specific segmental distribution and degree of segmental affection in AFD patients. Material and Methods Twenty-five patients with AFD, 14 patients with hypertrophic cardiomyopathy (HCM), and 21 controls were included. A Modified Look-Locker Inversion Recovery sequence (MOLLI) was used for non-enhanced T1 mapping at 1.5 T in addition to standard cardiac imaging in 10-12 short axis views. T1 values were evaluated with a mixed model ANOVA and regression analysis to determine the best diagnostic cutoff values for T1 for each myocardial segment. Results Regression analysis showed the best diagnostic cutoff compared to controls in cardiac segments 1-4, 8-9, and 14. Mean differences between T1 for AFD versus HCM were greatest in segment 3, 4, and 9 (99 ms, 103 ms, 86 ms, respectively). Overall T1 times were 888 ± 70 ms and 903 ± 14 ms (AFD with and without LVH); 1014 ± 17 ms and 1001 ± 22 ms (HCM and controls, P < 0.05). Conclusion Myocardial segments are affected by a varying degree of T1 shortening in AFD patients. Segment-specific cutoff values allow the most specific detection and quantification of the extent of myocardial affection.

Quantification of both the area-at-risk and acute myocardial infarct size in ST-segment elevation myocardial infarction using T1-mapping

BACKGROUND

A comprehensive cardiovascular magnetic resonance (CMR) in reperfused ST-segment myocardial infarction (STEMI) patients can be challenging to perform and can be time-consuming. We aimed to investigate whether native T1-mapping can accurately delineate the edema-based area-at-risk (AAR) and post-contrast T1-mapping and synthetic late gadolinium (LGE) images can quantify MI size at 1.5 T. Conventional LGE imaging and T2-mapping could then be omitted, thereby shortening the scan duration.

METHODS

Twenty-eight STEMI patients underwent a CMR scan at 1.5 T, 3 ± 1 days following primary percutaneous coronary intervention. The AAR was quantified using both native T1 and T2-mapping. MI size was quantified using conventional LGE, post-contrast T1-mapping and synthetic magnitude-reconstructed inversion recovery (MagIR) LGE and synthetic phase-sensitive inversion recovery (PSIR) LGE, derived from the post-contrast T1 maps.

RESULTS

Native T1-mapping performed as well as T2-mapping in delineating the AAR (41.6 ± 11.9% of the left ventricle [% LV] versus 41.7 ± 12.2% LV, P = 0.72; R2 0.97; ICC 0.986 (0.969-0.993); bias -0.1 ± 4.2% LV). There were excellent correlation and inter-method agreement with no bias, between MI size by conventional LGE, synthetic MagIR LGE (bias 0.2 ± 2.2%LV, P = 0.35), synthetic PSIR LGE (bias 0.4 ± 2.2% LV, P = 0.060) and post-contrast T1-mapping (bias 0.3 ± 1.8% LV, P = 0.10). The mean scan duration was 58 ± 4 min. Not performing T2 mapping (6 ± 1 min) and conventional LGE (10 ± 1 min) would shorten the CMR study by 15-20 min.

CONCLUSIONS

T1-mapping can accurately quantify both the edema-based AAR (using native T1 maps) and acute MI size (using post-contrast T1 maps) in STEMI patients without major cardiovascular risk factors. This approach would shorten the duration of a comprehensive CMR study without significantly compromising on data acquisition and would obviate the need to perform T2 maps and LGE imaging.

Blood correction reduces variability and gender differences in native myocardial T1 values at 1.5 T cardiovascular magnetic resonance - a derivation/validation approach

BACKGROUND

Myocardial native T1 measurements are likely influenced by intramyocardial blood. Since blood T1 is both variable and longer compared to myocardial T1, this will degrade the precision of myocardial T1 measurements. Precision could be improved by correction, but the amount of correction and the optimal blood T1 variables to correct with are unknown. We hypothesized that an appropriate correction would reduce the standard deviation (SD) of native myocardial T1.

METHODS

Consecutive patients (n = 400) referred for CMR with known or suspected heart disease were split into a derivation cohort for model construction (n = 200, age 51 ± 18 years, 50% male) and a validation cohort for assessing model performance (n = 200, age 48 ± 17 years, 50% male). Exclusion criteria included focal septal abnormalities. A Modified Look-Locker inversion recovery sequence (MOLLI, 1.5 T Siemens Aera) was used to acquire T1 and T1* maps. T1 and T1* maps were used to measure native myocardial T1, and blood T1 and T1*. A multivariate linear regression correction model was implemented using blood measurement of R1 (1/T1), R1* (1/T1*) or hematocrit. The correction model from the derivation cohort was applied to the validation cohort, and assessed for reduction in variability with the F-test.

RESULTS

Blood [LV + RV] mean R1, mean R1* and hematocrit correlated with myocardial T1 (Pearson's r, range 0.37 to 0.45, p < 0.05 for all) in both the derivation and validation cohorts respectively, suggesting that myocardial T1 measurements are influenced by intramyocardial blood. Mean myocardial native T1 did not differ between the derivation and validation cohorts (1030 ± 42.6 ms and 1023 ± 45.2 ms respectively, p = 0.07). In the derivation cohort, correction using blood mean R1 and mean R1* yielded a decrease in myocardial T1 SD (45.2 ms to 36.6 ms, p = 0.03). When the model from the derivation cohort was applied to the validation cohort, the SD reduction was maintained (39.3 ms, p = 0.049). This 13% reduction in measurement variability leads to a 23% reduction in sample size to detect a 50 ms difference in native myocardial T1.

CONCLUSIONS

Correcting native myocardial T1 for R1 and R1* of blood improves the precision of myocardial T1 measurement by ~13%, and could consequently improve disease detection and reduce sample size needs for clinical research.

Myocardial T1 and T2 Mapping: Techniques and Clinical Applications

Cardiac magnetic resonance (CMR) imaging is widely used in various medical fields related to cardiovascular diseases. Rapid technological innovations in magnetic resonance imaging in recent times have resulted in the development of new techniques for CMR imaging. T1 and T2 image mapping sequences enable the direct quantification of T1, T2, and extracellular volume fraction (ECV) values of the myocardium, leading to the progressive integration of these sequences into routine CMR settings. Currently, T1, T2, and ECV values are being recognized as not only robust biomarkers for diagnosis of cardiomyopathies, but also predictive factors for treatment monitoring and prognosis. In this study, we have reviewed various T1 and T2 mapping sequence techniques and their clinical applications.

Myocardial T1-mapping at 3T using saturation-recovery: reference values, precision and comparison with MOLLI

BACKGROUND

Myocardial T1-mapping recently emerged as a promising quantitative method for non-invasive tissue characterization in numerous cardiomyopathies. Commonly performed with an inversion-recovery (IR) magnetization preparation at 1.5T, the application at 3T has gained due to increased quantification precision. Alternatively, saturation-recovery (SR) T1-mapping has recently been introduced at 1.5T for improved accuracy. Thus, the purpose of this study is to investigate the robustness and precision of SR T1-mapping at 3T and to establish accurate reference values for native T1-times and extracellular volume fraction (ECV) of healthy myocardium.

METHODS

Balanced Steady-State Free-Precession (bSSFP) Saturation-Pulse Prepared Heart-rate independent Inversion-REcovery (SAPPHIRE) and Saturation-recovery Single-SHot Acquisition (SASHA) T1-mapping were compared with the Modified Look-Locker inversion recovery (MOLLI) sequence at 3T. Accuracy and precision were studied in phantom. Native and post-contrast T1-times and regional ECV were determined in 20 healthy subjects (10 men, 27 ± 5 years). Subjective image quality, susceptibility artifact rating, in-vivo precision and reproducibility were analyzed.

RESULTS

SR T1-mapping showed <4 % deviation from the spin-echo reference in phantom in the range of T1 = 100-2300 ms. The average quality and artifact scores of the T1-mapping methods were: MOLLI:3.4/3.6, SAPPHIRE:3.1/3.4, SASHA:2.9/3.2; (1: poor - 4: excellent/1: strong - 4: none). SAPPHIRE and SASHA yielded significantly higher T1-times (SAPPHIRE: 1578 ± 42 ms, SASHA: 1523 ± 46 ms), in-vivo T1-time variation (SAPPHIRE: 60.1 ± 8.7 ms, SASHA: 70.0 ± 9.3 ms) and lower ECV-values (SAPPHIRE: 0.20 ± 0.02, SASHA: 0.21 ± 0.03) compared with MOLLI (T1: 1181 ± 47 ms, ECV: 0.26 ± 0.03, Precision: 53.7 ± 8.1 ms). No significant difference was found in the inter-subject variability of T1-times or ECV-values (T1: p = 0.90, ECV: p = 0.78), the observer agreement (inter: p > 0.19; intra: p > 0.09) or consistency (inter: p > 0.07; intra: p > 0.17) between the three methods.

CONCLUSIONS

Saturation-recovery T1-mapping at 3T yields higher accuracy, comparable inter-subject, inter- and intra-observer variability and less than 30 % precision-loss compared to MOLLI.

Heart deformation analysis: the distribution of regional myocardial motion patterns at left ventricle

The aim of the present study was to test the hypothesis that heart deformation analysis (HDA) is able to discriminate regional myocardial motion patterns on the left ventricle (LV). Totally 21 healthy volunteers (15 men and 6 women) without documented cardiovascular diseases were recruited. Cine MRI was performed on those subjects at four-chamber, two-chamber, and short-axis views. The variations of segmental myocardial motion indices of the LV, which were measured with the HDA tool, were investigated. Regional displacement, velocity, strain and strain rate were compared between lateral wall and septal wall using t tests. There are significant variations (CoV = 18.0-72.4%) of myocardial motion indices (average over 21 subjects) among 16 myocardial segments. There are significant differences (p < 0.05) between displacement, velocity, strain and strain rate measured at lateral and septal areas of the LV. In conclusion, HDA is able to present different regional LV motion patterns from multiple aspects in healthy volunteers.

Myocardial Mapping With Cardiac Magnetic Resonance: The Diagnostic Value of Novel Sequences

Cardiac magnetic resonance has evolved into a crucial modality for the evaluation of cardiomyopathy due to its ability to characterize myocardial structure and function. In the last few years, interest has increased in the potential of "mapping" techniques that provide direct and objective quantification of myocardial properties such as T1, T2, and T2* times. These approaches enable the detection of abnormalities that affect the myocardium in a diffuse fashion and/or may be too subtle for visual recognition. This article reviews the current state of myocardial T1 and T2-mapping in both health and disease.

Native myocardial T1 mapping in pulmonary hypertension: correlations with cardiac function and hemodynamics

Objectives

To analyze alterations in left ventricular (LV) myocardial T1 times in patients with pulmonary hypertension (PH) and to investigate their associations with ventricular function, mass, geometry and hemodynamics.

Methods

Fifty-eight patients with suspected PH underwent right heart catheterization (RHC) and 3T cardiac magnetic resonance imaging. Ventricular function, geometry and mass were derived from cine real-time short-axis images. Myocardial T1 maps were acquired by a prototype modified Look-Locker inversion-recovery sequence in short-axis orientations. LV global, segmental and ventricular insertion point (VIP) T1 times were evaluated manually and corrected for blood T1.

Results

Septal, lateral, global and VIP T1 times were significantly higher in PH than in non-PH subjects (septal, 1249 ± 58 ms vs. 1186 ± 33 ms, p < 0.0001; lateral, 1190 ± 45 ms vs. 1150 ± 33 ms, p = 0.0003; global, 1220 ± 52 ms vs. 1171 ± 29 ms, p < 0.0001; VIP, 1298 ± 78 ms vs. 1193 ± 31 ms, p < 0.0001). In PH, LV eccentricity index was the strongest linear predictor of VIP T1 (r = 0.72). Septal, lateral and global T1 showed strong correlations with VIP T1 (r = 0.81, r = 0.59 and r = 0.75, respectively).

Conclusions

In patients with PH, T1 times in VIPs and in the entire LV myocardium are elevated. LV eccentricity strongly correlates with VIP T1 time, which in turn is strongly associated with T1 time changes in the entire LV myocardium.

Key Points

• Native T1 mapping detects left ventricular myocardial alterations in pulmonary hypertension

• In pulmonary hypertension, native T1 times at ventricular insertion points are increased

• These T1 times correlate strongly with left ventricular eccentricity

• In pulmonary hypertension, global and segmental myocardial T1 times are increased

• Global, segmental and ventricular insertion point T1 times are strongly correlated

Systolic MOLLI T1 mapping with heart-rate-dependent pulse sequence sampling scheme is feasible in patients with atrial fibrillation

BACKGROUND

T1 mapping enables assessment of myocardial characteristics. As the most common type of arrhythmia, atrial fibrillation (AF) is often accompanied by a variety of cardiac pathologies, whereby the irregular and usually rapid ventricle rate of AF may cause inaccurate T1 estimation due to mis-triggering and inadequate magnetization recovery. We hypothesized that systolic T1 mapping with a heart-rate-dependent (HRD) pulse sequence scheme may overcome this issue.

METHODS

30 patients with AF and 13 healthy volunteers were enrolled and underwent cardiovascular magnetic resonance (CMR) at 3 T. CMR was repeated for 3 patients after electric cardioversion and for 2 volunteers after lowering heart rate (HR). A Modified Look-Locker Inversion Recovery (MOLLI) sequence was acquired before and 15 min after administration of 0.1 mmol/kg gadopentetate dimeglumine. For AF patients, both the fixed 5(3)3/4(1)3(1)2 and the HRD sampling scheme were performed at diastole and systole, respectively. The HRD pulse sequence sampling scheme was 5(n)3/4(n)3(n)2, where n was determined by the heart rate to ensure adequate magnetization recovery. Image quality of T1 maps was assessed. T1 times were measured in myocardium and blood. Extracellular volume fraction (ECV) was calculated.

RESULTS

In volunteers with repeated T1 mapping, the myocardial native T1 and ECV generated from the 1st fixed sampling scheme were smaller than from the 1st HRD and 2nd fixed sampling scheme. In healthy volunteers, the overall native T1 times and ECV of the left ventricle (LV) in diastolic T1 maps were greater than in systolic T1 maps (P < 0.01, P < 0.05). In the 3 AF patients that had received electrical cardioversion therapy, the myocardial native T1 times and ECV generated from the fixed sampling scheme were smaller than in the 1st and 2nd HRD sampling scheme (all P < 0.05). In patients with AF (HR: 88 ± 20 bpm, HR fluctuation: 12 ± 9 bpm), more T1 maps with artifact were found in diastole than in systole (P < 0.01). The overall native T1 times and ECV of the left ventricle (LV) in diastolic T1 maps were greater than systolic T1 maps, either with fixed or HRD sampling scheme (all P < 0.05).

CONCLUSION

Systolic MOLLI T1 mapping with heart-rate-dependent pulse sequence scheme can improve image quality and avoid T1 underestimation. It is feasible and with further validation may extend clinical applicability of T1 mapping to patients with atrial fibrillation.