Skip to main content

ORIGINAL RESEARCH article

Front. Cardiovasc. Med., 04 June 2021
Sec. Cardiovascular Imaging
This article is part of the Research Topic Highlights in Cardiovascular Imaging: 2021 View all 11 articles

Age- and Sex-Specific Changes in CMR Feature Tracking-Based Right Atrial and Ventricular Functional Parameters in Healthy Asians

\nShuang Leng&#x;Shuang Leng1Jiajun Guo&#x;Jiajun Guo2Ru-San Tan,Ru-San Tan1,3Ping Chai,Ping Chai4,5Lynette Teo,Lynette Teo5,6Marielle V. Fortier,,Marielle V. Fortier3,7,8Chao GongChao Gong2Xiaodan ZhaoXiaodan Zhao1Ching Ching Ong,Ching Ching Ong5,6John C. AllenJohn C. Allen3Wen Ruan,Wen Ruan1,3Angela S. Koh,Angela S. Koh1,3Teng Hong Tan,Teng Hong Tan3,7James W. Yip,James W. Yip4,5Ju Le Tan,Ju Le Tan1,3Yucheng Chen
Yucheng Chen2*Liang Zhong,
Liang Zhong1,3*
  • 1National Heart Centre Singapore, Singapore, Singapore
  • 2Cardiology Division, Department of Medicine, West China Hospital, Sichuan University, Chengdu, China
  • 3Duke-NUS Medical School, Singapore, Singapore
  • 4Department of Cardiology, National University Heart Centre, Singapore, Singapore
  • 5Yong Loo Lin School of Medicine, National University of Singapore, Singapore, Singapore
  • 6Department of Diagnostic Imaging, National University Hospital, Singapore, Singapore
  • 7KK Women's and Children's Hospital, Singapore, Singapore
  • 8Singapore Institute for Clinical Sciences, A*STAR, Singapore, Singapore

Cardiovascular magnetic resonance (CMR) is the reference standard for non-invasive assessment of right-sided heart function. Recent advances in CMR post-processing facilitate quantification of tricuspid annular (TA) dynamics and longitudinal strains of the right ventricle (RV) and right atrium (RA). We aimed to determine age- and sex-specific changes in CMR-derived TA dynamics, and RV and RA functional parameters in healthy Asian adults. We studied 360 healthy subjects aged 21–79 years, with 30 men and 30 women in each of the six age groups. Functional parameters of RV and RA were measured on standard four-chamber cine CMR using fast feature tracking: (1) TA peak velocities (systolic velocity S′, early diastolic velocity E′, late diastolic velocity A′) and TA plane systolic excursion (TAPSE); (2) RV global longitudinal strain (GLS) and strain rates; and (3) RA phasic longitudinal strains and strain rates. S′ and TAPSE exhibited negative correlations with age. RV GLS was significantly higher in females than in males but not associated with age in both sexes. Females had similar E′, lower A′, and higher E′/A′ ratios compared to males. Positive associations of E′ and E′/A′, and negative association of A′ with age were observed in both sexes. Females had higher RA reservoir and conduit strains compared to males. There were significantly negative and positive associations between RA conduit and booster strains, respectively, with age. Age- and sex-specific reference ranges were established, and associations revealed, for fast CMR feature tracking parameters of right heart function in a large normal Asian population.

Introduction

The right side of the heart has traditionally received less attention than the left, yet there is a growing body of evidence showing that right heart size and function are perhaps of equal importance in providing diagnostic and prognostic information in a wide range of cardiovascular diseases (1).

Cardiovascular magnetic resonance (CMR) imaging is the current gold standard for quantitation of right ventricular (RV) geometry and function (2). Guidelines recommend CMR for assessment of global RV functional parameters such as RV volumes, ejection fraction (EF), and cardiac output (3). Recent advances in CMR post-processing have facilitated reliable quantitation of tricuspid annular (TA) dynamics and longitudinal strains of the RV and right atrium (RA) using feature tracking (48). CMR-based TA systolic velocity (S′), early diastolic velocity (E′), late diastolic velocity (A′), and TA plane systolic excursion (TAPSE) have been reported to be significantly lower in patients diagnosed with heart failure (4), hypertrophic cardiomyopathy (4), repaired tetralogy of Fallot (rTOF) (4, 5), unrepaired atrial septal defect (5), and severe tricuspid insufficiency (6) than in normal controls. RV and RA longitudinal strains can be quantified either by feature tracking of the whole RV and RA endocardium, respectively (9, 10), or more rapidly by tracking discrete anatomical points on standard cine CMR (7, 8). A prior study demonstrated the prognostic utility of impaired fast CMR feature tracking-derived RA strain parameters for predicting clinical and hemodynamic deterioration in patients with pulmonary arterial hypertension (PAH) (7). The quantitative RV fast strain and strain rate parameters assessed from CMR identify abnormalities of RV function in rTOF and PAH and are predictive of exercise capacity, RV decompensation, and clinical risks in these patients (8).

