Three-Dimensional Volumetric Assessment of Diastolic Function by Cardiac Magnetic Resonance Imaging: The Multi-Ethnic Study of Atherosclerosis (MESA)

Abstract

Background:

Cardiac Magnetic Resonance is in need of a simple and robust method for diastolic function assessment that can be done with routine protocol sequences.

Objective:

To develop and validate a three-dimensional (3D) model-based volumetric assessment of diastolic function using cardiac magnetic resonance (CMR) imaging and compare the results obtained with the model with those obtained by echocardiography.

Methods:

The study participants provided written informed consent and were included if having undergone both echocardiography and cine steady-state free precession (SSFP) CMR on the same day. Guide points at the septal and lateral mitral annulus were used to define the early longitudinal relaxation rate (E'), while a time-volume curve from the 3D model was used to assess diastolic filling parameters. We determined the correlation between 3D CMR and echocardiography and the accuracy of CMR in classifying the diastolic function grade.

Results:

The study included 102 subjects. The E/A ratio by CMR was positively associated with the E/A ratio by echocardiography (r = 0.71, p < 0.0001). The early diastolic relaxation velocity by tissue Doppler and longitudinal relaxation rate for the lateral mitral annulus displacement were positively associated (p = 0.007), as were the ratio between Doppler E/e' and CMR E/E' (p = 0.01). CMR-determined normalized peak E (NE) and deceleration time (DT) were able to predict diastolic dysfunction (areas under the curve [AUCs] = 0.70 and 0.72, respectively). In addition, the lateral E/E' ratio showed good utility in identifying diastolic dysfunction (AUC = 0.80). Overall, echocardiography and CMR interobserver and intraobserver agreements were excellent (intraclass correlation coefficient range 0.72 - 0.97).

Conclusion:

3D modeling of standard cine CMR images was able to identify study subjects with reduced diastolic function and showed good reproducibility, suggesting a potential for a routine diastolic function assessment by CMR.

Keywords:
Ventricular Function; Evaluation; Magnetic Resonance; Imaging Three Dimensional; Echocardiography, Three -Dimensional

Resumo

Fundamento:

A ressonância magnética cardíaca necessita de um método simples e robusto para a avaliação da função diastólica que pode ser feito com sequências protocolares de rotina.

Objetivo:

Desenvolver e validar a avaliação volumétrica da função diastólica através de um modelo tridimensional (3D) com utilização de imagens de ressonância magnética cardíaca (RMC) e comparar os resultados obtidos com este modelo com os obtidos por ecocardiografia.

Métodos:

Os participantes do estudo assinaram um termo de consentimento e foram incluídos se tivessem sido submetidos no mesmo dia tanto à ecocardiografia quanto à cine RMC com precessão livre no estado estacionário (steady-state free precession, SSFP). Pontos-guia foram utilizados no anel mitral septal e lateral para definir a velocidade de estiramento no início da diástole (E'), enquanto curvas de volume-tempo do modelo 3D foram utilizadas para avaliar os parâmetros de enchimento diastólico. Foram determinadas a correlação entre a RMC 3D e a ecocardiografia, além da acurácia da RMC em classificar o grau de função diastólica.

Resultados:

Ao todo, 102 sujeitos foram incluídos no estudo. A razão E/A pela RMC esteve positivamente associada com a razão E/A obtida pela ecocardiografia (r = 0,71, p < 0,0001). Estiveram positivamente associadas a velocidade de relaxamento diastólico inicial ao Doppler tecidual e a velocidade de relaxamento longitudinal de deslocamento do anel mitral lateral (p = 0,007), bem como a razão entre E/e' por Doppler e E/E' pela RMC (p = 0,01). A velocidade normalizada de pico de enchimento (EM) determinada pela RMC e o tempo de desaceleração (TD) foram capazes de predizer a disfunção diastólica (áreas sob a curva [AUCs] = 0,70 e 0,72, respectivamente). Além disso, a razão E/E' lateral mostrou boa utilidade para a identificação da disfunção diastólica (AUC = 0,80). No geral, a ecocardiografia e a RMC apresentaram excelente concordância interobservador e intraobservador (coeficiente de correlação intraclasse 0,72 - 0,97).

Conclusão:

Uma modelagem 3D de imagens padrões de cine RMC foi capaz de identificar os indivíduos do estudo com função diastólica reduzida e mostrou uma boa reprodutibilidade, sugerindo ter potencial na avaliação rotineira da função diastólica por RMC.

Palavras-chave:
Função Ventricular; Avaliação; Ressonância Magnética; Imagem Tridimensional; Ecocardiografia Tridimensional

Introduction

the prevalence and cost of treatment of heart failure (HF) in the United States are high. In 2008, this condition was estimated to affected 5.3 million adults and was associated with a total spending of 34.8 billion dollars.^1,2 Approximately 50% of the patients were reported to have diastolic HF.^1,2 Diastolic dysfunction is an increasingly recognized component of a variety of diseases of the myocardium,^3,4 and its recognition is necessary for patient management.⁵

Echocardiography is currently used as the standard of reference to evaluate diastolic dysfunction.^6-10 With cardiac magnetic resonance (CMR) imaging, diastolic function is assessed using special pulse sequences such as phase-contrast analysis or myocardial tissue tagging.^5,6,8,11-16 These assessments require additional time and software for acquisition and analysis. As a result, the diastolic assessment with CMR is not routinely applied.^5,17,18 Thus, CMR is in need of a simple and robust method for diastolic function assessment that can be done with routine protocol sequences.

A three-dimensional (3D) model of myocardial function has been developed to assess the myocardial function based on standard steady-state free precession (SSFP) CMR cine images.¹⁹ A model-based analysis of the systolic function is relatively fast (~15 minutes per CMR study) and allows extraction of time-varying function parameters that may characterize the diastolic function.^19-23

Thus, the purpose of this study was to perform an intraindividual analysis to develop and validate a 3D model-based volumetric assessment of diastolic function using CMR imaging and compare the results obtained with this model with those obtained by echocardiography.

Methods

Study population

The study included participants who underwent both echocardiography and CMR between 2008 and 2009 in a substudy of the Multi-Ethnic Study of Atherosclerosis (MESA) at the Johns Hopkins Hospital. Details of the MESA study have been previously described.²⁴ In brief, 1096 participants free of clinically apparent cardiovascular disease and aged 45-84 years, were enrolled at the Baltimore field center at baseline in 2000-2002. A total of 149 consecutive participants were invited to participate in the CMR-echocardiography substudy. Participants were excluded if they had not undergone both studies at the same day; if they had a heart rate variability of more than 15 beats per minute between both studies, severe mitral annular calcification or mitral valve regurgitation; or if the qualitative assessments of the left ventricular (LV) function was impaired by arrhythmias or poor image quality by either modality (Figure 1). The study was approved by the local ethics committee, and all subjects gave a written informed consent for participation.

Figure 1
Flowchart of the study population. Abbreviation: CMR: cardiac magnetic resonance.

Since this study included a correlation between echocardiography and CMR, not all variables were used in the analysis. We will describe the variables that can be acquired by echocardiography and the 3D model-based volumetric assessment of diastolic function using CMR.

Echocardiography

Echocardiograms were obtained by expert sonographers according to the recommendations of the American Society of Echocardiography (ASE).¹⁰ The examinations were reviewed offline by two readers. Readers 1 (A.L.C.A.) and 2 (A.C.A.) had 20 and 5 years of experience, respectively, in reading echocardiograms. Two-dimensional (2D) echocardiograms were recorded using an Aplio scanner (Toshiba Medical Systems Corp, Tochigi, Japan). The images were acquired from an LV apical four-chamber view. Image acquisitions were performed using B-mode harmonic images adjusting transducer frequencies (1.7-3.5 MHz), frame rate (40-80 frames per second), focus, sector width (as narrow as possible), sector depth (minimal), and gain, in order to optimize myocardial image quality. The images were digitally recorded, stored on compact discs, and transferred to a computer terminal for post processing.

Mitral inflow velocities: All Doppler measurements were assessed according to the ASE recommendations.²⁵ From the transmitral recordings, the following measurements were carried out: a) transmitral early peak filling velocity during diastole (early peak filling rate [E]), in centimeters per second; b) (transmitral late peak atrial filling velocity during diastole [peak atrial velocity [A]), in centimeters per second; c) time elapsed between E and the point where the extrapolation of the deceleration slope of the E velocity crosses the zero baseline (deceleration time [DT]), in milliseconds; d) time elapsed between the systolic peak to E (time to peak E [relative TPE]), in milliseconds; e) time elapsed between the systolic peak to A (time to peak A [relative TPA], in milliseconds.

Tissue Doppler measurement of mitral annular velocity: Pulsed wave tissue Doppler imaging (TDI) was performed in the apical views to acquire the mitral annular velocities according to the ASE recommendations.²⁵ The sample volume was placed in the ventricular myocardium immediately adjacent to the mitral annulus in the septal and lateral walls. With this method, the early diastolic myocardial relaxation velocity (e'; cm/s), as the annulus ascends away from the apex, was assessed in this study.

Cardiac magnetic resonance

Cine-CMR images were acquired on a 1.5 T scanner (Avanto, Siemens, Malvern, PA, USA) using a 2D SSFP acquisition in vertical long-axis, horizontal long-axis, and short-axis orientations with the following parameters: TE 1.16 ms, TR 3.2 ms, flip angle 60°, receiver bandwidth ±1220 kHz, FOV 36 cm, slice thickness 8 mm, slice gap 2 mm, acquisition matrix 205 × 256, number of averages = 1, number of frames = 30. The mean reconstructed temporal resolution (R-R interval/number of cardiac phases) was 30.43 ± 5.44 ms.

CMR images were analyzed using a research version of the CIM 6.2 program modified to assess diastolic function (Auckland MRI Research Group, University of Auckland, New Zealand).¹⁹ CMR image analyses were done by two readers accredited by the Auckland MRI Research Group. Readers 1 (M.S.N.) and 2 (E.Y.) had 7 years and 1 year of experience, respectively, in reading CMR.

Time-volume curve: All timing measurements were defined semiautomatically with manual correction with the observer using a slider on the time/rate curve (Figure 2). The following measurements were assessed: a) diastolic volume recovery (DVR), defined as the time from end-systole (ES) to the time at which the volume has filled to 80% of the stroke volume (msec); b) E (mL/sec), the first maximum filling rate detected after ES. Peak E was also divided by the end-diastolic volume (EDV) to generate a normalized peak E filling rate (NE). Additional measurements included: c) relative time to early peak filling rate (RTPE) (msec), the trigger time to peak E from the ES phase; d) A (mL/sec), the second peak filling rate after ES. Peak A was also indexed by EDV to generate a normalized peak A filling rate (NA); e) relative time of atrial peak filling rate (RTPA; msec), the trigger time to peak A from the ES phase; and f) DT (msec), or the time delay of E subtracted from the E wave downslope intersecting the baseline.

