Reference ranges for cardiac structure and function using cardiovascular magnetic resonance (CMR) in Caucasians from the UK Biobank population cohort

Abstract

Background

Cardiovascular magnetic resonance (CMR) is the gold standard method for the assessment of cardiac structure and function. Reference ranges permit differentiation between normal and pathological states. To date, this study is the largest to provide CMR specific reference ranges for left ventricular, right ventricular, left atrial and right atrial structure and function derived from truly healthy Caucasian adults aged 45–74.

Methods

Five thousand sixty-five UK Biobank participants underwent CMR using steady-state free precession imaging at 1.5 Tesla. Manual analysis was performed for all four cardiac chambers. Participants with non-Caucasian ethnicity, known cardiovascular disease and other conditions known to affect cardiac chamber size and function were excluded. Remaining participants formed the healthy reference cohort; reference ranges were calculated and were stratified by gender and age (45–54, 55–64, 65–74).

Results

After applying exclusion criteria, 804 (16.2%) participants were available for analysis. Left ventricular (LV) volumes were larger in males compared to females for absolute and indexed values. With advancing age, LV volumes were mostly smaller in both sexes. LV ejection fraction was significantly greater in females compared to males (mean ± standard deviation [SD] of 61 ± 5% vs 58 ± 5%) and remained static with age for both genders. In older age groups, LV mass was lower in men, but remained virtually unchanged in women. LV mass was significantly higher in males compared to females (mean ± SD of 53 ± 9 g/m² vs 42 ± 7 g/m²). Right ventricular (RV) volumes were significantly larger in males compared to females for absolute and indexed values and were smaller with advancing age. RV ejection fraction was higher with increasing age in females only. Left atrial (LA) maximal volume and stroke volume were significantly larger in males compared to females for absolute values but not for indexed values. LA ejection fraction was similar for both sexes. Right atrial (RA) maximal volume was significantly larger in males for both absolute and indexed values, while RA ejection fraction was significantly higher in females.

Conclusions

We describe age- and sex-specific reference ranges for the left ventricle, right ventricle and atria in the largest validated normal Caucasian population.

Background

Quantitative assessment of the cardiac chambers is vital for the determination of pathological states in cardiovascular disease. Intrinsic to this is knowledge of reference values for morphological and functional cardiovascular parameters specific to cardiovascular magnetic resonance (CMR), the most advanced tool for imaging the human heart. CMR has rapidly evolved towards faster and more detailed imaging methods limiting the generalisability of earlier results from relatively small studies [1–4]. More recent studies detailing “normal” ranges for CMR are limited by inclusion of individuals with cardiovascular risk factors such as obesity, diabetes and current smokers in their reference cohort [5, 6].

The UK Biobank is amongst the world’s largest population-based prospective studies, established to investigate the determinants of disease in middle and old age [7]. In addition to the collection of extensive baseline questionnaire data, biological samples and physical measurements, CMR is utilized to provide cardiovascular imaging-derived phenotypes [8].

Based on the UK Biobank participant demographics and health status in ~5000 consecutive participants from the early phase of CMR [8, 9], we aim to select validated normal healthy Caucasian participants in order to establish reference values for left ventricular, right ventricular, left atrial and right atrial structure and function.

Methods

Study population

CMR examinations of 5,065 consecutive UK Biobank participants were assessed. Participants with non-Caucasian ethnicity, known cardiovascular disease, hypertension, respiratory disease, diabetes mellitus, hyperlipidaemia, haematological disease, renal disease, rheumatological disease, malignancy, symptoms of chest pain or dyspnoea, current- or ex-tobacco smokers, those taking medication for diabetes, hyperlipidaemia or hypertension and those with BMI ≥30 kg/m² [10] were excluded from the analysis. In order to create evenly distributed age-decade groups (45–54, 55–64, 65–74), all participants older than 74 years were also excluded from the cohort. (See Appendix 1 for the full list of exclusions).

CMR protocol

The full CMR protocol in the UK Biobank has been described in detail elsewhere [9]. In brief, all CMR examinations were performed in Cheadle, United Kingdom, on a clinical wide bore 1.5 Tesla scanner (MAGNETOM Aera, Syngo Platform VD13A, Siemens Healthcare, Erlangen, Germany).

