Abstract
Right ventricular function has prognostic significance in patients with pulmonary hypertension. We evaluated whether cardiac magnetic resonance-derived strain and strain rate parameters could reliably reflect right ventricular systolic and diastolic function in precapillary pulmonary hypertension. End-systolic elastance and the time constant of right ventricular relaxation tau, both derived from invasive high-fidelity micromanometer catheter measurements, were used as gold standards for assessing systolic and diastolic right ventricular function, respectively. Nineteen consecutive precapillary pulmonary hypertension patients underwent cardiac magnetic resonance and right heart catheterization prospectively. Cardiac magnetic resonance data were compared with those of 19 control subjects. In pulmonary hypertension patients, associations between strain- and strain rate-related parameters and invasive hemodynamic parameters were evaluated. Longitudinal peak systolic strain, strain rate, and early diastolic strain rate were lower in PAH patients than in controls; peak atrial-diastolic strain rate was higher in pulmonary hypertension patients. Similarly, circumferential peak systolic strain rate was lower and peak atrial-diastolic strain rate was higher in pulmonary hypertension. In pulmonary hypertension, no correlations existed between cardiac magnetic resonance-derived and hemodynamically derived measures of systolic right ventricular function. Regarding diastolic parameters, tau was significantly correlated with peak longitudinal atrial-diastolic strain rate (r = −0.61), deceleration time (r = 0.75), longitudinal systolic to diastolic time ratio (r = 0.59), early diastolic strain rate (r = −0.5), circumferential peak atrial-diastolic strain rate (r = −0.52), and deceleration time (r = 0.62). Strain analysis of the right ventricular diastolic phase is a reliable non-invasive method for detecting right ventricular diastolic dysfunction in PAH.
Keywords
Introduction
In pulmonary arterial hypertension (PAH), a chronic elevation of pulmonary arterial pressure causes progressive right ventricular (RV) dysfunction leading to significant morbidity and premature death.1,2 Accurate evaluation of RV morphology and function is critical in the management of PAH considering that RV failure is the leading cause of death in this syndrome. While there is no consensus on the best method for assessing RV function in patients with PAH, recent studies have demonstrated the prognostic value of some non-invasive measures, including echocardiography-derived tricuspid annular plane systolic excursion 3 and cardiac magnetic resonance (CMR)-based assessment of RV ejection fraction. 4
RV strain and strain rate (SR), which can be obtained from CMR and calculated from feature tracking (or multimodality tissue tracking), have emerged as powerful tools to accurately quantify myocardial mechanics in various regional segments and directions. 5 Furthermore, RV strain analysis has been shown to be a useful tool for assessing disease severity and prognosis in patients with PAH.6–9
Invasive indices based on micromanometer catheter measurements remain the gold standard references in the evaluation of global RV systolic 10 and diastolic function. 11 However, these indices have only been used in a limited number of studies for RV assessment in PAH12,13 because of their invasiveness and need for special setups and expertise.
In this context, we aimed to ascertain whether CMR-based strain and SR imaging could provide accurate assessment of RV systolic and diastolic function in PAH patients. Accordingly, we first measured CMR-derived strain and SR parameters of PAH patients and compared them with those of age- and sex-matched subjects, using data from the multiethnic study of atherosclerosis subjects (MESA). 14 We then aimed to directly compare the strain and SR parameters from CMR (performed within 5 h of catheterization) with invasive indices obtained by high-fidelity micromanometer catheter measurements in PAH patients.
Methods
We enrolled patients with known or suspected PH referred for right heart catheterization (RHC) and same-day CMR from January 2013 to March 2016 at our Hospital. Patients with precapillary PH at RHC were eligible for inclusion in this analysis. Precapillary PH was defined as a mean pulmonary artery pressure ≥25 mmHg and pulmonary vascular resistance >3 Wood Units, with pulmonary artery wedge pressure ≤15 mmHg during right heart catheterization. 15 The research protocol was conducted in accordance with the Declaration of Helsinki and was approved by our Institutions Institutional Review Board. Informed consent was obtained from all patients. Age- and sex-matched subjects without cardiac and/or respiratory diseases from the MESA 4 study were included as controls for comparison. Baseline clinical characteristics were obtained from standardized clinic data prospectively obtained from our University Pulmonary Hypertension Registry.
CMR imaging
Volumetric cine images were acquired and analyzed to assess left ventricular (LV) and RV mass and function. All CMR images were read at our University by two investigators (C.P.C.-V. and S.L.Z.) blinded to the subjects and study timing. The cine imaging was steady-state free-precession cine imaging. A short-axis stack of 8–10 slices covering the entire ventricle, and 2 long-axis cine loops of 2- and 4-chamber views were obtained. Typical imaging parameters were slice thickness 6–8 mm, field of view 32–36 cm, minimum TR/TE, pixel size 1.4 x 1.4 mm, and temporal resolution <50 ms. CMR was performed within 5 h of catheterization.
Tissue-tracking CMR
We used offline semi-automated multimodality tissue-tracking software version 6.0 (Toshiba, Tokyo, Japan) to analyze the RV strain and SR. Multimodality tissue tracking reads characteristic pixel patterns in each 10 × 10 mm area as template blocks from the reference image. A search is performed in the next frame to find image blocks that best match the template block, such that the best match minimizes the mean squared error of the image pixel intensity (between current image block and template block). This procedure was repeated for all pixels in each image and for each frame to track moving pixels through the whole cardiac cycle.5,16,17 After uploading the CMR image, its brightness was optimized to ensure optimal endocardial/blood pool discrimination. The endocardial and epicardial borders of the ventricle were then manually traced on the end-diastolic frame, and the software automatically propagated the borders through the cardiac cycle. Adjustment of the tracked border was performed, as necessary, after visual assessment during cine loop playback by the reader.
RV longitudinal strain and SR data were directly obtained from cine images in the 4-chamber view (Fig. 1a). RV circumferential strain and SR data were obtained at mid-cavity level from short-axis images (Fig. 1b). A negative strain value indicates shortening with respect to the reference configuration at end-diastole, defined as the peak of the R wave on surface electrocardiography. Fig. 1c and d shows representative strain images of the control and PAH group, respectively. RV peak systolic strain was identified from the strain curve. Systolic time and diastolic time were also obtained from the strain analysis (Fig. 1c and d). Systolic to diastolic duration ratio (S:D ratio) was calculated as systolic time divided by diastolic time intervals from the strain curve. Fig. 1e and f shows representative SR images of the control and PAH subjects, respectively. Peak systolic strain rate (SRs), peak early diastolic strain rate (SRe), peak atrial-diastolic strain rate (SRa), and deceleration time (DT) of SRe were obtained from the SR curve.