Figure 2
Screenshots of different diastolic function examples using the program CIM. A) Normal, B) impaired, and C) reduced. The following measurements were assessed: end-systole (ES), end-diastole (ED), early peak filling rate (E), atrial peak filling rate (A), and deceleration time (DT). All timing measurements were defined semiautomatically with manual correction with the observer using a slider on the time/rate curve.

Guide points at the junction of the LV wall with the septal mitral annulus and at the junction of the LV wall with the lateral mitral annulus in the four-chamber view were used to define g) E' septal, and h) E' lateral, respectively. The ratio between E and E' was also calculated (Figure 3).

Figure 3
Three-dimensional displays of the model fits throughout the cardiac cycle in one R-R interval of 870 ms (outside images) for volume- and derivative-curve assessment (mL/s). Septal and lateral guide points motion can be evaluated through time, calculating the distance between the defined point and the model apex for myocardial longitudinal relaxation rate (mm/s). The endocardial surface is shaded in red and drawn with green lines.

Note that CMR rates are expressed as volume (mL) per unit of time, whereas echocardiographic parameters are expressed as distance (cm) per unit of time. However, CMR-derived E' is expressed as a linear velocity similar to its echocardiographic correlate.

Data and statistical analysis

The diastolic function classification used three echocardiographic parameters recommended by the ASE for this purpose: (1) septal e' < 8 cm/s, (2) lateral e'< 10 cm/s, and (3) ratio between average E and average e' ≥ 10.²⁵ If all three criteria were present, the diastolic function was rated as type II (reduced). If only two criteria were present or one criteria plus LV hypertrophy, the diastolic function was rated as type I (impaired). The LV mass was assessed by echocardiography and divided by the body surface area to define the LV mass index (LVMi). LV hypertrophy was defined as an LVMi > 115 g/m² for men and > 95 g/m² for women, as recommended by the ASE.²⁶

