Article Text

Original research
Cardiac MR modelling of systolic and diastolic blood pressure
  1. Hosamadin Assadi1,2,
  2. Gareth Matthews1,2,
  3. Xiaodan Zhao3,
  4. Rui Li1,2,
  5. Samer Alabed4,
  6. Ciaran Grafton-Clarke1,2,
  7. Zia Mehmood2,
  8. Bahman Kasmai1,2,
  9. Vaishali Limbachia2,
  10. Rebecca Gosling4,
  11. Gurung-Koney Yashoda2,
  12. Ian Halliday4,
  13. Peter Swoboda5,
  14. David Paul Ripley6,
  15. Liang Zhong3,7,
  16. Vassilios S Vassiliou1,2,
  17. Andrew J Swift4,
  18. Rob J van der Geest8 and
  19. Pankaj Garg1,2
  1. 1Department of Cardiovascular and Metabolic Health, University of East Anglia, Norwich, UK
  2. 2Norfolk and Norwich University Hospitals NHS Foundation Trust, Norwich, UK
  3. 3National Heart Research Institute, National Heart Centre, Singapore
  4. 4Department of Infection, Immunity & Cardiovascular Disease, University of Sheffield, Sheffield, UK
  5. 5LICAMM, University of Leeds, Leeds, UK
  6. 6Department of Cardiology, Northumbria Specialist Emergency Care Hospital, Cramlington, UK
  7. 7Cardiovascular Science Academic Program, Duke-NUS Medical School, Singapore
  8. 8Department of Radiology, Leiden University Medical Center, Leiden, The Netherlands
  1. Correspondence to Dr Hosamadin Assadi; hosamassadi{at}gmail.com

Abstract

Aims Blood pressure (BP) is a crucial factor in cardiovascular health and can affect cardiac imaging assessments. However, standard outpatient cardiovascular MR (CMR) imaging procedures do not typically include BP measurements prior to image acquisition. This study proposes that brachial systolic BP (SBP) and diastolic BP (DBP) can be modelled using patient characteristics and CMR data.

Methods In this multicentre study, 57 patients from the PREFER-CMR registry and 163 patients from other registries were used as the derivation cohort. All subjects had their brachial SBP and DBP measured using a sphygmomanometer. Multivariate linear regression analysis was applied to predict brachial BP. The model was subsequently validated in a cohort of 169 healthy individuals.

Results Age and left ventricular ejection fraction were associated with SBP. Aortic forward flow, body surface area and left ventricular mass index were associated with DBP. When applied to the validation cohort, the correlation coefficient between CMR-derived SBP and brachial SBP was (r=0.16, 95% CI 0.011 to 0.305, p=0.03), and CMR-derived DBP and brachial DBP was (r=0.27, 95% CI 0.122 to 0.403, p=0.0004). The area under the curve (AUC) for CMR-derived SBP to predict SBP>120 mmHg was 0.59, p=0.038. Moreover, CMR-derived DBP to predict DBP>80 mmHg had an AUC of 0.64, p=0.002.

Conclusion CMR-derived SBP and DBP models can estimate brachial SBP and DBP. Such models may allow efficient prospective collection, as well as retrospective estimation of BP, which should be incorporated into assessments due to its critical effect on load-dependent parameters.

  • Magnetic Resonance Imaging
  • Hypertension
  • Aortic Diseases
  • Cardiac Imaging Techniques

Data availability statement

Data are available on reasonable request. The datasets generated and analysed during the current study are not publicly available. Access to the raw images of patients is not permitted since specialised postprocessing imaging-based solutions can identify the study patients in the future. Data are available from the corresponding author on reasonable request.

https://creativecommons.org/licenses/by/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

WHAT IS ALREADY KNOWN ON THIS TOPIC

  • Prior research has established the significance of blood pressure (BP) in cardiovascular health and its potential influence on cardiac imaging assessments. It is recognised that accurate BP measurements are essential for a comprehensive understanding of cardiovascular function. However, conventional outpatient cardiovascular MR (CMR) imaging procedures typically do not involve BP measurements before image acquisition. This gap in the field highlights the need for a method to estimate BP using CMR data and patient characteristics, which is the focus of this study.

WHAT THIS STUDY ADDS

  • This study adds a novel approach to the field by proposing a model to estimate brachial systolic and diastolic BP using patient characteristics and CMR imaging data. By demonstrating associations between specific CMR parameters and BP, the study provides a potential means to collect BP data prospectively during CMR imaging efficiently and retrospectively estimate BP values. This approach fills a crucial gap in standard CMR procedures, allowing for a more comprehensive assessment of cardiovascular health by incorporating BP, which is known to impact load-dependent parameters significantly.

HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

  • This study’s novel approach to estimating could revolutionise cardiovascular research by streamlining data collection, enhancing the precision of cardiac imaging assessments and facilitating retrospective BP estimation. It may guide treatment decisions and risk assessment in clinical practice, potentially influencing healthcare policies to incorporate BP measurements into CMR protocols. Ultimately, it can potentially elevate the standard of cardiovascular care by integrating BP measurements into CMR protocols.

Introduction

Blood pressure (BP) plays a vital role in haemodynamic assessment. Hypertension is the most common preventable risk factor for cardiovascular disease and is a major contributor to all-cause mortality, given its pathogenic role in stroke and chronic kidney disease. Worldwide, an estimated 31.1% of adults have hypertension, with an increased prevalence in low-income and middle-income countries.1 Physiologically, BP is dependent on numerous factors, including heart rate, stroke volume, systemic vascular resistance, blood volume, arterial compliance and neuroendocrine axes.2 In most cases, hypertension is straightforward to address with lifestyle and pharmacological strategies. Hypotension can also cause morbidity and mortality, particularly in older adults and those with heart failure.3 4

From a cardiac imaging perspective, BP influences the evaluation of ventricular size, chamber function and severity of valve pathology during functional assessment.5–7 Hypertension is commonly present in heart failure,8 especially in the preserved ejection fraction phenotype, left-sided valvular disease,9 ischaemic heart disease and atrial fibrillation. Numerous imaging guidelines incorporate the assessment of loading conditions to accurately grade specific lesions, with changes in loading conditions recognised as potential causes for misclassification of disease severity.10 11 BP is dynamic and impacted by body positioning and activity. Standard outpatient guidelines recommend BP assessment after a patient is seated and relaxed with feet on the floor for more than 5 min. Furthermore, to minimise random error, 2–3 BP measurements should be obtained on 2–3 different occasions.12

In standard outpatient cardiovascular MR (CMR) imaging procedures, BP measurement is usually not taken prior to image acquisition. As a result, during routine CMR examinations, it is challenging to make a comprehensive afterload-dependent assessment of load-dependent imaging. Moreover, CMR is emerging as a prognostically relevant test for aortic and mitral valve diseases, particularly mitral regurgitation and aortic regurgitation—both dynamic valvular incompetencies with dependence on BP. The time efficiency of CMR protocols is improving; however, brachial BP recording would increase radiographers’ workload and workflow durations. Additionally, for retrospective analysis, there is presently no method to infer BP and, therefore, loading conditions in the absence of contemporaneous measurement.

We hypothesise that both systolic BP (SBP) and diastolic BP (DBP) can be modelled using both patient characteristics and CMR-derived parameters, which are already collected as part of the standard assessment.

Methods

Study cohort

This study included patients from several databases referred to our centre for further assessment of breathlessness. The derivation cohort included 57 patients from the PREFER-CMR registry (ClinicalTrials.gov: NCT05114785) in Norwich, UK and 163 from research registries in Sheffield, UK.13 These registries were developed to broadly characterise and subphenotype healthy populations and patients with heart disease by CMR. The inclusion criteria were individuals over the age of 18 undergoing CMR evaluation who provided written informed consent and had a brachial BP measured using a sphygmomanometer before the scan.

To externally validate our findings, we used a cohort of 169 healthy individuals from the INITIATE registry in Singapore (ClinicalTrials.gov: NCT03217240) that recruited both patients and healthy subjects without known cardiovascular disease or cardiovascular risk factors such as hypertension, diabetes and hyperlipidaemia. A flow chart illustrating the recruitment process and the steps taken for deriving predictive models for SBP and DBP is shown in figure 1.

Figure 1

Flow chart illustrating the process of deriving predictive models for systolic and diastolic blood pressures using CMR data. CMR, cardiovascular MR.

Cardiac MR protocol

CMR study was performed on a 3.0 Tesla Ingenia (Philips Healthcare, Best, the Netherlands) and 1.5 Tesla Magnetom Sola Siemens system with a superconducting magnet (Siemens Healthineers AG, Erlangen, Germany). All patients were examined in the supine position, head first, using a respiratory sensor and ECG gating. Additionally, the scanner was equipped with a biometric body with 18 coils.

