Background
Patients with functional single ventricles after the Fontan procedure have abnormal cardiac mechanics. The aims of this study were to determine factors that influence diastolic function and to describe associations of diastolic function with current clinical status.
Methods
Echocardiograms were obtained as part of the Pediatric Heart Network Fontan Cross-Sectional Study. Diastolic function grade (DFG) was assessed as normal (grade 0), impaired relaxation (grade 1), pseudonymization (grade 2), or restrictive (grade 3). Studies were also classified dichotomously (restrictive pattern present or absent). Relationships between DFG and pre-Fontan variables (e.g., ventricular morphology, age at Fontan, history of volume-unloading surgery) and current status (e.g., systolic function, valvar regurgitation, exercise performance) were explored.
Results
DFG was calculable in 326 of 546 subjects (60%) (mean age, 11.7 ± 3.3 years). Overall, 32% of patients had grade 0, 9% grade 1, 37% grade 2, and 22% grade 3 diastolic function. Although there was no association between ventricular morphology and DFG, there was an association between ventricular morphology and E′, which was lowest in those with right ventricular morphology ( P < .001); this association remained significant when using Z scores adjusted for age ( P < .001). DFG was associated with achieving maximal effort on exercise testing ( P = .004); the majority (64%) of those not achieving maximal effort had DFG 2 or 3. No additional significant associations of DFG with laboratory or clinical measures were identified.
Conclusions
Assessment of diastolic function by current algorithms results in a high percentage of patients with abnormal DFG, but few clinically or statistically significant associations were found. This may imply a lack of impact of abnormal diastolic function on clinical outcomes in this cohort, or it may indicate that the methodology may not be applicable to pediatric patients with functional single ventricles.
Although surgical palliation of patients with functional single ventricles by the Fontan procedure has resulted in improved survival during childhood, long-term morbidity and mortality remain major concerns. In particular, the relationships of pre-Fontan variables to diastolic function and of diastolic function to long-term morbidity and mortality in this cohort are still largely unknown. In the National Heart, Lung, and Blood Institute–sponsored Pediatric Heart Network Fontan Cross-Sectional Study of 546 Fontan survivors, the majority of patients (73%) had normal ejection fractions (EFs) by echocardiography, but only 28% had normal indices of diastolic function. Although patients with single-ventricle physiology have been described as having diastolic dysfunction with both reduced compliance and impaired relaxation, associations with anatomic, clinical, and historical factors have not been described. Additionally, because an acutely increased mass-to-volume ratio immediately after a Fontan operation has been postulated by some to be detrimental to diastolic function, we sought to assess correlations between diastolic function and current mass-to-volume ratio and to determine whether mass-to-volume ratio correlates with clinical characteristics.
Methods
The National Heart, Lung, and Blood Institute–sponsored Pediatric Heart Network Fontan Cross-Sectional Study characterized the health status of 546 Fontan survivors aged 6 to 18 years, enrolled by seven clinical centers in North America. Prospective data collection included two-dimensional and Doppler echocardiography, exercise testing, health status questionnaires, and laboratory tests. Full details of the study have been published. Written informed consent and assent were obtained according to local guidelines, and the study protocol was approved by the institutional review board or research ethics board of each center.
Echocardiograms were obtained according to a predetermined study protocol. None were obtained under sedation. The echocardiograms were submitted to the data-coordinating center for initial quality control and then forwarded to the echocardiography core laboratory for analysis. Ventricular morphology was characterized as left ventricular (LV) dominant (e.g., tricuspid atresia), right ventricular (RV) dominant (e.g., hypoplastic left heart syndrome), or mixed (e.g., unbalanced atrioventricular canal defect). A single representative beat was selected for all systolic and diastolic measures. To determine ventricular size and EF, the ventricle was analyzed from the apical transverse plane (ventricular long-axis) and the parasternal short-axis views. End-diastolic volume (EDV) and end-systolic volumes (ESV) were calculated using the biplane modified Simpson algorithm. Percentage EF was then determined as [(EDV − ESV)/EDV] × 100. Ventricular mass was calculated as myocardial EDV (epicardial volume − endocardial volume) × myocardial density (1.05 g/mL).
Assessment of diastolic function typically involves a multiparametric approach, predominantly using spectral Doppler velocities and tissue Doppler velocities. Lester et al . described a practical approach to diastolic function assessment that involved assessments of mitral inflow characteristics, tissue Doppler early diastolic velocities, atrial size, pulmonary vein flow, isovolumics relaxation time, and mitral inflow propagation velocity. Some of these parameters were not assessable in single-ventricle anatomic variants, particularly atrial size, because of the markedly abnormal and variable geometries of the atrial chambers in the various iterations of single-ventricle circulations.
Therefore, our diastolic assessment was based on the following measurements: Doppler parameters of atrioventricular valve (AVV) inflow (ratio of early mitral inflow velocity [E] to atrial wave inflow [A] velocity [E/A ratio] and deceleration time [DT] of early inflow), tissue Doppler assessment of peak AVV annular velocities in early diastole (E′), and mitral inflow propagation velocity using M-mode color Doppler across the AVV inflow. For E′ values, we used the average of the velocities of the two walls bounding the systemic ventricle (septum and AVV annulus for single ventricles, right and left lateral walls of unbalanced atrioventricular canal systemic ventricle). Two grading systems were used. First, a diastolic function grade (DFG) was assigned to each study on the basis of an algorithm proposed by others. For application to a pediatric cohort, an E/E′ Z score of 3 was used, the equivalent of an absolute value of E/E′ ratio of 10 in adults:
grade 0: normal = [(1 < E/A ≤ 2) and (DT ≥ 140 msec) and (E/E′ Z ≤ 3)];
grade 1: impaired relaxation = [E/A ≤ 1];
grade 2: pseudonormalization = {(1 < E/A ≤ 2) and [(DT < 140 msec) or (E/E′ Z > 3) or (mitral inflow propagation velocity < 55 cm/sec)]}; and
grade 3: restrictive = [E/A > 2].
Second, studies were classified as having restrictive pattern present or absent; a restrictive pattern was considered to be present if E/A > 2 or if 1 < E/A < 2 and DT < 140 msec. Because determination of both the DFG and presence or absence of a restrictive pattern required E/A ratios, 84 studies with E-wave and A-wave fusion (either partially or fully, as demonstrated in Figure 1 ) were excluded from analysis.
Historical factors were obtained from medical chart review, including pre-Fontan catheterization variables, history of coarctation requiring intervention at the time of the pre-Fontan catheterization, history of superior cavopulmonary anastomosis (e.g., bidirectional Glenn shunt, hemi-Fontan), and type of Fontan procedure (e.g., lateral tunnel connection, direct right atrial–to–pulmonary artery connection). Clinical data obtained at the time of echocardiography included B-type natriuretic peptide level, Child Health Questionnaire (CHQ) parental report physical and psychosocial summary scores, and medications, including history of and current angiotensin-converting enzyme (ACE) inhibitor use. Subjects also underwent exercise stress testing, and the following parameters were recorded: peak oxygen consumption (V o 2 ), ventilatory anaerobic threshold, percentage predicted peak V o 2 , percentage predicted ventilatory anaerobic threshold, maximal oxygen pulse/body surface area (BSA), maximal work rate, percentage predicted maximal oxygen pulse, percentage predicted maximal work rate, and percentage predicted maximal heart rate.
Statistical Analysis
Descriptive statistics are reported as mean ± SD or as median (interquartile range). Logistic regression was used to model the presence versus absence of a restrictive pattern as a function of age, and a generalized additive model was used to assess potential nonlinearity between age and restrictive pattern. Spearman correlation coefficients were used to examine the association between continuous diastolic function measures and age. Cumulative logistic regression was used to model the relationship between DFG and age. An interaction was fit between age and ventricular morphology to determine whether associations between DFG and age were dependent on anatomy. Similarly, tests of interaction were fit to assess whether associations of DFG and other clinical measurements depended on ventricular morphology. We used a χ 2 test to assess the association of DFG and categorical patient and clinical variables. Analysis of variance and Kruskal-Wallis tests were used to compare parameters of systolic function (EF, ESV, EDV, and mass/volume Z scores ), current AVV or semilunar valve regurgitation, current use of ACE inhibition, exercise performance, and functional status by DFG. Multivariate modeling was not conducted, because few variables were associated with DFG. Analysis of variance was used to compare the distributions of E/E′ ratio and other continuous diastolic function parameters and their Z scores indexed to age by age group and by ventricular morphology. To determine whether mass-to-volume ratio differed according to DFG and ACE inhibitor use, we used the Kruskal-Wallis test. A P value < .05 was considered to indicate statistical significance. All analyses were conducted in SAS version 9.2 (SAS Institute, Cary, NC) and S-Plus (Insightful, Seattle, WA).
Results
A summary of the parameters by DFG, with Z scores where applicable, is included in Table 1 , and baseline demographic characteristics for the cohort are presented in Table 2 . Of the 546 subjects who underwent echocardiography, 326 (60%) had echocardiograms adequate to assign DFGs, and in 343 (63%), it was possible to determine the presence or absence of a restrictive pattern. The 220 subjects for whom DFGs could not be calculated were similar to those with calculated grades with respect to time since Fontan procedure, gender, and ventricular morphology but differed slightly with regard to race (84% vs 77% white). The distribution of DFG in the overall cohort was as follows: grade 0, 32%; grade 1, 9%; grade 2, 37%; and grade 3, 22%.
Variable | Normal (grade 0) | Impaired relaxation (grade 1) | Pseudonormalization (grade 2) | Restrictive (grade 3) |
---|---|---|---|---|
n | 104 | 30 | 120 | 72 |
E/E′ ratio | 7.3 ± 2.2 | 6.8 ± 3.0 | 9.6 ± 5.0 | 9.0 ± 3.7 |
Z score | 0.6 ± 1.1 | 0.4 ± 1.4 | 1.7 ± 2.4 | 1.4 ± 1.8 |
E/A ratio | 1.5 ± 0.3 | 0.9 ± 0.1 | 1.4 ± 0.3 | 2.5 ± 0.4 |
Z score | −1.2 ± 0.4 | −2.1 ± 0.2 | −1.3 ± 0.5 | 0.4 ± 0.6 |
DT (msec) | 204 ± 56 | 125 ± 34 | 120 ± 39 | 193 ± 87 |
Z score | 1.6 ± 1.6 | −0.7 ± 1.0 | −0.8 ± 1.2 | 1.3 ± 2.6 |
V fp (cm/sec) ∗ | 63 ± 22 | 67 ± 17 | 64 ± 20 | 64 ± 17 |
Z score | −0.3 ± 1.0 | −0.1 ± 0.7 | −0.2 ± 0.9 | −0.2 ± 0.8 |
EF (%) † | 59 ± 10 | 56 ± 12 | 61 ± 11 | 60 ± 10 |
Z score | −1.0 ± 2.3 | −1.7 ± 2.7 | −0.6 ± 2.2 | −0.7 ± 2.1 |
Mass/volume ‡ | 1.1 ± 0.3 | 1.3 ± 0.6 | 1.2 ± 0.4 | 1.2 ± 0.4 |
Z score | 2.1 ± 2.1 | 3.2 ± 5.2 | 2.5 ± 3.4 | 2.9 ± 3.0 |
∗ V fp sample sizes by grade were 39, 7, 57, and 38.
† EF sample sizes by grade were 84, 25, 94, and 61.
‡ Mass/volume sample sizes by grade were 80, 25, 92, and 61.
Variable | No assignable DFG | All with assigned DFG | P 1 | Normal (grade 0) | Impaired relaxation (Grade 1) | Pseudonormalization (grade 2) | Restrictive (grade 3) | P 2 |
---|---|---|---|---|---|---|---|---|
n | 220 | 326 | 104 | 30 | 120 | 72 | ||
Time since Fontan (y) | 8.6 ± 3.6 | 8.5 ± 3.3 | .83 | 8.9 ± 3.3 | 7.8 ± 3.8 | 8.2 ± 3.2 | 8.9 ± 3.2 | .18 |
Age at echocardiography (y) | 12.1 ± 3.5 | 11.7 ± 3.3 | .31 | 12.3 ± 3.5 | 11.4 ± 3.4 | 11.2 ± 3.1 | 12.0 ± 3.4 | .09 |
Heart rate (beats/min) ∗ | 79 ± 15 | 71 ± 13 | <.001 | 68 ± 11 | 80 ± 9 | 77 ± 13 | 64 ± 12 | <.001 |
Male | 57% | 62% | .24 | 61% | 50% | 66% | 64% | .43 |
Ventricular type | .75 | .30 | ||||||
LV | 48% | 49% | 53% | 47% | 41% | 58% | ||
RV | 35% | 33% | 28% | 37% | 39% | 26% | ||
Mixed | 17% | 18% | 19% | 17% | 20% | 15% |
∗ For heart rate, n = 160 without DFGs, n = 277 with DFGs, n = 88 in DFG 0, n = 25 in DFG 1, n = 100 in DFG 2, and n = 64 in DFG 3.
Age-Related Patterns of Diastolic Function
DFG was not associated with time since Fontan or age at echocardiography ( Table 2 ). There were few associations with age among the individual variables used to determine DFG. Mitral valve DT was positively correlated with age (Spearman R = 0.15, P = .007), but this was not significant after adjustment for heart rate ( P = .14). None of the other variables had statistically significant associations with age after adjustment for heart rate.
DFG Associations with Pre-Fontan Variables
Ventricular Morphology
E′ was lowest in subjects with RV morphology (RV morphology, 8.0 ± 2.8; LV morphology, 10.0 ± 3.2; mixed morphology, 10.4 ± 3.9; P < .001) ( Table 3 ), and this finding remained statistically significant when using Z scores for age ( P < .001). The E/E′ ratio was higher in the RV group than in the other groups, because of the lower E′ ( Figure 2 ).
LV | RV | Mixed | P | ||||
---|---|---|---|---|---|---|---|
n | Mean ± SD | n | Mean ± SD | n | Mean ± SD | ||
E | 169 | 0.7 ± 0.2 | 108 | 0.8 ± 0.2 | 66 | 0.7 ± 0.1 | <.001 |
Z score ∗ | −1.4 ± 1.1 | −0.9 ± 1.1 | −1.3 ± 0.9 | <.001 | |||
E′ | 222 | 10.0 ± 3.2 | 155 | 8.0 ± 2.8 | 75 | 10.4 ± 3.9 | <.001 |
Z score ∗ | −1.3 ± 0.7 | −1.8 ± 0.6 | −1.2 ± 0.9 | <.001 | |||
E/E′ ratio | 147 | 7.2 ± 2.3 | 93 | 10.9 ± 5.1 | 56 | 7.4 ± 3.0 | <.001 |
Z score ∗ | 0.6 ± 1.1 | 2.4 ± 2.4 | 0.6 ± 1.4 | <.001 |
Age at Fontan Surgery
We stratified the cohort by age at Fontan surgery as <2, 2 to <3, 3 to <4, and ≥4 years. DFG did not differ significantly among the groups ( P = .50); nor did E/E′ Z score ( P = .45) ( Table 4 ).
Ventricular type ( n ) | Age at Fontan (y) | P ∗ | |||
---|---|---|---|---|---|
<2 ( n = 66) | 2 to <3 ( n = 109) | 3 to <4 ( n = 52) | ≥4 ( n = 69) | ||
LV (147) | 0.67 ± 1.07 (27) | 0.42 ± 1.18 (57) | 0.79 ± 1.17 (29) | 0.71 ± 1.00 (34) | .27 |
RV (93) | 3.58 ± 3.19 (25) | 2.03 ± 1.90 (36) | 1.91 ± 1.61 (12) | 1.82 ± 2.01 (20) | .15 |
Mixed (56) | 0.48 ± 1.57 (14) | 0.81 ± 1.56 (156) | 0.73 ± 1.28 (11) | 0.54 ± 1.45 (15) | .79 |
All (296) | 1.73 ± 2.61 (66) | 1.01 ± 1.66 (109) | 0.97 ± 1.41 (52) | 0.99 ± 1.53 (69) | .45 |
∗ P values were calculated by analysis of variance. Similar inferences were drawn using nonparametric testing (Kruskal-Wallis test).
Other Pre-Fontan Medical History Factors
Neither the individual diastolic function variables nor DFGs differed between subjects who underwent prior superior cavopulmonary anastomosis surgery compared with those who did not, even after adjusting for age ( Table 5 ). Among those who underwent volume-unloading surgery, the age at surgery was very mildly negatively correlated with E ( R = −0.13, P = .02) and E/E′ ( R = −0.15, P = .01) and positively correlated with E′ ( R = 0.15, P = .002). DFG also did not differ by type of Fontan procedure performed, presence or absence of a fenestration, presence of coarctation intervention at the time of the pre-Fontan catheterization, or pre-Fontan cardiac catheterization pressures ( Table 5 ).
Variable | Normal | Impaired relaxation | Pseudonormalization | Restrictive | chiao |
---|---|---|---|---|---|
N | 104 | 30 | 120 | 72 | |
SCPA performed | 68% | 73% | 78% | 69% | .43 |
LV morphology | 56% | 64% | 67% | 71% | .45 |
RV morphology | 86% | 91% | 91% | 84% | .81 |
Mixed morphology | 75% | 60% | 71% | 36% | .15 |
Type of Fontan | .60 | ||||
Aortopulmonary connection | 17% | 10% | 13% | 13% | |
Intracardiac lateral tunnel | 50% | 70% | 65% | 63% | |
Extracardiac conduit | 29% | 21% | 20% | 22% | |
Extracardiac lateral tunnel | 1% | 0% | 2% | 0% | |
Other | 3% | 0% | 1% | 3% | |
Fenestrated Fontan | 65% | 67% | 70% | 65% | .80 |
Pre-Fontan coarctation catheterization intervention | 3% | 3% | 8% | 7% | .42 |
n with catheterization | 104 | 30 | 120 | 72 | |
EDP (mm Hg), mean ± SD (median) | 7.1 ± 4.4 (8.0) | 7.4 ± 3.6 (8.0) | 7.0 ± 4.2 (7.0) | 7.4 ± 6.0 (7.0) | .95/.67 ∗ |
PA pressure (mm Hg), mean ± SD (median) | 9.6 ± 6.8 (11.0) | 11.0 ± 4.5 (11.0) | 8.9 ± 7.6 (11.0) | 8.7 ± 7.6 (10.0) | .65/.13 ∗ |
Oxygen saturation (%), mean ± SD (median) | 81 ± 21 (85) | 83 ± 5 (84) | 80 ± 19 (85) | 78 ± 22 (83) | .57/.75 ∗ |