Data are presented as mean ± standard deviation (SD) for continuous variables and as percentage for categorical variables. Multiple comparisons were tested by one-way analysis of variance (ANOVA) with post hoc Bonferroni correction. Fisher's exact test was used to examine the differences between proportions. As the variables were normally distributed, linear regression analysis was performed using Pearson's correlation coefficient (r) and setting echocardiography as the predictor variable and CMR as the dependent variable. We used Bland-Altman to compare variables with the same units. However, in many cases, the CMR surrogates for echocardiographic parameters were represented in different units, so the Bland-Altman analysis was inappropriate.

Receiver operating characteristic (ROC) curve analysis was used to identify the diagnostic performance of CMR in predicting diastolic dysfunction. This was achieved by using the group with reduced diastolic function assessed by echocardiography as the "true positive" surrogate marker for diastolic dysfunction in this population, compared with the group with normal function as the "true negative" (area under the curve [AUC] ≥ 0.5 to < 0.7 = poor fit, AUC ≥ 0.7 to < 0.9 = good fit, and AUC ≥ 0.9 to 1.0 = excellent fit).

Intraobserver and interobserver agreements were assessed using intraclass correlation coefficient (ICC) with a two-way random model (ICC < 0.40 = poor agreement, ICC ≥ 0.40 to 0.75 = fair to good agreement, ICC > 0.75 = excellent agreement).

The statistical analysis was performed using Stata, version 12.0 (StataCorp LP, College Station, Texas, USA). A p value < 0.05 was considered significant.

Results

A total of 102 participants met the inclusion criteria (Figure 1). On echocardiography, the diastolic function was classified as normal in 66 (64.7%) patients, impaired in 21 (20.6%), and reduced in 15 (14.7%) of them. The mean duration of the CMR analysis (systolic and diastolic function) was 18.3 ± 4.5 minutes. Note that the CMR analysis also yields parameters such as LV volume and mass, since the analysis is performed over the full cardiac cycle. The mean duration of the echocardiographic analysis (diastolic function only) was 4.6 ± 0.6 minutes (p < 0.0001 compared with the CMR analysis). Reduced diastolic function was more frequent in diabetic and hypertensive participants. Major variables, e.g., age, gender, body mass index (BMI), systolic blood pressure (SBP), LV mass, EDV, and heart rate, showed no significant variance between groups. The characteristics of the subjects and the clinical data related to their LV function are summarized in Table 1.

Table 1
Population characteristics by diastolic function grades

	Normal n = 66 (64.70%)	Type I n = 21 (20.60%)	Type II n = 15 (14.70%)	p value
