Quantitative Assessment of Left and Right Atrial Strains Using Cardiovascular Magnetic Resonance Based Tissue Tracking

Feasibility and reproducibility of cardiovascular magnetic resonance-feature tracking for quantitative right atrial function in dilated cardiomyopathy patients

Background

The importance of right heart assessment in dilated cardiomyopathy (DCM) is increasingly recognized. The development of cardiovascular magnetic resonance-feature tracking (CMR-FT) has provided a novel approach to quantify myocardial deformation and evaluate cardiac function. In this study, we aimed to evaluate the feasibility and reproducibility of CMR-FT for the quantitative derivation of right atrial (RA) strain and strain rate (SR) in patients with DCM.

Methods

A total of 68 DCM patients (84% male; aged 50.6±13.2 years) and 58 healthy controls (81% male; aged 48.4±11.2 years) were retrospectively enrolled from September 2018 to August 2022 at the First Affiliated Hospital of Zhejiang Chinese Medical University and Shenzhen Clinical Medical College of Guangzhou University of Chinese Medicine. RA reservoir, conduit, and booster strain (εs, εe, and εa) and peak positive, peak early negative, and peak late negative SR (SRs, SRe, and SRa) were measured using CMR-FT and compared between 2 groups using Student's t-test. Intra- and inter-observer reproducibility was evaluated using intraclass correlation coefficients (ICC) and Bland-Altman plots.

Results

Compared to healthy controls, DCM patients showed significantly lower RA strain (εs: 19.7%±9.0% vs. 44.4%±9.7%; εe: 7.9%±5.3% vs. 25.8%±8.6%; εa: 11.8%±6.2% vs. 18.6%±5.1%, all P<0.001) and SR (SRs: 1.17±0.48 vs. 1.92±0.62 s-1; SRe: -0.85±0.56 vs. -1.94±0.63 s-1; SRa: -1.39±0.71 vs. -2.01±0.65 s-1, all P<0.001). There was no significant difference in RA maximum volume index between the 2 groups. Simple linear regression analysis demonstrated a significant correlation between N-terminal B-type natriuretic peptide (NT-proBNP), RA emptying fraction passive (RAEF passive), and RA εe [(NT-proBNP and εe): r=-0.48, P<0.001, 95% confidence interval (CI): -0.64 to -0.26; and (RAEF passive and εe): r=0.41, P=0.001, 95% CI: 0.22 to 0.56, respectively] in DCM patients. Intra- and inter-observer reproducibility was excellent (all ICCs >0.85) for RA deformation measurements.

Conclusions

CMR-FT is a promising, noninvasive approach for the quantitative assessment of RA phasic function in patients with DCM. DCM patients exhibit impaired RA reservoir, conduit, and booster pump function prior to visible RA enlargement.

Myocardial Strain Measurements Derived From MR Feature-Tracking: Influence of Sex, Age, Field Strength, and Vendor

BACKGROUND

Cardiac magnetic resonance feature tracking (CMR-FT) is a novel technique for assessing myocardial deformation and dysfunction. However, a comprehensive assessment of normal values of strain parameters in all 4 cardiac chambers using different vendors is lacking.

OBJECTIVES

This study aimed to characterize the normal values for myocardial strain in all 4 cardiac chambers and identify factors that contribute to variations in FT strain through a systematic review and meta-analysis of the CMR-FT published reports.

METHODS

The investigators searched PubMed, Embase, and Scopus for myocardial strains of all 4 chambers measured by CMR-FT in healthy adults. The pooled means of all strain parameters were generated using a random-effects model. Subgroup analyses and meta-regressions were performed to identify the sources of variations.

RESULTS

