Learning associations between clinical information and motion-based descriptors using a large scale MR-derived cardiac motion atlas
Esther Puyol-Anton Bram Ruijsink Helene Langet Mathieu De Craene Paolo Piro Julia A. Schnabel Andrew P. King




empty  








E Puyol-Antón et al.   

School of Biomedical Engineering & Imaging Sciences, King's College London, UK. St Thomas' Hospital NHS Foundation Trust, London, UK. Philips Research, Medisys, Paris, France Sanofi R&D, Paris

Learning associations between clinical information and motion-based descriptors using a large scale MR-derived cardiac motion atlas
    Esther Puyol-Antón 1 Bram Ruijsink 1,2 Hélène Langet 3 Mathieu De Craene 3 Paolo Piro 4 Julia A. Schnabel 1 Andrew P. King 1
    December 30, 2023
===================================================================================================================================


 The availability of large scale databases containing imaging and non-imaging data, such as the UK Biobank, represents an opportunity to improve our understanding of healthy and diseased bodily function. Cardiac motion atlases provide a space of reference in which the motion fields of a cohort of subjects can be directly compared. In this work, a cardiac motion atlas is built from cine MR data from the UK Biobank (≈ 6000 subjects). Two automated quality control strategies are proposed to reject subjects with insufficient image quality. Based on the atlas, three dimensionality reduction algorithms are evaluated to learn data-driven cardiac motion descriptors, and statistical methods used to study the association between these descriptors and non-imaging data. Results show a positive correlation between the atlas motion descriptors and body fat percentage, basal metabolic rate, hypertension, smoking status and alcohol intake frequency. The proposed method outperforms the ability to identify changes in cardiac function due to these known cardiovascular risk factors compared to ejection fraction, the most commonly used descriptor of cardiac function. In conclusion, this work represents a framework for further investigation of the factors influencing cardiac health.


 Cardiac motion atlas; non-imaging data; dimensionality reduction; multivariate statistics







§ INTRODUCTION

Although much is known about the factors influencing healthy and diseased cardiac development, the recent  availability of large scale databases such as the UK Biobank represents an excellent opportunity to extend this knowledge. The wide range of clinical information contained in such databases can be used to learn new associations with indicators of cardiac function. In this paper we propose a framework for performing such an investigation. 

Left-ventricular (LV) function has been traditionally assessed using global indicators such as ejection fraction and stroke volume. The main drawback of such indicators is that they do not provide localised information about ventricular function, and so their diagnostic power is limited. Motion parameters such as displacements allow for a more spatially and temporally localised assessment of LV function.  Statistical cardiac motion atlases have been proposed as a way of making direct comparisons between motion parameters from multiple subjects <cit.>. We use the UK Biobank database to build a cardiac motion atlas from cine MR images using around 6000 subjects with no known cardiovascular disease. To ensure that the estimated motion parameters are accurate, two automated quality control strategies are proposed to reject subjects with insufficient image quality. We use motion descriptors extracted from the atlas to identify relationships between cardiac function and other clinical information.

Previously, several papers have tried to analyse the influence of clinical information such as smoking, alcohol and body weight on cardiovascular disease <cit.>. However, these studies were based on a limited number of factors and used relatively crude indicators of cardiac health, such as categorical variables indicating disease state. In contrast, we exploit the large scale atlas to make use of more localised motion parameters and also examine a wider range of clinical information. 

Contributions: Our scientific contributions in this paper are threefold. First, we present, for the first time, a large scale cardiac motion atlas built from the UK Biobank database. Second, we propose an automated quality control approach to ensure the robustness of the extracted motion descriptors. Third, we use the atlas to learn associations between descriptors of cardiac function and non-imaging clinical information. This automatic framework represents an ideal tool to further explore factors influencing cardiac health.





§ MATERIALS


The UK Biobank data set contains multiple imaging and non-imaging information from more than half a million 40-69 year-olds. From this data set, only participants with MR imaging data were considered for analysis. Furthermore, participants with non-Caucasian ethnicity, known cardiovascular disease, respiratory disease, diabetes, hyperlipidaemia, haematological disease, renal disease, rheumatological disease, malignancy, symptoms of chest pain or dyspnoea were excluded. 