The CMR protocol included baseline survey images and cines and gadolinium enhancement imaging acquisition methods previously described by our group.14–20 For standard cines, we acquired 30 phases throughout the cardiac cycle. Other cine acquisition parameters include TR: 2.71, TE: 1.13, field of view (FOV): 360 × 289.3 mm2 with phase FOV—80.4%, number of signal averages: 1, matrix: 224×180 (phase), bandwidth: 167.4 kHz, (930 Hz/Px), flip angle: 80, slice thickness: 8 mm and Grappa acceleration with a factor of 2. Additionally, two-dimensional phase contrast CMR data were acquired in the ascending aorta at an orthogonal plane just above the sinotubular junction.

Cardiac MR analysis

All image analyses were postprocessed with the in-house developed MASS research software (MASS; Version 2022-EXP, Leiden University Medical Center, Leiden, The Netherlands). For left ventricle (LV) and right ventricle (RV) volume analysis, we used the integrated AI-based segmentation tool for endocardial and epicardial borders using the short-axis cine stack of images.21 The trabeculae and papillary muscles were included in the analysis. The end-diastolic volum and end-systolic volume were defined as the maximum and minimum values on the volume curves, respectively. LV mass was recorded at the end-diastolic phase.

For flow analysis, we used semiautomated segmentation methods of the ascending and descending aorta throughout the cardiac cycle previously described by our group.20 Aortic forward and backward flows were acquired from the resultant flow curve (figure 2). The peak systolic phase was defined as the peak flow rate on the flow curve (figure 2). The calculation of flow displacement involved measuring the distance between the centre point of the vessel and the centre of the velocity of the forward flow. This distance was then adjusted to account for the overall size of the artery during each phase of the cardiac cycle.22 23 To calculate the rotational angle (RA), the angle change between the end-systolic point and the point where the flow angle became stable after reaching its peak during systole was measured on the RA curve. Rotational speed (RS) was calculated by dividing the total change in angle by the time interval between consecutive phases. Specifically, the average value of RS was determined from the time after peak systole until the end of systole. The maximal and minimal areas of the ascending aorta were defined by measuring its cross-sectional area throughout the cardiac cycle. The relative area change (RAC) was then calculated as the percentage change between the maximal and minimal areas, using the formula: (maximal area−minimal area)/minimal area×100%.

Figure 2

Central illustration demonstrating the demographics and multiparametric CMR volumetric and flow parameters associated with brachial systolic and diastolic blood pressure. CMR, cardiovascular MR.

Statistical analysis

Data analyses were performed using MedCalc (MedCalc Software, Ostend, Belgium, V.20.215). Normality was assessed using a visual assessment of histograms and the Shapiro-Wilk test. Continuous parametric variables were expressed as mean±SD. Variables that displayed significant correlations (p<0.05) were selected as independent variables for multivariate linear regression analysis, from which a regression model was developed. For comparing CMR-modelled with brachial-measured SBP and DBP, we used Pearson’s correlation coefficient to assess the correlation and Bland-Altman plots to assess agreement and bias. The discriminative capability was evaluated by calculating the area under the curve (AUC) using receiver operating characteristic curves. The significance threshold was set at p<0.05.

Results

The overall results of the study are summarised in figure 2.

Study population

Table 1 shows the baseline characteristics of each cohort. Some notable baseline demographic differences were observed, with the validation cohort being statistically younger with a lower body surface area (BSA), a higher proportion of females and a lower mean SBP. The mean DBP was similar between cohorts. Left ventricular CMR parameters were not significantly different, aside from a higher LV Mass Index (LVMI) in the derivation cohort and a slightly higher LV stroke volume index in the validation cohort. Aortic flow parameters were more divergent between the cohorts, with only the AO peak systolic flow and systolic retrograde flow being non-significant between the groups. The validation cohort generally had a greater forward and smaller backward aortic flow with less flow rotation. The derivation cohort had a lower RAC, likely a combination of reduced stroke volume and higher vascular stiffness. The observed differences reflect the focus of the original studies, with the derivation cohort predominantly focused on heart failure with preserved ejection fraction and the validation cohort being healthy, overall resulting in heterogeneity for subsequent modelling.

Table 1

Study demographics, CMR functional and aortic flow parameters of the derivation and validation cohorts and their significance

Derivation cohort

Parameters which were significantly associated with SBP and DBP are shown in table 2. CMR functional and aortic flow parameters that displayed significant correlations with SBP and DBP were selected as independent variables for multivariate linear regression analysis. For SBP, only age (r=0.27, p<0.001) and LV ejection fraction (LVEF) (r=0.24, p<0.001) were found to be significantly correlated (figure 3) (online supplemental table S1). For DBP, only BSA (r=0.16, p=0.01), LVMI (r=0.14, p=0.04) and AO forward flow (r=−0.15, p=0.02) were found to be significantly correlated (figure 4) (online supplemental table S2).

