Skip to main content

Association between left ventricular mechanics and diffuse myocardial fibrosis in patients with repaired Tetralogy of Fallot: a cross-sectional study



Patients with repaired tetralogy of Fallot (TOF) have progressive, adverse biventricular remodeling, leading to abnormal contractile mechanics. Defining the mechanisms underlying this dysfunction, such as diffuse myocardial fibrosis, may provide insights into poor long-term outcomes. We hypothesized that left ventricular (LV) diffuse fibrosis is related to impaired LV mechanics.


Patients with TOF were evaluated with cardiac magnetic resonance in which modified Look-Locker (MOLLI) T1-mapping and spiral cine Displacement encoding (DENSE) sequences were acquired at three LV short-axis positions. Linear mixed modeling was used to define the association between regional LV mechanics from DENSE based on regional T1-derived diffuse fibrosis measures, such as extracellular volume fraction (ECV).


Forty patients (26 ± 11 years) were included. LV ECV was generally within normal range (0.24 ± 0.05). For LV mechanics, peak circumferential strains (−15 ± 3%) and dyssynchrony indices (16 ± 8 ms) were moderately impaired, while peak radial strains (29 ± 8%) were generally normal. After adjusting for patient age, sex, and regional LV differences, ECV was associated with log-adjusted LV dyssynchrony index (β = 0.67) and peak LV radial strain (β = −0.36), but not LV circumferential strain. Moreover, post-contrast T1 was associated with log-adjusted LV diastolic circumferential strain rate (β = 0.37).


We observed several moderate associations between measures of fibrosis and impaired mechanics, particularly the LV dyssynchrony index and peak radial strain. Diffuse fibrosis may therefore be a causal factor in some ventricular dysfunction in TOF.


Early surgical repair of tetralogy of Fallot (TOF) results in 20-year survival rates over 90%, but mortality rates increase substantially 25 years after surgery, mostly from cardiac causes [1]. Most patients have right ventricular (RV) volume overload stemming from chronic pulmonary regurgitation, and adverse ventricular remodeling, such as RV dilation, is common. This remodeling may contribute to the development of systolic dysfunction of both the RV and left ventricle (LV), as quantified by ejection fraction (EF) [2, 3], or more complex descriptors of cardiac mechanics, such as strains and dyssynchrony [4,5,6]. Systolic dysfunction is associated with poor outcomes, such as sudden death, sustained ventricular tachycardia, or exercise intolerance [2, 3, 7]. Understanding of the mechanisms underlying the development of this dysfunction is incomplete.

Adverse ventricular remodeling may include diffuse myocardial fibrosis; for example, remodeling may cause a fibrotic expansion of the extracellular volume fraction (ECV) [8]. In fact, recent studies have identified increased myocardial ECV in patients with TOF using T1 mapping cardiovascular magnetic resonance imaging (CMR) [9, 10]. One of those studies additionally identified elevated LV ECV as a risk factor for ventricular arrhythmia [10].

The presence of both diffuse myocardial fibrosis and impaired ventricular mechanics in patients with repaired TOF may be associated with each other through physical interaction. That is, the stiffening effects of fibrosis on the myocardium may impair cardiac function [11], or the loss of myofibrils through processes leading to fibrosis may alter myocardial mechanics. In fact, recent studies have demonstrated weak-to-moderate associations between ECV and myocardial strains, strain rates, or dyssynchrony in hypertensive patients with hypertrophy, patients with heart failure, as well as subjects from the multi-ethnic study of atherosclerosis (MESA) cohort [12,13,14]. However, this association has not been studied in TOF, which represents a unique anatomic and physiologic setting. We hypothesized that LV diffuse fibrosis is associated with impaired LV mechanics in patients with repaired TOF.


Patient selection

The study protocol was approved by the institutional review board at The Children’s Hospital of Philadelphia, and informed consent for research was obtained from all patients. Between December 2012 and July 2015, a modified Look-Locker (MOLLI) T1 mapping sequence [15, 16] and spiral cine Displacement Encoding with Stimulated Echoes (DENSE) sequence [17,18,19] were both acquired in 43 patients with repaired TOF undergoing a clinical CMR evaluation.

Image acquisition and analysis

All images were acquired using a 1.5-T Siemens Avanto (Siemens Healthineers, Erlangen, Germany) CMR system. For each subject, the imaging protocol included cine balanced steady-state free precession (bSSFP) images, 2D phase contrast velocity, phase-sensitive inversion recovery, spiral cine DENSE, and MOLLI acquired both before and 15–20 min after administration of gadolinium contrast (gadopentetate dimeglumine, 0.2 mmole/kg, Bayer Health Care, Whippany New Jersey, USA). The MOLLI and DENSE acquisitions were co-localized at 3 ventricular short-axis locations (basal, mid-ventricular, and apical), such that cardiac mechanics and diffuse myocardial fibrosis could be interrogated over the same regions and directly compared. Details for the acquisition and post-processing for each specific sequence are as follows.