This meta-analysis included 44 studies with a total of 3,359 healthy subjects. The pooled means of left ventricular global longitudinal strain (LV-GLS), LV global radial strain, and LV global circumferential strain (GCS) were -18.4% (95% CI: -19.2% to -17.6%), 43.7% (95% CI: 40.0%-47.4%), and -21.4% (95% CI: -22.3% to -20.6%), respectively. The pooled means of left atrial (LA)-GLS (corresponding to total strain, passive strain, and active strain) were 34.9% (95% CI: 29.6%-40.2%), 21.3% (95% CI: 16.6%-26.1%) and 14.3% (95% CI: 11.8%-16.8%), respectively. The pooled means of right ventricular (RV)-GLS and right atrial global longitudinal total strain were -24.0% (95% CI: -25.8% to -22.1%) and 36.3% (95% CI: 15.5%-57.0%), respectively. Meta-regression identified field strength (P < 0.001; I2 = 98.6%) and FT vendor (P < 0.001; I2 = 98.5%) as significant confounders contributing to heterogeneity of LV-GLS. The variations of LA-GLSactive were associated with regional distribution (P < 0.001; I2 = 97.3%) and FT vendor (P < 0.001; I2 = 97.4%). Differences in FT vendor were attributed to variations of LV-GCS and RV-GLS (P = 0.02; I2 = 98.8% and P = 0.01; I2 = 93.8%).

CONCLUSIONS

This study demonstrated the normal values of CMR-FT strain parameters in all 4 cardiac chambers in healthy subjects. Differences in FT vendor contributed to the heterogeneity of LV-GLS, LV-GCS, LA-GLSactive, and RV-GLS, whereas sex, age, and MR vendor had no effect on the normal values of CMR-FT strain measurements.

Association of serum uric acid with right cardiac chamber remodeling assessed by cardiovascular magnetic resonance feature tracking in patients with connective tissue disease

Background

Right cardiac chamber remodeling is widespread in patients with connective tissue disease (CTD). Serum uric acid (SUA) is considered a potential independent risk factor for cardiovascular disease, and elevated SUA levels are often observed in patients with CTD. The correlation between SUA levels and right cardiac chamber remodeling remains unclear. This study investigated the association of SUA with right cardiac chamber remodeling as assessed by cardiac magnetic resonance feature-tracking (CMR-FT) in CTD patients.

Methods and results

In this cross-sectional study, a total of 104 CTD patients and 52 age- and sex-matched controls were consecutively recruited. All individuals underwent CMR imaging, and their SUA levels were recorded. The patients were divided into three subgroups based on the tertiles of SUA level in the present study. CMR-FT was used to evaluate the right atrial (RA) longitudinal strain and strain rate parameters as well as right ventricular (RV) global systolic peak strain and strain rate in longitudinal and circumferential directions for each subject. Univariable and multivariable linear regression analyses were used to explore the association of SUA with RV and RA strain parameters. Compared with the controls, the CTD patients showed significantly higher SUA levels but a lower RV global circumferential strain (GCS) and RA phasic strain parameters (all p < 0.05), except the RA booster strain rate. RV GCS remained impaired even in CTD patients with preserved RV ejection fraction. Among subgroups, the patients in the third tertile had significantly impaired RV longitudinal strain (GLS), RV GCS, and RA reservoir and conduit strain compared with those in the first tertile (all p < 0.05). The SUA levels were negatively correlated with RV GLS and RV GCS as well as with RA reservoir and conduit strain and strain rates (the absolute values of β were 0.250 to 0.293, all P < 0.05). In the multivariable linear regression analysis, the SUA level was still an independent determinant of RA conduit strain (β = -0.212, P = 0.035) and RV GCS (β = 0.207, P = 0.019).

Conclusion

SUA may be a potential risk factor of right cardiac chamber remodeling and is independently associated with impaired RA conduit strain and RV GCS in CTD patients.

Associations of Circulating Vascular Cell Adhesion Molecule-1 and Intercellular Adhesion Molecule-1 With Long-Term Cardiac Function

BACKGROUND

Although VCAM-1 (vascular cell adhesion molecule-1) and ICAM-1 (intercellular adhesion molecule-1) have been associated with incident heart failure with preserved ejection fraction (HFpEF) and atrial fibrillation (AF), the associations of VCAM-1 and ICAM-1 with sensitive measures of cardiac structure/function are unclear. The objective of this study is to evaluate associations between VCAM-1, ICAM-1, and measures of cardiac structure and function as potential pathways through which cellular adhesion molecules promote HFpEF and AF risk.

METHODS AND RESULTS