The emerging importance of these parameters motivates efforts to define normal ranges and distributions in a healthy population. In addition, most CMR normal ranges for right heart size and volume have been established for Caucasians and may not be generalizable to other ethnicities (11). Accordingly, the aim of this study is to determine age- and sex-specific changes in CMR-derived RV and RA dimensions and functional (systolic and diastolic) parameters among healthy Asian subjects.

Materials and Methods

Study Population

In this multicenter study, 360 subjects (aged 21–79 years, 180 males and 180 females) without known cardiovascular disease were enrolled from three hospitals between 2014 and 2019. The subjects were recruited from the (1) Cardiac Aging study (12), (2) INITIATE study, and (3) Database of healthy controls in West China Hospital. All subjects (1) had no known cardiovascular disease and demonstrated no signs of cardiovascular disease in prior electrocardiographic or echocardiographic investigations, (2) displayed no uncontrolled cardiovascular risk factors at the time of enrollment (systolic blood pressure ≤140 mmHg and diastolic blood pressure ≤90 mmHg without anti-hypertensive treatment, total cholesterol <6.2 mmol/L, fasting glucose <7 mmol/L, body mass index <30 kg/m2, and current non-smoker status), and (3) had no significant kidney or lung disease. Other exclusion criteria were cerebrovascular disease or nervous system disease, cancer, autoimmune diseases, recent systemic infection (within a month), recent surgery or severe trauma (within a month), any recent medications, and a history of implantation of a pacemaker or other metals that are contraindicated for CMR. The institutional review board of each hospital approved the study protocol. Informed consent was obtained from each participant.

CMR Acquisition

CMR acquisitions were performed using a 3T magnetic resonance scanner (Ingenia, Philips Healthcare, The Netherlands) at National Heart Centre Singapore, a 1.5T MAGNETOM Aera magnetic resonance scanner (Siemens Healthineers, Erlangen, Germany) at National University Hospital Singapore, and a 3T MAGNETOM Tim Trio magnetic resonance scanner (Siemens Healthineers, Erlangen, Germany) at West China Hospital (Sichuan China). End-expiratory breath-hold-balanced steady-state free precession cine images were acquired in standard short- and long-axis views. Typical parameters for the Philips scanner were as follows: repetition time (TR)/echo time (TE), 3/1 ms; matrix, 240 × 240; flip angle, 45°; field of view, 300 × 300 mm; pixel spacing, 1.25 × 1.25 mm; slice thickness, 8 mm; number of frames, 30/40 per cardiac cycle. Parameters for the 1.5T Siemens scanner were: TR/TE, 33/1 ms; matrix, 192 × 180; flip angle, 58°; field of view, 320 × 300 mm; pixel spacing, 1.67 × 1.67 mm; slice thickness, 6 mm; number of frames, 30 per cardiac cycle. Parameters for the 3T Siemens scanner were: TR/TE, 3.4/1.3 ms; matrix, 192 × 162; flip angle, 50°; field of view, 320 × 270 mm; pixel spacing, 1.67 × 1.67 mm; slice thickness, 8 mm; and number of frames, 25 per cardiac cycle.

Echocardiography

A sub-study was conducted to validate CMR-derived TA dynamics against those measured using echocardiography. Echocardiography was performed on the same day as CMR using a commercial ultrasound system (Aloka α10, Japan). TA velocities and displacement were measured using tissue Doppler imaging (TDI) and M-mode, respectively, in the apical four-chamber view.

CMR Data Analysis

The parameters measured were (1) RV volumes and RVEF (derived from standard CMR volumetric analysis), (2) Two-dimensional chamber dimensions (RV chamber diameters and RA area and diameters), (3) TA velocities (systolic velocity S′, early diastolic velocity E′, and late diastolic velocity A′) and displacement (TAPSE), (4) RV longitudinal strain and strain rates, and (5) RA phasic (reservoir, conduit, and booster) longitudinal strains and strain rates.

Volumetric Analysis and Chamber Dimensions

Endocardial contours were manually traced from the stack of short-axis cine images to obtain RV end-diastolic volume (EDV) and end-systolic volume (ESV), from which stroke volume (SV) and EF were derived (Figures 1A,B). Papillary muscles and trabeculae were included in the blood volume (3). RV size was measured from the four-chamber view at end-diastole. RV basal diameter was measured at the level of the tricuspid valve, and mid-cavity diameter was measured in the middle third of the RV at the level of the left ventricular (LV) papillary muscle (13) (Figure 1C).

FIGURE 1
www.frontiersin.org

Figure 1. Derivation of right ventricular (RV) volumes by contouring of RV at (A) end-diastole and (B) end-systole; (C) RV basal and mid-cavity diameters at end-diastole; (D) right atrial (RA) area, longitudinal and transverse diameters at RV end-systole; (E) tricuspid annular dynamics at RV medial (septum) and lateral (free wall); (F) longitudinal strains of RV and RA; (G) Example of RV strain and strain rate curves; (H) Example of RA strain and strain rate curves.

