Right-left ventricular shape variations in tetralogy of Fallot: associations with pulmonary regurgitation
Journal of Cardiovascular Magnetic Resonance volume 23, Article number: 105 (2021)
Relationships between right ventricular (RV) and left ventricular (LV) shape and function may be useful in determining optimal timing for pulmonary valve replacement in patients with repaired tetralogy of Fallot (rTOF). However, these are multivariate and difficult to quantify. We aimed to quantify variations in biventricular shape associated with pulmonary regurgitant volume (PRV) in rTOF using a biventricular atlas.
In this cross-sectional retrospective study, a biventricular shape model was customized to cardiovascular magnetic resonance (CMR) images from 88 rTOF patients (median age 16, inter-quartile range 11.8–24.3 years). Morphometric scores quantifying biventricular shape at end-diastole and end-systole were computed using principal component analysis. Multivariate linear regression was used to quantify biventricular shape associations with PRV, corrected for age, sex, height, and weight. Regional associations were confirmed by univariate correlations with distances and angles computed from the models, as well as global systolic strains computed from changes in arc length from end-diastole to end-systole.
PRV was significantly associated with 5 biventricular morphometric scores, independent of covariates, and accounted for 12.3% of total shape variation (p < 0.05). Increasing PRV was associated with RV dilation and basal bulging, in conjunction with decreased LV septal-lateral dimension (LV flattening) and systolic septal motion towards the RV (all p < 0.05). Increased global RV radial, longitudinal, circumferential and LV radial systolic strains were significantly associated with increased PRV (all p < 0.05).
A biventricular atlas of rTOF patients quantified multivariate relationships between left–right ventricular morphometry and wall motion with pulmonary regurgitation. Regional RV dilation, LV reduction, LV septal-lateral flattening and increased RV strain were all associated with increased pulmonary regurgitant volume. Morphometric scores provide simple metrics linking mechanisms for structural and functional alteration with important clinical indices.
The survival rate of repaired tetralogy of Fallot (rTOF) patients has greatly increased due to the improvement of surgical repair. Currently, the early mortality rate is below 2%  and the 25-year survival rate is above 94% . However, pulmonary regurgitation (PR) is a common consequence of the surgical relief of right ventricular (RV) outflow tract narrowing. rTOF patients are therefore monitored for PR, and its effects on RV dysfunction and long-term outcomes [3, 4]. Despite being well tolerated in childhood, PR has damaging effects on RV dysfunction and leads to RV dilation over an ill-defined period of time. In follow-up, rTOF patients may demonstrate progressive exercise intolerance, arrhythmia, right- or left-sided heart failure, RV and/or left ventricular (LV) dysfunction or sudden cardiac death [5,6,7]. Patients with a history of severe PR are at risk for irreversible RV functional impairment making timing for pulmonary valve replacement a critical issue [8, 9]. However, the timing of valve replacement based on RV size and systolic function remains a subject of debate, with wide inter-institutional, qualitative variations of clinical practice [3, 4, 10, 11].
Studies of RV remodeling in rTOF have demonstrated a strong relationship between RV dilation, RV bulging, and apical dilation with the presence of pulmonary regurgitation [12,13,14]. However, the relationships with LV function and the underlying mechanisms are less understood in rTOF population and few studies have focused on RV-and-LV interactions [15,16,17,18,19]. Increased RV dilation has been known to affect LV function in rTOF [19,20,21,22] and approximately 20% of adult rTOF patients develop LV dysfunction . Furthermore, the close relationship between LV ejection fraction (LVEF) and RV ejection fraction (RVEF) has been previously demonstrated [24, 25], indicating the importance of ventricular interactions in rTOF. This has led to the suggestion that LV function and inter-ventricular interactions should be considered in the timing of pulmonary valve replacement [26, 27]. Conventional clinical metrics of ventricular volume and function remain limited in their scope and do not fully explain the contribution of biventricular shape and interaction, thus limiting the value of imaging data in clinical evaluation and understanding of pathophysiology in relation to clinical outcomes.
Here, we use a biventricular atlas  to quantify relationships between RV and LV shape and function, and examine their associations with PR, in a retrospective cross-sectional study of 88 rTOF patients with no history of pulmonary valve replacement. We aimed to determine whether novel 3D shape features, calculated from standard 3D imaging exams, can be used to quantify specific shape and function alterations associated with clinically important metrics such as PR. If so, these methods could be readily applied to any cross-sectional imaging exam to provide shape scores in relation to factors of interest, thereby providing a mechanistic link between these factors and alterations in cardiac geometry and function and reducing inherent inaccuracies of qualitative clinical decision-making.
Cardiovascular magnetic resonance (CMR) examinations and clinical data from 88 rTOF patients were obtained from the Cardiac Atlas Project (CAP) congenital heart disease database .
CAP is a large-scale database of cardiac images and limited associated clinical data that facilitates data sharing for the development of new methods for image analysis and collaborative statistical analysis of heart shape and function across multiple cohorts (cardiacaltas.org). Datasets related to congenital heart disease were added to the CAP beginning in 2015.
In this paper, deidentified datasets were contributed from two clinical centers (Auckland, New Zealand and San Diego, California, United States) with approval from local institutional review boards compatible with data sharing. Demographic data are shown in Table 1. Patients with pulmonary valve replacement, or severe tricuspid regurgitation from either echocardiography or CMR were excluded.
CMR imaging protocol and image analysis
CMR images were acquired with either prospectively or retrospectively electrocardiogram (ECG) gated balanced steady-state free-precession cine sequence on 1.5 T CMR scanners (Avanto, Siemens Healthineers, Erlangen, Germany, or Discovery MR450, GE Healthcare Systems, Chicago, Illinois, United States) during breath-holding. The short axis slices were acquired parallel to the tricuspid annulus plane and spanned both ventricles. Long-axis slices were obtained through all valves in standard 4-chamber, 2-chamber, LV outflow tract and RV outflow tract views. Typical imaging parameters were: repetition time 24-32 ms; echo time 1.1–1.5 ms; flip angle 70–80°; pixel size reconstructed to in plane 0.59–1.75 \(\times\) 0.59–1.75 mm; slice thickness 4–6 mm; number of time frames 20–35; image matrix 180–224 \(\times\) 208–256, and field of view 200–300 mm.
Antegrade and retrograde pulmonary flow measurements were obtained from two-dimensional phase contrast (2D PC) imaging. PC analysis of antegrade and retrograde flows in the main artery was performed on a plane at a location just below the pulmonary artery (PA) bifurcation and perpendicular to the axis of the PA. Typical imaging parameters were: acceleration factor 3, echo time = 2.3–3.0 ms, repetition time = 4.8–5.0 ms, field of view = 169–315 mm × 300–420 mm, spatial resolution = (1.4–2.0) × (1.4–2.0) × (5–8) mm3, temporal resolution = 37–41 ms, flip angle of 15°-30°. Scouts were used to set the velocity encoding.
Contours were drawn manually on both long axis and short axis cine slices by one expert analyst using Segment . Contours from another independent expert analyst were also performed for 35 cases. Ventricular masses and volumes at end-diastole (ED) and end-systole (ES) for both ventricles were calculated using the volumetric summation of discs method. Papillary muscles were excluded from the masses and endocardial contours segmentation. Tricuspid and mitral valve hinge points were defined from the intersection of the left atrial and ventricular contours delineated on the 2-chamber and 4-chamber long axis images, and right atrial and ventricular contours on the 4-chamber long axis images. Aortic valve hinge points were extracted from the LV outflow images and the ventricular extent of the RV outflow tract was extracted from the RV outflow tract images. When aortic and pulmonary leaflets were not visible, boundary points were defined by the transition in appearance from myocardium to vessel wall.
A biventricular subdivision surface template mesh was constructed as described previously . An overview of the analysis pipeline is shown in Fig. 1. This included the LV, RV, and the four valves (aortic, mitral, pulmonary, and tricuspid). The template mesh was automatically customized to each patient and breath-hold slice misregistration was also automatically corrected using an iterative registration algorithm . Valve locations were customized to the manual landmarks by using landmark registration and surfaces were customized by using diffeomorphic non-rigid registration to the manual contours. LV and RV volume and mass were calculated by numerical integration of mesh volumes. To build the biventricular atlas, all the ED mesh points were first aligned to the mean mesh surface points by a rigid registration (translation and rotation). This transformation was then applied to the ES mesh. The ED and ES surface points were then concatenated to form a single combined shape. Principal component analysis (PCA) was used to evaluate the distribution of shape variation across the cohort . This results in a relatively small number of components (24 in this study) that describe the shape variation across all patients, while accounting for correlations between points in the model. The first component explains the most variance, the second component explains the most remaining variance, and so on. For each patient, morphometric scores were calculated, which quantified the amount of each component present (“Appendix”). Owing to the combination of ED and ES shapes in each mode, both shape and shape changes between ED and ES could be captured using this method.
Pulmonary regurgitant volume
Pulmonary regurgitation was quantified from CMR imaging using pulmonary regurgitant volume indexed to body surface area (PRVI) as this has been shown to be a more accurate reflection of regurgitation severity compared with pulmonary regurgitant fraction [12, 33,34,35]. PR was quantified using PC imaging, and forward and regurgitant flows in the main PA were quantified using commercially available software (Argus Flow, Siemens Healthineers, and cvi42, Circle Cardiovascular Imaging, Calgary, Alberta, Canada). The range of PRVI was 1.5–95.7 ml/m2 with median of 23.7 and IQR (14.2, 33.3) ml/m2 (Table 1).
Multivariate associations with pulmonary regurgitant volume
A linear regression model was constructed to estimate associations with biventricular heart shape, using the morphometric scores as the response (dependent) variables . PRVI, height, weight and age were included as continuous predictor variables, and sex and tricuspid regurgitation severity from the CMR report (none/mild/mild-to-moderate/moderate) were included as categorical predictors. Height, weight, sex, and age were added to the model to control for body habitus as done previously . LV and RV volumes and mass could be calculated from the resulting morphometric scores, as well as LV and RV dimensions.
Regional shape variations associated with PRVI were confirmed using univariate regression models for specific features derived from the 3D patient-specific biventricular models. These included: basal bulge, tricuspid tilting, apical dilation, RV and LV anterior–posterior dimensions and lateral to septal dimensions (Fig. 2). Basal bulge was calculated as the distance between the most basal point on the RV free-wall and its perpendicular projection onto a plane perpendicular to the LV long axis and going through the mitral valve centroid. Tricuspid tilt was calculated as the angle between a plane perpendicular to the LV long axis and going through the centroid of the mitral valve and a plane fitted to the tricuspid valve points, similar to . Apical dilation was quantified using the apical angle. Apical angle was defined as in , using two lines defined at the RV apex in a four-chamber view: one aligned with the endocardium of the septum and another line aligned parallel to the most linear portion of the RV free wall endocardium. LV and RV lateral to septal dimensions were defined as the length of the minor axis perpendicular to the interventricular septum at mid-ventricle in both LV and RV. LV and RV anterior–posterior dimensions were defined as the length of the minor axis parallel to the interventricular septum at mid-ventricle in both LV and RV. Figure 2 summarizes how those remodeling features were calculated.
To further study interventricular relationships, model-derived systolic deformations for both chambers were assessed. Geometric strain was defined as the change in geometric arc length from ED to ES. Longitudinal strain (LS), which represents the longitudinal shortening of the cardiac muscle, and circumferential strain (CS), shortening along the circular perimeter, and radial strain (RS), thickening of the wall during systole, were derived from length changes between ED and ES using the Cauchy strain formula:
where ε is the engineering strain, L the length at ES and L0 the length at ED. Changes in geometric arc length have previously shown good agreement with myocardial strain derived from tagging and displacement encoded CMR[38, 39]. For RS, myocardial wall thickness was calculated using a modified version of the centerline method  using the intersection between short axis slices and the biventricular model. For CS and LS, arc lengths were calculated from the model mesh. RV free wall and LV RS and CS were averaged over three different levels: base, mid-ventricle and apex. LV and RV LS were averaged over the apical 4-chamber and RV outflow tract views. The interventricular septum was included in the LV strain.
Univariate regression models were used to quantify the association between LV and RV strains and the amount of regurgitation. Although systolic CS and LS are conventionally negative, their absolute values were used in this study for simpler interpretation of the association with PRVI.
Statistical analysis was carried out with R  (R Foundation for Statistical Computing, Vienna, Austria). All data are reported as mean and standard deviation or median and interquartile range, depending on the distribution, for continuous variable and as frequency for categorical variables. Normality was tested using Shapiro–Wilk tests. All variables were standardized before regression. Statistical differences were presented by p-values using one-way ANOVA or Kruskal–Wallis test depending on the distribution. A p-value of 0.05 was considered significant for the overall effect. Parameter estimates (PE), which represent the change in the response variable associated with a 1-unit change of the predictor, were used to measure the strength of the association between strain measurements, morphometric scores, and pulmonary regurgitation.
RVLV model customization
All models were successfully customized to the manual contours and landmarks. Table 2 shows the mean volume and mass of LV and RV calculated from the biventricular models by numerical integration; they had good agreement with those calculated from the manual contours by short axis slice summation. Inter-observer errors in manual contour results in 35 participants are also reported in Table 2. The differences between model and manual estimates were similar to the inter-observer differences. The larger model RV volumes may be due to inaccuracies of short axis slice summation at the base and apex (in particular due to the basal bulge common in rTOF) since the biventricular model incorporates information from the long axis contours as well as the short axis contours. The larger LV mass in the model may be due to differences in the definition of LV vs RV myocardial partitions in the model relative to the contours.
The first four principal components of shape variation, accounting for the most variation in biventricular shape across the cohort (total 59%), are shown in Fig. 3 (anterior view). Animations of these shape variations can be found in Additional file 1. These represent the largest variations of biventricular shape as well as shape changes between ED to ES within the cohort. The first mode was associated with overall size and accounted for 37.5% of the total shape variance. The second mode accounted for 9.3% of the total shape variance and was associated with septal-free wall dimension (i.e., expansion and contraction of both ventricles toward and away from the interventricular septum) as well as the transition from a LV dominant apex shape to a RV dominant shape. The third mode (7.1%) was associated with basal vs apical bulging of the RV, with basal bulging associated with a tilt of the LV base away from the septum. The fourth mode (5.3%) was associated with a systolic septal deviation towards the RV, a common feature in rTOF remodeling due to volume overload [42, 43] in concert with LV diminution.
The first 24 PCA modes accounted for 90.2% of the total shape variance and were selected to perform the multivariate regression. Morphometric scores were used as outcome (dependent) variables and sex, height, weight, age, degree of tricuspid regurgitation, and PRVI were used as predictor (independent) variables. The relationships between each predictor variable and the associated shapes were visualized using the method described in “Appendix”.
Height was a significant predictor in modes 1, 4, 9 and 15 (p < 0.001, p < 0.001, p = 0.005, and p = 0.04 respectively), accounting for 27% of the total variation. Weight was a significant predictor in modes 1 and 4 (p < 0.001 and p < 0.001 respectively, accounting for 16.3% of the total variation), age in modes 6, 9 and 14 (p = 0.01, p = 0.003, and p = 0.04 respectively) and sex in modes 1, 8, 9 and 14 (p = 0.003, p = 0.002, p = 0.006, and p = 0.04 respectively).
Tricuspid regurgitation was significant predictor in modes 4 and 11 (p < 0.001 and p = 0.003 respectively). PRVI had significant effects on mode 1 (p < 0.001, PE = − 0.15), mode 2 (p = 0.007, PE = 0.3), mode 4 (p < 0.001, PE = − 0.4), mode 7 (p = 0.02, PE = − 0.27) and mode 8 (p = 0.002, PE = 0.3). Combining these scores, PRVI accounted for 12.3% of the total shape variation. A positive (negative) PE indicates increasing (decreasing) morphometric scores. For example, PRVI was associated with a negative coefficient for mode 1, meaning that as PRVI increased, patients had larger hearts. Figure 4a shows the regional shape differences associated with increasing PRVI computed with the regression model.
An animation of these shape differences can be found in Additional file 2. Figure 4b shows the quantitative differences in shape of the endocardial surfaces (expressed as a displacement in mm) relative to the mean PRVI of 25.5 ml/m2, at both ED and ES. Displacements outward from the LV are shown in red, displacements inwards towards the LV are shown in blue. The dilation of the RV outflow tract and apex are confirmed, as well as diminution of the LV. The associations between shape and PRVI were similar if height and weight were replaced by BSA in the predictor variables, and also if height and weight and age were replaced by BSA. The results were therefore robust to different methods of accounting for body habitus.
Results of the univariate regression models performed with specific shape features are summarized in Table 3. Significant correlations with PRVI are shown with corresponding PE (p < 0.05). In the RV, as PRVI increased, RV size increased, the RV apex dilated, and the tricuspid annulus tilted. In the LV, the distance between the LV free wall and LV septum reduced, resulting in a flattening of the LV. In addition, larger PRVI was associated with a systolic septal motion towards the RV (paradoxical septal motion), which is consistent with RV overload and decreased LV dimension at ED.
To illustrate and quantify the paradoxical septal motion, Fig. 5 shows changes in systolic motion of the endocardial surfaces, expressed as a displacement in mm from ED to ES, for a PRVI of 5 ml/m2, and a PRVI of 40 ml/m2. Displacements outward from the LV are shown in red, displacements inwards towards the LV are shown in blue.
LS, CS, and RS for LV and RV are shown in Table 4 and significant correlations with PRVI are shown with corresponding PE. Although all RV strains were correlated with PRVI, for the LV only RS was significantly correlated with PRVI.
Adults with rTOF are exposed to ongoing physiologic sequelae of surgical interventions, causing remodeling which affects both the LV and RV. This cross-sectional study used atlas-based shape analysis to enable precise quantification of complex shape differences associated with pulmonary regurgitation. We used linear regression methods to control for confounding variables, and quantified independent effects associated with PRVI. PRVI explained more shape variance than any principal component except for the first (size mode), highlighting its importance in biventricular remodeling. Increasing PRVI was associated with dilation of the RV apex and outflow tract, together with a diminution of the LV and flattening in which the lateral free wall to septum distance decreased. Increasing PRVI was also associated with increased RV RS, LS and CS, as well as an increase in LV RS.
In the INDICATOR cohort study [44, 45], preoperative risk factors for postoperative occurrence of death or sustained tachycardia included RV systolic dysfunction, age at pulmonary valve replacement > 28, and elevated RV mass-to-volume ratio. Current guidelines recommend pulmonary valve replacement in asymptomatic patients if progressive RV dilation to RV end-systolic volume index (ESVI) > 80 ml/m2 and/or RV end-diastolic volume index (EDVI) > 160 ml/m2 is observed [46, 47]. However, optimal timing of replacement remains a challenge particularly in asymptomatic patients and there is wide qualitative institutional bias in practice. Owing to the continuity of muscle fibers between the LV and RV, a confined pericardial space and a relatively non distensible pericardium, volumes and function of one ventricle directly affect the other ventricle. Increased LVESVI and decreased LVEF have been associated with adverse outcomes in the INDICATOR cohort [45, 48] and other studies [49, 50]. Smaller LV diameters were previously associated with RV dilation and the severity of pulmonary regurgitation , as in the current study. PR has been associated with decreasing contribution of longitudinal shortening to the RV ejection and increasing lateral pumping, which results in larger volume changes and septal motion towards the RV , as seen in our data (Additional file 2). RV shortening has also been associated with exercise capacity in adult rTOF . Fernandes et al.  also suggested that RV apical dilation alters LV geometry, thereby decreasing LV stroke volume. Alterations in LV geometry may also affect the mitral valve, contributing to the relatively high proportion of rTOF patients with mitral valve prolapse . Our methods enable precise quantification of these relationships, including new features of LV flattening and motion, which highlight the importance of understanding the contributions of both RV and LV shape and their interactions.
The atlas captured regional shape differences and also reproduced global mass and volume with acceptable bias and precision (Table 2). Similar errors were previously obtained for 4329 cases from UK Biobank . There was consistent overestimation of LV and RV mass in the model relative to manual contour estimations, likely due to differences in the apex and base. The biventricular model numerically integrated the volume contained within the 3D surfaces using information from both long and short axis images, whereas contour slice summation was performed only on short axis images, in which the contours are often difficult to estimate at the apex and base.
Shape models and atlas-based methods have previously shown RV shape changes associated with PR [12, 17, 36, 53]. These studies demonstrated RV dilation, basal bulging, tricuspid tilting, and apical dilation associated with pulmonary regurgitation, in agreement with the current study (Additional file 2, posterior view). This tilting may be a specific effect of RV dilation that contributes to tricuspid regurgitation, which develops in one third of rTOF patients . Our biventricular analysis enabled extension of these methods to quantify the associations between biventricular shape and function and amount of PR.
In the German Competence Network rTOF study, reduced LV CS and RV LS were independent predictors of adverse outcomes . Decreased LV LS has also been associated with adverse events . In our study, we found that LV RS as well as RV RS, CS and LS increased as pulmonary regurgitation increased. A compensatory increase in RV free wall RS in response to PR and volume overload was also observed in [51, 57], which may be due to increased RV SV and a higher preload. The increase in LV RS found the current study may be due to the inclusion of the septum into the LV RS measurement.
Our study also found a relationship between RV longitudinal shortening and the amount of regurgitation, in agreement with Ylitalo et al. . However, conclusions regarding the effects of PR on longitudinal shortening have varied between studies. In pediatric populations, Ylitalo et al.  found the PR led to an increase in RV LS, on the other hand, Eyskens et al.  and Ouyang et al.  found an inverse correlation between PR and RV LS measurements. Studies in adult rTOF patients did not show any variation in RV and LV strain in relation to PR [16, 22, 43]. However, reduced contractile function in relation to PR was demonstrated in  but did not have any influence on RV strain measurements. Our study contained mostly younger adult patients, but decreasing strain beyond a certain point may occur with more severe disease.
The type of primary repair has also been shown to affect RV shape . Zaidi et al.  found differences in volume and regional curvatures in patients with transannular repair vs those with pulmonary valve preserving repair. In our study, 78.4% had transannular patch primary repair. Also, our findings of reduced septal-lateral dimension are consisted with increased curvature in the free wall.
Although data sharing mechanisms provided through the CAP enable merging of data from different institutions, access to clinical information was limited. Data on ethnicity was not available. Indices of RV pressure loading were not available, and investigation of pressure effects requires further study, ideally using measurements from catheter recordings. The interpretation of strain measurements in the rTOF population is also difficult due to the competing effects of increased force of contraction with higher preload, and decompensating RV contractility. Another limitation of the study is the use of cross-sectional data. The temporal evolution of ventricular dysfunction in relation to regurgitation should be performed in longitudinal studies. Finally, our study did not quantify dyssynchrony (due to right bundle branch block) or late gadolinium enhancement patterns, which may be have important effects on biventricular function [61, 62].
Future application of these methods will examine quantitative shape changes in relation to interventions and clinically significant metrics, including pulmonary valve replacement and adverse outcomes. The distillation of complex shape features into a small number of morphometric scores allow precise quantification of mechanistic effects linking structural and functional alterations to a variety of clinical measures, including rhythm disturbances, functional score, restrictive RV physiology, and exercise capacity. Relationships between morphometric scores and future adverse outcomes will enable examination of mechanisms of developing risk. Fully automated analysis including quantification of contours, biventricular shape model customization, and computation of scores relative to reference populations, is now possible using machine learning and AI methods.
Biventricular morphometric relationships with important clinical features can be quantified in rTOF patients. Increasing PR is associated with LV diminution and septal flattening, in conjunction with RV dilation, especially at the base and apex, as well as an increase in RV systolic strain. Specific and quantifiable biventricular shape characteristics, revealed in biventricular models of rTOF, can be used to quantitatively evaluate RV and LV dysfunction in these patients. This study provides a basis for including quantitative shape analysis into the clinician’s toolbox, providing simple scores which express complex features that experienced clinicians may well recognize but which have been difficult to report until now.
Availability of data and materials
Data and shape models are available from www.cardiacatlas.org.
Body surface area
Cardiac atlas project
Cardiovascular magnetic resonance
End-diastolic volume index
End-systolic volume index
Left ventricle/left ventricular
Left ventricular ejection fraction
Left ventricular mass
Left ventricular mass index
Principal component analysis
Pulmonary regurgitant volume index
Repaired tetralogy of Fallot
Right ventricle/right ventricular
Right ventricular mass
Right ventricular mass index
Townsley MM, Windsor J, Briston D, Alegria J, Ramakrishna H. Tetralogy of Fallot: perioperative management and analysis of outcomes. J Cardiothorac Vasc Anesth. 2019;33(2):556–65.
Smith CA, McCracken C, Thomas AS, Spector LG, St Louis JD, Oster ME, Moller JH, Kochilas L. Long-term outcomes of tetralogy of Fallot: a study from the pediatric cardiac care consortium. JAMA Cardiol. 2019;4(1):34–41.
Frigiola A, Redington AN, Cullen S, Vogel M. Pulmonary regurgitation is an important determinant of right ventricular contractile dysfunction in patients with surgically repaired tetralogy of Fallot. Circulation. 2004;110(11 Suppl 1):II153-157.
Geva T. Repaired tetralogy of Fallot: the roles of cardiovascular magnetic resonance in evaluating pathophysiology and for pulmonary valve replacement decision support. J Cardiovasc Magn Reson. 2011;13:9.
Nieminen HP, Jokinen EV, Sairanen HI. Causes of late deaths after pediatric cardiac surgery: a population-based study. J Am Coll Cardiol. 2007;50(13):1263–71.
Probst J, Diller GP, Reinecke H, Leitz P, Frommeyer G, Orwat S, Vormbrock J, Radke R, de Torres AF, Kaleschke G, et al. Prevention of sudden cardiac death in patients with Tetralogy of Fallot: Risk assessment and long term outcome. Int J Cardiol. 2018;269:91–6.
Walsh EP. Sudden death in adult congenital heart disease: risk stratification in 2014. Heart Rhythm. 2014;11(10):1735–42.
Oechslin EN, Harrison DA, Harris L, Downar E, Webb GD, Siu SS, Williams WG. Reoperation in adults with repair of tetralogy of Fallot: indications and outcomes. J Thorac Cardiovasc Surg. 1999;118(2):245–51.
Therrien J, Siu SC, Harris L, Dore A, Niwa K, Janousek J, Williams WG, Webb G, Gatzoulis MA. Impact of pulmonary valve replacement on arrhythmia propensity late after repair of tetralogy of Fallot. Circulation. 2001;103(20):2489–94.
Quail MA, Frigiola A, Giardini A, Muthurangu V, Hughes M, Lurz P, Khambadkone S, Deanfield JE, Tsang V, Taylor AM. Impact of pulmonary valve replacement in tetralogy of Fallot with pulmonary regurgitation: a comparison of intervention and nonintervention. Ann Thorac Surg. 2012;94(5):1619–26.
Geva T. Indications for pulmonary valve replacement in repaired tetralogy of Fallot: the quest continues. Circulation. 2013;128(17):1855–7.
Mansi T, Voigt I, Leonardi B, Pennec X, Durrleman S, Sermesant M, Delingette H, Taylor AM, Boudjemline Y, Pongiglione G, et al. A statistical model for quantification and prediction of cardiac remodelling: application to tetralogy of Fallot. IEEE Trans Med Imaging. 2011;30(9):1605–16.
Sheehan FH, Kilner PJ, Sahn DJ, Vick GW 3rd, Stout KK, Ge S, Helbing WA, Lewin M, Shurman AJ, Buechel EV, et al. Accuracy of knowledge-based reconstruction for measurement of right ventricular volume and function in patients with tetralogy of Fallot. Am J Cardiol. 2010;105(7):993–9.
Nyns EC, Dragulescu A, Yoo SJ, Grosse-Wortmann L. Evaluation of knowledge-based reconstruction for magnetic resonance volumetry of the right ventricle in tetralogy of Fallot. Pediatr Radiol. 2014;44(12):1532–40.
Zervan K, Male C, Benesch T, Salzer-Muhar U. Ventricular interaction in children after repair of tetralogy of Fallot: a longitudinal echocardiographic study. Eur J Echocardiogr. 2009;10(5):641–6.
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(5):494–503.
Agarwal A, Harris IS, Mahadevan VS, Foster E. Coexistence of abnormal systolic motion of mitral valve in a consecutive group of 324 adult Tetralogy of Fallot patients assessed with echocardiography. Open Heart. 2016;3(2):000518.
Ylitalo P, Lehmonen L, Lauerma K, Holmstrom M, Pitkanen-Argillander O, Jokinen E. Severe pulmonary regurgitation in adolescents with tetralogy of Fallot leads to increased longitudinal strain. MAGMA. 2020;33(2):309–16.
Khalaf A, Tani D, Tadros S, Madan S. Right- and left-ventricular strain evaluation in repaired pediatric Tetralogy of Fallot patients using magnetic resonance tagging. Pediatr Cardiol. 2013;34(5):1206–11.
Ylitalo P, Jokinen E, Lauerma K, Holmstrom M, Pitkanen-Argillander OM. Additional mechanism for left ventricular dysfunction: chronic pulmonary regurgitation decreases left ventricular preload in patients with tetralogy of Fallot. Cardiol Young. 2018;28(2):208–13.
Takayasu H, Takahashi K, Takigiku K, Yasukochi S, Furukawa T, Akimoto K, Kishiro M, Shimizu T. Left ventricular torsion and strain in patients with repaired tetralogy of Fallot assessed by speckle tracking imaging. Echocardiography. 2011;28(7):720–9.
Tretter JT, Redington AN. The forgotten ventricle? The left ventricle in right-sided congenital heart disease. Circ Cardiovasc Imaging. 2018;11(3):007410.
Broberg CS, Aboulhosn J, Mongeon FP, Kay J, Valente AM, Khairy P, Earing MG, Opotowsky AR, Lui G, Gersony DR, et al. Prevalence of left ventricular systolic dysfunction in adults with repaired tetralogy of Fallot. Am J Cardiol. 2011;107(8):1215–20.
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(6):1068–74.
Davlouros PA, Kilner PJ, Hornung TS, Li W, Francis JM, Moon JC, Smith GC, Tat T, Pennell DJ, Gatzoulis MA. Right ventricular function in adults with repaired tetralogy of Fallot assessed with cardiovascular magnetic resonance imaging: detrimental role of right ventricular outflow aneurysms or akinesia and adverse right-to-left ventricular interaction. J Am Coll Cardiol. 2002;40(11):2044–52.
Tweddell JS, Simpson P, Li SH, Dunham-Ingle J, Bartz PJ, Earing MG, Pelech AN. Timing and technique of pulmonary valve replacement in the patient with tetralogy of Fallot. Semin Thorac Cardiovasc Surg Pediatr Card Surg Annu. 2012;15(1):27–33.
Tobler D, Crean AM, Redington AN, Van Arsdell GS, Caldarone CA, Nanthakumar K, Stambach D, Dos L, Wintersperger BJ, Oechslin EN, et al. The left heart after pulmonary valve replacement in adults late after tetralogy of Fallot repair. Int J Cardiol. 2012;160(3):165–70.
Mauger C, Gilbert K, Lee AM, Sanghvi MM, Aung N, Fung K, Carapella V, Piechnik SK, Neubauer S, Petersen SE, et al. Right ventricular shape and function: cardiovascular magnetic resonance reference morphology and biventricular risk factor morphometrics in UK Biobank. J Cardiovasc Magn Reson. 2019;21(1):41.
Fonseca CG, Backhaus M, Bluemke DA, Britten RD, Chung JD, Cowan BR, Dinov ID, Finn JP, Hunter PJ, Kadish AH, et al. The Cardiac Atlas Project–an imaging database for computational modeling and statistical atlases of the heart. Bioinformatics. 2011;27(16):2288–95.
Heiberg E, Sjogren J, Ugander M, Carlsson M, Engblom H, Arheden H. Design and validation of Segment–freely available software for cardiovascular image analysis. BMC Med Imaging. 2010;10:1.
Sinclair M, Bai W, Puyol-Anton E, Oktay O, Rueckert D, King AP. Fully automated segmentation-based respiratory motion correction of multiplanar cardiac magnetic resonance images for large-scale datasets. Med Image Comp Comp-Assisted Intervention. 2017:332–340.
Jolliffe I. Principal component analysis. John Wiley & Sons; 2005.
Gorter TM, van Melle JP, Freling HG, Ebels T, Bartelds B, Pieper PG, Berger RM, van Veldhuisen DJ, Willems TP. Pulmonary regurgitant volume is superior to fraction using background-corrected phase contrast MRI in determining the severity of regurgitation in repaired tetralogy of Fallot. Int J Cardiovasc Imaging. 2015;31(6):1169–77.
Wald RM, Redington AN, Pereira A, Provost YL, Paul NS, Oechslin EN, Silversides CK. Refining the assessment of pulmonary regurgitation in adults after tetralogy of Fallot repair: should we be measuring regurgitant fraction or regurgitant volume? Eur Heart J. 2009;30(3):356–61.
Kang IS, Redington AN, Benson LN, Macgowan C, Valsangiacomo ER, Roman K, Kellenberger CJ, Yoo SJ. Differential regurgitation in branch pulmonary arteries after repair of tetralogy of Fallot: a phase-contrast cine magnetic resonance study. Circulation. 2003;107(23):2938–43.
Sheehan FH, Ge S, Vick GW 3rd, Urnes K, Kerwin WS, Bolson EL, Chung T, Kovalchin JP, Sahn DJ, Jerosch-Herold M, et al. Three-dimensional shape analysis of right ventricular remodeling in repaired tetralogy of Fallot. Am J Cardiol. 2008;101(1):107–13.
Lopez-Candales A, Dohi K, Iliescu A, Peterson RC, Edelman K, Bazaz R. An abnormal right ventricular apical angle is indicative of global right ventricular impairment. Echocardiography. 2006;23(5):361–8.
Wehner GJ, Jing L, Haggerty CM, Suever JD, Chen J, Hamlet SM, Feindt JA, Dimitri Mojsejenko W, Fogel MA, Fornwalt BK. Comparison of left ventricular strains and torsion derived from feature tracking and DENSE CMR. J Cardiovasc Magn Reson. 2018;20(1):63.
Cowan BR, Peereboom SM, Greiser A, Guehring J, Young AA. Image feature determinants of global and segmental circumferential ventricular strain from cine CMR. JACC Cardiovasc Imaging. 2015;8(12):1465–6.
Bolson EL, Kliman S, Sheehan FH, Dodge HT. Left ventricular segmental wall motion: a new method using local direction information. Comp Cardiol. 1980:245–248.
R: A language and environment for statistical computing. http://www.R-project.org/.
Muzzarelli S, Ordovas KG, Cannavale G, Meadows AK, Higgins CB. Tetralogy of Fallot: impact of the excursion of the interventricular septum on left ventricular systolic function and fibrosis after surgical repair. Radiology. 2011;259(2):375–83.
Weyman AE, Wann S, Feigenbaum H, Dillon JC. Mechanism of abnormal septal motion in patients with right ventricular volume overload: a cross-sectional echocardiographic study. Circulation. 1976;54(2):179–86.
Valente AM, Gauvreau K, Assenza GE, Babu-Narayan SV, Evans SP, Gatzoulis M, Groenink M, Inuzuka R, Kilner PJ, Koyak Z, et al. Rationale and design of an International Multicenter Registry of patients with repaired tetralogy of Fallot to define risk factors for late adverse outcomes: the INDICATOR cohort. Pediatr Cardiol. 2013;34(1):95–104.
Geva T, Mulder B, Gauvreau K, Babu-Narayan SV, Wald RM, Hickey K, Powell AJ, Gatzoulis MA, Valente AM. Preoperative predictors of death and sustained ventricular tachycardia after pulmonary valve replacement in patients with repaired tetralogy of Fallot enrolled in the INDICATOR cohort. Circulation. 2018;138(19):2106–15.
Baumgartner H, De Backer J, Babu-Narayan SV, Budts W, Chessa M, Diller GP, Lung B, Kluin J, Lang IM, Meijboom F et al. 2020 ESC Guidelines for the management of adult congenital heart disease. Eur Heart J 2020.
Stout KK, Daniels CJ, Aboulhosn JA, Bozkurt B, Broberg CS, Colman JM, Crumb SR, Dearani JA, Fuller S, Gurvitz M, et al. 2018 AHA/ACC guideline for the management of adults with congenital heart disease: executive summary: a report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines. J Am Coll Cardiol. 2019;73(12):1494–563.
Valente AM, Gauvreau K, Assenza GE, Babu-Narayan SV, Schreier J, Gatzoulis MA, Groenink M, Inuzuka R, Kilner PJ, Koyak Z, et al. Contemporary predictors of death and sustained ventricular tachycardia in patients with repaired tetralogy of Fallot enrolled in the INDICATOR cohort. Heart. 2014;100(3):247–53.
Ghai A, Silversides C, Harris L, Webb GD, Siu SC, Therrien J. Left ventricular dysfunction is a risk factor for sudden cardiac death in adults late after repair of tetralogy of Fallot. J Am Coll Cardiol. 2002;40(9):1675–80.
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(2):211–6.
Stephensen S, Steding-Ehrenborg K, Munkhammar P, Heiberg E, Arheden H, Carlsson M. The relationship between longitudinal, lateral, and septal contribution to stroke volume in patients with pulmonary regurgitation and healthy volunteers. Am J Physiol-Heart Circ Physiol. 2014;306:H895-903.
Leng S, Tan RS, Guo J, Chai P, Zhang G, Teo L, et al. Cardiovascular magnetic resonance-assessed fast global longitudinal strain parameters add diagnostic and prognostic insights in right ventricular volume and pressure loading disease conditions. J Cardiovasc Magn Reson. 2021;23(1):1–13.
Zhong L, Gobeawan L, Su Y, Tan JL, Ghista D, Chua T, Tan RS, Kassab G. Right ventricular regional wall curvedness and area strain in patients with repaired tetralogy of Fallot. Am J Physiol Heart Circ Physiol. 2012;302(6):H1306-1316.
Mahle WT, Parks WJ, Fyfe DA, Sallee D. Tricuspid regurgitation in patients with repaired Tetralogy of Fallot and its relation to right ventricular dilatation. Am J Cardiol. 2003;92(5):643–5.
Orwat S, Diller GP, Kempny A, Radke R, Peters B, Kuhne T, Boethig D, Gutberlet M, Dubowy KO, Beerbaum P, et al. Myocardial deformation parameters predict outcome in patients with repaired tetralogy of Fallot. Heart. 2016;102(3):209–15.
Diller GP, Kempny A, Liodakis E, Alonso-Gonzalez R, Inuzuka R, Uebing A, Orwat S, Dimopoulos K, Swan L, Li W, et al. Left ventricular longitudinal function predicts life-threatening ventricular arrhythmia and death in adults with repaired tetralogy of fallot. Circulation. 2012;125(20):2440–6.
Kovacs A, Lakatos B, Tokodi M, Merkely B. Right ventricular mechanical pattern in health and disease: beyond longitudinal shortening. Heart Fail Rev. 2019;24(4):511–20.
Eyskens B, Brown SC, Claus P, Dymarkowski S, Gewillig M, Bogaert J, Mertens L. The influence of pulmonary regurgitation on regional right ventricular function in children after surgical repair of tetralogy of Fallot. Eur J Echocardiogr. 2010;11(4):341–5.
Ouyang R, Leng S, Sun A, Wang Q, Hu L, Zhao X, et al. (2021). Detection of persistent systolic and diastolic abnormalities in asymptomatic pediatric repaired tetralogy of Fallot patients with preserved ejection fraction: a CMR feature tracking study. Eur Radiol. 1–13.
Zaidi SJ, Cossor W, Singh A, Maffesanti F, Kawaji K, Woo J, Mor-Avi V, Roberson DA, Kutty S, Patel AR. Three-dimensional analysis of regional right ventricular shape and function in repaired tetralogy of Fallot using cardiovascular magnetic resonance. Clin Imaging. 2018;52:106–12.
Ortega M, Triedman JK, Geva T, Harrild DM. Relation of left ventricular dyssynchrony measured by cardiac magnetic resonance tissue tracking in repaired tetralogy of Fallot to ventricular tachycardia and death. Am J Cardiol. 2011;107(10):1535–40.
Ghonim S, Ernst S, Keegan J, Giannakidis A, Spadotto V, Voges I, Smith GC, Boutsikou M, Montanaro C, Wong T et al. 3D late gadolinium enhancement cardiovascular magnetic resonance predicts inducibility of ventricular tachycardia in adults with repaired tetralogy of Fallot. Circ Arrhythm Electrophysiol 2020.
Funding was provided by the National Institutes of Health (USA) R01HL121754 and AHA Grant 19AIML35120034. SG acknowledges National Institutes of Health NHLBI 1T32HL105373. AY acknowledges Health Research Council of New Zealand grant 17/234. CM acknowledges New Zealand Heart Foundation Grant 1695. RC and TH acknowledge the support of Inria-UTSW Medical Center Dallas Associated Team TOFMOD, and RC additionally the support of the Ministry of Health of the Czech Republic (project No. NV19-08-00071). LZ acknowledges National Medical Research Council of Singapore Grant (NMRC/OFIRG/0018/2016).
Ethics approval and consent to participate
De-identified data were contributed to the Cardiac Atlas Project from multiple clinical centers with approval from local institutional review boards.
Consent for publication
No identifiable individual information was included in this study.
ADM and JHO are co-founders of and have an equity interest in Insilicomed, Inc., and he serves on the scientific advisory board. Some of their research grants, including those acknowledged here, have been identified for conflict of interest management based on the overall scope of the project and its potential benefit to Insilicomed, Inc. The authors are required to disclose this relationship in publications acknowledging the grant support, however the research subject and findings reported here did not involve the company in any way and have no relationship whatsoever to the business activities or scientific interests of the company. The terms of this arrangement have been reviewed and approved by the University of California San Diego in accordance with its conflict of interest policies. The other authors have no competing interests to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Animation of the first four PCA components. Left: anterior view. Right: posterior view. “s” indicates the number of standard variations.
Shape changes due to pulmonary regurgitation. Left: anterior view. Right: posterior view. “s” indicates the number of standard variations.
Each shape model was subdivided into 5,810 points on the LV and RV endocardial and epicardial surfaces as well as the four valves. These points (x,y,z) coordinates at ED and ES were concatenated into a single shape vector with 34,860 elements. The shape vectors from each patient were assembled into a data matrix (one row per case). The PCA shape modes were calculated from the covariance matrix and morphometric scores for each patient were calculated by projection onto the shape modes and then standardized to provide z-scores. A multivariate linear regression was performed to estimate the contribution of predictor variables to each principal component score:
where T are the morphometric scores, A is a matrix of predictor variables, B is a matrix of regression coefficients and R is a residual matrix. The predictor variables comprised sex, age, height, weight, and PRVI. For each predictor variable, the associated shape changes were visualized by calculating the contribution to the morphometric scores for that predictor.
About this article
Cite this article
Mauger, C.A., Govil, S., Chabiniok, R. et al. Right-left ventricular shape variations in tetralogy of Fallot: associations with pulmonary regurgitation. J Cardiovasc Magn Reson 23, 105 (2021). https://doi.org/10.1186/s12968-021-00780-x
- Cardiovascular magnetic resonance
- Ventricular function
- Myocardial deformation
- Tetralogy of Fallot