Assessment of cardiac function was performed based on combination of several cine series: long axis cines (horizontal long axis – HLA, vertical long axis – VLA, and left ventricular outflow tract –LVOT cines, both sagittal and coronal) and a complete short axis stack covering the left ventricle (LV) and right ventricle (RV) were acquired at one slice per breath hold. All acquisitions used balanced steady-state free precession (bSSFP) with typical parameters (subject to standard radiographer changes to planning), as follows: TR/TE = 2.6.1.1 ms, flip angle 80°, Grappa factor 2, voxel size 1.8 mm × 1.8 mm × 8 mm (6 mm for long axis). The actual temporal resolution of 32 ms was interpolated to 50 phases per cardiac cycle (~20 ms). No signal or image filtering was applied besides distortion correction.

Image analysis

Manual analysis of LV, RV, LA and RA were performed across two core laboratories based in London and Oxford, respectively. Standard operating procedures for analysis of each chamber were developed and approved prior to study commencement. CMR scans were analysed using cvi⁴² post-processing software (Version 5.1.1, Circle Cardiovascular Imaging Inc., Calgary, Canada).

In each CMR examination, the end-diastolic phase was selected as the first phase of the acquisition. Observers selected the end-systolic phase by determining the phase in which the LV intra-cavity blood pool was at its smallest by visual assessment at the mid-ventricular level. LV endocardial and epicardial borders were manually traced in both the end-diastolic and end-systolic phases in the short-axis view. In both end-diastole and end-systole, the most basal slice for the LV was selected when at least 50% of the LV blood pool was surrounded by myocardium. In order to reduce observer variability, LV papillary muscles were included as part of LV end-diastolic volume and end-systolic volume, and excluded from LV mass. As an internal quality control measure, the LV mass values in both diastole and systole were checked to ensure they are almost identical. In cases with significant discrepancy, the contours were reviewed and corrected through consensus group approach.

For the RV, endocardial borders were manually traced in end-diastole and end-systole in the short axis view. Volumes below the pulmonary valve were included. At the inflow tract, thin-walled structures without trabeculations were not included as part of the RV. RV end-diastolic and end-systolic phases were denoted to be the same as those for the LV. LV and RV stroke volumes were checked to ensure they were similar.

LA and RA end-diastolic volume, end-systolic volume, stroke volume and ejection fraction were derived by manually tracing endocardial LA contours at end-systole (maximal LA area) and end-diastole (minimal LA area) in the HLA (4-chamber) view. For LA, the same measurements were also derived from the VLA (2-chamber) view and LA volumes were calculated according to the biplane area-length method. Example contours for all four cardiac chambers are provided in Fig. 1.

Fig. 1

Examples of ventricular and atrial contours. The above panels are representative of analysis undertaken on each CMR examination. a and b demonstrate contouring of the left and right ventricle from base to apex at end-diastole and end-systole, respectively. d and e demonstrate contouring of the left and right atrium in the four-chamber view. f and g demonstrate contouring of the left atrium in the two-chamber view

Inter-observer and inter-centre quality assurance aspects

Image analysis was undertaken by a team of eight observers under guidance of three principal investigators. For all cases, analysts filled in progress sheets to monitor any problems in evaluation of CMR data, with any problematic cases flagged, such as a significant discrepancy (defined as more than 10% difference). For such flagged cases all contours and images were reviewed looking for presence of artefacts or slice location problems, operator error or evidence of pathology, such as significant shunt or valve regurgitation. These cases were discussed in regular inter-centre meetings by teleconferencing with respective decisions closed by consensus of at least three team members with relevant knowledge. The team included two biomedical engineers, one radiologist, two career image analysts and six cardiologists. The quality assessment outputs were subject to formal ontological analysis [11]. Inter- and intra-observer variability between analysts for atrial and ventricular measurements was assessed by analysis of fifty, randomly-selected CMR examinations, repeated after a one-month interval.

Statistical analysis