Age (years)	66.8 ± 8.9	65.5 ± 7.5	64.4 ± 9.7	0.60
45 to 64 years	24 (36.3)	11 (52.3)	7 (46.6)	0.48*
65 to 84 years	42 (63.6)	10 (47.6)	8 (53.3)	0.48*
Gender (male)	26 (39.3)	7 (33.3)	6 (40.0)	0.91
Race
White, Caucasian	41 (62.0)	11 (55.0)	6 (40.0)	0.25*
Black, African-American	25 (38.0)	10 (45.0)	9 (60.0)	0.25*
Weight (kg)	77.5 ± 15.1	80.3 ± 19.4	80.3 ± 22.2	0.73
Height (cm)	168.0 ± 9.4	166.0 ± 11.2	166.1 ± 9.7	0.65
BMI (kg/m2)	28.0 ± 4.4	29.1 ± 5.7	28.2 ± 7.1	0.71
BSA	1.8 ± 0.2	1.8 ± 0.2	1.8 ± 0.2	0.89
Smoking status
Never	27 (40.9)	9 (42.8)	5 (33.3)	0.37*
Former	33 (50.0)	12 (57.1)	7 (46.6)	0.37*
Current	6 (9.0)	0 (0.0)	3 (20.0)	0.37*
Systolic blood pressure (mmHg)	121.8 ± 18.7	119.8 ± 14.6	121.3 ± 25.9	0.91
Diastolic blood pressure (mmHg)	71.2 ± 11.2	66.3 ± 10.6	69.1 ± 10.7	0.21
Hypertension (%)	33 (50.0)	7 (33.3)	9 (60.0)	0.28*
Any hypertension medication	31 (46.9)	6 (28.5)	9 (60.0)	0.20*
Diabetes (%)	3 (4.5)	2 (9.5)	3 (20.0)	0.11*
Triglycerides (mg/dL)	111.9 ± 60.5	100.3 ± 67.3	101.3 ± 54.8	0.68
LDL cholesterol (mg/dL)	111.5 ± 32.3	109.1 ± 34.2	112.4 ± 42.9	0.95
HDL cholesterol (mg/dL)	58.9 ± 18.4	62.4 ± 24.9	52.6 ± 12.0	0.33
Total cholesterol (mg/dL)	192.8 ± 38.5	191.5 ± 38.8	185.2 ± 52.7	0.81
Metabolic syndrome	21 (31.8)	4 (19.0)	2 (13.0)	0.26*
Echocardiographic measurements
Heart rate (beats/min)	64.8 ± 9.6	65.0 ± 9.4	62.6 ± 5.7	0.66
End-diastolic diameter (mm)	4.4 ± 0.5	4.5 ± 0.5	4.6 ± 0.4	0.44
Diastolic septal thickness (mm)	1.0 ± 0.2	1.0 ± 0.1	0.9 ± 0.1	0.18
Diastolic inferolateral thickness (mm)	0.9 ± 0.1	0.9 ± 0.1	0.9 ± 0.1	0.66
CMR measurements
Heart rate (beats/min)	65.2 ± 10.4	66.4 ± 9.7	61.6 ± 5.5	0.31
Ejection fraction (%)	69.0 ± 7.3	70.7 ± 7.0	70.6 ± 10.2	0.51
End-diastolic volume (mL)	106.8 ± 24.4	110.6 ± 28.7	99.6 ± 22.1	0.43
End-systolic volume (mL)	33.8 ± 13.6	33.5 ± 14.2	28.8 ± 10.6	0.42
LV mass (g)	124.8 ± 34.4	132.5 ± 38.2	121.8 ± 26.2	0.59
Stroke volume (mL)	73.0 ± 15.1	76.5 ± 18.0	73.1 ± 17.5	0.67

• BMI: body mass index; BSA: body surface area; CMR: cardiac magnetic resonance; LV: left ventricular. Note:
• ^*
  Fisher’s exact test was used to compare proportions between diastolic severity grades.

Echocardiographic parameters showed increasing mean values in association with diastolic dysfunction severity (p < 0.05, Table 2). However, A alone showed no statistically significant difference between groups. NE and DT obtained from derivative volume-curves by CMR showed trends toward diastolic dysfunction severity similar to those obtained by echocardiography (p < 0.05). The E/A ratio by CMR was 1.10 ± 0.38 in the normal group, and was lower in the impaired group (1.01 ± 0.26) and higher in the reduced diastolic function group (1.33 ± 0.45, p = 0.03). All other variables showed no difference between groups (Table 2).

Table 2
Diastolic measurements by echocardiography and cardiac magnetic resonance

	Normal n = 66 (64.70%)	Type I n = 21 (20.60%)	Type II n = 15 (14.70%)	p value
Echocardiography
Mitral inflow velocities
E (cm/s)	74.53 ± 16.43	74.89 ± 20.76	87.68 ± 20.94	0.03
DT (ms)	220.14 ± 45.19	247.95 ± 76.70	258.4 ± 77.69	0.03
A (cm/s)	77.66 ± 19.59	75.40 ± 19.55	77.18 ± 23.09	0.90
E/A	0.99 ± 0.24	1.02 ± 0.22	1.23 ± 0.47	0.01
Tissue Doppler velocities
Septal
e’ (cm/s)	9.38 ± 1.69	8.33 ± 2.12	6.00 ± 1.26	<0.0001
E/e’	8.20 ± 2.24	9.23 ± 2.22	15.37 ± 5.91	<0.0001
Lateral
e’(cm/s)	11.61 ± 2.45	8.23 ± 1.68	6.97 ± 1.80	<0.0001
E/e’	6.65 ± 1.82	9.37 ± 3.22	13.36 ± 5.21	<0.0001
Mean
e’ (cm/s)	10.47 ± 1.59	8.27 ± 1.45	6.48 ± 1.44	<0.0001
E/e’	7.25 ± 1.78	9.12 ± 2.06	14.19 ± 5.40	<0.0001
CMR
Volume-curves
E (mL/s)	189.30 ± 66.39	206.30 ± 62.58	213.60 ± 71.67	0.33
NE (s-1)	1.77 ± 0.46	1.89 ± 0.50	2.11 ± 0.43	0.03
DT (ms)	186.61 ± 43.94	211.08 ± 43.75	218.37 ± 42.59	0.01
TPE (ms)	504.86 ± 82.41	493.46 ± 68.75	517.54 ± 37.80	0.63
A (mL/s)	181.13 ± 72.08	211.09 ± 75.17	164.73 ± 43.96	0.11
NA(s-1)	1.70 ± 0.53	1.98 ± 0.76	1.71 ± 0.58	0.16
TPA (ms)	837.27 ± 193.40	861.57 ± 155.17	866.00 ± 115.64	0.78
E/A	1.10 ± 0.38	1.01 ± 0.26	1.33 ± 0.45	0.03
DVR (ms)	535.32 ± 117.96	542.44 ± 122.45	516.08 ± 78.16	0.80
Longitudinal relaxation rate
Septal
E’ (mm/s)	75.35 ± 24.49	66.49 ± 25.31	58.22 ± 24.11	0.03
E/E’ (mL/mm)	2.64 ± 0.96	3.45 ± 1.60	4.65 ± 3.38	0.0002
Lateral
E’ (mm/s)	82.36 ± 26.14	70.88 ± 28.45	61.06 ± 27.73	0.01
E/E’ (mL/mm)	2.40 ± 0.83	3.32 ± 1.80	4.52 ± 3.54	0.0001
Mean
E’ (mm/s)	78.86 ± 24.85	68.69 ± 26.26	59.64 ± 25.45	0.02
E/E’ (mL/mm)	2.50 ± 0.87	3.33 ± 1.53	4.55 ± 3.44	0.0001