For MOLLI, the pre-contrast acquisition scheme was 5(3)2 or 5(3)3, while the post-contrast scheme was 4(1)3(1)2 or 4(1)2(1)2. The pulse sequence parameters included: Initial TI 90–120 ms, TE 1.0–1.3 ms, TR 3.6–12.9 ms, flip angle 35°, resolution 1.0–1.9 × 1.0–1.9 mm, and slice thickness 6–8 mm. Motion correction and curve fitting for T1 computation were performed online automatically during image reconstruction.

LV borders were manually contoured on MOLLI T1 images for fibrosis analysis, using custom MATLAB software (The Mathworks, Natick, Massachusetts, USA). Care was taken to restrict the region of interest to the mid-myocardial layer to minimize contamination of T1 signal at blood or air interfaces. The anterior RV insertion point was manually marked to annotate LV segments based on the standard 16-segment model. Native and post-contrast T1 values were quantified for the myocardium and blood pool for each region (slice) and each segment. From these measures, gadolinium partition coefficient (λGd) was computed:

$$ {\lambda}_{Gd}=\frac{\frac{1}{T_{1\ myo\ with\ Gd}}-\frac{1}{T_{1\ myo\ with out\ Gd}}}{\frac{1}{T_{1\ blood with\ Gd}}-\frac{1}{T_{1\ blood with out\ Gd}}}. $$

Hematocrit measurements were not available; however, hematocrit was estimated based on a validated relationship with the native T1 of blood [20]. Hence, the ‘synthetic’ ECV was computed as:

$$ ECV={\uplambda}_{Gd}\left(1- Hematocrit\right). $$


DENSE acquisition parameters included: 6 spiral interleaves (2 interleaves acquired per temporal frame); 0.10 cycles/mm in-plane encoding frequency, with simple encoding [18]; 0.08 cycles/mm through-plane de-phasing frequency [21]; Complementary spatial modulation of magnetization (CSPAMM) echo suppression [19]; TR 15 ms (30 ms temporal resolution without view sharing); TE 1.08 ms; flip angle 20°; resolution 1.8–2.8 × 1.8–2.8 mm; slice thickness 8 mm.

DENSE images were analyzed in DENSEanalysis [22], an open-source MATLAB program. Endocardial and epicardial boundaries were manually delineated on the DENSE magnitude images, followed by automatic phase unwrapping, spatial smoothing, and temporal fitting of phase displacement data [23], and calculation of LV mechanics. These measures included peak myocardial strains (radial and circumferential), peak circumferential strain rates (systolic and diastolic), and a dyssynchrony index (DI), each of which were quantified for each region (base, mid-ventricle, apex), in addition to a global mean value. To quantify dyssynchrony, the mean contraction delay was computed for each segment by comparing the circumferential strain curves for each element to a patient-specific reference curve with a cross-correlation analysis [24]. The segmental delays were circumferentially smoothed with a cubic spline interpolant, and DI for each region was computed as the standard deviation of the segmental delays [4]. Data were partitioned to the 16-segment model, based on the manual definition of the anterior RV insertion point, so peak strain results were also defined at the segmental level, matching the T1 data. The reproducibility of segmental strain rates has not been demonstrated, so those measures were only used to derive regional and global strain rates, and segmental strain rates were not included.


A contiguous short-axis stack of bSSFP images spanning the ventricles was also acquired. The pulse sequence parameters included: TR 2.7–9.0 ms, TE 1.2–4.5 ms, flip angle 15–90°, resolution 1.0–1.9 × 1.0–1.9 mm, slice thickness 6.5–10 mm. These images were manually segmented to quantify LV and RV volumes and ejection fraction and LV mass based on a sum of slices technique. Volumes and mass were indexed to body surface area [25].

2D phase contrast (velocity)

Pulmonary regurgitation was quantified using a single through-plane encoded phase contrast velocity encoded image positioned in the main pulmonary artery. The velocity encoding for these acquisitions ranged from 150 to 350 cm/s (adjusted to avoid phase aliasing). Other parameters included TR 11.4–13.3 ms, TE 2.5–3.0 ms, flip angle 25°, resolution 1.0–1.8 × 1.0–1.8 mm, slice thickness 4–5 mm. These images were manually segmented using custom MATLAB software (The Mathworks).

Phase-sensitive inversion recovery

Late gadolinium enhancement (LGE) was clinically evaluated using phase-sensitive inversion recovery images. For each subject, LGE was qualified as “present” or “absent” based on the clinical interpretation.

Intra- and inter-observer reproducibility

Segmentation of the MOLLI T1 images was repeated 2 months following the initial analysis by the primary analyst, and independently by a second analyst to quantify intra- and inter-observer reproducibility, respectively. The mean coefficient of variation (CoV) and the 95% limits of agreement [26] were quantified. The mean CoV was computed as:

$$ \frac{\sum_{i=1}^N\left[ St. Dev.{\left({X}_{Obs.1}\ {X}_{Obs.2}\right)}_i\right]}{\left|{\sum}_{i=1}^N\left[{\left(\frac{X_{Obs.1}+{X}_{Obs.2}}{2}\right)}_i\right]\right|} $$

Statistical analysis