Figure 3

Scatter diagrams with heat maps demonstrating correlations between systolic blood pressure with age (A) and left ventricular (LV) ejection fraction (B).

Figure 4

Scatter diagrams with heat maps demonstrating correlations between diastolic blood pressure with aortic forward flow (A), body surface area (B) and Left Ventricular Mass Index (C).

Table 2

Correlation coefficient of the study demographics, CMR functional and aortic flow variables with brachial measured systolic and diastolic blood pressures in the derivation cohort

A multivariate linear regression model was constructed to predict SBP and DBP. The derived equations were as follows:

CMR-modelled SBP=83.845+(0.4225×age)+(0.4187×LVEF).

CMR-modelled DBP=58.8591+(−0.1229×AO forward flow)+(8.2279×BSA)+(0.1738×LVMI).

Validation cohort

We applied the above equations to the validation cohort (n=169) to predict brachial SBP and DBP. The correlation coefficient of the CMR-derived SBP composite model and brachial SBP was (r=0.16, 95% CI 0.011 to 0.305, p=0.03). The diagnostic accuracy (figure 5A) of CMR-derived SBP to predict SBP>120 mmHg had a sensitivity of 54% and a specificity of 71% with an AUC of 0.59 (95% CI 0.74 to 0.89, p=0.038). In Bland-Altman’s analysis, the mean difference between CMR-derived SBP and brachial-measured SBP was −3 mmHg, with limits of agreement (LOA) of −35.8 to 29.7 mmHg (figure 6A).

Figure 5

(A) Receiver operator curves demonstrate satisfactory agreement between CMR-modelled and brachial-measured systolic blood pressures (BP) in estimating pressures >120 mmHg. (B) Receiver operator curves show satisfactory agreement between CMR-modelled and brachial-measured diastolic BPs in estimating pressures >80 mmHg. AUC, area under the curve; CMR, cardiovascular MR.

Figure 6

Bland-Altman plots demonstrating the degree of agreement between brachial measured and CMR-modelled systolic blood pressures (BPs) (A) and between brachial measured and CMR-modelled diastolic BP (B). CMR, cardiovascular MR.

Similarly, with DBP, the correlation coefficient of the CMR-derived DBP composite model and brachial DBP were (r=0.27, 95% CI 0.122 to 0.403, p=0.0004). The diagnostic accuracy (figure 5B) of CMR-derived DBP to predict brachial DBP>80 mmHg had a sensitivity of 47% and a specificity of 81% with an AUC of 0.64 (95% CI 0.56 to 0.71, p=0.002). Moreover, the mean bias between CMR-derived DBP and brachial-measured DBP was 4 mmHg, with LOA of −18.2 to 26.2 mmHg (figure 6B).

Discussion

This study aimed to develop and test a physiological model to estimate SBP and DBPs using patient characteristics and CMR functional and flow data. Our findings suggest that a model incorporating patient demographics, LVEF, left ventricular mass index and aortic forward flow can accurately predict SBP and DBP. CMR-modelled SBP and DBP exhibit modest correlation with brachial measurements, especially with systolic pressures of >120 mmHg and diastolic pressures of >80 mmHg. Furthermore, our model performed well in an external validation cohort from a different centre, indicating its potential robustness in heterogeneous patient populations.

Developing a reliable and accurate model to estimate BP using non-invasive variables, including patient demographics and CMR imaging data, has significant clinical advantages. It could provide haemodynamic indices for CMR imaging, which would aid the assessment of afterload-dependent pathologies such as valvular heart disease and heart failure. We have previously reported on non-invasive CMR-derived pulmonary capillary wedge pressure, itself an index of LV end-diastolic pressure, and therefore, providing an index of preload.13 Therefore, CMR can provide both preload and afterload assessment to interpret contractility.

The present model could be retrospectively applied where brachial BP measurements were not obtained prior to the study. Additionally, the ability to estimate BP without additional measurements would save time, reduce patient discomfort and minimise errors from transient BP fluctuation.

Elevated BP is a leading risk factor for cardiovascular diseases and a significant cause of morbidity and mortality worldwide. The ability to estimate pressure using imaging data has been investigated in several studies using different imaging modalities, including MRI,24–26 ultrasound27 and CT.28–30 CMR imaging is particularly well suited for estimating BP due to its high spatial and temporal resolution and the ability to measure both flow and function.