MR imaging was carried out on a 1.5 Tesla scanner (Siemens Healthcare, Erlangen, Germany) <cit.>. Short-axis (SA) stacks covering the full heart, and two orthogonal long-axis (LA) planes (2-chamber (2Ch) and 4-chamber (4Ch) views) were available for each subject (TR/TE = 2.6/1.10, flip angle = 80). In-plane resolution of the SA stack and LA images was 1.8, with slice thickness of  8 and 6 for SA and LA respectively. 50 frames were acquired per cardiac cycle  (temporal resolution ≈14-24 ms/frame). In addition, the following non-imaging parameters were used for the statistical analysis: age, smoking status (smoker/non-smoker), body mass index (BMI), body fat percentage (BFP), basal metabolic rate (BMR), hypertension and alcohol intake frequency (occasionally/regularly). Demographics of the analysed cohort are presented in Table <ref>.






§ METHODS


In the following Sections we describe the estimation and analysis of the atlas-based motion descriptors. Details of the atlas formation procedure are reported in Sect. <ref>. Sect. <ref> details the quality control (QC) strategies proposed.  Sect. <ref> describes the different dimensionality reduction algorithms used to characterise cardiac function. Sect. <ref> presents the statistical methods used to identify relationships between the motion descriptors and non-motion clinical information.



 §.§ Motion Atlas Formation


This section outlines the procedure used to build the motion atlas, which is based on the work published in <cit.> and <cit.>.

LV Geometry Definition. A fully-convolutional network (FCN) with a 17 convolutional layer VGG-like architecture was used for automatic segmentation of the LV myocardium and blood-pool at end-diastole (ED) for SA and LA slices <cit.>. Each convolutional layer of the network is followed by batch normalisation and ReLU, except the last one, which is followed by the softmax function. In the case of the SA stack, each slice is segmented independently, i.e. in 2D. From the segmentations, a bounding box was generated and used to crop the image to only include the desired field of view, improving pipeline speed and reducing errors in motion tracking. The MR SA and LA segmentations were used to correct for breath-hold induced motion artefacts using the iterative registration algorithm proposed in <cit.>. The motion-corrected LA/SA segmentations were fused to form a 3D smooth myocardial mask. An open-source statistical shape model (SSM) was optimised to fit to this LV binary mask following an initial rigid registration using anatomical landmarks <cit.>. The use of landmarks ensured that the LV mesh was aligned to the same anatomical features for each patient, and the endocardial and epicardial surfaces of the mesh were smoothly fitted to the myocardial mask, providing surface meshes with point correspondence for each of the subjects and a 17 AHA regions segmentation. 

Motion Tracking. A 3D GPU-based B-spline free-form deformation (FFD) was used <cit.> to estimate LV motion between consecutive frames. Subsequently, the inter-frame transformations were estimated and then composed using a 3D+t B-spline to estimate a full cycle  3D+t transformation.

Spatial Normalisation. Similar to <cit.>, we transform each mesh to an unbiased atlas coordinate system using a combination of Procrustes alignment and Thin Plate Spline transformation. We denote the transformation for each subject n from its subject-specific coordinate system to the atlas by  ϕ_n.

Medial Surface Generation. A medial surface mesh with regularly sampled vertices (≈ 1000) was generated from the SSM epicardial and endocardial meshes using ray-casting and homogeneous downsampling followed by cell subdivision. 

Motion Reorientation. The displacements (u_n) for each subject n were reoriented into the atlas coordinate system under a small deformation assumption using a push-forward action: u^atlas_n = J(ϕ_n(r_n))u_n <cit.>, where J(·) refers to the Jacobian matrix, and r_n are the positions of the vertices in the ED frame. 