All data is presented as mean ± standard deviation unless stated otherwise. Continuous variables were visually assessed for normality using histograms and Q-Q plots. Independent sample Student’s t-test was used to compare the mean values of CMR parameters between men and women. Outliers were defined a priori as CMR measurements more than three interquartile ranges below the first quartile or above the third quartile and removed from analysis. Mean values for all cardiac parameters are presented by gender and decade (45–54, 55–64, 65–74). Reference ranges for measured (volume, mass) and derived (ejection fraction) data are defined as the 95% prediction interval which is calculated by mean ± t_{0.975, n-1} (√(n + 1)/n) (standard deviation) [12]. Absolute values were indexed to body surface area (BSA) using the DuBois and DuBois formula [13].

The normal ranges for the whole cohort (aged 45–74) were defined as the range where the measured value fell within the 95% prediction interval for the whole cohort regardless of age decade. The borderline zone was defined as the upper and lower ranges where the measured value lay outside the 95% prediction interval for at least one age group. The abnormal zone was defined as the upper and lower ranges where the measured values were outside the 95% prediction interval for any age group.

Pearson’s correlation coefficient was used to assess the impact of age on ventricular and atrial volumes and function. Intra-class correlation coefficients (ICC) were calculated to assess inter- and intra-observer variability, and were visually assessed using Bland-Altman plots [14]. Two-way ICC (2,1) was computed for inter-observer ICCs, to reflect the fact that a sample of cases and a sample of raters were observed, whilst a one-way ICC (1,1) was computed for intra-observer ICC [15]. A p-value <0.05 was considered statistically significant for all tests performed. Statistical analysis was performed using R (version 3.3.0) Statistical Software [16].

Results

A total of 5,065 CMR examinations underwent manual image analysis. 90 subjects were excluded as either the CMR data was of insufficient quality or the CMR identifier did not match the participant identifier. Of the remaining 4,975, 804 (16.2%) met the inclusion criteria. The breakdown of the number of participants meeting individual exclusion criterion is available in Appendix 1. The mean age of the cohort was 59 ± 7 (range 45–74) years. Upon removing outliers, a total of 800 participants (368 males, 432 females) were included in the ventricular analysis and 795 participants (363 male, 432 female) in the atrial analysis (Fig. 2). Baseline characteristics for all participants are provided in Table 1. A summary of CMR parameters stratified by gender is presented in Appendix 2, Tables 13 and 14. The association between CMR parameters and age stratified by gender is included in Appendix 2, Tables 14 and 15.

Fig. 2

Case selection flowchart

Table 1 Baseline Characteristics

CMR left ventricular, right ventricular, left atrial and right atrial reference ranges are provided in a traffic light format for males and females for the whole cohort regardless of their age groups for both absolute and indexed values in numerical format (Tables 2, 3, 4 and 5). These tables are also presented together in a user-friendly poster format for clinical use which is available in Additional file 1. Age-specific reference ranges are also provided in ‘look-up’ tables for those measured CMR values in the borderline (yellow) zone. (Tables 6, 7, 8, 9)

Table 2 Ventricular reference range for Caucausian men

Table 3 Ventricular reference range for Caucausian women

Table 4 Atrial reference range for Caucausian men

Table 5 Atrial reference range for Caucausian women

Table 6 Age-specific ventricular reference ranges for Caucausian men

Table 7 Age-specific ventricular reference ranges for Caucausian women

Table 8 Age-specific atrial reference ranges for Caucausian men

Table 9 Age-specific atrial reference ranges for Caucausian women

Left ventricle

LV end-diastolic volume and LV end-systolic volume were significantly larger in males (LV EDV: absolute = 166 ± 32 ml, indexed = 85 ± 15 ml; LV ESV: absolute = 69 ± 16 ml, indexed = 36 ± 8 ml) compared to females (LV EDV: absolute = 124 ± 21 ml, indexed = 74 ± 12 ml; LV ESV: absolute = 49 ± 11 ml, indexed = 29 ± 6 ml) for both absolute and indexed values. (Appendix 2, Table 12) In men, LV end-diastolic volumes and stroke volumes were lower with older age for both absolute and indexed values. (Appendix 2, Table 14) In women, LV end-diastolic volume, end-systolic volume and stroke volume were smaller with advancing age for absolute and indexed values. LV ejection fraction was significantly greater in females (61 ± 5%) compared to males (58 ± 5%). LV ejection fraction demonstrated no correlation with age in neither males nor females. LV mass was significantly higher in males (103 ± 21 g) compared to females (70 ± 13 g). Upon normalization for body surface area, LV mass did not change significantly with age in either gender. In females, LV mass to end-diastolic volume ratio, a measure of distinct patterns of anatomical adaptations [17], increased significantly (r = 0.14, p <0.01) with age; this was not demonstrated in males.

Right ventricle

RV end-diastolic volume and RV end-systolic volume were significantly larger in males (RV EDV: absolute = 182 ± 36 ml, indexed = 93 ± 17 ml; RV ESV: absolute = 85 ± 22 ml, indexed = 43 ± 11 ml) compared to females (RV EDV: absolute = 130 ± 24 ml, indexed = 77 ± 13 ml; RV ESV: absolute = 55 ± 15 ml, indexed = 33 ± 9 ml) for both absolute and indexed values. Both RV end-diastolic volume and end-systolic volume were lower in older age groups in males and females for absolute and indexed values. RV ejection fraction was significantly higher in females (58 ± 6%) compared to males (54 ± 6%). RV ejection fraction demonstrated a weak but significant positive correlation with advancing age in females only (r = 0.1, p < 0.05).

Left and right atria

Left and right atrial reference ranges are presented in Tables 4, 5, 8 and 9. LA maximal volume and stroke volume, as determined by the biplane method, were significantly larger in males compared to females for absolute values (71 ± 19 vs 62 ± 17 ml) but not for BSA-indexed values (36 ± 9 vs 37 ± 10 ml). LA ejection fraction was almost identical (60% vs 61%) in males and females. Upon normalization for BSA, there was no change in left atrial volumes or function with age in men. In women, indexed LA stroke volume was significantly lower (r = −0.2, p < 0.001) with advancing age.

RA maximal volume and stroke volume were significantly larger in males (RA absolute maximal volume = 93 ± 27 ml, RA absolute stroke volume = 38 ± 14 ml) compared to females (RA absolute maximal volume = 69 ± 17 ml, RA absolute stroke volume = 32 ± 10 ml) for absolute values; upon indexing for BSA, this effect was seen for RA maximal volume only (48 ± 14 vs 41 ± 10 ml). RA ejection fraction was significantly higher (46% vs 41%, p < 0.001) in females compared to males. Upon normalization for BSA, there was no change in right atrial volumes or function with age in males or females.

Intra- and inter-observer variability

Intra and inter-observer variability data is presented in Table 10 and as Bland-Altman plots (representative examples of all observers) in Appendix 3, Figures 3, 4 and 5. Good to excellent intra- and inter-observer variability was achieved for LV and RV end-diastolic volume, end-systolic volume and stroke volume and LA and RA maximal volume and stroke volume.

Table 10 Inter- and intra-observer variability

Discussion

The present study provides clinically relevant age- and gender-specific CMR reference ranges in a traffic light system for the left ventricular, right ventricular, left atrial and right atrial chambers derived from a cohort of 804 Caucasian adults aged 45–74 strictly free from pathophysiological or environmental risk factors affecting cardiac structure or function at 1.5 Tesla.

Whilst determination of reference ranges for CMR has been performed by several previous studies, this work is novel for a number or reasons. Firstly, the substantially larger cohort with strict evidence to ensure participants are free of biological or environmental factors known to impact upon cardiac structure or function differentiates this study from its predecessors. Secondly, reference ranges for CMR parameters are detailed not only by gender but also by age decade, thereby providing increased granularity and clinical utility. Thirdly, previously described findings are reinforced, particularly with respect to age- and gender-related differences in ventricular and atrial parameters. Fourthly, in-depth data surrounding intra- and inter-observer variability is provided.

The validity of a reference range is dependent on a number of factors, including the number of observations available in order to determine the reference interval [12]. This study utilises 800 participants for derivation of left and right ventricular reference ranges. This is a substantial increase compared to the majority of previous studies describing ventricular reference ranges using the SSFP technique: Alfakih et al. [3] (n = 60), Hudsmith et al. [2] (n = 108), Maceira et al. [1] (n = 120) and similar to those published by the Framingham Heart Study group. Similarly, 795 participants are included for derivation of left and right atrial reference ranges. Although previous studies outlining atrial reference ranges have used differing techniques, again, all utilise substantially fewer participants: Sievers et al. [18] (n = 111), Hudsmith et al. [2] (n = 108), Maceira et al. [19, 20] (n = 120). Even a recent systematic review and meta-analysis of normal values for CMR in adults and children is based on smaller numbers than the normal reference ranges presented here [4]. A recently published paper by Gandy and colleagues presents LV reference ranges for 1,515 UK individuals scanned at 3 Tesla [21]. However, their study population includes participants with high plasma B type natriuretic peptide (BNP) levels and blood pressure >149/95 mmHg by design, thus, could not be considered strictly healthy. Le Van et al. describes ventricular and atrial reference values derived from 434 Caucasian adults with similar exclusion criteria to the present study [22]. However, their study examines a much younger cohort, aged 18 to 35 years, and thus the present study complements their findings by investigating an older age range.

Furthermore, this study complied with approved statistical recommendations on derivation of reference limits [12]. Data has been partitioned – dividing reference values by age and sex – in order to reduce variation. The distribution of the reference values was inspected and assessed for normality and values identified as outliers discarded as per our a priori definition.

A total of 5,065 CMR examinations of UK Biobank participants were analysed for this study. Utilising this large population sample permitted a posteriori (retrospective) selection of the reference sample, the preferred method when compiling reference values from healthy individuals [23]. Indeed, only 16% of the original sample were included in this study, with rule-out criteria extending beyond known cardiovascular disease to include traditional cardiovascular risk factors (diabetes mellitus, hypercholesterolaemia, hypertension, current- and ex-tobacco smokers, obesity), cardiovascular symptoms, current or previous cancer, stroke, respiratory, renal or haematological disease and use of certain pharmacological agents. In doing so, a robust definition of what constitutes “health” was created, permitting confidence that reference ranges for cardiovascular structure and function in CMR have been derived from an appropriately selected cohort. This contrasts to the LV reference values published from the Framingham Heart Study Offspring Cohort where the healthy reference group consisted of 47.5% of the total cohort, and exclusion criteria were a history of hypertension, history of use of antihypertensive medication, previous myocardial infarction and heart failure only. Similarly, in the RV reference values study published by the same group, the “healthy reference” cohort included participants with hypertension, diabetes, hypercholesterolaemia and those who were current tobacco smokers [6].

For the left ventricle, our findings that men demonstrated greater volumes and mass compared to females is consistent with both the CMR literature [4] and that derived from other imaging modalities [24, 25]. Our demonstration of decreasing LV end-diastolic and end-systolic volumes with advancing age is also consistent with previous findings. Values for LV end-diastolic volumes are similar to those described by Hudsmith [2], Kawel-Boehm [4] and the Framingham Offspring Cohort group. LV end-systolic volumes were larger, reflecting this study’s methodology of including papillary muscles as part of the LV cavity – the technique most commonly employed when analysing clinical CMR examinations. Consequently, LV ejection fraction mean values and reference intervals were lower than previously reported. Despite this, the finding of a marginally, but significantly, lower LV ejection fraction in men compared to women is consistent with other large cohorts, including the Framingham Offspring Cohort, the Dallas Heart Study cohort [26] and the Multi-Ethnic Study of Atherosclerosis (MESA) cohort [27], although the latter two studies utilised the older gradient-recalled echo sequences. Our study demonstrated no change in LV ejection fraction across age groups, this is consistent with studies across imaging modalities [28, 29]. LV mass, upon normalization for BSA, did not change significantly across age groups in either gender. This is consistent with findings from the MESA cohort, but differs from the Framingham Offspring cohort which demonstrated a significant decrease in BSA-normalised LV mass with age. Autopsy-derived data concerning LV mass in individuals free from hypertension and coronary artery disease and corrected for BSA corroborate findings from our study, suggesting no change in cardiac mass with ageing [30].

For the right ventricle, our findings that males exhibited greater absolute and indexed volumes than females and that volumes were lower with advancing age in both genders are consistent with previously published literature. We demonstrated a larger RV ejection fraction in women compared to men, this is corroborated by Alfakih [3] using both SSFP and gradient-recalled echo sequences and by Foppa and Arora in the Framingham Offspring cohort [6].