Analysis was performed in R (version 3.4.0; R Foundation for Statistical Computing, Vienna, Austria) [27]. Data are presented as mean ± standard deviation. Z-scores for ventricular volumes, mass, and EF were computed based on population data assumed from Alfakih et al. [28] Pearson correlation was used to detect simple associations between strain/T1 derived measures and ventricular volumes, LV mass, and LV/RV EF. Linear mixed models [29] were fit to capture the mean association of regional (base, mid-ventricle, or apex) or segmental LV mechanics with T1-derived measures, while accounting for effects of sex, age, measurement region or segment, and within-subjects repeated measures (e.g., multiple regions or segments per patient). Repeated measures for the primary regional analysis were modeled using an unstructured correlation matrix. For the segmental analysis, a compound symmetric correlation matrix was assumed to ensure model convergence. Continuous outcome and independent variables were standardized for model fitting to report standardized β coefficients. Outcome variables were log transformed, as appropriate, to adjust for heteroscedasticity or non-normality. The statistical significance level was set at 0.05, with a Benjamini-Hochberg adjustment for multiple comparisons of each outcome variable [30]. Preliminary power analysis [31], assuming a linear correlation coefficient of 0.30, demonstrated that 84 samples were required to achieve 80% power.


Study group

From the original 43 patients, two were excluded for either poor DENSE image quality (n = 1) or errors with the reconstruction of T1 images (n = 1). An additional patient was excluded because of DENSE and MOLLI image slice mis-registration, leaving a total of 40 subjects in the final study group. An additional eight DENSE slices (4 base, 4 apex) from seven different subjects were also excluded for poor quality. Because of the need to compute a ventricular reference curve for the dyssynchrony analysis, the dyssynchrony index was not computed for these seven subjects. However, the remaining data for these subjects were included. Demographic details as well as ventricular volume indices, EF, and LV mass indices for these subjects are provided in Table 1. The RVs on average were severely dilated (Z = 4.6) with reduced EF (Z = −3.1), while the LVs had mostly normal chamber volume (Z = −0.6) and myocardial mass (Z = −0.7), but impaired EF (Z = −2.0). These structural and global functional characteristics are generally consistent with contemporary findings from the INDICATOR study [3].

Table 1 Demographic and ventricular volume details for 40 patients with repaired TOF

Most patients in this study (62.5%) were negative for LGE (Table 1). Of those with LGE reported, the predominant finding (in 73%) was a small region of enhancement at the inferior RV hinge point.

Diffuse fibrosis and strain findings

Representative T1 and DENSE images are shown in Fig. 1, and summary measures are presented in Table 2. Figure 2 represents the measured regional values of ECV, DI, and circumferential and radial strains, as well as global (“mean”) values. LV ECV (Fig. 2a) was generally within normal ranges, although mean values for seven patients (18%) were greater than 0.28, which has previously been reported as the upper bound of the normal range [10]. LV ECV was moderately associated with age (r = 0.35, p = 0.03), but not with LV or RV EF, LV mass index, LV or RV end-diastolic volume indices, or pulmonary regurgitant fraction (data not shown). Additionally, mean ECV was not different between patients based on the presence or absence of LGE (data not shown).

Fig. 1
figure 1

Example MOLLI T1 and DENSE images. Top row shows Pre- and Post-Contrast MOLLI T1 images for a representative patient. Bottom row shows DENSE magnitude and phase images for the same patient and same slice near end-systole, as well as the resulting displacement vectors at the same instant. The derived measures from each imaging sequence are listed for clarity

Table 2 Primary imaging findings
Fig. 2
figure 2

Box Plots of Circumferential and Radial Strains and ECV Results. Measures of a) extracellular volume fraction (ECV), b) peak circumferential strain, c) dyssynchrony index (DI), and d) peak radial strain represented as boxplots. Data are shown for each region (Base, Mid-ventricle, and Apex) as well as the global average (“Mean”) across slices

Peak LV circumferential strain (Fig. 2b) demonstrated a gradient of increasing magnitude from base to apex, which is consistent with normal LV function [32]. However, mean values were at the low (magnitude) end of the normal range, with 12 patients (31%) having strain magnitudes below the normal range (i.e., > −14.4%). Only two of these 12 patients also had an elevated mean LV ECV. This moderate impairment in LV circumferential strain mirrored the impairment in LV EF, as the two measures were significantly correlated (r = −0.54, p < 0.001), but mean LV circumferential strain was not associated with LV mass, LV/RV volumes, RV EF, or pulmonary regurgitation, and did not differ between patients with or without LGE. LV dyssynchrony was moderately elevated from the normal mean (DI = 16.1 ± 8 ms; Fig. 2c), and seven patients (21% with measured DI) had a value above the normal range (>21.1 ms) [4]. Only three of these seven patients also had reduced mean LV circumferential strain. Finally, mean LV radial strains were generally within the normal range (Fig. 2d), and did not differ between patients with or without LGE (data not shown).