In MESA (Multi-Ethnic Study of Atherosclerosis), we evaluated the associations of circulating VCAM-1 and ICAM-1 at examination 2 (2002-2004) with measures of cardiac structure/function on cardiac magnetic resonance imaging at examination 5 (2010-2011) after multivariable adjustment. Mediation analysis of left atrial (LA) strain on the association between VCAM-1 or ICAM-1 and AF or HFpEF was also performed. Overall, 2304 individuals (63±10 years; 47% men) with VCAM-1 or ICAM-1, cardiac magnetic resonance imaging, and covariate data were included in analysis. Higher VCAM-1 and ICAM-1 were associated with lower LA peak longitudinal strain and worse global circumferential left ventricular strain but were not associated with left ventricular myocardial scar or interstitial fibrosis. Lower LA peak longitudinal strain mediated 8% (95% CI, 2-30) of the relationship between VCAM-1 and HFpEF and 9% (95% CI, 2-21) of the relationship between VCAM-1 and AF.

CONCLUSIONS

Higher VCAM-1 and ICAM-1 were associated with lower LA function and left ventricular systolic function but were not associated with myocardial scar or interstitial fibrosis. VCAM-1 and ICAM-1 may promote HFpEF and AF risk through impaired LA reservoir function.

Heart Failure Preserved Ejection Fraction in Women: Insights Learned from Imaging

While the prevalence of heart failure, in general, is similar in men and women, women experience a higher rate of HFpEF compared to HFrEF. Cardiovascular risk factors, parity, estrogen levels, cardiac physiology, and altered response to the immune system may be at the root of this difference. Studies have found that in response to increasing age and hypertension, women experience more concentric left ventricle remodeling, more ventricular and arterial stiffness, and less ventricular dilation compared to men, which predisposes women to developing more diastolic dysfunction. A multi-modality imaging approach is recommended to identify patients with HFpEF. Particularly, appreciation of sex-based differences as described in this review is important in optimizing the evaluation and care of women with HFpEF.

Personalized biomechanical insights in atrial fibrillation: opportunities & challenges

INTRODUCTION

Atrial fibrillation (AF) is an increasingly prevalent and significant worldwide health problem. Manifested as an irregular atrial electrophysiological activation, it is associated with many serious health complications. AF affects the biomechanical function of the heart as contraction follows the electrical activation, subsequently leading to reduced blood flow. The underlying mechanisms behind AF are not fully understood, but it is known that AF is highly correlated with the presence of atrial fibrosis, and with a manifold increase in risk of stroke.

AREAS COVERED

In this review, we focus on biomechanical aspects in atrial fibrillation, current and emerging use of clinical images, and personalized computational models. We also discuss how these can be used to provide patient-specific care.

EXPERT OPINION

Understanding the connection betweenatrial fibrillation and atrial remodeling might lead to valuable understanding of stroke and heart failure pathophysiology. Established and emerging imaging modalities can bring us closer to this understanding, especially with continued advancements in processing accuracy, reproducibility, and clinical relevance of the associated technologies. Computational models of cardiac electromechanics can be used to glean additional insights on the roles of AF and remodeling in heart function.

Differing Impact of Preterm Birth on the Right and Left Atria in Adulthood

Background Preterm birth affects 10% of live births and is associated with an altered left ventricular and right ventricular phenotype and increased cardiovascular disease risk in young adulthood. Because left atrial (LA) and right atrial (RA) volume and function are known independent predictors of cardiovascular outcomes, we investigated whether these were altered in preterm-born young adults. Methods and Results Preterm-born (n=200) and term-born (n=266) adults aged 18 to 39 years underwent cardiovascular magnetic resonance imaging. LA and RA maximal and minimal volumes (absolute, indexed to body surface area, and as a ratio to ventricular volumes) were obtained to study atrial morphology, while LA and RA stroke volume, strain, and strain rate were used to assess atrial function. Secondary analyses consisted of between-group comparisons based on degree of prematurity. Absolute RA volumes and RA volumes indexed to right ventricular volumes were significantly smaller in preterm-born compared with term-born adults. In addition, RA reservoir and booster strain were higher in preterm-born adults, possibly indicating functional compensation for the smaller RA volumes. LA volumes indexed to left ventricular volumes were significantly greater in preterm-born adults as compared with term-born adults, although absolute LA volumes were similar between groups. LA and RA changes were observed across gestational ages in the preterm group but were greatest in those born very-to-extremely preterm. Conclusions Preterm-born adults show changes in LA and RA structure and function, which may indicate subclinical cardiovascular disease. Further research into underlying mechanisms, opportunities for interventions, and their prognostic value is warranted.