(a) Representative image of right ventricular (RV) longitudinal peak systolic strain and peak strain rate in the 4-chamber view. (b) Representative image of RV circumferential peak strain and strain rate measurement in the short axis view at the mid-level. RV and left ventricular free wall were traced as a whole after the confirmation of six segments. Thereafter, three left ventricular free wall segments were excluded. (c) Representative strain curve image of control subjects. (d) Representative strain curve image of patients with pulmonary arterial hypertension. Peak systolic strain, systolic time, and diastolic time were identified from the strain curve. (e) Representative strain rate curve image of control subjects. (f) Representative strain rate image of image of patients with pulmonary arterial hypertension. Peak systolic strain rate, peak early diastolic strain rate (SRe), peak atrial-diastolic strain rate, and deceleration time of SRe were also analyzed from the strain rate curve.
Pressure volume analysis
The RV volume signal was calibrated to resting RV volume measurements obtained by same-day CMR. End-systolic pressure volume points were determined from a set of loops with varying preload volumes during phase 2 of the Valsalva maneuver, and fit by perpendicular regression using an iterative algorithm to derive the slope (end-systolic elastance, Ees, which is a gold standard for evaluating systolic function) of the end-systolic pressure volume relationship. This slope was then applied to the resting pressure volume loop to determine the x intercept (V0) of the end-systolic pressure volume relationship, as described. 12 Effective arterial elastance (Ea), which is a gold standard for assessing afterload, was determined by dividing end-systolic pressure by stroke volume. RV-pulmonary artery coupling was assessed by the Ees/Ea ratio. Tau, which is the time constant of RV relaxation, was also calculated. Representative images of micromanometer catheter measurements of Ees, Ea, and tau are shown in Supplemental Figure 1.
Data analysis
CMR measurements were compared between control subjects and PH patients. To identify CMR parameters most reflective of systolic function, correlation analyses were performed between hemodynamically derived Ees and global RV measures (RV ejection fraction, RV end-diastolic volume, and RV mass), as well as measures of RV systolic function (peak longitudinal systolic strain, peak longitudinal SRs, peak circumferential systolic strain, peak circumferential SRs, and S:D ratio). We also correlated Ea and Ees/Ea with these CMR parameters. To identify CMR parameters most reflective of diastolic function, correlation analyses were performed between hemodynamically derived tau, global RV measures, and measures of RV diastolic function (peak longitudinal SRe, peak longitudinal SRa, longitudinal DT of SRe, longitudinal S:D ratio, peak circumferential SRe, peak circumferential SRa, circumferential DT of SRe, and circumferential S:D ratio).
Statistical methods
Categorical variables are expressed as percentages, and continuous variables are expressed as mean ± SD for normally distributed variables or as medians with interquartile ranges for non-normal distributions. Departures from normality were detected with the Shapiro–Wilk statistic. In case of a difference in distribution between the control and PH groups, we added the interquartile range data to adequately compare these two groups. Differences in measurements between the control and PH groups were assessed with Fisher’s exact test, Student’s t test, or the Wilcoxon test as appropriate. Pearson correlation coefficients were used to determine correlations between continuous variables. Log transformations were performed as necessary to normalize distributions before correlation analysis. A two-sided P-value of <0.05 was considered significant throughout the analyses.
Given estimations based on previous studies 18 that the correlation coefficient between CMR-derived strain parameters and pressure volume loop-derived systolic and diastolic parameters was 0.6 or greater, using a two-sided test, 5% significance level (α = 0.05) with 80% power (β = 0.2), the required sample size was estimated to be approximately 19 individuals. 19
The intra-observer agreement for measurements was assessed by comparing the measurements of repeated analysis in 10 randomly chosen subjects (T.S.). The inter-observer agreement was assessed using the same patients by comparing the results measured by T.S. and those obtained by an experienced reader (T.F.). The second reader was not aware of the CMR measurements of the first examiner. Bland–Altman analysis and intraclass correlations between the two measurements were used to assess reproducibility.
All statistical analyses were performed using JMP® Version 11 (SAS Institute Inc., Cary, NC).
Results
Study population and demographics
Sixteen PAH patients and three patients with precapillary PH due to systemic sclerosis-related interstitial lung disease (SSc-ILD) met the entry criteria during the study period. Of the 16 PAH patients, six were diagnosed with idiopathic PAH, and 10 were diagnosed with SSc-associated PAH. Patients with SSc-associated PH met the American College of Rheumatology/European League Against Rheumatism diagnostic criteria for systemic sclerosis. 20 The diagnosis of PH due to SSc-ILD was based on previously reported criteria. 21 Nineteen age- and sex-matched healthy controls were included from the MESA 4 study. There were no significant demographic differences between the PH and control groups (Table 1).
Demographics of participants.
PH: pulmonary hypertension.
Among patients with PH, 18 (95%) were categorized as World Health Organization functional class II or III (Table 2). The mean pulmonary artery pressure was 33 mmHg, and pulmonary vascular resistance was 5.5 Wood Units (Table 2).
Characteristics of patients with pulmonary arterial hypertension.
IPAH: idiopathic PAH; SSc-PAH: scleroderma-associated PAH; SSc-ILD-PH: scleroderma with interstitial lung disease and PH; WHO: World Health Organization; PDE: phosphodiesterase; RV: right ventricular.
Comparison of CMR parameters between the control and PH groups
Table 1 shows the comparison of CMR parameters between the control and PH groups. PH patients exhibited significantly greater RV end-diastolic volume, RV end-systolic volume, and RV mass compared to controls. RV ejection fraction and LV mass were lower in PH patients than in control subjects.
For strain-related parameters, RV peak longitudinal systolic strain, RV peak longitudinal SRs, RV peak longitudinal SRe, and RV peak circumferential SRs were significantly decreased in PH patients than in controls. In contrast, RV peak longitudinal SRa and peak circumferential SRa were higher in PH patients than in controls. Both longitudinal and circumferential S:D ratios were higher in PH patients than in controls.
Correlations between invasive parameters and CMR indices of the right ventricle
There was no significant correlation between invasively measured Ees or Ees/Ea, and longitudinal CMR measures of systolic function; likewise, there was no correlation between Ees or Ees/Ea and circumferential CMR systolic measures. Log Ea was significantly correlated with RV longitudinal systolic strain (r = 0.46, P = 0.0498) alone (Supplemental Figure 2). Fig. 2 shows relationships between tau and CMR-derived strain, SR, and S:D ratio. For RV longitudinal analysis, tau was significantly correlated with RV peak longitudinal SRa (r = −0.61, P = 0.006), RV peak longitudinal SRe of DT (r = 0.75, P = 0.0002), and longitudinal S:D ratio (r = 0.59, P = 0.0083). In circumferential analysis, tau was significantly correlated with RV peak circumferential SRe (r = −0.5, P = 0.027), RV peak circumferential SRa (r = −0.52, P = 0.0257), and RV peak circumferential SRe of DT (r = 0.62, P = 0.0058).