At the segmental level, Fig. 3 presents the values for LV ECV, LV circumferential strain, and LV radial strain for the standard 16-segment LV model. LV ECV values tended to be higher in anterior and septal segments than the posterior and free wall, although the overall spatial variation was low (standard deviation of the segment means, SD = 0.02). LV circumferential strains in the mid-ventricular and apical free wall segments were preserved; however, strain in other segments were considerably lower (SD = 2.0%). LV radial strain also demonstrated segmental heterogeneity (SD = 4.4%), with higher strains in the basal and mid-ventricular free wall compared with the apex and mid-ventricular septum.

Fig. 3
figure 3

Bullseye Plots for Segmental Means. ECV (left), circumferential (middle) and radial strain (right) data presented with respect to the standard 16-segment left ventricular (LV) model. The mean of each segment is overlaid on the segment and graphically represented by the color relative to the other segments

Association between fibrosis and mechanics

From mixed model analyses, T1-derived measures had the strongest association with regional LV DI (Table 3), as λGd, ECV, and post-contrast T1 time were all significantly associated with increased regional LV DI after adjusting for age, sex, and regional differences. In each case, the β coefficients of these associations ranged between 0.33 (absolute) and 0.67, with ECV having the largest individual coefficient (Fig. 4a). Moreover, four of the seven patients with elevated mean ECV also had an abnormal mean LV DI. All T1-derived measures were also significantly associated with reduced regional LV radial strain, with ECV again having the largest magnitude coefficient (β = −0.36; Fig. 4b). The other statistically significant models were post-contrast T1 being related with regional LV systolic (β = −0.30) and (log-transformed) diastolic circumferential strain rates (β = 0.37; Fig. 4c). There was also a trend with post-contrast T1 and regional LV circumferential strain (p < 0.05; Table 3), but this did not meet statistical significance thresholds after adjustment for multiple hypothesis testing. No associations were observed at the segmental level (Table 4).

Table 3 Results of mixed model analyses for LV regional measures
Fig. 4
figure 4

T1-Based Measures are Associated with LV Mechanics. Scatter plots demonstrating the regional multivariate models in which ECV is associated with a) the log-adjusted LV dyssynchrony index (DI) and b) peak LV radial strain, and c) post-contrast T1 predicts the log-adjusted LV diastolic circumferential strain rate. The standardized β coefficient of the primary independent variable is reported in each case

Table 4 Results of mixed model analyses for LV segmental measures


Consistent with previous studies [10], T1-derived measures were highly reproducible at a regional level on both intra- and inter-observer bases (Fig. 5). In all cases, the mean CoV was less than or equal to 3%, denoting tight limits of agreement relative to measurement means. The reproducibility of segmental-level values was similarly strong, with mean CoVs ≤5% (not shown).

Fig. 5
figure 5

Reproducibility of T1 mapping. Bland-Altman comparisons of intra- (top) and inter-observer (bottom) analyses of T1 mapping data with respect to the quantification of post-contrast T1 (left), λGd (middle), and ECV (right). In addition to the bias and limits of agreement, the mean coefficient of variation (COV) is reported for each plot


Diffuse myocardial fibrosis and cardiac dysfunction of the LV have both been implicated as risk factors for adverse cardiac events, such as arrhythmia, in patients with repaired TOF [3, 10]. However, the extent to which these factors are linked—as has been shown in acquired heart disease [12,13,14, 33]—is unclear. Therefore, using MOLLI T1 mapping and DENSE CMR, we quantified markers of diffuse fibrosis and cardiac mechanics throughout the LV to elucidate the potential interplay between these characteristics for the first time in patients with repaired TOF. As hypothesized, we found associations at a regional level, indicating that T1-derived markers of fibrosis were related with LV dyssynchrony, peak LV radial strains, and LV circumferential strain rates. However, there was also evidence that components of LV function/dysfunction were independent of diffuse fibrosis, such as peak circumferential strains and ejection fraction. Collectively, these findings suggest that diffuse fibrosis is one of multiple drivers of cardiac dysfunction in TOF. The prognostic significance of these different factors will require further longitudinal study.

LV fibrosis

From our findings, the overall LV myocardial burden of diffuse fibrosis in patients with repaired TOF was generally low, based on data from healthy controls in other studies [34, 35]. This finding is consistent with the results of Chen et al., who similarly reported a mean LV ECV of 0.24 ± 0.03 in a cohort of 84 TOF patients of similar mean age [10]. However, seven patients (18%) had an elevated ECV (>0.28), which was previously associated with increased arrhythmic risk [10]. Continued follow-up of our patients is needed to determine if a similar risk association is observed. Our findings do conflict with those of Broberg et al., who reported a slightly higher LV ECV (0.29 ± 0.03) in a smaller (n = 17) patient group [9]. The patients in that study were older (34 ± 12 years), which may explain this difference given the positive association we observed between ECV and age.

LV function and fibrosis

Findings of reduced systolic LV function, including impaired myocardial strain, are common in patients with repaired TOF [4, 5, 36,37,38,39]; the fact that mean LV circumferential strain was moderately impaired in this study was not surprising. However, none of the T1-derived measures were associated with circumferential strain at either the regional or segmental levels, which was unexpected given that such associations were reported by Kuruvilla et al. [12] and Donekal et al. [14] in acquired heart diseases. This difference suggests a distinct etiology of impaired circumferential strain in the setting of repaired TOF, which is not fibrosis-related. Instead, these changes likely result from other mediating factors, such as surgical techniques [2] or septal shifting secondary to abnormal ventricular-ventricular interactions [36, 38]. Alternatively, the pattern of fibrosis may be such that myofibril contraction is affected to a much greater extent in the radial direction than in the circumferential direction.