Cardiovascular Magnetic Resonance Imaging-Based Right Atrial Strain Analysis of Cardiac Amyloidosis

Background: Cardiac amyloidosis (CA) manifests in a hypertrophic phenotype with a poor prognosis, making differentiation from hypertrophic cardiomyopathy (HCM) challenging and delaying early treatment. The extent to which magnetic resonance imaging (MRI) quantifies the right atrial strain (RAS) and strain rate (RASR), providing valuable diagnostic information, is not yet clinically established. Aims: This study assesses diagnostic differences in the longitudinal RAS and RASR between CA and HCM patients, control subjects (CTRL) and CA subtypes in addition to the impact of atrial fibrillation (AF) on the right atrial function in CA patients. The RAS and RASR of tricuspid regurgitation (TR) patients are used to assess the potential for diagnostic overlap. Methods: RAS and RASR quantification was conducted via MRI feature-tracking for biopsy-confirmed CA patients with subtypes identified. Strain parameters were compared for CTRL, HCM and TR patients. Post hoc testing identified intergroup differences. Results: In total, 41 CA patients were compared to 47 CTRL, 20 HCM and 31 TR patients. Reservoir (R), conduit and booster RAS and RASRs allow for significant differentiation (p < 0.001) between CA and HCM patients (R: 10.6 ± 14.3% vs. R: 33.5 ± 16.3%) and CTRL (R: 44.6 ± 15.7%). Booster and reservoir RAS and RASRs qualified as reliable diagnostic tests (AUC > 0.8). CA patients with AF, in contrast to sinus rhythm, demonstrated a significantly impaired reservoir RAS and RASR and booster RASR. The discriminative power of RAS for CA vs. TR was insufficient (R: 10.6% ± 14.3% vs. 7.0% ± 6.0%, p = 0.069). Differentiation between 21 transthyretin and 20 light-chain amyloidosis subtypes was not achievable (R: 0.7% ± 1.0% vs. 0.7% ± 1.0%, p = 0.827). Conclusion: The MRI-derived RAS and RASR are impaired in CA patients and may support noninvasive differentiation between CA, HCM and CTRL.

Multi-parametric analyses to investigate dependencies of normal left atrial strain by cardiovascular magnetic resonance feature tracking

Left-atrial (LA) strain is the result of complex hemodynamics, which may be better characterized using a multiparametric approach. Cardiovascular magnetic resonance (CMR) feature tracking was used to perform a comprehensive LA strain assessment of 183 enrolled healthy volunteers (11–70 years, 97 females, median 32.9 ± 28.3 years). Novel strain dependencies were assessed using multi-parametric regression (MPR) analyses. LA volumetric data, left ventricular strain, transmitral and pulmonary venous blood flow parameters were utilized to create clusters for MPR of all subjects and a heart rate controlled subgroup (pulse: 60–75/min, N = 106). The LA reservoir(r) and conduit(c) strains of the total cohort were significantly elevated (p ≤ 0.001) in women (r: 49.7 ± 12.9%, c: 32.0 ± 11.0%) compared to men (r: 42.9 ± 11.4%, c: 26.1 IQ 10.5%). In contrast, there were no gender-specific differences (p > 0.05) for subgroup LA reservoir, conduit and booster(b) strains (all, r: 47.3 ± 12.7%; c: 29.0 IQ 15.5%; b: 17.6 ± 5.4%) and strain rates (all, 2.1 IQ 1.0 s⁻¹; − 2.9 IQ 1.5 s⁻¹; − 2.3 IQ 1.0 s⁻¹). MPR found large effect sizes (|R²|≥ 0.26) for correlations between strain and various cardiac functional parameters. Largest effect size was found for the association between LA conduit strain and LA indexed booster volume, LA total ejection fraction, left ventricular global radial strain and E-wave (|R²|= 0.437). In addition to providing normal values for sex-dependent LA strain and strain rate, no gender differences were found with modified heart rate. MPR analyses of LA strain/strain rate and various cardiac functional parameters revealed that heart rate control improved goodness-of-fit for the overall model.

Cardiac Computed Tomography-Derived Left Atrial Strain and Volume in Pediatric Patients With Congenital Heart Disease: A Comparative Analysis With Transthoracic Echocardiography

Purpose

This study aimed at exploring the feasibility and reproducibility of CCT for the measurement of Left Atrial (LA) strain and volume compared with transthoracic echocardiography (TTE) in pediatric patients with congenital heart disease (CHD).

