Skip to main content

ORIGINAL RESEARCH article

Front. Med., 14 March 2022
Sec. Pulmonary Medicine
Volume 9 - 2022 | https://doi.org/10.3389/fmed.2022.840196

CMR Measures of Left Atrial Volume Index and Right Ventricular Function Have Prognostic Value in Chronic Thromboembolic Pulmonary Hypertension

Yousef Shahin1,2* Samer Alabed1,2 Syed Rehan Quadery3 Robert A. Lewis3 Christopher Johns2 Dheyaa Alkhanfar1 Maria Sukhanenko2 Faisal Alandejani1 Pankaj Garg1 Charlie A. Elliot3 Abdul Hameed1,3 Athaniosis Charalampopoulos3 James M. Wild1,4 Robin Condliffe3 Andrew J. Swift1,2,4 David G. Kiely1,3,4
  • 1Department of Infection, Immunity and Cardiovascular Disease, University of Sheffield, Sheffield, United Kingdom
  • 2Department of Clinical Radiology, Sheffield Teaching Hospitals NHS FT, Sheffield, United Kingdom
  • 3Sheffield Pulmonary Vascular Disease Unit, Royal Hallamshire Hospital, Sheffield Teaching Hospitals NHS FT, Sheffield, United Kingdom
  • 4INSIGNEO, Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom

Providing prognostic information is important when counseling patients and planning treatment strategies in chronic thromboembolic pulmonary hypertension (CTEPH). The aim of this study was to assess the prognostic value of gold standard imaging of cardiac structure and function using cardiac magnetic resonance imaging (CMR) in CTEPH. Consecutive treatment-naive patients with CTEPH who underwent right heart catheterization and CMR between 2011 and 2017 were identified from the ASPIRE (Assessing-the-Specturm-of-Pulmonary-hypertensIon-at-a-REferral-center) registry. CMR metrics were corrected for age and sex where appropriate. Univariate and multivariate regression models were generated to assess the prognostic ability of CMR metrics in CTEPH. Three hundred and seventy-five patients (mean+/-standard deviation: age 64+/-14 years, 49% female) were identified and 181 (48%) had pulmonary endarterectomy (PEA). For all patients with CTEPH, left-ventricular-stroke-volume-index-%predicted (LVSVI%predicted) (p = 0.040), left-atrial-volume-index (LAVI) (p = 0.030), the presence of comorbidities, incremental shuttle walking test distance (ISWD), mixed venous oxygen saturation and undergoing PEA were independent predictors of mortality at multivariate analysis. In patients undergoing PEA, LAVI (p < 0.010), ISWD and comorbidities and in patients not undergoing surgery, right-ventricular-ejection-fraction-%predicted (RVEF%pred) (p = 0.040), age and ISWD were independent predictors of mortality. CMR metrics reflecting cardiac function and left heart disease have prognostic value in CTEPH. In those undergoing PEA, LAVI predicts outcome whereas in patients not undergoing PEA RVEF%pred predicts outcome. This study highlights the prognostic value of imaging cardiac structure and function in CTEPH and the importance of considering left heart disease in patients considered for PEA.

Introduction

Pulmonary hypertension (PH) is heterogeneous and treatment depends on the underlying cause (1). Chronic thromboembolic pulmonary hypertension (CTEPH) is a potentially curable form of PH and a recent meta-analysis has identified a cumulative incidence of 2.9% in patients surviving an acute pulmonary embolism (2). It can also present as PH with no previous evidence of venous thromboembolism (3). It is characterized by non-resolution of thrombus and remodeling of the pulmonary arteries resulting in PH and right ventricular (RV) dysfunction (4) and without treatment a poor prognosis. Pulmonary endarterectomy (PEA), however, provides a potentially curative treatment for selected patients with CTEPH (5) with other options including pulmonary arterial hypertension (PAH) therapies and balloon pulmonary angioplasty (6). Increasingly patients with CTEPH are presenting with comorbidities and additional information that could aid decision making would be helpful.

Several clinical and haemodynamic measurements have been used to assess disease severity and risk stratify patients with PAH. These include an assessment of symptoms (WHO Functional class), exercise capacity (6-min walk test) and measures reflecting RV function including blood based biomarkers (N-terminal pro brain natriuretic peptide) and measures from cardiac catheterization (right atrial pressure, cardiac index and mixed venous oxygen saturation) (79). However, some of these measurements are limited by their subjectivity and invasive nature. Nonetheless, a multiparameter risk assessment incorporating a number of these measurements is now recommended in patients with PAH (9) and there is evidence that this approach can also be used in patients with CTEPH (10, 11).

Magnetic resonance imaging using CMR has been shown to have diagnostic value in suspected PH (1214), prognostic value in PAH (1517) and in screening for CTEPH using MR perfusion maps (18, 19) and aiding the surgical assessment of CTEPH using MR pulmonary angiography (12, 20). Septal angle, pulmonary artery area and ventricular mass index have additive value in a model to estimate pulmonary artery pressure (14, 21). Measures of RV function including stroke volume, right ventricular end systolic volume and right ventricular ejection fraction predict clinical worsening and mortality and left ventricular measures such as left ventricular end systolic volume predict mortality in PAH (1517, 22). However, there is limited data on the utility of CMR measures of right ventricular function to aid mortality prediction in CTEPH and to our knowledge no data on the use of left ventricular or atrial measurements.

The aim of this study was, therefore, to assess the prognostic value of gold standard imaging of cardiac structure and function using CMR in a large cohort of patients with CTEPH.

Materials and Methods

Patients

Consecutive treatment-naive patients diagnosed with CTEPH who had CMR and right heart catheterization (RHC) up to February 2017 were prospectively recorded in hospital databases as part of the ASPIRE (Assessing the Spectrum of Pulmonary Hypertension Identified at a Referral Center) registry as previously described (23, 24) (ClinicalTrials.gov identifier NCT02565030). Patients' demographics, imaging, and clinical metrics with follow-up data were prospectively collected using a census date of April 29th 2019. CMR and RHC metrics, exercise test, pulmonary function, and treatment were included.