Our study used a multiparametric approach to estimate BP, including demographic variables (age and body mass index) and CMR-derived functional and flow parameters, including LVEF, LVMI and AO forward flow. Several studies have investigated the relationship between CMR-derived parameters and BP in different patient populations. For example, studies have shown that LVMI, left atrial volume index, and aortic distensibility are strongly associated with BP in patients with hypertension.31–33 Our study similarly shows that LVMI correlated with SBP. BP has previously been shown to affect aortic and mitral regurgitation severity, which are dynamic valvular incompetencies.34 35 These studies highlight the importance of accurate BP estimation in assessing cardiovascular diseases.

Limitations

Although we included patients from multiple databases, a limitation of our study is the relatively small sample size of the derivation cohort, which could limit the generalisability of our results. Furthermore, patients included in the study were either healthy volunteers or those with cardiovascular diseases. The variations in demographics and clinical parameters, such as age, gender distribution and SBP, between the two cohorts might introduce some level of heterogeneity that could potentially impact the generalisability of our findings. It is unclear whether our model would perform similarly in other patient populations, such as those with renal or pulmonary diseases. Therefore, the transferability of our results to broader populations should be interpreted with caution, as our findings may have specific relevance to cohorts with similar characteristics. Future studies should aim to address this limitation by using more homogeneous populations, where possible, to validate further and extend the applicability of our results. Brachial BP measurements were taken before the CMR scan in the resting state so that they reflected the haemodynamics during imaging assessment and would be akin to a clinic BP recording. They do not necessarily reflect the average long-term or ambulatory BP. They might also be subject to fluctuations such as those due to white-coat syndrome or anxiety regarding CMR assessment on the day.

Conclusion

In conclusion, our study demonstrates that a CMR-modelled SBP and DBP can estimate BP using patient demographics and CMR-derived functional and flow data. This model has the potential to provide afterload assessment during CMR imaging and could aid in the management of complex cardiovascular diseases. Furthermore, the ability to estimate BP without additional measurements would save time, could be retrospectively applied, reduce patient discomfort and minimise errors. Future studies are needed to validate our findings in larger patient populations and also with ambulatory BP assessments.

Data availability statement

Data are available on reasonable request. The datasets generated and analysed during the current study are not publicly available. Access to the raw images of patients is not permitted since specialised postprocessing imaging-based solutions can identify the study patients in the future. Data are available from the corresponding author on reasonable request.

Ethics statements

Patient consent for publication

Ethics approval

This study involves human participants and the Norfolk and Norwich University Hospital and the University of East Anglia approved the study for the Norwich cohort as an audit and observational retrospective study (2020/21-075). Because this was an observational retrospective study, patient consent was waived. The study for the Sheffield cohort was approved by Sheffield Teaching Hospitals and the National Research Ethics Service (16/YH/0352) in the UK. Written informed consent was obtained. The institutional review boards approved the study for the Singapore cohort, and every subject provided written informed consent. The study complied with the Declaration of Helsinki. Participants gave informed consent to participate in the study before taking part.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Twitter @HosamAssadi, @NA, @Vass_Vassiliou

  • Contributors HA: data curation, data collection, formal analysis, methodology, statistical analysis, visualisation, writing-original draft preparation, writing-review and editing. GM: formal analysis, visualisation, writing-original draft preparation, writing-review and editing. XZ: data curation, data collection, formal analysis, writing review, editing and manuscript revision. RL and BK: data anonymisation, writing review, editing and manuscript revision. SA: data collection, independent review, editing and manuscript revision. CG-C, ZM and VL: contextual input, writing review, editing and manuscript revision. RG, G-KY, IH, PS, DPR, LZ, VV and AJS: critical contextual review, independent review, editing and manuscript revision. RvdG: conceptualisation, software solution, methodology, critical contextual review, independent review, editing and manuscript revision. PG: conceptualisation and conceived the study design, supervision, data curation, funding acquisition, resources, formal analysis, methodology, statistical analysis, critical contextual review, visualisation, writing-original draft preparation, writing review and editing. PG guarantees the overall content of the article on behalf of all authors. All authors listed have contributed sufficiently to the project to be included, and all those qualified to be authors are listed in the author byline. All authors read and approved the final manuscript.

  • Funding PG and AJS are funded by Wellcome Trust Clinical Research Career Development Fellowships (220703/Z/20/Z & 205188/Z/16/Z). GM is funded by the National Institute of Health Research (NIHR). For the purpose of Open Access, the authors have applied a CC BY public copyright licence to any Author Accepted Manuscript version arising from this submission.

  • Disclaimer The funders had no role in study design, data collection, analysis, publishing decisions, or manuscript preparation.

  • Competing interests PG is a clinical advisor for Pie Medical Imaging and Medis Medical Imaging. PG is consultant for Anteris and Edwards Lifesciences. All other authors have no competing interests to declare.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.