Materials and Methods

The present study included 43 postoperative patients with CHD (7.39 ± 3.64 years, 56% male) who underwent clinically indicated CCT, and all patients underwent additional TTE on the same day. LA strain and volume parameters were measured by dedicated software. The correlation and agreement of LA strain and volume parameters were assessed using Pearson’s correlation coefficient and Bland-Altman analysis. Intra-class correlation coefficients (ICC) were used to assess CCT intra-observer and inter-observer reproducibility.

Results

All strain parameters of CCT were lower compared to TTE (reservoir strain: 28.37 ± 6.92 vs. 32.15 ± 8.15, respectively; conduit strain: 21.33 ± 6.46 vs. 24.23 ± 7.75, respectively; booster strain: 7.04 ± 2.74 vs. 7.92 ± 3.56). While the volume parameters of CCT were higher compared to TTE (LAV: 29.60 ± 19.01 vs. 25.66 ± 17.60, respectively; LAVi: 30.36 ± 22.31 vs. 28.63 ± 19.25, respectively). Both LA strain and volume measurements showed good correlation and agreement between the two modalities (r = 0.63–0.87, p < 0.001). CT-derived LA strain and volume measurements showed good intra- and inter-observer reproducibility using prototype software (ICC = 0.78–0.96).

Conclusions

CCT was feasible for measuring LA strain and volume with good correlation and high reproducibility as compared with TTE. As a complementary modality, CCT can regard as an accepted method in the evaluation of LA function in pediatric patients with CHD

Training and clinical testing of artificial intelligence derived right atrial cardiovascular magnetic resonance measurements

BACKGROUND

Right atrial (RA) area predicts mortality in patients with pulmonary hypertension, and is recommended by the European Society of Cardiology/European Respiratory Society pulmonary hypertension guidelines. The advent of deep learning may allow more reliable measurement of RA areas to improve clinical assessments. The aim of this study was to automate cardiovascular magnetic resonance (CMR) RA area measurements and evaluate the clinical utility by assessing repeatability, correlation with invasive haemodynamics and prognostic value.

METHODS

A deep learning RA area CMR contouring model was trained in a multicentre cohort of 365 patients with pulmonary hypertension, left ventricular pathology and healthy subjects. Inter-study repeatability (intraclass correlation coefficient (ICC)) and agreement of contours (DICE similarity coefficient (DSC)) were assessed in a prospective cohort (n = 36). Clinical testing and mortality prediction was performed in n = 400 patients that were not used in the training nor prospective cohort, and the correlation of automatic and manual RA measurements with invasive haemodynamics assessed in n = 212/400. Radiologist quality control (QC) was performed in the ASPIRE registry, n = 3795 patients. The primary QC observer evaluated all the segmentations and recorded them as satisfactory, suboptimal or failure. A second QC observer analysed a random subcohort to assess QC agreement (n = 1018).

RESULTS

All deep learning RA measurements showed higher interstudy repeatability (ICC 0.91 to 0.95) compared to manual RA measurements (1st observer ICC 0.82 to 0.88, 2nd observer ICC 0.88 to 0.91). DSC showed high agreement comparing automatic artificial intelligence and manual CMR readers. Maximal RA area mean and standard deviation (SD) DSC metric for observer 1 vs observer 2, automatic measurements vs observer 1 and automatic measurements vs observer 2 is 92.4 ± 3.5 cm2, 91.2 ± 4.5 cm2 and 93.2 ± 3.2 cm2, respectively. Minimal RA area mean and SD DSC metric for observer 1 vs observer 2, automatic measurements vs observer 1 and automatic measurements vs observer 2 was 89.8 ± 3.9 cm2, 87.0 ± 5.8 cm2 and 91.8 ± 4.8 cm2. Automatic RA area measurements all showed moderate correlation with invasive parameters (r = 0.45 to 0.66), manual (r = 0.36 to 0.57). Maximal RA area could accurately predict elevated mean RA pressure low and high-risk thresholds (area under the receiver operating characteristic curve artificial intelligence = 0.82/0.87 vs manual = 0.78/0.83), and predicted mortality similar to manual measurements, both p < 0.01. In the QC evaluation, artificial intelligence segmentations were suboptimal at 108/3795 and a low failure rate of 16/3795. In a subcohort (n = 1018), agreement by two QC observers was excellent, kappa 0.84.