• E: early peak filling rate; DT: deceleration time; A: atrial peak filling rate; E/A:E/A ratio; e’: early diastolic myocardial relaxation velocity; E/e’: E/e’ ratio; NE: normalized peak E filling rate; NA: normalized peak A filling rate; DVR: diastolic volume recovery; E’: early longitudinal relaxation rate; CMR: cardiac magnetic resonance, TPE: time to peak E; TPA: time to peak A.

Tissue Doppler velocities by echocardiography assessed e' and the E/e' ratio. In all regions (septal and lateral mitral annulus), e' showed significantly decreased mean values in the normal diastolic function group (e' lateral = 11.6 ± 2.4 cm/s) and in the reduced diastolic function group (e' lateral = 6.9 ± 1.8 cm/s, p < 0.05). Also, E/e' increased from the group with a normal diastolic function to the one with reduced diastolic function (6.65 ± 1.8 and 13.3 ± 5.2, respectively, p < 0.0001). Compared with CMR, E' and E/E' showed similar trends toward worse diastolic function for both septal and lateral walls (p < 0.05 and p < 0.001, respectively) (Table 2).

Table 3 highlights the associations between the diastolic function measured by echocardiography and CMR. E/A ratios on echocardiography were positively associated with E/A ratios on CMR (r = 0.71, p < 0.0001). The 95% limits of agreement between the two methods were -0.45% to +0.62%. A small bias (0.081%) toward a higher E/A ratio by CMR was detected (Figure 4).

Table 3
Associations between measures of diastolic function by echocardiography and cardiac magnetic resonance (n = 102)

Echocardiography	CMR	Pearson’s correlation coefficient (r)	p value
Mitral inflow velocities	Volume-curves
E (cm/s)	E (mL/s)	0.06	0.51
E (cm/s)	NE (s-1)	0.1	0.18
A (cm/s)	A (mL/s)	0.22	0.01
A (cm/s)	NA (s-1)	0.28	0.003
E/A	E/A	0.71	< 0.0001
Tissue Doppler	Longitudinal relaxation rate
Septal	Septal
e’ (cm/s)	E’ (mm/s)	0.11	0.26
E/e’	E/E’ (mL/mm)	0.11	0.30
Lateral	Lateral
e’ (cm/s)	E’ (mm/s)	0.26	0.007
E/e’	E/E’ (mL/mm)	0.24	0.01
Mean	Mean
e’ (cm/s)	E’ (mm/s)	0.22	0.02
E/e’	E/E’ (mL/mm)	0.17	0.07

• CMR: cardiac magnetic resonance; early peak filling rate; A: atrial peak filling rate; e’: early diastolic myocardial relaxation velocity; NE: normalized peak E filling rate; NA: normalized peak A filling rate; E/A: E/A ratio; E’: early longitudinal relaxation rate. Echocardiography corresponds to Doppler echocardiography.

Figure 4
Results obtained using cardiac magnetic resonance (CMR) three-dimensional volume-curve and echocardiography Doppler mitral valve inflow. The ratio between the early peak filling (E) and atrial peak filling rate (A) using velocity (cm/s) by echocardiography and flow (mL/s) by CMR. (A) Linear regression and Pearson's correlation; (B) Bland-Altman analysis.

Values of e' by tissue Doppler and E' for the lateral mitral annulus displacement were positively correlated (r = 0.26, p = 0.007), as were E/e' by CMR and echocardiography (r = 0.24, p = 0.01). However, both septal measurements were not correlated (p > 0.05).

Prediction of reduced diastolic function by cardiac magnetic resonance

Table 4 shows the ROC curve analysis for reduced diastolic function for all CMR parameters. CMR-determined NE and DT were able to predict diastolic dysfunction (AUCs = 0.70 and 0.72, respectively). In addition, the lateral E/E' ratio appeared to be useful in the classification of diastolic dysfunction (AUC = 0.80) (Table 4).

Table 4
Prediction of reduced diastolic function by cardiac magnetic resonance (n = 81)

CMR	Area under the ROC curve	p value
Volume-curves
E (mL/s)	0.60	0.21
NE (s-1)	0.70	0.008
DT (ms)	0.72	0.01
A (mL/s)	0.53	0.37
NA (s-1)	0.48	0.92
DVR (ms)	0.51	0.57
E/A	0.66	0.05
Longitudinal relaxation rate
Septal
E’ (mm/s)	0.67	0.01
E/E’ (mL/mm)	0.76	0.0003
Lateral
E’ (mm/s)	0.70	0.0004
E/E’ (mL/mm)	0.80	< 0.0001
Mean
E’ (mm/s)	0.69	0.006
E/E’ (mL/mm)	0.78	0.0001

• CMR: cardiac magnetic resonance; ROC: receiver operating characteristic; E: early peak filling rate; NE: normalized peak E filling rate; DT: deceleration time; A: atrial peak filling rate; NA: normalized peak A filling rate; DVR: diastolic volume recovery; E/A: E/A ratio; E’: early longitudinal relaxation rate; E/E’: E/E’ ratio.

Diastolic time periods and cardiac cycle duration

No significant differences were detected in relative TPE and RTPA values obtained by CMR compared with those obtained by echocardiography (mean RTPA: 183.3 ± 47.32 ms versus 181.5 ± 27.45 ms, respectively, p = 0.90; mean TPE: 544.32 ± 145.62 ms versus 550.77 ± 196.19 ms, respectively, p = 0.91). The cardiac cycle duration (R-R interval) was also not significantly different by CMR versus echocardiography (mean 943.65 ± 135.11 ms versus 944.77 ± 135.42 ms, respectively, p = 0.95).

Interobserver and intraobserver agreements

Overall, echocardiography and CMR interobserver and intraobserver agreements were excellent (Table 5). The mean ICC for measurements by echocardiography was excellent (0.89) and slightly higher than those obtained by CMR (0.86).

Table 5
Intraobserver and interobserver agreement (n = 20)

	Intraclass correlation coefficient (ICC)	Bias	95% limits of agreement
Echocardiography R1 versus R2
Mitral inflow velocities
E (cm/s)	0.93	-1.66	-11.70 to 8.36
DT (ms)	0.84	9.84	-38.67 to 58.36
A (cm/s)	0.95	-1.12	-14.84 to 12.59
Tissue Doppler velocities
Septal
e’ (cm/s)	0.85	0.42	-1.84 to 2.68
Lateral
e’ (cm/s)	0.89	-0.37	-1.76 to 2.49
Echocardiography R1 versus R1
Mitral inflow velocities
E (cm/s)	0.95	-1.39	-9.22 to 6.44
DT (ms)	0.72	6.77	-62.40 to 75.96
A (cm/s)	0.96	-0.22	-12.85 to 12.39
Tissue Doppler velocities
Septal
e’ (cm/s)	0.89	0.28	-1.78 to 2.35
Lateral
e’ (cm/s)	0.92	-0.59	-2.06 to 0.86
CMR R1 versus R2
Volume-curves
E (mL/s)	0.84	2.54	-79.77 to 84.86
DT (ms)	0.77	-21.52	-81.75 to 38.70
A (mL/s)	0.82	22.89	-51.20 to 97.00
Longitudinal relaxation rate
Septal
E’ (mm/s)	0.75	-4.90	-32.59 to 22.63
Lateral
E’ (mm/s)	0.89	-5.48	-25.24 to 14.27
CMR R1 versus R1
Volume-curves
E (mL/s)	0.97	-1.36	-33.73 to 31.00
DT (ms)	0.84	12.93	-28.46 to 28.46
A (mL/s)	0.96	-15.51	-64.82 to 33.79
Longitudinal relaxation rate
Septal
E’ (mm/s)	0.85	-2.93	-23.91 to 18.03
Lateral
E’ (mm/s)	0.94	-4.11	-18.93 to 10.70

• Note: R1: reader 1 and R2: reader 2. E: early peak filling rate; DT: deceleration time; A: atrial peak filling rate; e’: early diastolic myocardial relaxation velocity; E’: early longitudinal relaxation rate.

Discussion

The purpose of this study was to evaluate the role of cine CMR for diastolic function assessment and compare values obtained with this method with those obtained with echocardiography. Using a relatively fast and reproducible method, CMR-derived parameters were shown to be comparable to those obtained by echocardiography, with good correlations. Importantly, this study demonstrated that CMR was capable of identifying diastolic dysfunction in most patients with diastolic dysfunction detected by echocardiography. This suggests a role for CMR in the assessment of LV diastolic function in the general population.

Echocardiography has long been used to evaluate diastolic dysfunction. The combination of mitral inflow velocity curves and tissue Doppler velocities of the mitral annulus are known to provide better estimates of LV filling pressures than other methods.²⁷ Although routinely reported by echocardiography, diastolic function by CMR is usually not routinely assessed due to the requirement of additional phase contrast or tagged sequences, as well as separate post processing. Automated segmentation of LV volumes for all temporal phases holds the potential to rapidly assess diastolic filling patterns;²⁸ however, this method alone only provides partial information regarding the diastolic physiology needed to differentiate all degrees of diastolic dysfunction severity.

Recently, CMR software innovations^19,29,30 have allowed the assessment of similar parameters using SSFP cine CMR with 3D post processing. HF with preserved ejection fraction is increasing in incidence and has a high clinical relevance,⁸ although a clear consensus for its diagnosis has yet to be established.³¹ In this study, we decided to follow the ASE recommendations²⁵ to delineate normal versus reduced diastolic function groups.

CMR is considered a reference standard for ventricular systolic function, including the analysis of regional wall motion, mass, and volumes, and estimation of ejection fraction.³² The assessment of diastolic function by CMR is usually not routinely performed in our clinical practice. CMR diastolic assessment typically requires an increased scan time for image acquisition (e.g., additional phase-contrast sequences), as well as a tedious imaging post-processing analysis. Automated segmentation of LV volumes for all temporal phases holds the potential to assess rapidly diastolic filling patterns;²⁸ however, this method relies on a sequential cross-sectional analysis that only provides partial information regarding the diastolic physiology needed to differentiate all degrees of diastolic dysfunction severity.

In our study, we were able to overcome several limitations of CMR diastolic function analysis using a new 3D method with an average analysis time of fewer than 20 minutes with no need to add more sequences in our routine protocol. In our experience, this analysis time is comparable to that obtained for full 3D volumetric assessment of systole alone. For CMR, it appears that E/E', NE, DT, and E/A were the most useful parameters obtained from time-volume curves. When the longitudinal shortening was measured, both septal and lateral measurements were able to categorize diastolic dysfunction. However, E' at the lateral wall was shown to be more reproducible and easily measured by CMR. By comparison, septal E' had lower reader reproducibility.

This study had several limitations. Echocardiography was used as the reference standard, but reader variability and diastolic classification are known to be imperfect with this method.³¹ The correlation between LV mitral valve inflow velocities and time-volume curves from CMR represents different physiological processes. The time-volume curves from CMR should not be adversely affected by mitral valve disease or angle of acquisition.^9,25 In this study population, hemodynamic data were not available. In addition, our study population did not include subjects with restrictive cardiomyopathy. CMR time-volume curves represent the average of several cardiac cycles, whereas echocardiography shows peak values for each cardiac cycle. Finally, although a good correlation of echocardiographic and CMR data appears to be present, outcome data is still needed to validate further the CMR approach.

Conclusion

The 3D CMR method was relatively fast, reproducible, and successfully applied to routine SSFP cine CMR data. CMR was able to identify most patients with reduced diastolic function identified by echocardiography. This suggests a role for CMR in the assessment of LV diastolic function in the general population and in patients with mild and moderate diastolic dysfunction.

• Sources of Funding
  This study was supported by a National Institutes of Health (NIH) intramural research program, a National Heart, Lung, and Blood Institute grant (RO1-HL66075-01), and by the MESA study contracts NO1-HC-9808, NO1-HC-95168, and NO1-HC-95169.
• Study Association
  This article is part of the thesis of Post Doctoral submitted by Marcelo Souto Nacif, from Johns Hopkins School of Medicine and National Institutes of Health Clinical Center.

Acknowledge

The authors thank the other investigators, the staff, and the participants of the MESA study for their valuable contributions. A full list of participating MESA investigators and institutions can be found at http://www.mesa-nhlbi.org.

References

• ¹
  Rosamond W, Flegal K, Furie K, Go A, Greenlund K, Haase N, et al; American Heart Association Statistics Committee and Stroke Statistics Subcommittee. Heart disease and stroke statistics--2008 update: a report from the American Heart Association Statistics Committee and Stroke Statistics Subcommittee. Circulation. 2008;117(4):e25-146. Erratum in: Circulation. 2010;122(1):e10.
• ²
  Sharma V, Zehtabchi S, Rojas N, Birkhahn R. Ethnic variations in quality of life and depressive symptoms among black Americans with acute decompensated heart failure. J Natl Med Assoc. 2009;101(10):985-91.
• ³
  Sanders D, Dudley M, Groban L. Diastolic dysfunction, cardiovascular aging, and the anesthesiologist. Anesthesiol Clin. 2009;27(3):497-517.
• ⁴
  Schertel ER. Assessment of left-ventricular function. Thorac Cardiovasc Surg. 1998;46 Suppl 2:248-54.
• ⁵
  Bollache E, Redheuil A, Clement-Guinaudeau S, Defrance C, Perdrix L, Ladouceur M, et al. Automated left ventricular diastolic function evaluation from phase-contrast cardiovascular magnetic resonance and comparison with Doppler echocardiography. J Cardiovasc Magn Reson. 2010;12:63.
• ⁶
  Daneshvar D, Wei J, Tolstrup K, Thomson LE, Shufelt C, Merz CN. Diastolic dysfunction: improved understanding using emerging imaging techniques. Am Heart J. 2010;160(3):394-404.
• ⁷
  Paelinck BP, de Roos A, Bax JJ, Bosmans JM, van Der Geest RJ, Dhondt D, et al. Feasibility of tissue magnetic resonance imaging: a pilot study in comparison with tissue Doppler imaging and invasive measurement. J Am Coll Cardiol. 2005;45(7):1109-16. Erratum J Am Coll Cardiol. 2005;45(10):1737.
• ⁸
  Rubinshtein R, Glockner JF, Feng D, Araoz PA, Kirsch J, Syed IS, et al. Comparison of magnetic resonance imaging versus Doppler echocardiography for the evaluation of left ventricular diastolic function in patients with cardiac amyloidosis. Am J Cardiol. 2009;103(5):718-23.
• ⁹
  Evangelista A, Flachskampf F, Lancellotti P, Badano L, Aguilar R, Monaghan M, et al; European Association of Echocardiography. European Association of Echocardiography recommendations for standardization of performance, digital storage and reporting of echocardiographic studies. Eur J Echocardiogr. 2008;9(4):438-48.
• ¹⁰
  Lang RM, Bierig M, Devereux RB, Flachskampf FA, Foster E, Pellikka PA, et al. Recommendations for chamber quantification. Eur J Echocardiogr. 2006;7(2):79-108.
• ¹¹
  Marsan NA, Westenberg JJ, Tops LF, Ypenburg C, Holman ER, Reiber JH, et al. Comparison between tissue Doppler imaging and velocity-encoded magnetic resonance imaging for measurement of myocardial velocities, assessment of left ventricular dyssynchrony, and estimation of left ventricular filling pressures in patients with ischemic cardiomyopathy. Am J Cardiol. 2008;102(10):1366-72.
• ¹²
  Gatehouse PD, Rolf MP, Graves MJ, Hofman MB, Totman J, Werner B, et al. Flow measurement by cardiovascular magnetic resonance: a multi-centre multi-vendor study of background phase offset errors that can compromise the accuracy of derived regurgitant or shunt flow measurements. J Cardiovasc Magn Reson. 2010;12:5.
• ¹³
  Codreanu I, Robson MD, Golding SJ, Jung BA, Clarke K, Holloway CJ. Longitudinally and circumferentially directed movements of the left ventricle studied by cardiovascular magnetic resonance phase contrast velocity mapping. J Cardiovasc Magn Reson. 2010;12:48.
• ¹⁴
  Shehata ML, Cheng S, Osman NF, Bluemke DA, Lima JA. Myocardial tissue tagging with cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2009;11:55.
• ¹⁵
  Chirinos JA, Segers P. Noninvasive evaluation of left ventricular afterload: part 1: pressure and flow measurements and basic principles of wave conduction and reflection. Hypertension. 2010;56(4):555-62.
• ¹⁶
  Chirinos JA, Segers P. Noninvasive evaluation of left ventricular afterload: part 2: arterial pressure-flow and pressure-volume relations in humans. Hypertension. 2010;56(4):563-70.
• ¹⁷
  Hartiala JJ, Mostbeck GH, Foster E, Fujita N, Dulce MC, Chazouilleres AF, et al. Velocity-encoded cine MRI in the evaluation of left ventricular diastolic function: measurement of mitral valve and pulmonary vein flow velocities and flow volume across the mitral valve. Am Heart J. 1993;125(4):1054-66.
• ¹⁸
  Soldo SJ, Norris SL, Gober JR, Haywood LJ, Colletti PM, Terk M. MRI-derived ventricular volume curves for the assessment of left ventricular function. Magn Reson Imaging. 1994;12(5):711-7.
• ¹⁹
  Young AA, Cowan BR, Thrupp SF, Hedley WJ, Dell'Italia LJ. Left ventricular mass and volume: fast calculation with guide-point modeling on MR images. Radiology. 2000;216(2):597-602.
• ²⁰
  Hung J, Francois C, Nelson NA, Young A, Cowan BR, Jerecic R, et al. Cardiac image modeling tool for quantitative analysis of global and regional cardiac wall motion. Invest Radiol. 2009;44(5):271-8.
• ²¹
  Cowan BR, Young AA, Anderson C, Doughty RN, Krittayaphong R, Lonn E, et al; ONTARGET Investigators. Left ventricular mass and volume with telmisartan, ramipril, or combination in patients with previous atherosclerotic events or with diabetes mellitus (from the ONgoing Telmisartan Alone and in Combination With Ramipril Global Endpoint Trial [ONTARGET]). Am J Cardiol. 2009;104(11):1484-9.
• ²²
  Boudoulas H. Systolic time intervals. Eur Heart J. 1990;11 Suppl I:93-104.
• ²³
  Nasir K, Katz R, Mao S, Takasu J, Bomma C, Lima JA, et al. Comparison of left ventricular size by computed tomography with magnetic resonance imaging measures of left ventricle mass and volumes: the multi-ethnic study of atherosclerosis. J Cardiovasc Comput Tomogr. 2008;2(3):141-8.
• ²⁴
  Bild DE, Bluemke DA, Burke GL, Detrano R, Diez Roux AV, Folsom AR, et al. Multi-ethnic study of atherosclerosis: objectives and design. Am J Epidemiol. 2002;156(9):871-81.
• ²⁵
  Nagueh SF, Appleton CP, Gillebert TC, Marino PN, Oh JK, Smiseth OA, et al. Recommendations for the evaluation of left ventricular diastolic function by echocardiography. J Am Soc Echocardiogr. 2009;22(2):107-33.
• ²⁶
  Lang RM, Bierig M, Devereux RB, Flachskampf FA, Foster E, Pellikka PA, et al; Chamber Quantification Writing Group; American Society of Echocardiography's Guidelines and Standards Committee; European Association of Echocardiography. Recommendations for chamber quantification: a report from the American Society of Echocardiography's Guidelines and Standards Committee and the Chamber Quantification Writing Group, developed in conjunction with the European Association of Echocardiography, a branch of the European Society of Cardiology. J Am Soc Echocardiogr. 2005;18(12):1440-63.
• ²⁷
  Ommen SR, Nishimura RA, Appleton CP, Miller FA, Oh JK, Redfield MM, et al. Clinical utility of Doppler echocardiography and tissue Doppler imaging in the estimation of left ventricular filling pressures: a comparative simultaneous Doppler-catheterization study. Circulation. 2000;102(15):1788-94.
• ²⁸
  Kawaji K, Codella NC, Prince MR, Chu CW, Shakoor A, LaBounty TM, et al. Automated segmentation of routine clinical cardiac magnetic resonance imaging for assessment of left ventricular diastolic dysfunction. Circ Cardiovasc Imaging. 2009;2(6):476-84.
• ²⁹
  Fonseca CG, Dissanayake AM, Doughty RN, Whalley GA, Gamble GD, Cowan BR, et al. Three-dimensional assessment of left ventricular systolic strain in patients with type 2 diabetes mellitus, diastolic dysfunction, and normal ejection fraction. Am J Cardiol. 2004;94(11):1391-5.
• ³⁰
  Fonseca CG, Oxenham HC, Cowan BR, Occleshaw CJ, Young AA. Aging alters patterns of regional nonuniformity in LV strain relaxation: a 3-D MR tissue tagging study. Am J Physiol Heart Circ Physiol. 2003;285(5):H621-30.
• ³¹
  Unzek S, Popovic ZB, Marwick TH; Diastolic Guidelines Concordance Investigators. Effect of recommendations on interobserver consistency of diastolic function evaluation. JACC Cardiovasc Imaging. 2011;4(5):460-7.
• ³²
  Hundley WG, Bluemke DA, Finn JP, Flamm SD, Fogel MA, Friedrich MG, et al; American College of Cardiology Foundation Task Force on Expert Consensus Documents. ACCF/ACR/AHA/NASCI/SCMR 2010 expert consensus document on cardiovascular magnetic resonance: a report of the American College of Cardiology Foundation Task Force on Expert Consensus Documents. J Am Coll Cardiol. 2010;55(23):2614-62.

Publication Dates

• Publication in this collection
  25 May 2017
• Date of issue
  2017

History

• Received
  19 Apr 2016
• Reviewed
  19 Dec 2016
• Accepted
  19 Dec 2016