The functional factor that demonstrated the strongest relationship with T1-derived measures of fibrosis was the LV dyssynchrony index. Findings of dyssynchrony in TOF are also common [4, 36, 40], and have been thought to be evidence of inter-ventricular interaction [4]. Our results suggest that diffuse fibrosis may be a primary mediator of the development of LV dyssynchrony in these patients. This linkage may help to explain the association between elevated LV ECV and increased risk of arrhythmia [10] through mechanisms such as electromechanical coupling [36, 41], but further study is warranted to definitively establish such connections. More generally, the prognostic significance of ventricular dyssynchrony in TOF has not been fully evaluated. A recent study found that dyssynchrony derived from CMR feature tracking did not predict longitudinal changes in LV or RV EF [5], although endpoints such as arrhythmias and mortality were not considered.

The remaining regional mechanical factors found to depend on T1-derived measures of fibrosis were LV radial strains and LV circumferential strain rates. Interestingly, the mean values for these measures were comparable to data for healthy controls from other studies [42,43,44], suggesting that they may serve as a compensatory mechanism to prevent more severe reductions in EF. The dependence on ECV, particularly for radial strain, is thus a potential cause for concern if the fibrotic expansion or cell atrophy [10] weakens that compensatory ability. Longitudinal follow-up, particularly with continued increase in ECV expansion, is needed to further evaluate these relationships.


This was a small cross-sectional study, so the ability to firmly establish causal relationships or long-term significance was limited. However, the use of CMR to measure all endpoints did ensure excellent reproducibility and reduced variance to enhance statistical power [45].

We did not quantify fibrosis or cardiac mechanics for the RV, which would be of obvious interest in repaired TOF, because of the generally insufficient resolution of the images for the thin wall of the RV. Future studies will leverage newer, higher resolution T1 mapping techniques [46] and higher resolution DENSE to better resolve the RV [44]. However, given the fact that LV diffuse fibrosis [10] and LV systolic dysfunction [3] have been implicated as predictors of poor outcomes in TOF, our focus on the LV is highly relevant.

Hematocrit data were not systematically available for all patients, necessitating the use of a “synthetic” hematocrit estimation to define ECV. This model may have produced inaccurate estimates for some patients, particularly younger individuals. While not the current gold standard for ECV estimation, it is important to note that the synthetic hematocrit model was derived and validated using a MOLLI sequence on 1.5 T Siemens scanners (inclusive of the Avanto model used for the present study), and demonstrated excellent agreement with the conventional measurement (95% limits of agreement of approximately ±4%) [20]. Given these findings, we feel its use, in the absence of hematocrit data, is justified. Moreover, associations of ECV and mechanics were generally similar to those of other T1-based endpoints, which did not require the synthetic hematocrit assumption, demonstrating the appropriateness of the ECV estimations.

No normal reference range for ECV was available from our laboratory, as is recommended for clinical reporting [47]. Differences in acquisition and analysis techniques may thus confound comparisons of absolute values to published reference ranges. However, as the primary objective of the study was to quantify associations with cardiac mechanics, this limitation had minimal effect on our primary findings.

Finally, no long-axis DENSE or T1 data were acquired to allow for the quantification of LV longitudinal strain. Previous studies have shown impairments in longitudinal strain in patients with TOF [36], so future studies should additionally explore the relationship of such impairments with fibrosis.


The factors governing LV mechanics in patients with repaired TOF are complex. We found that elevated diffuse LV myocardial fibrosis was associated with regional LV dyssynchrony and peak radial strain, while decreasing post-contrast T1 was associated with abnormal circumferential strain rates. Extracellular matrix expansion and fibrotic remodeling thus potentially represent one of the mediating factors of LV function. However, other important components of LV function, such as regional circumferential strain, appeared to be impaired independent of fibrosis. Hence, while it is tempting to speculate that diffuse fibrosis may precede frank systolic dysfunction and provide early prognostic insight, it appears more likely that detailed characterization of both fibrosis and cardiac mechanics as part of routine surveillance of these patients is needed to fully capture their independent contributions.



Balanced steady state free precession


Cardiovascular magnetic resonance


Coefficient of Variation


Displacement encoding with stimulated echoes


Dyssynchrony index


Extracellular volume fraction


Ejection fraction


Late gadolinium enhancement


Left ventricle/left ventricular


Modified Look-Locker


Right ventricle/right ventricular


Echo time


Tetralogy of Fallot


Repetition time

λGd :