CONCLUSION

Automatic artificial intelligence CMR derived RA size and function are accurate, have excellent repeatability, moderate associations with invasive haemodynamics and predict mortality.

The Prognostic Value of Right Atrial and Right Ventricular Functional Parameters in Systemic Sclerosis

Introduction

Right ventricular (RV) function is of particular importance in systemic sclerosis (SSc), since common SSc complications, such as interstitial lung disease and pulmonary hypertension may affect RV afterload. Cardiovascular magnetic resonance (CMR) is the gold standard for measuring RV function. CMR-derived RV and right atrial (RA) strain is a promising tool to detect subtle changes in RV function, and might have incremental value, however, prognostic data is lacking. Therefore, the aim of this study was to evaluate the prognostic value of RA and RV strain in SSc.

Methods

In this retrospective study, performed at two Dutch hospitals, consecutive SSc patients who underwent CMR were included. RV longitudinal strain (LS) and RA strain were measured. Unadjusted cox proportional hazard regression analysis and likelihood ratio tests were used to evaluate the association and incremental value of strain parameters with all-cause mortality.

Results

A total of 100 patients (median age 54 [46–64] years, 42% male) were included. Twenty-four patients (24%) died during a follow-up of 3.1 [1.8–5.2] years. RA reservoir [Hazard Ratio (HR) = 0.95, 95% CI 0.91–0.99, p = 0.009] and conduit strain (HR = 0.93, 95% CI 0.88–0.98, p = 0.008) were univariable predictors of all-cause mortality, while RV LS and RA booster strain were not. RA conduit strain proved to be of incremental value to sex, atrial fibrillation, NYHA class, RA maximum volume indexed, and late gadolinium enhancement (p < 0.05 for all).

Conclusion

RA reservoir and conduit strain are predictors of all-cause mortality in SSc patients, whereas RV LS is not. In addition, RA conduit strain showed incremental prognostic value to all evaluated clinical and imaging parameters. Therefore, RA conduit strain may be a useful prognostic marker in SSc patients.

Normal Ranges of Right Atrial Strain and Strain Rate by Two-Dimensional Speckle-Tracking Echocardiography: A Systematic Review and Meta-Analysis

Background: Normal range values of right atrial (RA) phasic function markers are essential for the identification of normal and abnormal values, comparison with reference values, and the clinical meaning of obtained values. Accordingly, we aimed to define the normal range values of RA phasic function markers obtained by 2D speckle-tracking echocardiography through a meta-analysis and determine the main sources of heterogeneity among reported values. Methods: PUBMED, SCOPUS, and EMBASE databases were searched for the following keywords: "right atrial/right atrium" and "strain/speckle/deformation" and "echocardiography." Studies were selected that included a human healthy adult group without any cardiovascular diseases or risk factors and that were written in the English language. For the calculation of each marker of RA phasic functions, a random-effect model was used. Meta-regression was employed to define the major sources of variabilities among reported values. Results: Fifteen studies that included 2,469 healthy subjects were selected for analysis. The normal range values for RA strain and strain rate were 42.7% (95% CI, 39.4 to 45.9%) and 2.1 s^-1 (95% CI, 2.0 to 2.1 s^-1) during the reservoir phase, respectively, 23.6% (95% CI, 20.7 to 26.6%) and -1.9 s^-1 (95% CI, -2.2 to -1.7 s^-1) during the conduit phase, correspondingly, and 16.1% (95% CI, 13.6 to 18.6%) and -1.8 s^-1 (95% CI, -2.0 to -1.5 s^-1) during the contraction phase, respectively. The sources of heterogeneity for the normal range of these markers were the number of participants, the type of software, the method of global value calculation, the right ventricular fractional area change, the left ventricular (LV) ejection fraction, the RA volume index, sex, the heart rate, the diastolic blood pressure, the body mass index, and the body surface area. Conclusions: Using 2D speckle-tracking echocardiography, we defined normal values for RA phasic function markers and identified the sources of heterogeneity as demographic, anthropometric, hemodynamic, and echocardiography factors. Systematic Review Registration: https://www.crd.york.ac.uk/prospero/display_record.php?ID=CRD42021236578, identifier: CRD42021236578.