RA area and diameters were measured from the four-chamber view at RV end-systole. RA longitudinal diameter was measured as the distance between the midpoint of the line joining the medial and lateral (free wall) insertions of the tricuspid valve and the posterior wall (roof) of the RA. RA transverse diameter was measured as the orthogonal line bisecting the longitudinal diameter (14) (Figure 1D).

All chamber dimension measurements were reported as both absolute values and indexed values normalized to body surface area (BSA).

Tricuspid Annular Velocities and Displacement

Custom software developed in the MATLAB environment (MathWorks Inc., MA, USA) and validated in prior studies (4, 7, 8, 1519) was used to perform the semi-automatic tracking of medial and lateral tricuspid valve insertions in the four-chamber view (Figure 1E). A video of the TA tracking and resultant TA velocity and displacement curves are shown in Supplementary Materials. Peak TA velocities were read off the velocity curve during systole (S′), early diastole (E′), late diastole (A′); TAPSE—TA displacement at end-systole—was read off the displacement curve (4, 7). We reported mean values for S′, E′, A′, and TAPSE averaged from medial and lateral TA measurements.

Longitudinal Strains of RV and RA

Using the same custom software, the RV epicardial apex and the RA roof point were tracked in the four-chamber view over the cardiac cycle (Figure 1F). We followed the convention of determining RV apex location by its proximity to the LV apex. The RV apex is close to but separated from the LV apex by the interventricular septal wall, which can be seen on the four-chamber view (13). The RA roof was localized to the intersection of the RA posterior wall and the RA longitudinal diameter.

The distance (L) from either the medial or the lateral tricuspid valve insertion to the RV epicardial apex on the four-chamber view was measured at any time point (t). RV strain value at time t with respect to RV end-diastole (time 0) was calculated using the formula for Lagrangian strain: (L(t)−L(0)) × 100/L(0) (8). RV global longitudinal strain (GLS)—strain value at end-systole—was read off the strain curve, and peak global longitudinal strain rates at RV systole (GLSRs), early diastole (GLSRe), and late diastole (GLSRa) were calculated as first time-derivatives of the strain curve at the corresponding cardiac phases (8) (Figure 1G).

Similarly, RA longitudinal strain (ε) was derived from the time-varying distances between either the medial or the lateral tricuspid valve insertion and the RA roof point (7). RA reservoir strain (εs), conduit strain (εe) and booster strain (εa) were read off the generated strain curve at RV end-systole, diastasis, and pre-RA systole, respectively, and the corresponding peak strain rates (SRs, SRe, and SRa) calculated as first time-derivatives of the strain curve at the corresponding cardiac phases (7) (Figure 1H).

A video showing the derivation of RV and RA strain is given in Supplementary Materials. Mean RV and RA strain and strain rate parameters averaged from medial and lateral TA measurements were used for all analyses. RV strain and strain rate parameters were presented as absolute values.

Conventional Feature Tracking Strain

A sub-study was conducted to validate the fast feature tracking-derived RA and RV strain measurements against those derived from conventional CMR RA and RV endocardial feature tracking using dedicated QStrain software (Version 2.0, Medis BV, Leiden, The Netherlands).

Statistical Analysis

The distribution normality of continuous variables was assessed using the Shapiro-Wilk test. Data were expressed as mean ± standard deviation (SD), and reported across six pre-specified age groups stratified by sex. Age-specific reference limits were defined as mean ± 1.96·SD. Values obtained in females and males were compared using Student's t-test. For either sex, linear regression was used to characterize the relationship of right heart measurements with age. Intra- and inter-observer variability of study parameters was studied in a selected subgroup of 20 cases using Bland-Altman analysis and coefficient of variation. Statistical analyses were performed using SPSS software (Version 17.0, IBM, Chicago, IL, USA). P ≤ 0.05 indicated statistical significance.

Results

Baseline Characteristics and Summary Results

We recruited 360 healthy individuals (male:female 180:180, age range 21–79 years). Baseline demographics and right heart CMR measurements stratified by sex are given in Table 1. Age group-specific CMR measurements stratified by sex are presented in Tables 2, 3. Figures 2, 3 show RV volume and systolic and diastolic RV function parameters plotted against age with 5th, 50th, and 95th percentile values in males and females, respectively. Figures 4, 5 show RA dimensions and phasic function parameters vs. age with 5th, 50th, and 95th percentile values in males and females, respectively. Measurements of right heart dimensions and volumes stratified by sex and age group are presented in Supplementary Tables 1–3.

TABLE 1
www.frontiersin.org

Table 1. Demographics and right heart function parameters of the control population.

TABLE 2
www.frontiersin.org

Table 2. Males: right heart function parameters by age group [mean ± SD (reference range, lower/upper limits calculated as mean ± 1.96·SD)].

TABLE 3
www.frontiersin.org

Table 3. Females: right heart function parameters by age group [mean ± SD (reference range, lower/upper limits calculated as mean ± 1.96·SD)].