Gadolinium partition coefficient


  1. Nollert G, Fischlein T, Bouterwek S, Bohmer C, Klinner W, Reichart B. Long-term survival in patients with repair of tetralogy of Fallot: 36-year follow-up of 490 survivors of the first year after surgical repair. J Am Coll Cardiol. 1997;30:1374–83.

    Article  CAS  PubMed  Google Scholar 

  2. Geva T, Sandweiss BM, Gauvreau K, Lock JE, Powell AJ. Factors associated with impaired clinical status in long-term survivors of tetralogy of Fallot repair evaluated by magnetic resonance imaging. J Am Coll Cardiol. 2004;43:1068–74.

    Article  PubMed  Google Scholar 

  3. Valente AM, Gauvreau K, Assenza GE, Babu-Narayan SV, Schreier J, Gatzoulis MA, Groenink M, Inuzuka R, Kilner PJ, Koyak Z, Landzberg MJ, Mulder B, Powell AJ, Wald R, Geva T. Contemporary predictors of death and sustained ventricular tachycardia in patients with repaired tetralogy of Fallot enrolled in the INDICATOR cohort. Heart. 2014;100:247–53.

    Article  PubMed  Google Scholar 

  4. Jing L, Haggerty CM, Suever JD, Alhadad S, Prakash A, Cecchin F, Skrinjar O, Geva T, Powell AJ, Fornwalt BK. Patients with repaired tetralogy of Fallot suffer from intra- and inter-ventricular cardiac dyssynchrony: a cardiac magnetic resonance study. Eur Heart J Cardiovasc Imaging. 2014;15:1333–43.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Jing L, Wehner GJ, Suever JD, Charnigo RJ, Alhadad S, Stearns E, Mojsejenko D, Haggerty CM, Hickey K, Valente AM, Geva T, Powell AJ, Fornwalt BK. Left and right ventricular dyssynchrony and strains from cardiovascular magnetic resonance feature tracking do not predict deterioration of ventricular function in patients with repaired tetralogy of Fallot. J Cardiovasc Magn Reson. 2016;18:49.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Menteer J, Weinberg P, Fogel M. Quantifying regional right ventricular function in Tetralogy of Fallot. J Cardiovasc Magn Reson. 2005;7:753–61.

    Article  PubMed  Google Scholar 

  7. Knauth AL, Gauvreau K, Powell AJ, Landzberg MJ, Walsh EP, Lock JE, del Nido PJ, Geva T. Ventricular size and function assessed by cardiac MRI predict major adverse clinical outcomes late after tetralogy of Fallot repair. Heart. 2008;94:211–6.

    Article  CAS  PubMed  Google Scholar 

  8. Weber KT, Brilla CG. Pathological hypertrophy and cardiac interstitium. Fibrosis and renin-angiotensin-aldosterone system. Circulation. 1991;83:1849–65.

    Article  CAS  PubMed  Google Scholar 

  9. Broberg CS, Chugh SS, Conklin C, Sahn DJ, Jerosch-Herold M. Quantification of diffuse myocardial fibrosis and its association with myocardial dysfunction in congenital heart disease. Circ Cardiovasc Imaging. 2010;3:727–34.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Chen C-A, Dusenbery SM, Valente AM, Powell AJ, Geva T. Myocardial ECV fraction assessed by CMR is associated with type of hemodynamic load and arrhythmia in repaired Tetralogy of Fallot. JACC Cardiovasc Imaging. 2016;9:1–10.

    Article  PubMed  Google Scholar 

  11. Díez J, Querejeta R, López B, González A, Larman M, Martínez Ubago JL. Losartan-dependent regression of myocardial fibrosis is associated with reduction of left ventricular chamber stiffness in hypertensive patients. Circulation. 2002;105:2512–7.

    Article  PubMed  Google Scholar 

  12. Kuruvilla S, Janardhanan R, Antkowiak P, Keeley EC, Adenaw N, Brooks J, Epstein FH, Kramer CM, Salerno M. Increased extracellular volume and altered mechanics are associated with LVH in hypertensive heart disease, not hypertension alone. JACC Cardiovasc Imaging. 2015;8:172–80.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Lin L-Y, Wu C-K, Juang J-MJ, Wang Y-C, Su M-YM, Lai L-P, Hwang J-J, Chiang F-T, Tseng W-YI, Lin J-L. Myocardial regional interstitial fibrosis is associated with left intra-ventricular Dyssynchrony in patients with heart failure: a cardiovascular magnetic resonance study. Sci Rep. 2016;6(2015):20711.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Donekal S, Venkatesh BA, Liu YC, Liu CY, Yoneyama K, Wu CO, Nacif M, Gomes AS, Hundley WG, Bluemke DA, Lima JA. Interstitial fibrosis, left ventricular remodeling, and myocardial mechanical behavior in a population-based multiethnic cohort: the multi-ethnic study of atherosclerosis (MESA) study. Circ Cardiovasc Imaging. 2014;7:292–302.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Look DC, Locker DR. Time saving in measurement of NMR and EPR relaxation times. Rev Sci Instrum. 1970;41:250–1.

    Article  CAS  Google Scholar 

  16. Messroghli DR, Radjenovic A, Kozerke S, Higgins DM, Sivananthan MU, Ridgway JP. Modified look-locker inversion recovery (MOLLI) for high-resolution T1 mapping of the heart. Magn Reson Med. 2004;52:141–6.

    Article  PubMed  Google Scholar 

  17. Aletras AH, Ding S, Balaban RS, Wen H. DENSE: displacement encoding with stimulated echoes in cardiac functional MRI. J Magn Reson. 1999;137:247–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Zhong X, Helm PA, Epstein FH. Balanced multipoint displacement encoding for DENSE MRI. Magn Reson Med. 2009;61:981–8.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Kim D, Gilson WD, Kramer CM, Epstein FH. Myocardial tissue tracking with two-dimensional cine displacement-encoded MR imaging: development and initial evaluation. Radiology. 2004;230:862–71.

    Article  PubMed  Google Scholar 

  20. Treibel TA, Fontana M, Maestrini V, Castelletti S, Rosmini S, Simpson J, Nasis A, Bhuva AN, Bulluck H, Abdel-Gadir A, White SK, Manisty C, Spottiswoode BS, Wong TC, Piechnik SK, Kellman P, Robson MD, Schelbert EB, Moon JC. Automatic measurement of the myocardial Interstitium synthetic extracellular volume quantification without Hematocrit sampling. JACC Cardiovasc Imaging. 2016;9:54–63.

    Article  PubMed  Google Scholar 

  21. Zhong X, Spottiswoode BS, Cowart EA, Gilson WD, Epstein FH. Selective suppression of artifact-generating echoes in cine DENSE using through-plane dephasing. Magn Reson Med. 2006;56:1126–31.

    Article  PubMed  Google Scholar 

  22. Gilliam AD, Suever JD, and contributors: DENSEanalysis. Version 0.4.0. 2016. [].

    Google Scholar 

  23. Spottiswoode BS, Zhong X, Hess AT, Kramer CM, Meintjes EM, Mayosi BM, Epstein FH. Tracking myocardial motion from cine DENSE images using spatiotemporal phase unwrapping and temporal fitting. IEEE Trans Med Imaging. 2007;26:15–30.

    Article  CAS  PubMed  Google Scholar 

  24. Suever JD, Fornwalt BK, Neuman LR, Delfino JG, Lloyd MS, Oshinski JN. Method to create regional mechanical dyssynchrony maps from short-axis cine steady-state free-precession images. J Magn Reson Imaging. 2014;39:958–65.

    Article  PubMed  Google Scholar 

  25. Mosteller R. Simplified calculation of body-surface area. N Engl J Med. 1987;317:1098.

    CAS  PubMed  Google Scholar 

  26. Bland JM, Altman D. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986;1:307–10.

    Article  CAS  PubMed  Google Scholar 

  27. R Development Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2017. []

    Google Scholar 

  28. Alfakih K, Plein S, Thiele H, Jones T, Ridgway JP, Sivananthan MU. Normal human left and right ventricular dimensions for MRI as assessed by turbo gradient echo and steady-state free precession imaging sequences. J Magn Reson Imaging. 2003;17:323–9.

    Article  PubMed  Google Scholar 

  29. Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team: {nlme}: linear and nonlinear mixed effects models. R package version 3.1–131. 2017. [].

    Google Scholar 

  30. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995;57:289–300.

    Google Scholar 

  31. Champely S: Pwr: basic functions for power analysis. R package version 1.2–1. 2017. [].

    Google Scholar 

  32. Wehner GJ, Suever JD, Haggerty CM, Jing L, Powell DK, Hamlet SM, Grabau JD, Mojsejenko WD, Zhong X, Epstein FH, Fornwalt BK. Validation of in vivo 2D displacements from spiral cine DENSE at 3T. J Cardiovasc Magn Reson. 2015;17:5.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Collins J, Sommerville C, Magrath P, Spottiswoode B, Freed BH, Benzuly KH, Gordon R, Vidula H, Lee DC, Yancy C, Carr J, Markl M. Extracellular volume fraction is more closely associated with altered regional left ventricular velocities than left ventricular ejection fraction in nonischemic Cardiomyopathy. Circ Cardiovasc Imaging. 2014;8:e001998.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Dusenbery SM, Jerosch-Herold M, Rickers C, Colan SD, Geva T, Newburger JW, Powell AJ. Myocardial extracellular remodeling is associated with ventricular diastolic dysfunction in children and young adults with congenital aortic stenosis. J Am Coll Cardiol. 2014;63:1778–85.

    Article  PubMed  Google Scholar 

  35. Jerosch-Herold M. Cardiac magnetic resonance imaging of myocardial contrast uptake and blood flow in patients affected with idiopathic or familial dilated cardiomyopathy. Am J Hear Circ Physiol. 2008;2115:1234–42.

    Article  Google Scholar 

  36. Tzemos N, Harris L, Carasso S, Subira LD, Greutmann M, Provost Y, Redington AN, Rakowski H, Siu SC, Silversides CK. Adverse left ventricular mechanics in adults with repaired tetralogy of Fallot. Am J Cardiol. 2009;103:420–5.

    Article  PubMed  Google Scholar 

  37. Diller GP, Kempny A, Liodakis E, Alonso-Gonzalez R, Inuzuka R, Uebing A, Orwat S, Dimopoulos K, Swan L, Li W, Gatzoulis MA, Baumgartner H. Left ventricular longitudinal function predicts life-threatening ventricular arrhythmia and death in adults with repaired tetralogy of fallot. Circulation. 2012;125:2440–6.

    Article  PubMed  Google Scholar 

  38. Fernandes FP, Manlhiot C, Roche SL, Grosse-Wortmann L, Slorach C, McCrindle BW, Mertens L, Kantor PF, Friedberg MK. Impaired left ventricular myocardial mechanics and their relation to pulmonary regurgitation, right ventricular enlargement and exercise capacity in asymptomatic children after repair of tetralogy of Fallot. J Am Soc Echocardiogr. 2012;25:494–503.

    Article  PubMed  Google Scholar 

  39. Li S-N, Yu W, Lai CT-M, Wong SJ, Cheung Y-F. Left ventricular mechanics in repaired tetralogy of Fallot with and without pulmonary valve replacement: analysis by three-dimensional speckle tracking echocardiography. PLoS One. 2013;8:e78826.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Abd El Rahman MY, Hui W, Yigitbasi M, Dsebissowa F, Schubert S, Hetzer R, Lange PE, Abdul-Khaliq H. Detection of left ventricular asynchrony in patients with right bundle branch block after repair of tetralogy of Fallot using tissue-Doppler imaging-derived strain. J Am Coll Cardiol. 2005;45:915–21.

    Article  PubMed  Google Scholar 

  41. Gatzoulis MA, Till JA, Somerville J, Redington AN. Mechanoelectrical interaction in Tetralogy of Fallot. Circulation. 1995;92:231–7.

    Article  CAS  PubMed  Google Scholar 

  42. Moody WE, Taylor RJ, Edwards NC, Chue CD, Umar F, Taylor TJ, Ferro CJ, Young AA, Townend JN, Leyva F, Steeds RP. Comparison of magnetic resonance feature tracking for systolic and diastolic strain and strain rate calculation with spatial modulation of magnetization imaging analysis. J Magn Reson Imaging. 2015;41:1000–12.

    Article  PubMed  Google Scholar 

  43. Jing L, Binkley CM, Suever JD, Umasankar N, Haggerty CM, Rich J, Nevius CD, Wehner GJ, Hamlet SM, Powell DK, Radulescu A, Kirchner HL, Epstein FH, Fornwalt BK. Cardiac remodeling and dysfunction in childhood obesity: a cardiovascular magnetic resonance study. J Cardiovasc Magn Reson. 2016;18:28.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Suever JD, Wehner GJ, Jing L, Powell DK, Hamlet SM, Grabau JD, Mojsejenko D, Andres KN, Haggerty CM, Fornwalt BK. Right ventricular strain, torsion, and Dyssynchrony in healthy subjects using 3D spiral cine DENSE magnetic resonance imaging. IEEE Trans Med Imaging. 2017;36:1076–85.

    Article  PubMed  Google Scholar 

  45. Grothues F, Smith GC, Moon JCC, Bellenger NG, Collins P, Klein HU, Pennell DJ. Comparison of interstudy reproducibility of cardiovascular magnetic resonance with two-dimensional echocardiography in normal subjects and in patients with heart failure or left ventricular hypertrophy. Am J Cardiol. 2002;90:29–34.

    Article  PubMed  Google Scholar 

  46. Mehta BB, Chen X, Bilchick KC, Salerno M, Epstein FH: Accelerated and navigator-gated look-locker imaging for cardiac t1 estimation (ANGIE): development and application to T1 mapping of the right ventricle. Magn Reson Med 2015, 160:150–160.

  47. Messroghli DR, Moon JC, Ferreira VM, Grosse-Wortmann L, He T, Kellman P, Mascherbauer J, Nezafat R, Salerno M, Schelbert EB, Taylor AJ, Thompson R, Ugander M, Van Heeswijk RB, Friedrich MG. Clinical recommendations for cardiovascular magnetic resonance mapping of T1, T2, T2* and extracellular volume: a consensus statement by the Society for Cardiovascular Magnetic Resonance (SCMR) endorsed by the European Association for Cardiovascular Imagi. J Cardiovasc Magn Reson. 2017;19:1–24.

    Article  Google Scholar 

Download references


Not applicable.


This work was supported by a National Institutes of Health (NIH) Director’s Early Independence Award (DP5 OD-012132). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available on reasonable request with approval of the corresponding author.

Author information

Authors and Affiliations



CH conceived the study, participated in the study design, analyzed the data, and drafted the manuscript. JS, AMS, AP, GW, and LJ assisted with data analysis and provided critical revisions to the manuscript. RC helped with the design of the statistical methods and provided critical revision of the manuscript. BF assisted with the study design and analysis, and provided critical revision of the manuscript. MF collected the data, conceived the study, participated in the study design, and provided critical revision of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Christopher M. Haggerty.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board at the Children’s Hospital of Philadelphia. All subjects provided informed consent for research participation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Haggerty, C.M., Suever, J.D., Pulenthiran, A. et al. Association between left ventricular mechanics and diffuse myocardial fibrosis in patients with repaired Tetralogy of Fallot: a cross-sectional study. J Cardiovasc Magn Reson 19, 100 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: