Abstract
Right atrial (RA) mechanics have been studied infrequently in children in the past due to technical constraints. With the advent of strain imaging, RA physiology can now be studied in greater detail. The principal aim of this study was to describe functional changes in right heart mechanics of children with idiopathic pulmonary arterial hypertension (PAH), by using new applications of RA strain. In this retrospective study, we evaluated RA mechanics of 20 patients (age range = 3–23 years) with PAH and 18 control patients. RA longitudinal strain (RALS) and longitudinal displacement (LD) were calculated by speckle-tracking echocardiography. RALS was plotted against LD, producing a characteristic strain-displacement (S-D) loop. Standard indices of right heart function and right heart catheterization data were obtained. Patients were clinically subdivided into “compensated” and “decompensated” PAH. A chart review was performed to identify patients who subsequently developed adverse outcomes, including death, awaiting or received lung and/or heart transplantation. RALS was significantly lower in decompensated PAH compared with both controls and compensated PAH. Area enclosed by S-D loops differed significantly between the compensated and decompensated PAH subgroups (5.33 [3.90–9.44] versus 1.83 [1.17–2.36],
Keywords
Introduction
Standard and advanced echocardiographic techniques continue to represent essential elements in the evaluation of patients with pulmonary hypertension (PAH). This disease process is characterized by increased pulmonary vascular resistance resulting in chronic right ventricular (RV) pressure overload, myocardial remodeling, leading to RV hypertrophy, dilation, and eventual RV failure.1,2 The echocardiographic evaluation of a patient with suspected or established PAH may reveal characteristic findings such as tricuspid regurgitation (TR), right atrial (RA) and RV dilation, RV hypertrophy, abnormal shape and function of the interventricular septum, and dilation of the main pulmonary artery. 3 However, many of these findings may not be present until quite late in the course of the disease, leading to the search for indices that may identify decompensation at earlier stages.
With the advent of speckle-tracking echocardiography (STE), it is now possible to quantify RA deformation over the course of the entire cardiac cycle from a single beat. Strain imaging has been applied in patients with PAH and right heart dysfunction. 4 In adults with chronic systolic heart failure, RA peak longitudinal strain (RALS) showed a strong inverse correlation with invasively measured systolic pulmonary artery pressure (PAP) and may represent an adjunct to invasive means of hemodynamic assessment. 5 In adults with PAH, RALS was noted to be low6–8 and it negatively correlated with RA pressure. 7
However, in the pediatric population the RA represents the “forgotten” chamber of the heart and studies on RA physiology and function in children are quite limited. This is especially true in pediatric PAH. STE can generate both strain and displacement data from the same heartbeat, which can simultaneously assess regional and global myocardial deformation of a cardiac chamber, including mechanical dyssynchrony. Simultaneous measurement of strain and displacement have been studied in other chambers of the heart, including the pressure-overloaded RV in adults. Notably, these studies in adults have shown that measurement of RV displacement and strain may facilitate reliable and quantitative estimation of RV function.9,10 However, when the RV is significantly dilated, it may not be possible to include all the walls of the RV in the same echocardiographic acquisition window. The apex of the RV may often be truncated in acquisition windows, particularly in retrospective studies. In contrast, the RA is a simpler geometric structure and invariably well-visualized from an apical four-chamber view. Therefore, RA mechanics may provide a simpler and more reliable measure of RV decompensation. Based on these findings, we were enticed by the feasibility of describing a strain-displacement (S-D) relationship in the RA of children with PAH.
The objective of our study is to characterize RA function throughout a single cardiac cycle in healthy children and those with PAH, specifically looking at differences in RA strain and the strain-displacement relationship. We hypothesize that evaluation of RA strain in combination with RA displacement may provide useful insight into the function of both the normal and failing RA in PAH. These indices may identify patients at risk for future adverse outcomes.
Methods
Patient population
Patients followed at The Children’s Hospital of Philadelphia with a diagnosis of PAH without post-tricuspid valve shunt were recruited for this study. Twenty children classified as having PAH due to idiopathic etiology (Nice Classification type 1.1), drug or toxin (Nice Classification type 1.3), congenital heart disease but without current shunt (Nice Classification type 1.4.4), or hematologic diagnosis (Nice Classification type 5.1) 3 were included in this study. We excluded PAH patients with significant defects in the atrial septum and patients with left ventricular (LV) dysfunction or pulmonary vein obstruction. Control individuals in this study were otherwise healthy pediatric patients referred for echocardiography for indications such as non-cardiac chest pain, vasodepressor syncope, and functional heart murmurs. All control participants were found to have structurally normal hearts without evidence of PAH by echocardiogram. This retrospective, descriptive study was approved by the Institutional Review Board of the Children’s Hospital of Philadelphia.
Strain imaging
Echocardiograms were obtained using Philips iE33 ultrasound machines (Philips Medical Systems, Andover, MA, USA). Two-dimensional apical four-chamber views were obtained and uploaded for analysis by a vendor-independent STE software (Tomtec Cardiac Performance version 1.2.4.2, TomTec Imaging Systems, Munich, Germany). The frame rate of the acquired images exceeded 60 frames/s. Manual tracing along the RA endocardial border was performed with termination of the tracing 0.5 cm above the atrioventricular junction, to limit influence of tricuspid annular motion. The software divided the right atrium into six segments, with longitudinal strain and displacement obtained for each segment over a complete cardiac cycle. The global RA strain and longitudinal displacement (LD) for each patient were averages of the values obtained from these six segments. By convention, displacement away from the centroid of the cardiac chamber is indicated by negative values and may be considered as a surrogate for increase in RA volume.
Right-heart function
To better assess the utility of strain imaging in evaluation of right-sided function in individuals with PAH, the patients were subdivided into two groups: decompensated and compensated. “Decompensated PAH” was defined as the subgroup of patients with functional and laboratory evidence of significant right heart dysfunction at time of the echocardiogram. Specifically, the participants in this subgroup were in World Health Organization (WHO) functional class III or IV. 3 Based on adult and pediatric studies on brain natriuretic peptide (BNP), which showed strikingly decreased survival in idiopathic PAH, we chose a cut-off value for BNP of 180 pg/mL.11,12 In our study, all decompensated PAH patients exhibited a BNP value of > 180 pg/mL. “Compensated PAH” was defined as individuals in WHO functional class I or II or a BNP < 180 pg/mL, or both.
In addition to RA mechanics and strain values, more conventional indices such as RV peak strain, RV area change (RVAC), and tricuspid annular plane systolic excursion (TAPSE) were recorded and used as indices of RV function. To compare TAPSE in participants of different ages and sizes, the values were indexed to the body surface area (BSA). We also compared percent-predicted TAPSE values using established reference ranges.13–15 In individuals with PAH, BNP values obtained on or as close as possible to the date of the echocardiogram were evaluated for this study: BNP levels were obtained on 16 patients on the same day as the echocardiogram, the remaining four were obtained within three months.
Six-minute walk distances (6MWD) were obtained on or as close as possible to the date of the echocardiogram. Of the 20 PAH participants, 19 had 6MWD measurements recorded. Twelve measurements took place on the same day as the echocardiogram and four underwent the 6-min walking text (6MWT) within three months of their echocardiograms. For the remaining three individuals, the 6MWT was performed approximately six months after their echocardiograms. To determine a percent-predicted 6MWT value for each participant, the observed 6MWD was divided by the distance expected for their sex and age. Several reference ranges are available for pediatric 6MWT times based on their ages.16–19 Using the appropriate reference ranges, we calculated predicted values for 6MWTs in our cohort.
Strain-displacement loops
Construction of S-D loops was undertaken by plotting RALS along the y-axis and LD along the x-axis. A typical S-D loop is illustrated in Fig. 2. The area enclosed by the S-D loop was calculated using the trapezoidal area of polygon formula, where the loop was broken into mini-trapezoids and the area of the loop was equal to the sum of all the trapezoids. The total area of the loop was the sum of [x(i+1)–x(i)][y(i)+y(i+1))/2 for i points of a polygon, where x(1), y(1) = x(i + 1), y(i + 1). 20 The unit for this measurement is %.mm/cm and is simply referred to as unit (U) in this paper. In order to compare the right atria of children of different ages, the LD was indexed to the BSA. We used an additional method of indexing LD by dividing it with the largest RA dimension with the tricuspid valve closed. 21
Cardiac catheterization data
We evaluated right heart catheterization data that were obtained within six months of the echocardiogram. Mean pulmonary artery pressure (mPAP) and indexed pulmonary vascular resistance (PVRi) were recorded for each individual.
Adverse outcomes
Based on > 2-year follow-up data, we identified participants who later developed adverse outcomes associated with their PAH. Adverse outcomes included death, awaiting or receiving lung and/or heart transplantation.
Statistical analysis
Continuous variables were non-normally distributed and are reported as median with interquartile ranges (IQR) noted parenthetically. Comparisons between the variables were performed using the Mann–Whitney U test. Non-parametric receiver operating characteristic (ROC) analysis was performed to test the sensitivity and specificity of the different echocardiographic indices, for the outcome described above. ROC curves were plotted and the areas under the curves were measured. For all significance testing, a difference was considered significant if
Results
Demographics
Demographic, clinical, and echocardiographic characteristics.
Data are expressed as mean ± standard deviation.
BSA, body surface area; FC, functional class; WHO, World Health Organization.
Comparison of control and PAH population
The control group differed significantly from the PAH group in several RV parameters, including percent-predicted TAPSE, and RV to LV ratio in diastole. RV longitudinal strain was lower in the PAH group versus controls (–15.7 ± 3.6% vs –23.3 ± 3.8%, Mean PAP shows an inverse correlation with RA peak longitudinal strain. Right heart data. Data are expressed as median with first and third interquartile range in parenthesis. RALS, right atrial longitudinal strain; Loop areaBSA, S-D loop area indexed for body surface area; Loop arealength, S-D loop area indexed for right atrial length; TAPSEBSA, tricuspid annular plane systolic excursion, indexed for body surface area; %p TAPSE, percent-predicted TAPSE; RV:LVdiastole, ratio of right ventricular to left ventricular diameter at end-diastole; BNP, brain natriuretic peptide; 6MWD, 6-min walk distance; %p 6MWD, percent-predicted 6-min walk distance; mPAP, mean pulmonary artery pressure. Typical RA S-D loop. The S-D loop is generated by plotting the RA LD along the x-axis and the longitudinal strain along the y-axis. This results in a characteristic elliptical loop in the counter-clockwise direction, starting with tricuspid valve closure (TVC). The area of the S-D loop was constituted by reservoir, conduit, and contractile phases of RA function. The reservoir phase continues from TVC to tricuspid valve opening (TVO), after which starts the conduit phase, followed by a brief contractile phase. The contractile phase in children is not as pronounced as in adults and therefore a distinct incisura is not well-detected in the waveforms of pediatric patients.

The S-D loop had a characteristic elliptical shape and counterclockwise rotation, and was reflective of the reservoir, conduit, and contractile phases of RA function (Fig. 2). The reservoir phase began with tricuspid valve closure (TVC) and continued until the tricuspid valve opened (TVO). The conduit phase began at TVO and was followed by a brief contractile (booster pump) phase. The contractile phase in children is not as pronounced as in adults and manifested only a subtle incisura on the strain curves. Therefore, the contractile (booster pump) phase was not calculated separately in our study. When indexed for RA lengths, the areas enclosed by S-D loops were significantly decreased in the group consisting of all PAH patients compared to control participants (3.8 [1.9–8.5] versus 7.9 [5.9–9.9],
Comparison of compensated and decompensated PAH groups
The entire PAH population was divided into compensated and decompensated groups to determine if RA indices such as RALS and S-D loops varied between these subgroups and whether a point of decline in RA function could be detected in these individuals. The data are summarized in Table 2. Regarding the more conventional data, it is noteworthy that the BNP values exceeded 190 pg/mL in all patients in the decompensated group (median =434.5 pg/mL). The decompensated PAH participants differed significantly from the compensated PAH group with regards to several measures of RA mechanics (RALS, S-D loop area indexed for either BSA or RA length), RV mechanics (RV to LV ratio in diastole), laboratory markers, and clinical function (percent-predicted 6MWT) as outlined in Table 2. Both the compensated and decompensated PAH populations exhibited elevated mPAPs. Despite a trend of observing higher mPAP and PVRi in the decompensated group, there was no statistical difference in mPAP or PVRi between the two subgroups.
We wish to point out that the area enclosed by the S-D loop was similar between the control and compensated groups, but significantly smaller in the decompensated PAH patients compared with both control and compensated PAH individuals (Fig. 3, Table 2).
Characteristic appearance of S-D loops in a normal subject, compensated PAH, and decompensated PAH. The area enclosed by the S-D loop is notably smaller in decompensated PAH, but not in compensated PAH. Here LD is indexed for body surface area. Gradually increasing negative values of LD indicate filling of RA.
Outcomes in PAH populations
ROC curve data for measures of right heart function and adverse outcomes.
AUC, area under ROC curve; Loop arealength, S-D loop area indexed for right atrial length; Loop areaBSA, S-D loop area indexed for body surface area; RALS, right atrial longitudinal strain; BNP, brain natriuretic peptide; CI, confidence interval.

ROC curves comparing (a) S-D loop area indexed for RA length vs outcome, (b) S-D loop area indexed for BSA vs outcome, (c) RALS vs outcome, (d) BNP vs outcome, (e) percent-predicted TAPSE vs outcome, and (f) TAPSE indexed for BSA vs outcome. Additional data for each ROC curve (cut-point, sensitivity, specificity, area under curve [AUC]) can be found in Table 3. In all ROC curves, cut-points are delineated by triangles.
Intra-observer and inter-observer variability
We found excellent reproducibility for RALS and LD measurements. ICC values for inter-observer variability were 0.88 for RALS and 0.94 for LD, whereas ICC values for intra-observer variability were 0.9 for RALS and 0.92 for LD, respectively.
Discussion
Although the RA is often regarded as the “forgotten chamber” of the heart, RA mechanics may serve as an important barometer of right heart dysfunction. Strain imaging represents an exciting new tool for evaluating RA mechanics. This single-center, pilot study represents one of the first efforts to use this tool in the evaluation of pediatric PAH. Additionally, in this study we have proposed a new concept, the S-D loop, that may provide useful discriminatory information regarding compensated and decompensated states of PAH in children.
Right atrial longitudinal strain
Measurement of RA strain may represent an exciting new tool in the assessment of RA physiology that was not available in the past. Previous research into pediatric RA strain have demonstrated abnormal indices in single ventricle patients 22 and in patients with repaired tetralogy of Fallot. 23 Adult studies have demonstrated a decreased RALS in participants with PAH compared with controls6,7 and this finding is reflected in our pediatric data. Of interest, RALS was also significantly different between compensated and decompensated subgroups of PAH.
Patients with chronic right heart pressure overload demonstrate diastolic dysfunction before systolic dysfunction. 24 The gradual increase in RA strain is reflected by the ascending limb and peak of the RA strain curve (Fig. 2). This phase occurs during atrial filling, referred to as reservoir function of the RA. Atrial filling, in turn, occurs during ventricular diastole. Therefore, the lower RALS in patients with PAH may reflect the abnormal myocardial relaxation and diastolic dysfunction of the RV.
Atrial strain-displacement loops
To gain insights into RA function and mechanics over the entire cardiac cycle, we propose a novel approach, the S-D loop. Plotting atrial longitudinal strain against LD produced a characteristic elliptical loop, with turning points corresponding to tricuspid valve closing and opening (Fig. 2). Strain represents local deformation of a small myocardial segment, whereas displacement evaluates the movement of the whole chamber. Moreover, displacement may serve as a surrogate for change in RA volume. The S-D loops may serve as a helpful index, comparing atrial segmental deformation against the overall movement of the atrial wall and may provide a useful “composite view” of atrial filling (reservoir phase) and emptying (conduit + contractile phases). The concept of simultaneous measurement of longitudinal strain and LD has been applied in adults with constrictive pericarditis, 25 RV pressure overload,9,10 LV pressure overload in aortic stenosis, 26 and in the evaluation of carotid artery stiffness. 27 Longitudinal strain and LD have been used to describe the elasticity of the dilated aorta. Elasticity of the aorta incorporates both the property of dilating and the property of recoiling to its initial shape. 28 Similarly, the atrial S-D relationship may provide indirect evidence of the elasticity of the RA tissue. We also speculate that during the conduit phase (when tricuspid valve remains open), the RA S-D loop may also reflect RV diastolic stiffness. The concept of the S-D loop is quite new. It provides simultaneous information regarding the regional shortening and global movement of the atrial wall during the entire cardiac cycle. The concept was derived from torsion-displacement loop described by us and others.28,29 To describe the area enclosed by torsion-volume loops, Notomi et al. coined the term “Filling Work.” They proposed that the area enclosed by the loops gives a measure of “work” being performed by the chamber and may provide a measure of the energy stored in the cardiac chamber. 30 Although a theoretical analogy is being drawn, we cannot extrapolate these theoretical concepts to our study, as ours is a pilot study that aimed to show the feasibility of constructing strain versus displacement loops in normal and diseased RA. Clearly, further work is needed to determine the physiologic significance of the area enclosed by the S-D loop. Based on other studies that have used strain and displacement simultaneously, we can merely speculate that the reduced area of the S-D loop may provide a visual display of the decreased elasticity of the atrial tissue. It is a first step in visually describing the RA mechanics during an entire cardiac cycle. This was nicely demonstrated in patients with decompensated PAH, where the S-D loops were notably “skinnier” (Fig. 3, lower panel).
Comparison with current means of assessing right heart mechanics in PAH
Non-invasive identification of impending RV decompensation is challenging in children with PAH, since the 6MWD and the degree of TR do not deteriorate with increase in disease severity at the same rate as in adult PAH. Moreover, in the absence of sufficient TR, current echocardiographic indices are poor predictors of disease severity. In this scenario, RA mechanics may provide a better insight into the state of RV decompensation. Both RALS and S-D loops may serve as upgrades to the clinician’s toolkit. It is important to recognize the best use of each tool available for the evaluation of PAH in children. It may be that a combination of tools will yield the most complete picture of right heart mechanics at a given point in time. This study and others in adults6,7 have demonstrated that the decrease in RALS with increasing disease severity makes RALS a useful measure of RV diastolic dysfunction. Since diastolic dysfunction of the RV precedes systolic dysfunction in chronic right heart pressure overload, 22 RALS may be a useful index to monitor serial diastolic changes in patients with PAH, similar to findings in serum BNP. In contrast, TAPSE and RV strain provide serial evaluation of RV during systole, while the RA S-D loop may provide a composite insight into both systolic and diastolic function of RV. Our study showed that the RA indices had greater sensitivity and specificity (and larger AUCs) compared to the more conventional indices like BNP and TAPSE. In our study, out of all the indices evaluated, the S-D loop indexed for RA length showed the best sensitivity (100%) and specificity (92.3%). This points to great potential value of this index but also warrants further validation. In particular, serial monitoring of these indices after initiating therapy for PAH would be of utmost interest to researchers and clinicians alike.
Clinical implications
Established clinical criteria used in the evaluation of adult PAH may be less useful in the pediatric population. For example, the 6MWT can be difficult to carry out according to established protocols in younger children and does not always reflect the degree of right heart dysfunction. Children often present with lower RA pressures, better compensation of the RV to the increased afterload, and higher mixed venous oxygen saturations.2,13 The study of RA mechanics may therefore be of value in this population. The S-D loop provides us with a composite view of the reservoir, conduit, and contractile phases of RA function. The mathematical calculations for measuring these new indices can be incorporated into calculation packages of commercial echocardiography machines. Automation of RALS and S-D loops into commercial ultrasound machines may encourage more use of these indices in busy clinical settings and not limit them to the research arena only.
Study limitations
Due to small sample size, further subgroup analyses were not feasible. Furthermore, statistical evaluations determining degree of association between the laboratory variables and adverse outcomes were not performed because of small sample size. We excluded individuals with ongoing lung parenchymal lesions, significant atrial level shunts, pulmonary vein stenoses, and LV dysfunction. In this fashion, we present a relatively homogeneous population which should have narrow limits for test results and clinical outcomes. Moreover, even in a tertiary referral center, idiopathic PAH remains an uncommon disease. As our primary aim in this study was to perform research on a new echocardiographic index, we pursued a statistically sufficient number of patients needed to provide meaningful results, at the same time limiting heterogeneity in the sample. This study is a single-center, pilot study that may pave the way for larger, multicenter studies. The primary focus of this study was on the RA physiology. Therefore, we did not focus significantly on the RV physiology.
Conclusions
RA mechanics may serve as novel metrics in the assessment of right heart dysfunction in children with PAH. Identification of impending RV decompensation is difficult in pediatrics, where sub-maximal exercise testing and the degree of TR are often better than expected for the degree of elevated mPAP. RA strain and S-D loops may serve as useful tools for identifying patients at risk for future adverse outcomes from progressive RV dysfunction.
Footnotes
Conflict of interest
The author(s) declare that there is no conflict of interest.
Funding
This research was funded by the Madeline Schwamm Memorial Fund for Pediatric Pulmonary Hypertension Research.