Transform to Local Coordinate System. For each subject n the displacements in the atlas coordinate system u^atlas_n  were projected onto a local cylindrical coordinate system x^atlas_n, providing radial, longitudinal, and circumferential information. The long axis of the LV ED atlas medial mesh was used as the longitudinal direction.


 §.§ Quality control


To ensure that the atlas was formed only from subjects with high quality cine MR images and robust motion descriptors, two QC methods were implemented to automatically reject subjects with insufficient quality or incorrect segmentations and/or motion tracking. These two methods are:

QC1. SA images with insufficient image quality were excluded using two deep learning based approaches. First the 3D convolutional neural network (CNN) described in  <cit.> was used to automatically identify images with motion-related artefacts such as mistriggering, arrhythmia and breathing artefacts. Second, subjects with incorrectly planned 4Ch LA images were automatically excluded using the CNN described in <cit.>. This pipeline identified images containing the left ventricular outflow tract, which is a common feature of poorly planned 4Ch images. A total of 806 subjects were rejected using the QC1 step.

QC2. Our second QC measure aimed to detect segmentation and/or motion tracking errors. LV volume curves over the cardiac cycle were computed based on the LV geometry and motion fields. Based on these curves, an automatic pipeline was used to identify curves with unrealistic properties. The LV curve was divided into diastole and systole based on the ED and end-systolic (ED) points, where ED corresponds to the point with maximum volume, and ES corresponds to the point with minimum volume. These two curves were treated independently. In the systole part of the curve, the downslope should be smooth and continuous and therefore the first derivative should only contain one peak. Cases with more than one (automatically identified using the first and the second derivative of the LV curve) peak were therefore excluded. In the diastole part the curves were divided into rapid inflow, diastasis and atrial systole phases. In the rapid inflow phase the upslope should be smooth and continuous and the first derivative should only contain one peak. Cases with more than one peak were therefore excluded. In the diastasis phase the curve should be flat and with no variation, so the first and second derivatives were used to exclude any cases with changes of volume higher than 10%. Finally, in the atrial systole phase the upslope should be smooth and continuous and the first derivative should only contain one peak. Cases with more than one peak were therefore excluded. The QC2 step excluded a total of 64 cases.  


 §.§ Dimensionality reduction


In this work, we propose to learn descriptors of cardiac function from the motion atlas using dimensionality reduction techniques. More specifically, three techniques were compared: Principal Component Analysis (PCA) <cit.>, Locally Linear Embedding (LLE) <cit.>, and stacked denoising autoencoders (SDA) <cit.>. To form the inputs to these three algorithms, the local displacements per AHA region were averaged and these were concatenated into a column vector such that for subject n, x^atlas∈ℝ^L, where L= 3AT, T is the number of cardiac phases and A is the number of AHA regions (17) in the atlas medial surface mesh. The column vectors for each subject were combined to produce a matrix X=[x^atlas_1,...,x^atlas_N], where N is the number of subjects. The matrix X is the input for all of the proposed dimensionality reduction techniques, and the output of all of them is a matrix D ∈ℝ^(d × N), where d is the size of the low dimensional space.


 §.§ Statistical analysis


Multivariate statistics were used to find relationships between the low-dimensional motion descriptors and non-imaging clinical information. For categorical variables (i.e. smoking status, hypertension and alcohol intake) a multivariate Hotelling's T^2 test (assuming unequal variance) <cit.> was used to test for differences between groups, and the p-value was recorded. For continuous variables (age, BMI, BFP, BMR) a multivariate regression was used to compute the R statistic that measures the degree of relationship between the data. The R statistic and associated p-values were recorded. All the recorded p-values were adjusted using Bonferroni's correction for multiple tests.  

Furthermore, we analysed the impact of using a cardiac motion atlas compared to standard clinical global measures such as ejection fraction (EF), the current most used clinical global measure of cardiac function. To this end, the EF for each subject was computed from the tracked meshes. For categorical variables, similar to the previous case, a paired t-test (assuming unequal variance) was used to test for differences between groups, and the p-value was recorded. For continuous variables a linear regression was used to compute the R statistic and the associated p-value. As before, all the recorded p-values were adjusted using Bonferroni's correction for multiple tests. 