FIGURE 2
www.frontiersin.org

Figure 2. Right ventricular parameters plotted against age in male subjects. (Central line) The 50th percentile; (Top line) The 95th percentile; (Bottom line) The 5th percentile.

FIGURE 3
www.frontiersin.org

Figure 3. Right ventricular parameters plotted against age in female subjects. (Central line) The 50th percentile; (Top line) The 95th percentile; (Bottom line) The 5th percentile.

FIGURE 4
www.frontiersin.org

Figure 4. Right atrial parameters plotted against age in male subjects. (Central line) The 50th percentile; (Top line) The 95th percentile; (Bottom line) The 5th percentile.

FIGURE 5
www.frontiersin.org

Figure 5. Right atrial parameters plotted against age in female subjects. (Central line) The 50th percentile; (Top line) The 95th percentile; (Bottom line) The 5th percentile.

Validation of Fast CMR Feature Tracking-Derived Measurements

Results from a sub-study of 60 subjects comprising 10 randomly selected from each of the six age groups showed that fast CMR feature tracking-derived TA dynamic measurements, RA and RV strain measurements were in good agreement with echocardiographic (r = 0.76–0.83, P < 0.0001, Supplementary Figure 1) and conventional CMR endocardial tracking-derived results (r = 0.75–0.94, P < 0.0001, Supplementary Figures 2–4).

Influence of Sex and Age on Right Ventricular Function

Among RV systolic functional parameters, RVEF and RV GLS were significantly higher in females than in males but there were no significant sex differences in S′, TAPSE, and RV GLSRs (Table 1). RVEF correlated positively with age in both males (r = 0.31, P < 0.0001) and females (r = 0.30, P < 0.0001). S′ and TAPSE exhibited negative correlations with age in both males (S′: r = −0.36, P < 0.0001; TAPSE: r = −0.28, P < 0.0001) and females (S′: r = −0.30, P < 0.0001; TAPSE: r = −0.31, P < 0.0001). RV GLS and RV GLSRs were not associated with age in both sexes (Figures 2, 3).

With respect to RV diastolic function, females had significantly lower A′ and higher E′/A′ ratios and RV GLSRe compared to males, but no significant sex differences in E′ and RV GLSRa were observed (Table 1). Among male subjects, E′ (r = −0.75, P < 0.0001), E′/A′ ratio (r = −0.80, P < 0.0001) and RV GLSRe (r = −0.58, P < 0.0001) decreased significantly with age, whereas A′ (r = 0.49, P < 0.0001) and RV GLSRa (r = 0.57, P < 0.0001) increased with age. Similar positive associations of E′, E′/A′, and RV GLSRe, and negative associations of A′ and RV GLSRa with age were observed among females (all P < 0.0001) (Figures 2, 3).

Influence of Sex and Age on Right Atrial Phasic Function

Females had higher RA reservoir strain, conduit strain, reservoir strain rate, and conduit strain rate compared to males, while RA booster strain and strain rate were similar (Table 1). RA reservoir strain, conduit strain, and strain rate correlated negatively with age among both males (reservoir strain: r = −0.26, P < 0.0001; conduit strain: r = −0.61, P < 0.0001; conduit strain rate: r = −0.70, P < 0.0001) and females (reservoir strain: r = −0.38, P < 0.0001; conduit strain: r = −0.70, P < 0.0001; conduit strain rate: r = −0.75, P < 0.0001). Booster strain and strain rate exhibited positive correlations with age among both males (booster strain: r = 0.46, P < 0.0001; booster strain rate: r = 0.42, P < 0.0001) and females (booster strain: r = 0.51, P < 0.0001; booster strain rate: r = 0.34, P < 0.0001) (Figures 4, 5).

Influence of Sex and Age on Right Ventricular Dimensions

Absolute and indexed RV EDV, ESV, and SV were higher in males compared to females (P < 0.0001 for all). Males had larger absolute RV basal and mid-cavity diameters than females (both P < 0.0001) but not when normalized to BSA (Supplementary Table 1). RV volumes correlated negatively with age among both males (EDV: r = −0.50, P < 0.0001; ESV: r = −0.49, P < 0.0001) and females (EDV: r = −0.39, P < 0.0001; ESV: r = −0.41, P < 0.0001), with similar correlations after normalization to BSA (males, EDV/BSA: r = −0.43, P < 0.0001; ESV/BSA: r = −0.44, P < 0.0001; females, EDV/BSA: r = −0.45, P < 0.0001; ESV/BSA: r = −0.44, P < 0.0001). Absolute and indexed RV basal and mid-cavity diameters were not significantly associated with age in both sexes.

Influence of Sex and Age on Right Atrial Dimensions

Males had larger absolute RA diameters and areas compared to females. After normalization to BSA, indexed RA longitudinal and transverse diameters were significantly smaller in males compared to females, while the indexed RA area was not significantly different (Supplementary Table 1). Among male subjects, age was positively correlated with indexed RA longitudinal diameter (r = 0.32, P < 0.0001). Among female subjects, age was positively correlated with indexed RA longitudinal diameter (r = 0.37, P < 0.0001) and indexed RA area (r = 0.30, P < 0.0001).