The diagnosis of CTEPH required patients to have undergone RHC with a measured mean pulmonary artery pressure (mPAP) ≥ 25 mmHg at rest and at least one segmental perfusion defect on perfusion lung Q scan or pulmonary artery (PA) obstruction seen by multidetector computed tomography pulmonary angiography (CTPA) or conventional pulmonary angiography with other causes of PH excluded (9).

Following multi-modality imaging patients were discussed at a multi-professional meeting. Patients with CTEPH were subsequently discussed at the national pulmonary endarterectomy MDT at Papworth Hospital, Cambridge where surgical accessibility and suitability for surgery were assessed.

All patients were followed up until the date of death or census date. Ethical approval for this study was granted by our institutional review board (ref c06/Q2308/8). This study was funded by grants from the Wellcome Trust (A.J.S.). The funding body was not involved in the study design or data interpretation.

Cardiac MRI Acquisition

CMR was performed using an eight channel cardiac coil on a GE HDx (GE Healthcare, Milwaukee, WI) whole body scanner at 1.5 T, as previously described (15). Short-axis cine images were acquired using a cardiac gated multislice balanced SSFP sequence (20 frames per cardiac cycle; slice thickness, 8 mm; field of view, 48 cm; matrix, 2563256; BW, 125 kHz/pixel; TR/TE, 3.7/1.6 ms). A stack of images in the short-axis plane with slice thickness of 8 mm (2-mm inter-slice gap) were acquired fully covering both ventricles from base to apex. End-systole was considered to be the smallest cavity area. End-diastole was defined as the first cine phase of the R-wave triggered acquisition or largest volume. Through plane phase contrast imaging was performed orthogonal to the main pulmonary trunk. Phase contrast imaging parameters were as follows: repetition time, TR 5.6 ms; echo time, TE 2.7 ms; slice thickness, 10 mm; field of view, 48 cm, bandwidth, 62.5 kHz; matrix, 256 3128; 20 reconstructed cardiac phases; and velocity encoding of flow, 150 cm/s. Patients were in the supine position with a surface coil and with retrospective ECG gating.

Image Analysis

Image analysis was performed on a GE Advantage Workstation 4.1 with the observer blinded to the patient clinical information, and cardiac catheter parameters. Right and left endocardial and epicardial surfaces were manually traced from the stack of short-axis cine images, using proprietary MR workstation software to obtain RV end-diastolic volume (RVEDV) and RV end-systolic volume (RVESV), and left ventricular (LV) end-diastolic volume (LVEDV) and LV end-systolic volume (LVESV). From end-diastolic and end-systolic volumes, RVEF and LV ejection fraction (LVEF) and RV and LV stroke volumes (SV) were calculated. With the exception of RVEF and LVEF, these measurements were all corrected for age and sex (%pred) based on previously published referenced data (25) and they were indexed for body surface area based on Mostellar formula (26). Based on previous work, SV was considered to be the most accurate from LV volumetry (27) and was used for MRI estimation of RV–PA coupling. For calculation of ventricular mass, the interventricular septum was considered as part of the LV. RV end-diastolic mass (RVEDM) and LV end-diastolic mass (LVEDM) were derived. Ventricular mass index (VMI) was defined as RV mass divided by LV mass, as previously described (28). Maximal and minimal PA areas were measured, and relative area change was defined by the following equation: relative area change = (maximum area – minimum area)/minimum area (29). Reproducibility for CMR measurements in our center has been described previously (15).

Right Heart Catheterization

RHC was performed using a balloon-tipped 7.5F thermodilution catheter (Becton-Dickinson, Franklin Lakes, NJ). RHC was usually performed via the internal jugular vein using a Swan-Ganz catheter. Measurements of right atrial pressure, pulmonary arterial wedge pressure (PAWP), cardiac output (CO) using thermodilution and mixed venous oxygen saturation were also made. Pulmonary vascular resistance (PVR) was calculated as [(mPAP-PAWP) / CO] / 80 and expressed as dyne.s.cm−5.

Pulmonary Arterial Stiffness and Coupling Measurements

As previously described (15, 3032), RV elastance (Ees) was estimated as mPAP divided by RVESV. PA elastance (Ea) was estimated using mPAP-PAWP divided by LVSV. Ees/Ea by a combined RHC and CMR approach was defined as follows: (mPAP/RVESV)/[(mPAP-PAWP)/LVSV]. CMR estimated Ees/Ea was defined by LVSV/RVESV. Distensibility, a measurement of PA stiffness, was defined as PA relative area change divided by pulse pressure.

Statistical Analysis