§ EXPERIMENTS AND RESULTS

 
For all of the proposed dimensionality reduction algorithms the number of dimensions d retained was optimised. In the case of PCA, d was selected to retain 99% of the variance; in the case of LLE a grid search algorithm was used to vary d between 2 and 256 and the number of neighbours s was fixed to 10. Then, d and s were selected to minimise the reconstruction error ϵ. Finally, SDA was applied with 5 layers with size [L, 2000, 1000, 500, d], a learning rate of 0.001, a corruption level of 50%, and number of epochs 500. The value of d was optimised using a grid search between 2 and 256, and d was selected to minimise the reconstruction error ϵ. For the three dimensionality reduction algorithms, the optimised d were: d_PCA = 932; d_LLE = 32 and ϵ_LLE = 0.37 mm; d_SDA = 50 and ϵ_SDE = 0.15mm.


Table <ref> shows the computed p-values with Bonferroni's correction and R statistics for the different clinical data for the different dimensionality reduction algorithms, and also using only EF as a motion descriptor. Bold numbers show p-values less than 0.05 (i.e. that showed statistical significance at 95% confidence). The results show that for all of the dimensionality reduction methods there was no statistical significance for age and BMI, while BFP, BMR, smoking status, hypertension and alcohol consumption are shown to be associated with cardiac function. Furthermore, the R statistics for BFP and BMR for all of the non-linear dimensionality reduction algorithms are around 60%, which indicates a clear linear relationship with the data. Using only EF we can see that there is no statistical significance with any of the clinical information, and the R values are lower compared to the use of the cardiac motion atlas. Non-linear dimensionality reduction algorithms, i.e. LLE and SDA show lower p-values and higher R statistics compared to linear PCA for all of the clinical information analysed. For BFP, BMR, smoking status, hypertension and alcohol consumption, LLE shows a better linear fitting, with slightly higher R statistics compared to SDA. However, the difference between the two methods is not significant. 







§ DISCUSSION


In this work, we have demonstrated the derivation of motion atlas-based descriptors of cardiac function and investigated their association with a range of clinical information. This is the first time that a cardiac motion atlas has been formed from a large number of subjects from the UK Biobank, and this has enabled us to produce a framework for investigating factors influencing cardiac health. We show that our method was able to identify factors known to impact cardiac health (directly or indirectly), such as hypertension, smoking, alcohol intake and body composition. Contrary to previous smaller studies, age was not related to large changes in cardiac function. We were unable to assess whether the apparent lack of association was influenced by a selection bias due to anticipated better survival of healthy versus relatively unhealthy individuals, the relatively small spread in ages (45-70) compared to other studies, the clustering of most participants in the middle age range (55-65), or was a consequence of the dimensionality reduction algorithms. We will investigate these possibilities in future work. In addition, we demonstrated that the use of a cardiac motion atlas has added value over using only global indicators such as EF. The proposed method represents an important contribution to furthering our understanding of the influences on cardiac function. We used the selected categories as an illustration of the capabilities of the atlas method. However, the univariate analysis that we performed does not capture the complex interrelation of health factors, age and cardiac function. Future work will focus on a deeper analysis of the impact of these parameters and will also include further clinical characteristics. In this work, displacements were used as a representation of cardiac motion. However, other representations such as velocity and strain can characterise motion in different ways. Our future work also includes extension of the motion atlas to incorporate strain parameters. In conclusion, this work represents an important contribution to furthering our understanding of the influences on cardiac function and opens up the possibility of learning the relationships between motion descriptors and further non-imaging information.




§ ACKNOWLEDGEMENTS

This work is funded by the Kings College London & Imperial College London EPSRC Centre for Doctoral Training in Medical Imaging (EP/L015226/1) and supported by the Wellcome EPSRC Centre for Medical Engineering at Kings College London (WT 203148/Z/16/Z). This research has been conducted using the UK Biobank Resource under Application Number 17806.




splncs03