Reproducibility

Good to excellent intra- and inter-observer reproducibility was observed for RV and RA dimensions and volumes, TA velocities and displacement, RV longitudinal strain and strain rates, and RA phasic strains and strain rates (Table 4).

TABLE 4
www.frontiersin.org

Table 4. Intra- and inter-observer reproducibility for right heart dimension and function parameters.

Discussion

In this study, we investigated age- and sex-specific changes in CMR-based right heart dimensions and functional measurements in a large, healthy Asian population. Increased age was associated with impairment of S′ and E′ velocities, TAPSE, and RV GLSRe and increases in A′ velocity and RV GLSRa, in both sexes. RV GLS and RV GLSRs were not affected by age in both sexes. As age increased, RA reservoir and conduit strains and strain rates decreased, while RA booster strain and strain rate increased.

CMR is recommended for many patients at the time of the transition from pediatric to adult congenital heart disease programmes, and this is the gold standard for RV volume, EF, flow quantification, and assessment of extracardiac anatomy. Additionally, CMR is recommended in the presence of clinical deterioration, non-diagnostic echo findings, and prior to surgical or transcatheter intervention (20). Data on age and sex-specific CMR reference ranges for right heart function in an Asian population are lacking at present. This information is crucial to CMR clinical practice in Asia. Moreover, our CMR-based method is simple, reproducible, and easily implemented for efficient RA and RV function assessment.

This work is novel for a number of reasons. First, it comprises the largest Asian CMR cohort reported to date, enrolling participants from multiple centers to provide broad generalizability of results and high precision in sample estimates of mean and SD. Second, CMR parameters were stratified not only by sex but also by age categories across a broad age range of 21–79 years. The increased granularity of results is key to clinical implementation and interpretation in light of the sex- and age-related differences observed in many of the studied parameters. Third, this study is the first to report age- and sex-specific reference ranges for various recently published fast CMR feature tracking parameters used to assess RV systolic and diastolic function as well as phasic RA function.

RV Functional Parameters

Asian subjects in our study have numerically smaller right heart size compared to Caucasians in another study (21) (RVEDV: 121 ± 34 vs. 154 ± 40 ml; indexed RVEDV: 72 ± 15 vs. 85 ± 17 ml/m2), which underscores the need for ethnic-specific reference values. Our study showed a negative correlation between RV volume and age and a positive correlation between RVEF and age, corroborating similar observations in Caucasian subjects (21). Beyond conventional RV volumes and EF, we quantitated TA dynamics and RV longitudinal strain and strain rates—measurements for assessing RV function in both systole and diastole—and reported their age- and sex-specific reference ranges. Increasing age was associated with reductions in TAPSE, S′, E′, E′/A′, and RV GLSRe and increases in A′ and RV GLSRa in both sexes. These results agree with prior studies involving echocardiographic TDI (22) and conventional feature tracking CMR with endocardial contour tracing (23). Of note, we found that RV GLS and RV GLSRs were not associated with age, which suggests that they are suited to be applied clinically as markers of disease progression unaffected by chronological aging. Our study replicated prior findings of higher systolic functional performance in females than in males as determined by RVEF (21) and RV GLS (24).

RA Functional Parameters

We previously demonstrated that RA function as assessed using our fast RA feature tracking method had important diagnostic and prognostic implications in patients with PAH (7). RA volume and function (i.e., emptying fraction) assessment using volumetric analysis requires additional cine CMR acquisitions of atrial short-axis images and post-acquisitional analysis, which are not routinely performed. The area-length method, through geometric assumptions, shortens analysis time slightly (calculation still has to be repeated at all time phases) but is dogged by accuracy and reproducibility issues (25). In our study, we characterized the RA reservoir, conduit, and booster pump functions using strain and strain rate parameters that can be reproducibly and expeditiously measured. These parameters have been shown to be useful for detecting RA dysfunction, RV decompensation, and monitoring disease progression in patients with PAH (7). Among the healthy subjects in our study, females had higher reservoir and conduit strains but similar booster strain compared to males. This is concordant with a previous study of Caucasian subjects (25). In addition, we found age-related changes in RA function manifested as a significant decrease in conduit strain and an increase in booster strain. We surmise the latter is necessary in order to maintain reservoir strain, which increased slightly with age. These results are in agreement with previous reports using speckle tracking echocardiography (26) and feature tracking CMR with endocardial contour tracing (27).

Advantages of Fast CMR Feature Tracking Parameters

This is the first study to establish age- and sex-specific reference ranges for CMR-derived TA dynamic parameters, which are a prerequisite for clinical applicability and adoption. The CMR-derived TA velocities correlated significantly with TDI-derived measurements, with no angle dependence and excellent reproducibility. This corroborates our previous work which demonstrated good correlation and agreement between CMR-derived mitral annular velocities and TDI measurements (18).