Correlations between tau and cardiac magnetic resonance-derived strain and strain rate parameters. S:D: systolic to diastolic duration
Reproducibility
Limits of agreement analyzed from Bland–Altman analysis of intra-observer variability of strain and SR indices and intraclass correlations are shown in Supplemental Table 1. Intraclass correlations were greater than 0.9, and the limit of agreement was within 5% in all of the 20 indices. Limits of agreement and intraclass correlation coefficients for inter-observer variability are shown in Supplemental Table 2. Intraclass correlations were similarly high, with coefficients >0.9 in all indices except three (RV longitudinal SRe, 0.87, RV circumferential SRe, 0.87, and RV circumferential DT, 0.89). Limits of agreement were within 5% in all 20 indices.
Discussion
In the present study, we compared CMR measures of RV function, including strain parameters, between 19 PH patients and 19 controls, and we examined relationships between CMR-measured RV systolic and diastolic function and invasively measured gold standard parameters of RV systolic and diastolic function in PH patients. The main findings of this study were (1) CMR-derived strain-related parameters were substantially different between PH and control groups, both longitudinally and circumferentially; (2) there was no significant correlation between CMR parameters of RV systolic function and invasively measured Ees in PH patients; (3) CMR parameters of RV diastolic function were correlated with tau, with peak atrial-diastolic SR and DT being the strongest correlates; and (4) S:D ratio derived from CMR strain was correlated with tau, but not with Ees.
In PH patients, RV function is a major determinant of survival. Both systolic and diastolic functions play important roles in RV performance, and both should be considered as important individual determinants of overall function. However, most studies to date have focused on RV systolic function,4,8,22,23 with few studies concentrating on RV diastolic dysfunction in PH. Imaging measures of RV systolic performance, however, have many clinical pitfalls. RV ejection fraction, for example, is known to be an insensitive measure of RV contractile performance, and correlates poorly with intrinsic, load-independent measurements of RV contractile function.12,24 This is because preload and afterload significantly impact RV contractility as assessed by non-invasive imaging. Similar to these findings, the present study shows no association between CMR measures of systolic function and invasively derived measures of load-independent contractile function.
Although CMR measures of systolic function fall short, CMR-derived RV diastolic measurements may in fact prove more clinically useful. In the present study, we demonstrated that diastolic RV SR, as measured by CMR, was quite different between controls and PH patients. Moreover, multiple diastolic indices were significantly correlated with invasively measured tau. These findings indicate that CMR-derived strain analysis could allow for non-invasive evaluation of RV diastolic function, particularly relaxation represented by tau, in PH.
Among RV diastolic strain or SR examined in this study, peak longitudinal SRa and peak circumferential SRa were higher in PH patients than in controls and were most strongly correlated with tau. There are several possible reasons why both longitudinal and circumferential RV SRa were elevated in our PH patients. For example, RV SRa reflects RV stiffness and right atrial contraction. Thus, RV SRa can theoretically increase either by RV stiffening or by augmented right atrial contraction. Indeed, recent CMR studies on PAH patients have reported a stiffened right ventricle 25 and increased right atrial contractility. 26 Thus, sustained afterload to the right heart might have caused RV hypertrophy and stiffening and a compensatory increase in right atrial contraction, with a subsequent increase in SRa in our PH cohort.
The correlation analysis indicated a negative association between RV SRa and tau in both longitudinal and circumferential analysis. Our results are consistent with a previous study by Okumura et al., 18 where tau and SRa showed an inverse relationship. Notably, a recent study reported that right atrial contraction was increased in patients with mild PH compared to controls, yet decreased in patients with advanced PH. 26 This “non-linear” relationship between right atrial contraction and the progression of PH might have caused an inverse relationship between RV SRa and tau. Another explanation could be that the two indices reflect different phases of RV diastole; tau is an early diastolic index, whereas RV SRa is a late diastolic parameter. This may have caused the lack of a linear positive relationship between the two indices seen in our study.
Previous echocardiographic studies have shown a mild correlation between longitudinal early diastolic SR and tau in children with congenital heart disease-associated PAH. 18 In contrast, in the present study, there was no significant association between the two indices. This may be due to differences in age or PH type of the participants, or to the different modalities used, i.e. echocardiography versus CMR, in the two studies. However, there was a significant correlation between circumferential peak SRe and tau in our CMR study, suggesting that circumferential strain indices may reflect RV diastolic dysfunction in a different manner compared with longitudinal parameters. CMR has become an established modality and the gold standard for evaluating right heart in patients with PH, but in recent years, 3D echocardiography has made it possible to evaluate strain in the entire right ventricle and in various directions,27–29 which make assessment of right ventricular contractility and diastolic function possible.
Tello et al. 30 have recently documented significant associations of RV longitudinal or radial peak strain with Ees/Ea, Ea, and Tau in 38 precapillary PH patients. Similar to our study, these findings indicate potential reliability of CMR-derived strain analysis for the assessment of RV function in PH patients. However, only peak strain was analyzed in this study by Tello et al., whereas strain rate was not. Thus, further evaluation for more appropriate use of different types of strain and strain rate indices is needed.
S:D ratio is an indicator of global ventricular performance. Recently, high S:D ratio of the left ventricle was found associated with poor prognosis in patients with dilated cardiomyopathy, and a similar predictive value of RV S:D ratio has also been demonstrated in PAH.31,32 In our study, a high longitudinal S:D ratio was correlated with an increased tau, consistent with the recent echocardiographic report by Okumura et al. 18 A similar positive correlation with tau was observed for DT in the present study. These findings suggest that time-related indices of CMR strain analysis could be used to evaluate RV diastolic function and prognosis in PAH.
Studies on patients with left heart disease 33 and rodent experimental models34,35 have demonstrated significant correlations between pressure volume loop analysis-derived Ees and echocardiography-derived peak systolic strain or SR of the left ventricle. However, in our study, significant correlations between Ees and CMR-derived parameters of RV systolic function, such as peak systolic strain and SR, were not found. Tello et al. 30 also failed to find a correlation between Ees and peak RV strain rate. The inconsistency between previous studies and ours may be due to the difference in the ventricle examined, i.e. left versus right ventricle, or to vendor- 36 or modality-specific variability, i.e. echocardiography versus CMR. In addition, there is a possibility that indices of strain analysis may reflect subtle or certain specific changes in RV function that cannot be detected by Ees. Importantly, a recent report by Kraigher-Krainer et al. 37 revealed decreased left ventricular longitudinal and circumferential strain indices in patients with heart failure whose LV ejection fraction was preserved. In addition, recent studies have reported significant associations between RV peak systolic strain and disease severity6,38 or mortality7,8 in PH. Thus, the clinical relevance of systolic strain and SR needs to be further evaluated in diverse PH types and in larger prospective studies.
In the present study, RV longitudinal peak systolic strain was significantly correlated with Ea. This is consistent with a recent echocardiographic study demonstrating a close correlation between RV free wall strain and pulmonary arterial pressure in PAH. 39 Taken together, these studies suggest that RV strain during systole is increased in response to an elevated RV afterload, and could be a unique marker of RV response to an increased afterload in PH. Alternatively, strain-related indices were not correlated with Ees/Ea in our study. Ees/Ea is a powerful load-independent indicator of RV-pulmonary artery coupling.12,24 Thus, our data indicated a limited sensitivity of RV strain-related CMR parameters for the assessment of RV-pulmonary artery coupling in PH.
While this was a prospective study for which a power analysis was performed to calculate the sample size needed to establish correlations between CMR parameters and invasive parameters, it included a small number of patients from a single center. In addition, for obvious reasons, pressure volume loops were not obtained in control subjects. Furthermore, CMR and right heart catheterization were not performed contemporaneously, which potentially hampered an accurate comparison of the techniques. However, all patients were in stable clinical condition, and the two studies were performed within 5 h. This study included a combination of patients with idiopathic PAH, SSc-associated PAH, and precapillary PH due to SSc-ILD, which might have also affected the results. Finally, variations in treatment regimens might have also influenced the results, as PAH-specific vasodilators may affect cardiac function. 40
In conclusion, the present study demonstrated differences in various CMR parameters including strain and SR between control subjects and PH patients. Notably, there was no significant correlation between Ees, the gold standard for the assessment of systolic function, and CMR strain parameters of RV systolic function. In contrast, there were significant correlations between tau and CMR-derived strain measures of the right ventricle, such as longitudinal SRe, longitudinal and circumferential SRa, DT of SRe, and longitudinal S:D ratio. These suggest that strain analysis using CMR is most appropriate to assess diastolic function in PAH.
Supplemental Material
sj-pdf-1-pul-10.1177_20458940211032529 - Supplemental material for Right ventricular function as assessed by cardiac magnetic resonance imaging-derived strain parameters compared to high-fidelity micromanometer catheter measurements
Supplemental material, sj-pdf-1-pul-10.1177_20458940211032529 for Right ventricular function as assessed by cardiac magnetic resonance imaging-derived strain parameters compared to high-fidelity micromanometer catheter measurements by Takahiro Sato, Bharath Ambale-Venkatesh, Stefan L. Zimmerman, Ryan J. Tedford, Steven Hsu, Ela Chamera, Tomoki Fujii, Christopher J. Mullin, Valentina Mercurio, Rubina Khair, Celia P. Corona-Villalobos, Catherine E. Simpson, Rachel L. Damico, Todd M. Kolb, Stephen C. Mathai, Joao A.C. Lima, David A. Kass, Ichizo Tsujino and Paul M. Hassoun in Pulmonary Circulation
Supplemental Material
sj-pdf-2-pul-10.1177_20458940211032529 - Supplemental material for Right ventricular function as assessed by cardiac magnetic resonance imaging-derived strain parameters compared to high-fidelity micromanometer catheter measurements
Supplemental material, sj-pdf-2-pul-10.1177_20458940211032529 for Right ventricular function as assessed by cardiac magnetic resonance imaging-derived strain parameters compared to high-fidelity micromanometer catheter measurements by Takahiro Sato, Bharath Ambale-Venkatesh, Stefan L. Zimmerman, Ryan J. Tedford, Steven Hsu, Ela Chamera, Tomoki Fujii, Christopher J. Mullin, Valentina Mercurio, Rubina Khair, Celia P. Corona-Villalobos, Catherine E. Simpson, Rachel L. Damico, Todd M. Kolb, Stephen C. Mathai, Joao A.C. Lima, David A. Kass, Ichizo Tsujino and Paul M. Hassoun in Pulmonary Circulation
Footnotes
Conflict of interest
The author(s) declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: PMH has served on a Scientific Advisory Board for Merck & Co..
Authors’ contributions
TS, BA-V, RJT, SH, VM, CJM, VM, JACL, PMH, DAK, and IC intellectually contributed to the analysis and interpretation of data; SLJ, EC, and TF managed the overall CMR imaging and assessment; RK, CES, RLD, TMK, and SCM contributed to the design data and acquisition; PMH gave final approval of the manuscript submitted.
Ethical approval
The research protocol was approved by the Johns Hopkins Medical Institutions Institutional Review Board, and informed consent was obtained from all patients.
Guarantor
Paul M. Hassoun.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship and/or publication of this article: The study was supported by The National Institutes of Health/National Heart, Lung, and Blood Institute awards P50 HL084946 (PMH), HL114910 (PMH), and U01HL125175-01 (PMH), and by a fellowship grant from the Uehara Memorial Foundation, the Mochida Memorial Foundation for Medical and Pharmaceutical Research, the Ito Foundation for the Promotion of Medical Science, and the Japanese Society for the Promotion of Science (TS).
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