For the atrial chambers, no consensus exists regarding the measurement of atrial volumes [4]. In this study, the LA was contoured in the 4-chamber and 2-chamber views and volumes calculated according to the biplane area-length method. Only Hudsmith presented LA reference ranges utilising a similar method with values for LA ejection fraction being almost identical to those described in this study. For the RA, the most recent work regarding reference ranges has been produced by Maceira et al. [20] using three-dimensional modelling which has not been undertaken in this study. Despite different methodology, general findings regarding absolute values being greater in males compared to females and no significant effect of age on RA volumes were replicated in our larger study.

Clinical utility

CMR measurements only provide meaningful information when compared to relevant reference values. However, comparison may be misleading if the CMR examination being considered does not adequately match the reference sample, particularly with regards to age and gender. It is known that cardiovascular disease predominantly affects individuals in middle- and old-age, and it is individuals in these age groups who most commonly undergo CMR examinations. Furthermore, atrial and ventricular structure and function do not remain static over time and undergo changes with age, even in those without evidence of cardiovascular disease. It is in this context that this study presents absolute and BSA-indexed CMR reference values for men and women at three different age groups: 45–54, 55–64 and 65–74.

Intra- and inter-observer variability

For LV and RV end-diastolic volume, end-systolic volume and stroke volume and LA and RA maximal volume and stroke volume, excellent inter- and intra-observer variability was achieved. It is notable, but perhaps not unsurprising, that ICC for derived parameters (i.e. ejection fraction) fell in comparison to those values for directly measured parameters. This is consistent with previous studies examining variability in CMR analysis, such as Margossian et al. [31] and Teo et al. [32], which reported very high inter-observer ICC’s for measured parameters which fell markedly when assessing the ejection fraction.

Study limitations

The reference intervals described were derived from a population of 45–74 year olds of Caucasian ethnicity and therefore may not be generalisable to other ethnic and age groups. As the UK Biobank Imaging project accumulates CMR imaging in up to 100,000 individuals in coming years, analysis of ethnicity effects will become feasible in due course. We included overweight participants with a BMI between 25 and 30 kg/m² in our reference range analysis, even though previous CMR publications, including our own, have shown that obesity affects cardiac structure and function even in an otherwise healthy population [33, 34]. Our rationale for this inclusion was two-fold: firstly, we aligned our inclusion criteria related to BMI with the “Recommendations for Cardiac Chamber Quantification by Echocardiography in Adults: An Update from the American Society of Echocardiography and the European Association of Cardiovascular Imaging” [10]; secondly, given that 2013 data from the UK demonstrates that only 32.9% of men and 42.8% of women had a BMI less than 25 kg/m², arguably our reference ranges represent the “new” normal range and are thus more applicable to the general population [35].

CMR examinations were not performed repeatedly on the same individuals over time, therefore the associations described between age and CMR parameters are not longitudinal, but rather cross-sectional.

Conclusions

This study provides normal reference ranges for all four cardiac chambers derived from the largest healthy cohort of Caucasian adults and will provide utility in the analysis of CMR examinations in both clinical and research settings.

Appendices

Appendix 1

Table 11 Exclusion criteria

Appendix 2

Table 12 Ventricular parameters stratified by gender

Table 13 Atrial parameters stratified by gender

Table 14 Correlation table for ventricular parameters with age

Table 15 Correlation table for atrial parameters with age

Appendix 3

Fig. 3

Exemplar Bland-Altman plots for inter- and intra-observer variability of left ventricular parameters

Fig. 4

Exemplar Bland-Altman plots for inter- and intra-observer variability of right ventricular parameters

Fig. 5

Exemplar Bland-Altman plots for inter- and intra-observer variability of atrial parameters

Appendix 4

UK Biobank data

UK Biobank data in a codified tabular format, received through our access application, was used to select the healthy cohort. Data was translated, using the data dictionary provided as part of the application and the coding tables available through the UK Biobank website, into a self-contained table which we used to perform the analysis. The data derived from the analysis of CMR studies were tested for gross errors such as non-physiological values (e.g., end-systolic volume larger than end-diastolic volume) and were removed from the final dataset.