In addition to TA dynamics, we presented simplified strain indices relating time-varying distances between TA medial and lateral points and a fixed anatomical point at the RV apex or RA roof. Good correlation between the simplified longitudinal strains and those obtained from conventional endocardial contour-based feature tracking CMR was demonstrated, which is consistent with the results of our prior studies involving patients with rTOF, PAH, and age- and sex-matched controls (7, 8). In fact, the accuracy of conventional endocardial feature tracking-derived longitudinal strain is degraded in subjects with vigorous TA motion, as contour tracking of the RV free wall segment adjacent to the tricuspid valve in the long-axis view becomes adversely affected (28). Compared to conventional contour-based strain measurements, the fast feature tracking strain measurements presented in this study are less dependent on RV and RA geometry, and more closely approximate the motion and reflect the function of longitudinal fibers in the RV, which are the greatest contributors to RV contraction (29). Our prior studies have shown that the simplified approach is not only faster but also more reproducible compared to the conventional CMR feature tracking approach (7, 8, 15, 16).

Limitations

There were some limitations in the present study. First, as CMR examinations were not performed repeatedly on the same subjects over time, the associations between age and CMR parameters are cross-sectional, not longitudinal. Nonetheless, the cross-sectional study design is commonly used in studies to establish reference ranges and make inferences regarding relationships in support of further research and clinical studies. We need a longitudinal study with repeated CMR scans on the same individual, which can help us investigate how age affects the heart structural and functional changes. Second, RV mass was not assessed in this study as RV mass is usually not quantified in a routine assessment because of the thin RV myocardium in healthy subjects.

Conclusions

We investigated age- and sex-related CMR measurements for right heart dimensions and function in a large Asian cohort that are of significant clinical and research utility in Asia. Using CMR, knowledge of age-, sex-, and ethnicity-specific distributions of right heart measurements should aid in the correct interpretation of disease states.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding authors.

Ethics Statement

The studies involving human participants were reviewed and approved by SingHealth Centralized Institutional Review Board. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

R-ST, PC, LT, MF, CO, WR, AK, TT, JY, JT, YC, and LZ conceived the study design. SL, JG, CG, and XZ analyzed data. SL, R-ST, PC, LT, MF, CO, WR, AK, TT, JY, JT, YC, and LZ interpreted results. JA performed statistical analysis. SL drafted manuscript. JG, R-ST, PC, LT, MF, CG, XZ, CO, JA, WR, AK, TT, JY, JT, YC, and LZ edited and revised manuscript. All authors read and approved the final manuscript.

Funding

This study was partially supported by the National Medical Research Council of Singapore (NMRC/OFIRG/0018/2016, MOH-000358, MOH-000351, NMRC/TA/0031/2015, and MOH-000153), SingHealth Duke-NUS Academic Medicine Research Grant [AM/TP032/2020 (SRDUKAMR2032)]. The funder had no role in the design and conduct of the study, collection, management, analysis, interpretation of the data, preparation, and review or approval of the manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2021.664431/full#supplementary-material

Abbreviations

A′, peak tricuspid annular late diastolic velocity; BSA, body surface area; CMR, cardiovascular magnetic resonance; εa, right atrial booster strain; εe, right atrial conduit strain; εs, right atrial reservoir strain; E′, peak tricuspid annular early diastolic velocity; EDV, end-diastolic volume; EF, ejection fraction; ESV, end-systolic volume; GLS, global longitudinal strain; GLSRa, global longitudinal strain rate during late diastole, GLSRe, global longitudinal strain rate during early diastole; GLSRs, global longitudinal strain rate during systole; LV, left ventricle; PAH, pulmonary arterial hypertension; RA, right atrium; rTOF, repaired tetralogy of Fallot; RV, right ventricle; S′, peak tricuspid annular systolic velocity; SD, standard deviation; SRa, right atrial booster strain rate; SRe, right atrial conduit strain rate; SRs, right atrial reservoir strain rate; SV, stroke volume; TA, tricuspid annular; TAPSE, tricuspid annular plane systolic excursion; TDI, tissue Doppler imaging.

References