Continuous variables were expressed as mean (SD) for parametric variables and median (interquartile range) for non-parametric variables. Categorical data was presented as the number of subjects and percentage. Continuous variables were compared using independent sample t-test. Categorical variables were compared using Pearson's Chi-square test. Survival analysis was conducted using Kaplan-Meier plots and survival was compared using the log-rank test. Survival was calculated from date of diagnosis to date of death or census date to compare survival between the PEA and non-PEA groups. To assess the prognostic value of CMR metrics, survival was calculated from date of CMR to date of death or census date. Univariate Cox proportional hazard regression analysis was used to assess the prognostic value of CMR metrics in terms of biventricular volume, function, mass and PA stiffness metrics in the whole cohort, PEA and non-PEA groups. Hazard ratios (HR) generated from univariate analysis were scaled by dividing the actual individual value of the variable by SD. This is to allow direct comparison of HR. Multivariate Cox proportional hazard regression analysis was performed for predicted CMR variables / MRI coupling measurement, clinical variables and combined CMR and clinical variables. Selected variables were entered into multivariate models if they were reported in literature as predictors of mortality (33, 34), ≤ 10% missing values and a p < 0.200 at univariate analysis. Multivariate models for the whole cohort, PEA and non-PEA groups were generated. To overcome multicollinearity, highly correlated variables (r > 0.80) were entered separately in the models. Receiver operating characteristics (ROC) analysis was used to assess the prognostic significance of the CMR, clinical and combined CMR and clinical models generated on the whole, PEA and non-PEA cohorts and presented with area under the curve (AUC). In order to evaluate the stability (internal validation) of the CMR prognostic model of the whole cohort, PEA group and non-PEA groups, the bootstrap approach with 1,000 bootstrap samples (default settings) was performed (35). In the bootstrap, simple sampling method with bias-corrected and accelerated 95% confidence interval (95% C.I.) type were selected. Locally estimated scatterplot smoothing regression analysis (LOESS) was performed for 1 year mortality, where significance was demonstrated at multivariate analysis, for patients undergoing and not undergoing PEA. All statistical tests were two-sided and a p-value of < 0.050 was considered statistically significant. A Statistical Package for the Social Sciences Program (SPSS) version 26 for Windows (SPSS Inc. Chicago, IL) was used for statistical analysis and for presentation of data GraphPad Prism 8.3.0 (GraphPad Software, San Diego, CA) was used.

Results

Baseline Demographics and Measurements

A total of 375 patients (mean+/-SD age 64+/- 14 years, 185 (49%) females) with CTEPH were included in the analysis. One hundred and eighty-one (48%) patients had PEA. Table 1 shows demographic, RHC, and lung function results and Table 2 CMR, RV-PA coupling metrics and pulmonary stiffness data, for the whole cohort and for the PEA and non-PEA groups. The mean+/- SD time delay between RHC and CMR was 45+/-15 days.

TABLE 1
www.frontiersin.org

Table 1. Patients demographics and results of baseline investigations for the whole CTEPH cohort, patients undergoing pulmonary endarterectomy, and not undergoing pulmonary endarterectomy.

TABLE 2
www.frontiersin.org

Table 2. Cardiac MRI imaging parameters for the whole CTEPH cohort, patients undergoing pulmonary endarterectomy, and not undergoing pulmonary endarterectomy.

Compared with the non-PEA group, patients in the PEA group were younger (p < 0.001), had a higher incremental shuttle walking distance (ISWD) (p = 0.001), lower RVESVI%pred (p=0.003), lower RVEDMI%pred (p = 0.046), higher forced expiratory volume at 1 min (FEV1) (p < 0.001) and higher forced vital capacity (FVC) (p = 0.003).

Survival Analysis

During the follow up period, 104 (28%) patients died. The median overall survival for the whole cohort was 143 months. Median survival of patients undergoing PEA was higher than patients not undergoing surgery (146 vs. 97 months, 95% C.I. (121–162); p < 0.001). Patients' demographics, CMR, RHC, lung function tests and pulmonary arterial stiffness data between survivors and non-survivors is summarized in Supplementary Table 1. Survivors were younger (p < 0.001) with a higher percentage of females (p < 0.001), had better lung function (p < 0.001), exercise capacity and less likely to have a history of malignancy, coronary artery disease, chronic kidney disease or COPD.

Univariate Analysis of Predictors of Mortality

Univariate Cox proportional hazard regression analysis for CMR metrics and clinical variables for the whole cohort is presented in Supplementary Table 2. Only variables with p < 0.200 are shown. CMR measures of RV size and function [RVEDVI%pred (p = 0.002), RVESVI%pred (p < 0.001), and RVEF%pred (p < 0.001)], LV function [LVSVI%pred (p = 0.001) and LVEF%pred (p = 0.011)], LA volume [LAVI (p = 0.010)] and invasive (p = 0.027) and non-invasive CMR-derived Ea/Ees ratio (p = 0.001), were significant predictors of mortality at univariate Cox regression analysis for the whole cohort. Separate univariate Cox regression analysis was performed for PEA and non-PEA groups and is shown in Tables 3, 4, respectively.

TABLE 3
www.frontiersin.org

Table 3. Univariate Cox proportional hazards regression analysis in patients with CTEPH undergoing pulmonary endarterectomy (metrics shown where p < 0.20).

TABLE 4
www.frontiersin.org

Table 4. Univariate Cox proportional hazards regression analysis in patients with CTEPH not undergoing pulmonary endarterectomy (metrics shown where p < 0.20).

Multivariate Analysis of Predictors of Mortality