1. Sanz J, Sánchez-Quintana D, Bossone E, Bogaard HJ, Naeije R. Anatomy, function, and dysfunction of the right ventricle: JACC state-of-the-art review. J Am Coll Cardiol. (2019) 73:1463–82. doi: 10.1016/j.jacc.2018.12.076

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Geva T. Is MRI the preferred method for evaluating right ventricular size and function in patients with congenital heart disease? MRI is the preferred method for evaluating right ventricular size and function in patients with congenital heart disease. Circ Cardiovasc Imaging. (2014) 7:190–7. doi: 10.1161/CIRCIMAGING.113.000553

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Schulz-Menger J, Bluemke DA, Bremerich J, Flamm SD, Fogel MA, Friedrich MG, et al. Standardized image interpretation and post processing in cardiovascular magnetic resonance: Society for Cardiovascular Magnetic Resonance (SCMR) Board of Trustees Task Force on Standardized Post Processing. J Cardiovasc Magn Reson. (2013) 15:35. doi: 10.1186/1532-429X-15-35

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Leng S, Jiang M, Zhao XD, Allen JC, Kassab GS, Ouyang RZ, et al. Three-dimensional tricuspid annular motion analysis from cardiac magnetic resonance feature-tracking. Ann Biomed Eng. (2016) 44:3522–38. doi: 10.1007/s10439-016-1695-2

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Ito S, McElhinney DB, Adams R, Bhatla P, Chung S, Axel L. Preliminary assessment of tricuspid valve annular velocity parameters by cardiac magnetic resonance imaging in adults with a volume-overloaded right ventricle: comparison of unrepaired atrial septal defect and repaired tetralogy of Fallot. Pediatr Cardiol. (2015) 36:1294–300. doi: 10.1007/s00246-015-1160-2

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Maffessanti F, Gripari P, Pontone G, Andreini D, Bertella E, Mushtaq S, et al. Three-dimensional dynamic assessment of tricuspid and mitral annuli using cardiovascular magnetic resonance. Eur Heart J Cardiovasc Imaging. (2013) 14:986–95. doi: 10.1093/ehjci/jet004

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Leng S, Dong Y, Wu Y, Zhao XD, Ruan W, Zhang GC, et al. Impaired CMR-derived rapid semi-automated right atrial longitudinal strain is associated with decompensated hemodynamics in pulmonary arterial hypertension. Circ Cardiovasc Imaging. (2019) 12:e008582. doi: 10.1161/CIRCIMAGING.118.008582

CrossRef Full Text | Google Scholar

8. Leng S, Tan RS, Guo JJ, Chai P, Zhang GC, 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:38. doi: 10.1186/s12968-021-00724-5

PubMed Abstract | CrossRef Full Text | Google Scholar

9. De Siqueira MEM, Pozo E, Fernandes VR, Sengupta PP, Modesto K, Gupta SS, et al. Characterization and clinical significance of right ventricular mechanics in pulmonary hypertension evaluated with cardiovascular magnetic resonance feature tracking. J Cardiovasc Magn Reson. (2016) 18:39. doi: 10.1186/s12968-016-0258-x

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Xie E, Yu R, Ambale-Venkatesh B, Bakhshi H, Heckbert SR, Soliman EZ, et al. Association of right atrial structure with incident atrial fibrillation: a longitudinal cohort cardiovascular magnetic resonance study from the Multi-Ethnic Study of Atherosclerosis (MESA). J Cardiovasc Magn Reson. (2020) 22:36. doi: 10.1186/s12968-020-00631-1

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Kawut SM, Lima JA, Barr RG, Chahal H, Jain A, Tandri H, et al. Sex and race differences in right ventricular structure and function: the multi-ethnic study of atherosclerosis-right ventricle study. Circulation. (2011) 123:2542–51. doi: 10.1161/CIRCULATIONAHA.110.985515

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Koh AS, Gao F, Leng S, Kovalik JP, Zhao XD, Tan RS, et al. Dissecting clinical and metabolomics associations of left atrial phasic function by cardiac magnetic resonance feature tracking. Sci Rep. (2018) 8:8138. doi: 10.1038/s41598-018-26456-8

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Valsangiacomo Buechel ER, Mertens LL. Imaging the right heart: the use of integrated multimodality imaging. Eur Heart J. (2012) 33:949–60. doi: 10.1093/eurheartj/ehr490

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Maceira AM, Cosín-Sales J, Roughton M, Prasad SK, Pennell DJ. Reference right atrial dimensions and volume estimation by steady state free precession cardiovascular magnetic resonance. J Cardiovasc Magn Reson. (2013) 15:29. doi: 10.1186/1532-429X-15-29

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Leng S, Tan RS, Zhao XD, Allen JC, Koh AS, Zhong L. Fast long-axis strain: a simple, automatic approach for assessing left ventricular longitudinal function with cine cardiovascular magnetic resonance. Eur Radiol. (2020) 30:3672–83. doi: 10.1007/s00330-020-06744-6

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Leng S, Tan RS, Zhao XD, Allen JC, Koh AS, Zhong L. Validation of a rapid semi-automated method to assess left atrial longitudinal phasic strains on cine cardiovascular magnetic resonance imaging. J Cardiovasc Magn Reson. (2018) 20:71. doi: 10.1186/s12968-018-0496-1

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Leng S, Ge H, He J, Kong LC, Yang YN, Yan FH, et al. Long-term prognostic value of cardiac MRI left atrial strain in ST-segment elevation myocardial infarction. Radiology. (2020) 296:299–309. doi: 10.1148/radiol.2020200176

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Leng S, Zhao XD, Huang FQ, Wong JI, Su BY, Allen JC, et al. Automated quantitative assessment of cardiovascular magnetic resonance-derived atrioventricular junction velocities. Am J Physiol Heart Circ Physiol. (2015) 309:H1923–35. doi: 10.1152/ajpheart.00284.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Ouyang RZ, Leng S, Sun AM, Wang Q, Hu LW, Zhao XD, et al. 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. (2021). doi: 10.1007/s00330-020-07643-6

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Di Salvo G, Miller O, Babu Narayan S, Li W, Budts W, Valsangiacomo Buechel ER, et al. Imaging the adult with congenital heart disease: a multimodality imaging approach-position paper from the EACVI. Eur Heart J Cardiovasc Imaging. (2018) 19:1077–98. doi: 10.1093/ehjci/jey102

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Petersen SE, Aung N, Sanghvi MM, Zemrak F, Fung K, Paiva JM, et al. Reference ranges for cardiac structure and function using cardiovascular magnetic resonance (CMR) in Caucasians from the UK Biobank population cohort. J Cardiovasc Magn Reson. (2017) 19:18. doi: 10.1186/s12968-017-0327-9

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Dalen H, Thorstensen A, Vatten LJ, Aase SA, Stoylen A. Reference values and distribution of conventional echocardiographic Doppler measures and longitudinal tissue Doppler velocities in a population free from cardiovascular disease. Circ Cardiovasc Imaging. (2010) 3:614–22. doi: 10.1161/CIRCIMAGING.109.926022

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Peng JP, Zhao XD, Zhao L, Fan ZM, Wang Z, Chen H, et al. Normal values of myocardial deformation assessed by cardiovascular magnetic resonance feature tracking in a healthy Chinese population: a multicenter study. Front Physiol. (2018) 9:1181. doi: 10.3389/fphys.2018.01181

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Qu YY, Li H, Rottbauer W, Ma GS, Buckert D, Rasche V. Right ventricular free wall longitudinal strain and strain rate quantification with cardiovascular magnetic resonance based tissue tracking. Int J Cardiovasc Imaging. (2020) 36:1985–96. doi: 10.1007/s10554-020-01895-5

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Maceira AM, Cosin-Sales J, Prasad SK, Pennell DJ. Characterization of left and right atrial function in healthy volunteers by cardiovascular magnetic resonance. J Cardiovasc Magn Reson. (2016) 18:64. doi: 10.1186/s12968-016-0284-8

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Peluso D, Badano LP, Muraru D, Dal Bianco L, Cucchini U, Kocabay G, et al. Right atrial size and function assessed with three-dimensional and speckle-tracking echocardiography in 200 healthy volunteers. Eur Heart J Cardiovasc Imaging. (2013) 14:1106–14. doi: 10.1093/ehjci/jet024

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Truong VT, Palmer C, Young M, Wolking S, Ngo TNM, Sheets B, et al. Right atrial deformation using cardiovascular magnetic resonance myocardial feature tracking compared with two-dimensional speckle tracking echocardiography in healthy volunteers. Sci Rep. (2020) 10:5237. doi: 10.1038/s41598-020-62105-9

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Bhave NM, Visovatti SH, Kulick B, Kolias TJ, McLaughlin VV. Right atrial strain is predictive of clinical outcomes and invasive hemodynamic data in group 1 pulmonary arterial hypertension. Int J Cardiovasc Imaging. (2017) 33:847–55. doi: 10.1007/s10554-017-1081-7

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Brown SB, Raina A, Katz D, Szerlip M, Wiegers SE, Forfia PR. Longitudinal shortening accounts for the majority of right ventricular contraction and improves after pulmonary vasodilator therapy in normal subjects and patients with pulmonary arterial hypertension. Chest. (2011) 140:27–33. doi: 10.1378/chest.10-1136

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: cardiovascular magnetic resonance, feature tracking, right ventricular function, right atrial function, age, sex

Citation: Leng S, Guo J, Tan R-S, Chai P, Teo L, Fortier MV, Gong C, Zhao X, Ong CC, Allen JC, Ruan W, Koh AS, Tan TH, Yip JW, Tan JL, Chen Y and Zhong L (2021) Age- and Sex-Specific Changes in CMR Feature Tracking-Based Right Atrial and Ventricular Functional Parameters in Healthy Asians. Front. Cardiovasc. Med. 8:664431. doi: 10.3389/fcvm.2021.664431

Received: 09 February 2021; Accepted: 26 April 2021;
Published: 04 June 2021.

Edited by:

Antonios Karanasos, Hippokration General Hospital, Greece

Reviewed by:

Alexios Antonopoulos, National and Kapodistrian University of Athens, Greece
Emmanuel Androulakis, Royal Brompton & Harefield NHS Foundation Trust, United Kingdom

Copyright © 2021 Leng, Guo, Tan, Chai, Teo, Fortier, Gong, Zhao, Ong, Allen, Ruan, Koh, Tan, Yip, Tan, Chen and Zhong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yucheng Chen, chenyucheng2003@126.com; Liang Zhong, zhong.liang@nhcs.com.sg

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.