Multivariate Cox proportional hazard regression analysis models of only CMR metrics, clinical variables and CMR metrics and clinical variables (p < 0.200) at univariate analysis are presented in Table 5. In the whole cohort, higher RVESVI%pred (p = 0.007), lower LVSVI%pred (p = 0.001) and higher LAVI (p = 0.001) were independent predictors of increased mortality in the CMR metrics model. The regression equation for the multivariate Cox regression model in the whole cohort is as follows: Expected hazard = (RVESV%pred*0.002)-(LVSVI%pred*0.019) + (LAVI*0.013). In the PEA group, lower LVSVI%pred (p = 0.035) and higher LAVI (p < 0.001) (Expected hazard = (LAVI*0.030) - (LVSVI%pred*0.018) and in the non-PEA group, lower RVEF%pred (p = 0.015) and lower LVSVI%pred (p = 0.024) were independent predictors of increased mortality [Expected hazard = (RVEF%pred*-0.167) - (LVSVI%pred*0.014)]. In the combined model, lower LVSVI%pred (p = 0.040) and higher LAVI (p = 0.033) in the whole cohort, higher LAVI (p = 0.009) in the PEA group and lower RVEF%pred (p = 0.040) in the non-PEA group were predictors of increased mortality.

TABLE 5
www.frontiersin.org

Table 5. Multivariate cardiac MR metrics, clinical variables, and combined cardiac MR and clinical variables Cox regression model in the whole CTEPH cohort, endarterectomy and non-endarterectomy groups.

Figure 1 shows Kaplan-Meier survival analysis of the PEA group (LVSVI%pred and LAVI) and non-PEA group (LVSVI%pred and RVEF%pred). LAVI is plotted above and below a value of 41 ml/m2 based on previous work and for other values based on median (36). Figure 2 shows LOESS regression analysis for risk of 1 year mortality for LVSVI%pred and LAVI for patients undergoing PEA and for LVSVI%pred and RVEF%pred for the Non-PEA group.

FIGURE 1
www.frontiersin.org

Figure 1. Kaplan-Meier survival analysis from date of cardiac MR showing the outcome of patients undergoing pulmonary endarterectomy (top) based on left ventricular stroke volume index %predicted (A) and left atrial volume index (B) and for patients not undergoing pulmonary endarterectomy (bottom) based on left ventricular stroke volume index %predicted (C) and right ventricular ejection fraction %predicted (D). Numbers at risk for each group are presented below the plot.

FIGURE 2
www.frontiersin.org

Figure 2. Locally estimated scatterplot smoothing regression analysis (LOESS) for risk of 1 year mortality for patients undergoing pulmonary endarterectomy (top) for left ventricular stroke volume index %predicted (A) and left atrial volume index (B) and patients not undergoing pulmonary endarterectomy for left ventricular stroke volume index %predicted (C) and for right ventricular ejection fraction %predicted (D). Black dots represent status (alive at the bottom and dead at the top). The blue line correspond to a smoothing function fitted by loess regression.

Accuracy, Sensitivity, and Specificity of Prognostic Models

The accuracy of the multivariate Cox models for the whole, PEA and non-PEA cohorts was tested using ROC analysis and presented as AUC. The prognostic accuracy of the CMR model in the whole cohort (AUC, 0.69, p = 0.001), clinical model (AUC, 0.62, p < 0.001) and in the combined CMR and clinical model (AUC, 0.65, p < 0.001). Sensitivity and specificity of the CMR whole cohort model was 71% (95% CI, 65–76) and 62% (95% CI, 52–72) respectively, clinical model 78% (95% CI, 69–86) and 43% (95% CI, 37–50) and for the combined clinical and CMR model 70% (95% CI, 63–76) and 55% (95% CI, 45–65), respectively (Figure 3A). In the PEA cohort, the prognostic accuracy of the CMR model (AUC, 0.73, p < 0.001), clinical model (AUC, 0.69, p = 0.002) and in the combined clinical and CMR model (AUC, 0.76, p < 0.001). Sensitivity and specificity of the CMR PEA model 80% (95% CI, 61–91) and 53% (95% CI, 46–62), respectively, clinical model 72% (95% CI, 65–79) and 65% (95% CI, 46–79) and for the combined model 80% (95% CI, 61–91) and 53% (95% CI, 46–62), respectively (Figure 3B). The prognostic accuracy of the CMR model in the non-PEA group (AUC, 0.63, p = 0.003), clinical model (AUC, 0.63, p = 0.003) and for the combined model (AUC, 0.56, p = 0.030). The sensitivity and specificity were 65% (95% CI, 56–73) and 60% (95% CI, 50–71), respectively for the CMR model, 60% (95% CI, 48–70) and 58% (95% CI, 50–67), respectively for the clinical model and 60% (95% CI, 48–70) and 58% (95% CI, 50–67), respectively for the combined model (Figure 3C).

FIGURE 3
www.frontiersin.org

Figure 3. Receiver operating curves for CMR model (green dashes), clinical model (blue dashes) and combined clinical and CMR models (red dashes) in the whole cohort (A), PEA group (B), and in the non-PEA (C) showing area under the curve (AUC) for overall mortality.

Internal Validation and Bootstrapping

Internal validation was performed using the bootstrap method as previously described. Supplementary Table 3 summarizes the results of the multivariate Cox regression model using the bootstrap method for the whole cohort CMR model. Variables remained significant with no bias after performing 1,000 bootstraps using bias corrected and accelerated confidence intervals. Similarly, variables for CMR prognostic model for the PEA and non-PEA groups remained significant after the bootstrap method (All p < 0.050).

Discussion

To our knowledge this is the first study to demonstrate that CMR has prognostic value in patients with CTEPH. We have also highlighted the importance of considering left heart disease in patients considered for PEA, by showing the prognostic value of LAVI in patients undergoing PEA, whilst confirming the value of established prognostic markers of RV function in patients with CTEPH who do not undergo PEA.

Despite extensive study of the potential value of CMR in patients with PAH, where measures of RV function are strongly prognostic (15, 17, 25, 3743), there is only limited data on the use of CMR in patients with CTEPH with large cohort studies in CTEPH primarily reporting on the utility of clinical, exercise and haemodynamic measures to predict mortality (24, 34, 44). Studies using CMR in CTEPH (45, 46) have focussed on changes in CMR metrics following PEA. Schoenfeld et al. (45) showed an improvement in RV mass and function 12 days post PEA and Mauritz et al. (46) showed an improvement in RV systolic wall stress 6 months post PEA. Non cardiac focused MRI based techniques such as MR pulmonary perfusion maps and MR pulmonary angiography for the screening of CTEPH and surgical planning have also shown to have diagnostic value (12, 1820, 47). Combining CMR with MR pulmonary perfusion maps and MR pulmonary angiography provides a potential one-stop shop for comprehensive cardiopulmonary of evaluation of suspected PH with MR imaging, providing information on the likelihood of PH, its severity and potential causes.

In this study we have shown that a number of measures reflecting RV function, including RVESVI%pred which is strongly prognostic in PAH (15, 25) and RVEF (%pred), were predictors of mortality in the whole cohort (Supplementary Table 1) and in those not undergoing PEA (Table 5), at univariate analysis. At multivariate analysis LVSVI %predicted, a measure of global cardiac function was an independent predictor of outcome in the whole cohort whilst RVEF%pred a measure of RV function was an independent predictor of outcome in patients not undergoing PEA. Increases in RV afterload in idiopathic PAH impair RV function resulting in reductions in LVSVI which has been shown to have prognostic value (39) and aid monitoring of treatment (40). In contrast, neither RV size nor RV functional metrics prior to PEA were predictors of mortality in patients undergoing surgery at univariate analysis. In carefully selected patients with surgical disease, PEA can dramatically reduce afterload, normalize pulmonary hemodynamics and improve RV function (45, 4850) and may explain in significant part why pre-operative measures of RV function are not prognostic in these patients.

A novel finding of this study relates to prognostic value of LA size in patients with CTEPH. LAVI was a significant independent predictor of mortality in the whole cohort and in patients undergoing PEA, highlighting the importance of considering left heart disease when considering the risk and benefits of PEA in CTEPH. Although a LAVI of >41 ml/m2 is used as a threshold to identify an enlarged left atrium on transthoracic echocardiography a normal LA size in patients undergoing PEA cannot be used to exclude left heart disease, and this may be unmasked following surgery. PEA results in increased LA and LV filling post-surgery (46) and increases in LAVI of ~20% post PEA have previously been noted and correlated with changes in PVR (51).

Markers of PA stiffness based on using a CMR only measurement method or a combination of CMR and RHC measurement methods were also predictors of mortality at univariate analysis but were not shown to be independent predictors of mortality in the whole, PEA or non-PEA cohorts. Vanderpool et al. (30) reported superior prognostic significance of CMR-derived estimate of RV-arterial coupling Ees/Ea over other invasive and non-invasive measures of RV function in patients referred with pulmonary hypertension. However, no additional prognostic value of a CMR-only measurement of RV–PA coupling over volumetric indices was demonstrated in the present study in patients with CTEPH.

This study also confirms findings from the International CTEPH Registry (52) and our previous work from the ASPIRE Registry regarding the prognostic impact of demographics such as age, comorbidities including malignancy, coronary artery disease and chronic kidney disease, pulmonary hemodynamics and measures of exercise capacity and lung function on outcomes in patients with CTEPH. However, despite inclusion of these prognostic factors MRI metrics still had independent predictive value (Table 5). The prognostic accuracy of the CMR only model was higher compared to the accuracy of the clinical only model and combining both models did not increase the prognostic ability of the model only for the PEA group (Figure 3).

Study Limitations

This was a single center study with retrospective analysis of prospectively collected data which might have introduced selection or misclassification bias. CMR metrics were corrected for age and sex but not for race, which impacts on CMR metrics (53). Despite identifying LAVI as an independent predictor of outcome in the whole cohort and patients undergoing PEA further work is required to identify thresholds to aid clinical decision making. This study also used CMR as a gold standard technique to assess cardiac volumetric and functional data rather than echocardiography and echocardiography may be able to provide comparative data although this is not answered by this study. Balloon pulmonary angioplasty has only recently been available in the UK and we were unable therefore to provide any data regarding the potential for CMR to assess for outcomes in patients undergoing this intervention.

Conclusion

CMR metrics reflecting RV function and the presence of left heart disease are of prognostic value in patients with CTEPH. In those undergoing PEA an elevated left atrial volume index predicts a worse outcome highlighting the importance of considering left heart disease in patients considered for PEA. This study also demonstrates the prognostic value of CMR imaging metrics in determining outcome in patients with CTEPH not-undergoing surgery. Whether CMR has a role in serial monitoring of such patients requires further evaluation.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by University of Sheffield, ref c06/Q2308/8. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

YS, AS, and DK conceived the idea for the study and prepared tables and figures. YS, SA, MS, DA, and RL collected and analyzed data for study. YS and AS did statistical analysis. SA, PG, CJ, and SR assisted with the demographic and MRI data. JW, RC, CE, AC, AH, AS, and DK supported this study management and assisted in the writing of the manuscript. AS and DK edited the final manuscript. All authors read and approved the manuscript.

Funding

This study was funded in part by the Wellcome Trust to AS, grant number: AJS 205188/Z/16/Z.

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.

Publisher's Note

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.

Supplementary Material

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

Glossary

CTEPH, chronic thromboembolic pulmonary hypertension; PAH, pulmonary arterial hypertension; mPAP, mean pulmonary artery pressure; RHC, right heart catheter; DLCO, diffusing capacity of the lungs for carbon monoxide; FEV1, forced expiratory volume; FVC, forced vital capacity; ISWD, incremental shuttle walking distance; mPAP, mean pulmonary arterial pressure; mRAP, mean right atrial pressure; PAWP, pulmonary artery wedge pressure; PEA, pulmonary endarterectomy; PVR, pulmonary vascular resistance; SvO2, mixed venous oxygen saturation; %pred, percentage predicted; Ea, arterial load; Ees, right ventricle elastance; LAVI, left atrium volume index; LVEDMI, left ventricle end diastolic mass index; LVEDVI, left ventricle end diastolic volume index; LVEF, left ventricle ejection fraction; LVESVI, left ventricle end systolic volume index; LVSVI, left ventricle stroke volume index; MRI, magnetic resonance imaging; PA, pulmonary artery; PA RAC, pulmonary artery relative area change; RVEDMI, right ventricle end diastolic mass index; RVEDVI, right ventricle end diastolic volume index; RVEF, right ventricle ejection fraction; RVSP, right ventricle systolic pressure; RVESVI, right ventricle end systolic volume index; RVSVI, right ventricle stroke volume index; VMI, ventricular mass index.

References

1. Kiely DG, Elliot CA, Sabroe I, Condliffe R. Pulmonary hypertension: diagnosis and management. BMJ. (2013) 346:f2028. doi: 10.1136/bmj.f2028

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Ende-Verhaar YM, Cannegieter SC, Vonk Noordegraaf A, Delcroix M, Pruszczyk P, Mairuhu AT, et al. Incidence of chronic thromboembolic pulmonary hypertension after acute pulmonary embolism: a contemporary view of the published literature. Eur Respir J. (2017) 49:1601792. doi: 10.1183/13993003.01792-2016

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Hoeper MM, Mayer E, Simonneau G, Rubin LJ. Chronic thromboembolic pulmonary hypertension. Circulation. (2006) 113:2011–20. doi: 10.1161/CIRCULATIONAHA.105.602565

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Moser KM, Braunwald NS. Successful surgical intervention in severe chronic thromboembolic pulmonary hypertension. Chest. (1973) 64:29–35. doi: 10.1378/chest.64.1.29

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Bonderman D, Skoro-Sajer N, Jakowitsch J, Adlbrecht C, Dunkler D, Taghavi S, et al. Predictors of outcome in chronic thromboembolic pulmonary hypertension. Circulation. (2007) 115:2153–8. doi: 10.1161/CIRCULATIONAHA.106.661041

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Delcroix M, Torbicki A, Gopalan D, Sitbon O, Klok FA, Lang I, et al. ERS statement on chronic thromboembolic pulmonary hypertension. Eur Respir J. (2021) 57:2002828. doi: 10.1183/13993003.02828-2020

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Benza RL, Miller DP, Gomberg-Maitland M, Frantz RP, Foreman AJ, Coffey CS, et al. Predicting survival in pulmonary arterial hypertension: insights from the Registry to Evaluate Early and Long-Term Pulmonary Arterial Hypertension Disease Management (REVEAL). Circulation. (2010) 122:164–72. doi: 10.1161/CIRCULATIONAHA.109.898122

PubMed Abstract | CrossRef Full Text | Google Scholar

8. D'Alonzo GE, Barst RJ, Ayres SM, Bergofsky EH, Brundage BH, Detre KM, et al. Survival in patients with primary pulmonary hypertension. Results from a national prospective registry. Ann Intern Med. (1991) 115:343–9. doi: 10.7326/0003-4819-115-5-343

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Galie N, Humbert M, Vachiery JL, Gibbs S, Lang I, Torbicki A, et al. 2015 ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension: the Joint Task Force for the Diagnosis and Treatment of Pulmonary Hypertension of the European Society of Cardiology (ESC) and the European Respiratory Society (ERS): endorsed by: association for European Paediatric and Congenital Cardiology (AEPC), International Society for Heart and Lung Transplantation (ISHLT). Eur Respir J. (2015) 46:903–75. doi: 10.1183/13993003.01032-2015

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Benza RL, Farber HW, Frost A, Grünig E, Hoeper MM, Busse D, et al. REVEAL risk score in patients with chronic thromboembolic pulmonary hypertension receiving riociguat. J Heart Lung Transplant. (2018) 37:836–43. doi: 10.1016/j.healun.2018.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Sandqvist A, Kylhammar D, Bartfay SE, Hesselstrand R, Hjalmarsson C, Kavianipour M, et al. Risk stratification in chronic thromboembolic pulmonary hypertension predicts survival. Scand Cardiovasc J. (2021) 55:43–9. doi: 10.1080/14017431.2020.1783456

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Swift AJ, Wild JM, Nagle SK, Roldán-Alzate A, François CJ, Fain S, et al. Quantitative magnetic resonance imaging of pulmonary hypertension: a practical approach to the current state of the art. J Thorac Imaging. (2014) 29:68–79. doi: 10.1097/RTI.0000000000000079

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Johns CS, Kiely DG, Rajaram S, Hill C, Thomas S, Karunasaagarar K, et al. Diagnosis of pulmonary hypertension with cardiac MRI: derivation and validation of regression models. Radiology. (2019) 290:61–8. doi: 10.1148/radiol.2018180603

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Whitfield AJ, Solanki R, Johns CS, Kiely D, Wild J, Swift AJ, et al. Prediction of precapillary pulmonary hypertension according to the sixth world symposium on pulmonary hypertension. Radiology. (2020) 294:482. doi: 10.1148/radiol.2019192078

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Swift AJ, Capener D, Johns C, Hamilton N, Rothman A, Elliot C, et al. Magnetic resonance imaging in the prognostic evaluation of patients with pulmonary arterial hypertension. Am J Respir Crit Care Med. (2017) 196:228–39. doi: 10.1164/rccm.201611-2365OC

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Kiely DG, Levin D, Hassoun P, Ivy DD, Jone PN, Bwika J, et al. EXPRESS: statement on imaging and pulmonary hypertension from the Pulmonary Vascular Research Institute (PVRI). Pulm Circ. (2019) 2019:2045894019841990. doi: 10.1177/2045894019841990

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Alabed S, Shahin Y, Garg P, Alandejani F, Johns CS, Lewis RA, et al. Cardiac-MRI predicts clinical worsening and mortality in pulmonary arterial hypertension: a systematic review and meta-analysis. JACC Cardiovasc Imaging. (2021) 14:931–42. doi: 10.1016/j.jcmg.2020.08.013

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Johns CS, Swift AJ, Rajaram S, Hughes PJC, Capener DJ, Kiely DG, et al. Lung perfusion: MRI vs. SPECT for screening in suspected chronic thromboembolic pulmonary hypertension. J Magn Reson Imaging. (2017) 46:1693–7. doi: 10.1002/jmri.25714

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Rajaram S, Swift AJ, Telfer A, Hurdman J, Marshall H, Lorenz E, et al. 3D contrast-enhanced lung perfusion MRI is an effective screening tool for chronic thromboembolic pulmonary hypertension: results from the ASPIRE Registry. Thorax. (2013) 68:677–8. doi: 10.1136/thoraxjnl-2012-203020

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Rajaram S, Swift AJ, Capener D, Telfer A, Davies C, Hill C, et al. Diagnostic accuracy of contrast-enhanced MR angiography and unenhanced proton MR imaging compared with CT pulmonary angiography in chronic thromboembolic pulmonary hypertension. Eur Radiol. (2012) 22:310–7. doi: 10.1007/s00330-011-2252-x

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Johns CS, Wild JM, Rajaram S, Tubman E, Capener D, Elliot C, et al. Identifying at-risk patients with combined pre- and postcapillary pulmonary hypertension using interventricular septal angle at cardiac MRI. Radiology. (2018) 289:61–8. doi: 10.1148/radiol.2018180120

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Lewis RA, Johns CS, Cogliano M, Capener D, Tubman E, Elliot CA, et al. Identification of cardiac magnetic resonance imaging thresholds for risk stratification in pulmonary arterial hypertension. Am J Respir Crit Care Med. (2020) 201:458–68. doi: 10.1164/rccm.201909-1771OC

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Hurdman J, Condliffe R, Elliot CA, Davies C, Hill C, Wild JM, et al. ASPIRE registry: assessing the Spectrum of Pulmonary hypertension Identified at a REferral centre. Eur Respir J. (2012) 39:945–55. doi: 10.1183/09031936.00078411

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Quadery SR, Swift AJ, Billings CG, Thompson AAR, Elliot CA, Hurdman J, et al. The impact of patient choice on survival in chronic thromboembolic pulmonary hypertension. Eur Respir J. (2018) 52:589. doi: 10.1183/13993003.00589-2018

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Swift AJ, Rajaram S, Campbell MJ, Hurdman J, Thomas S, Capener D, et al. Prognostic value of cardiovascular magnetic resonance imaging measurements corrected for age and sex in idiopathic pulmonary arterial hypertension. Circ Cardiovasc Imaging. (2014) 7:100–6. doi: 10.1161/CIRCIMAGING.113.000338

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Mosteller RD. Simplified calculation of body-surface area. N Engl J Med. (1987) 317:1098. doi: 10.1056/NEJM198710223171717

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Mauritz GJ, Marcus JT, Boonstra A, Postmus PE, Westerhof N, Vonk-Noordegraaf A. Non-invasive stroke volume assessment in patients with pulmonary arterial hypertension: left-sided data mandatory. J Cardiovasc Magn Reson. (2008) 10:51. doi: 10.1186/1532-429X-10-51

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Saba TS, Foster J, Cockburn M, Cowan M, Peacock AJ. Ventricular mass index using magnetic resonance imaging accurately estimates pulmonary artery pressure. Eur Respir J. (2002) 20:1519–24. doi: 10.1183/09031936.02.00014602

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Sanz J, Kariisa M, Dellegrottaglie S, Prat-González S, Garcia MJ, Fuster V, et al. Evaluation of pulmonary artery stiffness in pulmonary hypertension with cardiac magnetic resonance. JACC Cardiovasc Imaging. (2009) 2:286–95. doi: 10.1016/j.jcmg.2008.08.007

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Vanderpool RR, Pinsky MR, Naeije R, Deible C, Kosaraju V, Bunner C, et al. RV-pulmonary arterial coupling predicts outcome in patients referred for pulmonary hypertension. Heart. (2015) 101:37–43. doi: 10.1136/heartjnl-2014-306142

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Sanz J, Garcia-Alvarez A, Fernandez-Friera L, Nair A, Mirelis JG, Sawit ST, et al. Right ventriculo-arterial coupling in pulmonary hypertension: a magnetic resonance study. Heart. (2012) 98:238–43. doi: 10.1136/heartjnl-2011-300462

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Najjar SS, Schulman SP, Gerstenblith G, Fleg JL, Kass DA, O'Connor F, et al. Age and gender affect ventricular-vascular coupling during aerobic exercise. J Am Coll Cardiol. (2004) 44:611–7. doi: 10.1016/j.jacc.2004.04.041

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Escribano-Subías P, Del Pozo R, Román-Broto A, Domingo Morera JA, Lara-Padrón A, Elías Hernández T, et al. Management and outcomes in chronic thromboembolic pulmonary hypertension: From expert centers to a nationwide perspective. Int J Cardiol. (2016) 203:938–44. doi: 10.1016/j.ijcard.2015.11.039

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Delcroix M, Lang I, Pepke-Zaba J, Jansa P, D'Armini AM, Snijder R, et al. Long-term outcome of patients with chronic thromboembolic pulmonary hypertension: results from an International Prospective Registry. Circulation. (2016) 133:859–71. doi: 10.1161/CIRCULATIONAHA.115.016522

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Altman DG, Andersen PK. Bootstrap investigation of the stability of a cox regression model. Stat Med. (1989) 8:771–83. doi: 10.1002/sim.4780080702

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Aune E, Baekkevar M, Roislien J, Rodevand O, Otterstad JE. Normal reference ranges for left and right atrial volume indexes and ejection fractions obtained with real-time three-dimensional echocardiography. Eur J Echocardiogr. (2009) 10:738–44. doi: 10.1093/ejechocard/jep054

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Swift AJ, Telfer A, Rajaram S, Condliffe R, Marshall H, Capener D, et al. Dynamic contrast enhanced magnetic resonance imaging in patients with pulmonary arterial hypertension. Pulm Circ. (2014) 4:61–70. doi: 10.1086/674882

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Gan CTJ, Lankhaar JW, Westerhof N, Marcus JT, Becker A, Twisk JW, et al. Noninvasively assessed pulmonary artery stiffness predicts mortality in pulmonary arterial hypertension. Chest. (2007) 132:1906–12. doi: 10.1378/chest.07-1246

PubMed Abstract | CrossRef Full Text | Google Scholar

39. van Wolferen SA, Marcus JT, Boonstra A, Marques KM, Bronzwaer JG, Spreeuwenberg MD, et al. Prognostic value of right ventricular mass, volume, and function in idiopathic pulmonary arterial hypertension. Eur Heart J. (2007) 28:1250–7. doi: 10.1093/eurheartj/ehl477

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Van de Veerdonk MC, Kind T, Marcus JT, Mauritz GJ, Heymans MW, Bogaard HJ, et al. Progressive right ventricular dysfunction in patients with pulmonary arterial hypertension responding to therapy. J Am Coll Cardiol. (2011) 58:2511–9. doi: 10.1016/j.jacc.2011.06.068

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Freed BH, Gomberg-Maitland M, Chandra S, Mor-Avi V, Rich S, Archer SL, et al. Late gadolinium enhancement cardiovascular magnetic resonance predicts clinical worsening in patients with pulmonary hypertension. J Cardiovasc Magn Reson. (2012) 14:11. doi: 10.1186/1532-429X-14-11

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Yamada Y, Okuda S, Kataoka M, Tanimoto A, Tamura Y, Abe T, et al. Prognostic value of cardiac magnetic resonance imaging for idiopathic pulmonary arterial hypertension before initiating intravenous prostacyclin therapy. Circ J. (2012) 76:1737–43. doi: 10.1253/circj.CJ-11-1237

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Cho IJ, Oh J, Chang HJ, Park J, Kang KW, Kim YJ, et al. Tricuspid regurgitation duration correlates with cardiovascular magnetic resonance-derived right ventricular ejection fraction and predict prognosis in patients with pulmonary arterial hypertension. Eur Heart J Cardiovasc Imaging. (2014) 15:18–23. doi: 10.1093/ehjci/jet094

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Riedel M, Stanek V, Widimsky J, Prerovsky I. Longterm follow-up of patients with pulmonary thromboembolism. Late prognosis and evolution of hemodynamic and respiratory data. Chest. (1982) 81:151–8. doi: 10.1378/chest.81.2.151

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Schoenfeld C, Cebotari S, Hinrichs J, Renne J, Kaireit T, Olsson KM, et al. MR imaging-derived regional pulmonary parenchymal perfusion and cardiac function for monitoring patients with chronic thromboembolic pulmonary hypertension before and after pulmonary endarterectomy. Radiology. (2016) 279:925–34. doi: 10.1148/radiol.2015150765

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Mauritz GJ, Vonk-Noordegraaf A, Kind T, Surie S, Kloek JJ, Bresser P, et al. Pulmonary endarterectomy normalizes interventricular dyssynchrony and right ventricular systolic wall stress. J Cardiovasc Magn Reson. (2012) 14:5. doi: 10.1186/1532-429X-14-5

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Tsuchiya N, van Beek EJ, Ohno Y, Hatabu H, Kauczor HU, Swift A, et al. Magnetic resonance angiography for the primary diagnosis of pulmonary embolism: a review from the international workshop for pulmonary functional imaging. World J Radiol. (2018) 10:52–64. doi: 10.4329/wjr.v10.i6.52

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Berman M, Gopalan D, Sharples L, Screaton N, Maccan C, Sheares K, et al. Right ventricular reverse remodeling after pulmonary endarterectomy: magnetic resonance imaging and clinical and right heart catheterization assessment. Pulm Circ. (2014) 4:36–44. doi: 10.1086/674884

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Kreitner KF, Ley S, Kauczor HU, Mayer E, Kramm T, Pitton MB, et al. Chronic thromboembolic pulmonary hypertension: pre- and postoperative assessment with breath-hold MR imaging techniques. Radiology. (2004) 232:535–43. doi: 10.1148/radiol.2322030945

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Rolf A, Rixe J, Kim WK, Börgel J, Möllmann H, Nef HM, et al. Right ventricular adaptation to pulmonary pressure load in patients with chronic thromboembolic pulmonary hypertension before and after successful pulmonary endarterectomy—a cardiovascular magnetic resonance study. J Cardiovasc Magn Reson. (2014) 16:96. doi: 10.1186/s12968-014-0096-7

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Marston NA, Auger WR, Madani MM, Kimura BJ, Strachan GM, Raisinghani AB, et al. Assessment of left atrial volume before and after pulmonary thromboendarterectomy in chronic thromboembolic pulmonary hypertension. Cardiovasc Ultrasound. (2014) 12:32. doi: 10.1186/1476-7120-12-32

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Pepke-Zaba J, Delcroix M, Lang I, Mayer E, Jansa P, Ambroz D, et al. Chronic thromboembolic pulmonary hypertension (CTEPH): results from an international prospective registry. Circulation. (2011) 124:1973–81. doi: 10.1161/CIRCULATIONAHA.110.015008

PubMed Abstract | CrossRef Full Text | Google Scholar

53. 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

Keywords: cardiac MRI, chronic thromboembolic pulmonary hypertension, left atrium, survival, left heart disease

Citation: Shahin Y, Alabed S, Rehan Quadery S, Lewis RA, Johns C, Alkhanfar D, Sukhanenko M, Alandejani F, Garg P, Elliot CA, Hameed A, Charalampopoulos A, Wild JM, Condliffe R, Swift AJ and Kiely DG (2022) CMR Measures of Left Atrial Volume Index and Right Ventricular Function Have Prognostic Value in Chronic Thromboembolic Pulmonary Hypertension. Front. Med. 9:840196. doi: 10.3389/fmed.2022.840196

Received: 20 December 2021; Accepted: 04 February 2022;
Published: 14 March 2022.

Edited by:

Elena Goncharova, University of California, Davis, United States

Reviewed by:

Roberto J. Bernardo, University of Oklahoma Health Sciences Center, United States
Mark O. Wielpütz, Heidelberg University, Germany

Copyright © 2022 Shahin, Alabed, Rehan Quadery, Lewis, Johns, Alkhanfar, Sukhanenko, Alandejani, Garg, Elliot, Hameed, Charalampopoulos, Wild, Condliffe, Swift and Kiely. 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: Yousef Shahin, y.shahin@sheffield.ac.uk

Download