ABSTRACT
Aim To study the association between the life-course accumulation of health deficits and later-life heart size and function using data from the 1946 National Survey of Heath and Development (NSHD) British birth cohort, the world’s longest running birth cohort with continuous follow-up.
Methods and Results A multidimensional health deficit index (DI) looking at 45-health deficits was serially calculated at 4 time periods of the life-course in NSHD participants (0-16, 19-44, 45-54 and 60 to 64 years), and from these the mean and total DI for the life-course was derived (DImean, DIsum). The step change in deficit accumulation from one time period to another was also calculated. Echocardiographic data at 60-64 years provided: ejection fraction (EF), left ventricular mass indexed to body surface area (LVmassi, BSA), myocardial contraction fraction indexed to BSA (MCFi) and E/e’. Generalized linear models assessed the association between DIs and echocardiographic parameters after adjustment for sex, socio-economic position and body mass index.
1,375 NSHD participants were included (46.47% male). For each single new deficit accumulated at any one of the 4 time periods of the life-course, LVmassi increased by 0.91-1.44% (p<0.013), while MCFi decreased by 0.6-1.02% (p<0.05 except at 45-54 years). One unit increase in DI at age 45-54 and 60-64 decreased LV EF by 11-12% (p<0.013). A single deficit step change occurring between 60-64 years and one of the earlier time periods, translated into significantly higher odds (2.1—78.5, p<0.020) of elevated LV filling pressure defined as E/e’>13.
Conclusion The accumulation of health deficits at any time period of the life-course associates with a maladaptive cardiac phenotype in older age, dominated by myocardial hypertrophy and poorer function. The burden of health deficits appears to strain the myocardium potentially leading to future cardiac dysfunction.
INTRODUCTION
Aging is an inevitable process of biological life. In general, as humans age they accumulate more health deficits which are translated into adverse health outcomes, complex medical needs and high financial burden on society1. However, people do not age and accumulate health deficits at the same rate, highlighting the difference between chronological and biological age.
The accumulation of new health deficits throughout life results in a multidimensional loss of reserves spanning the physical and cognitive domains. This health deficit burden in an individual can be quantified and there are mainly three schools of thought for how to do this: 1) the rules-based operational definition approach relies on a set of phenotypical rules of performance against which the individual is scored2; 2) the clinical frailty score grades persons along a spectrum from very fit to very frail based on clinical judgment; and 3) the health deficit index (DI, also known as frailty index) counts health deficit accumulation and it is the ratio between the number of deficits in an individual, and the total number of deficits appraised, with sub-unitary scores for partial deficits3. The DI is most widely used in research because it is easy to score, it has the mathematical advantage of accommodating partial deficits, and it correlates with the clinical frailty score4. The DI should not be confused with a physical frailty phenotypic score that only appraises domains such as shrinking, exhaustion, low physical activity, slow gait speed and weak grip strength. The DI has been extensively validated and shown to predict mortality5 and other negative health outcomes6.
Cardiovascular diseases (CVDs) are an important component of the multi-morbidity syndrome accounting for more than 30% of all global deaths. Besides the traditional risk factors such as diabetes and hypertension, it is not known whether multi-morbidity and its quantifiable equivalent, the DI, are related to the appearance of CVD. It is also unknown whether a greater burden of health deficits in early-life predicts later-life CVD. To answer these questions, we analyzed life-course data of participants from the 1946 Medical Research Council (MRC) British National Survey of Health and Development (NSHD)-the world’s longest running birth cohort with continuous follow-up. We calculated DIs at different periods of the life-course and examined their association with cardiac size and function at age 60-64 years. The importance of the 7th decade stems from the current 80-year human life-expectancy-having cardiac dysfunction at 60 translates into either premature death or living two further decades but with a lower quality of life.
METHODS
Study design
Participants were from the MRC NSHD, a birth cohort study comprised of 5,362 individuals born in 1 week in 1946 in Britain. The cohort has been extensively followed up with periodic assessments which have been described elswehere7. The 2006-2010 NSHD data collection was granted ethical approval from the Greater Manchester Local Research Ethics Committee and the Scotland Research Ethics Committee. Participants’ weight and height were measured at 60-64 years and used to compute body mass index (BMI) and body surface area (BSA). Participants’ socioeconomic position (SEP) was evaluated at the time of echocardiography (60-64 years) or at 53 years where the former was not available, according to the UK Office of Population Censuses and Surveys Registrar General’s social class.
Construction of the health deficit index
The definition of a health deficit is broad and includes symptoms, signs, medical conditions, disabilities, laboratory and imaging abnormalities8. In line with the Rockwood approach3 deficits included in the DI (Figure 1, Supplementary Tables S1-S3) had to meet the following criteria: health related; prevalence between 1% to 80%; span a wide variety of medical systems; contain less than 5% missing data4. To sequentially apply the same index on cohort participants at different time points, we appraised the same health deficits at each period of the life-course. For reliability, we aimed for an index with at least 30 health deficits in line with previous works3.
The 45 health deficit categories, spanning a broad set of physical and cognitive domains, appraised in MRC NSHD 1946 British birth cohort participants.
MRC NSHD, Medical Research Council National Survey of Health and Development; GI, gastrointestinal.
Four time periods were examined: 0 to 16 years representing early-life (DI0_16), 19 to 44 years representing young adulthood (DI 19_44), 45 to 54 representing early rmiddle age (DI45_54) and 60 to 64 representing older age (DI60_64).
Health deficit score was “0” for complete absence, “1” for absolute presence, and “0.5” for partial deficits. Irreversible medical events (e.g. stroke) carried their deficit score onto all subsequent time periods, while reversible events (e.g. bone fracture) were scored de novo at each time point. To calculate the DI, the sum of health deficits present in an individual was divided by the total number of possible deficits in our index (that is 45) minus the number of missing deficits for that participant (Equation 1). The resultant DI can be a number from 0 (least) to 1 (most).
Equation 1: Formula for calculating the DI.
To assess whole-of-life health deficit burden the 4 individual DIs were averaged (DImean) and summed (DIsum). To assess health deficit change between time periods, the difference between DIs at the 4 time points were computed as follows: DI2–1=(DI19_44)-(DI0_16); DI3_1 = (DI45_54)-(DI0_16); DI4_1 = (DI60_64)-(DI0_16); DI3_2 = (DI45_54)-(DI19_44); D14_2 = (DI60_64)-(DI19_44); and DI4–3=(DI60_64)-(DI45_54).
Echocardiography
Between 2006-2010 when study members were 60-64 years, British-based NSHD participants who had not been lost to follow-up or withdrawn, were invited to attend a clinic-based assessment that included resting transthoracic echocardiography using General Electric (GE) Vivid I machines. Briefly, the echocardiographic protocol7 included long and short axis parasternal views, apical 5-, 4, 3- and 2-chamber views, aortic short axis view, and conventional/tissue Doppler in the apical 4-chamber. Reproducibility between readers was previously reported (intraclass correlation coefficients >0.80)9. Ejection Fraction (EF) was calculated by the biplane Simpson’s method. Left ventricular (LV) mass was indexed by BSA to obtain the LV mass index (LVmassi) as per the American Society of Echocardiography recommendations. Myocardial contraction fraction index (MCFi) was calculated as the ratio between stroke volume and myocardial volume. Myocardial volume was calculated by dividing LVmassi by mean myocardial density (1.04 g/ml)10. E/e’ ratio was calculated by dividing peak mitral valve velocity (E) by the average of lateral and septal mitral annular early diastolic velocity (e’)11. LV end-diastolic volume index (LVEDVI) was derived by dividing LVEDV by BSA.
Statistics
Statistical analysis was performed in R (version-3.6.0). Distribution of data were assessed on histograms and using Shapiro-Wilk test. Continuous variables are expressed as mean ± 1 standard deviation (SD); categorical variables, as counts and percent. Regression models were developed using DIs or their derivatives as exposures to predict LV EF, LVmassi, MCFi and E/e’ by echocardiography at 60-64 years. Model 1 was adjusted for sex, Model 2 for sex and SEP, and Model 3 for sex, SEP and BMI. For LVmassi and MCFi, anthropometry was already accounted for through BSA indexation, so BMI adjustment was not additionally pursued. This decision during study design was informed through the use of directed acrylic graphs (visual representations of causal assumptions that can identify sources of confounding) that helped identify potential colliders (in this case BMI) that should be left uncontrolled, versus non-colliders (or confounders) that should be controlled in the regression models12.
As a result of the skewed distribution of echocardiographic parameters, generalized linear models (glm) with gamma distribution and log link were used to investigate the association of DIs with the following continuous outcome variables: LV EF, LVmassi and MCFi; while glm with binomial distribution and logit link (equivalent to logistic regression) was employed for the binary outcome variable E/e’ using a cut-off of >13. A two-tailed p-value <0.05 was considered statistically significant. Model assumptions were verified with regression diagnostics to ascertain linearity, normality, homoscedasticity and independence of the residuals.
We ran sensitivity analyses in which we: (1) assessed whether the associations of DI scores with echocardiographic parameters remained unchanged if all cardiovascular-related health deficits were excluded from the DIs (we removed hypertension, ischaemic heart disease, myocardial infarction, heart failure, heart rhythm abnormality, congenital heart disease, rheumatic heart disease, and other cardiovascular diseases); (2) assessed whether additional adjustment for childhood SEP (over and above adjustment for sex, BMI and late-adulthood SEP) was pursued; (3) re-analyzed associations between DIs and EF/LVmassi as binary rather than continuous outcome variables using an EF cut-off of 55% and LVmassi cut-offs of 102 and 88g/m2 for men and women respectively, as per British Society of Echocardiography Guidelines for Chamber Quantification13; (4) compared the echocardiographic parameters at age 60-64 of excluded participants with >20% missing data (for which DIs were not calculated) vs. included participants who had <20% of missing data; (5) compared participants with complete echocardiographic variables and those with one or more missing echocardiographic parameter; (6) simulated a complete case analysis in terms of the outcome variables, through multiple imputation. To verify the reliability of the health deficit index, we performed a principal component analysis (PCA) of the included health deficits. We constructed a PCA-weighted health deficit index (wDI) via the weighted average of the health deficits using the first principal component weights (Supplementary Table S2). This procedure was done for two age-interval exemplars: young adulthood (19 to 44 years) and older age (60 to 64 years). For each of these wDIs, we then re-ran the fully adjusted regression models using wDI19_44, and wDI60_64 respectively. To check for different profiles in the scores, we performed a multiple correspondence analysis (MCA) based on the indicator matrix. A scree plot was used to extract the optimal number of dimensions and to determine the percentages of inertia. From the geometric space created in the MCA, participants were then classified into clusters according to proximity criteria using the k-means algorithm. This separated participants into homogenous groups while maximizing heterogeneity across groups.
RESULTS
Participant characteristics and life-course burden of health deficits
Ninety health deficits were initially considered (Supplementary Table S1). Forty-five health deficits met the criteria for inclusion in the DI (Figure 1, Supplementary Table S2-S3). Of the 2,856 participants invited to attend the clinic visit at 60-64 years, 1,690 attended and 1,653 had echocardiography of which 1,617 had acceptable image quality. Of these, 1,375 study members had <20% missing data and at least one outcome parameter of interest by echocardiography. Out of the 1375 NSHD participants studied, 1 had an DI=0 plus missing health deficit data of up to 20%, while 3 had DI=0 in the absence of missing data.
Participant characteristics and DI results are summarized in Table 1. Generally, higher DIsum and DImean were associated with female sex, lower SEP and higher BMI.
Characteristics of the MRC NSHD sample considering participants with less than <20% missing health deficit data and having at least one echocardiographic parameter of interest (EF, MCFi, LVmassi or E/e’).
For the vast majority of participants DI increased across the life-course, but in 86 participants (0.05%) DI decreased in spite of aging, largely due to the reversible nature of some health deficits. There was no significant difference in echocardiographic parameters between participants whose life-course DIs decreased when compared to the rest of the cohort (LV EF p=0.30, LVmassi p=0.99, MCFi p=0.08, E/e’ p=0.57).
There was no difference in Di between the ages 0-16 and 19-44 years, but increases were observed between young adulthood and early middle age, and between early middle age and older age. Raw health deficit index scores of NSHD participants are displayed in Supplementary Figure S1 and the life-course burden of health deficits in Figure 2. Health deficit contributors at the 4 time periods are shown in Supplementary Figure S2. Women consistently had higher DIs than men at every age interval (Figure 3).
Bar chart illustrating the life-course burden of health deficits in NSHD participants. Health deficit categories have been defined as follows: fit, DI<0.12; low burden, 0.12≤DI<0.24; moderate burden, 0.24≤DI<0.36; high burden, 0.36≤DI.
DI, health deficit index.
Sex differences in life-course burden of health deficits in the NSHD cohort. Females (pink) had higher DIs (y-axis) than males (blue) at every time period (x-axis). Vertical bars represent standard error of the means.
Abbreviations as in Figure 2.
Ejection fraction
The association between DIs and echocardiographic parameters are presented in Table 2. In fully adjusted models, a unit increase in DI (from 0 to 1) in early middle age (DI45_54) or older age (DI60_64) translated into an 11-12% decrease in EF (p=0.013 and 0.002). A unit increase in DIsum and DImean translated into a 4-15% decrease in EF (p=0.011 both). When comparing the DIs of participants with the highest and lowest deciles of EF, those with the lowest EF had higher DIs throughout life except in early life (0-16 years) (Figure 4A).
Associations between health deficit burden at the 4 time periods, step-change in deficit accumulation and whole-of-life health deficit burden, with echocardiographic parameters (EF, LVmassi, MCFi and E/e’) at 60-64 years.
Comparing trajectories of mean life-course DI (y-axis) at each age interval (x-axis) in terms of the older age echocardiographic phenotype. DI trajectories leading to favourable phenotypes shown in green and unfavourable trajectories in red as follows: (A) highest (green) vs. lowest (red) LV EF deciles; (B) highest (red) vs. lowest (green) LVmassi deciles; (C) highest (green) vs. lowest (red) MCFi deciles; and (D) highest (red) vs. lowest (green) E/e’ deciles. Vertical bars represent standard error of the means. Wilcoxon signed-rank test p-values (not shown) for median inter-decile differences were <0.05 at every time point except for the 0-16 years age interval.
EF, ejection fraction; LVmassi, left ventricular mass index; MCFi, myocardial contraction fraction index. Other abbreviations as in Figure 2.
Left ventricular mass index
In fully adjusted models, a unit increase in DIs translated into a 16-82% increase in LVmassi (all p<0.013). The accumulation of a single new health deficit in early-life, young adulthood, early middle age or older age, led to a 1.44%, 1.13%, 0.91% and 0.96% increase in LVmassi respectively. The group of participants with the highest decile of LVmassi had significantly higher DIs throughout life when compared to the group with the lowest decile, except in early-life (Figure 4B).
Myocardial contraction fraction index
In fully adjusted models, a unit increase in DI translated into a 15%-47% decrease in MCFi (all p<0.05). The accumulation of a single new health deficit in early-life, young adulthood, early middle age or older age, led to a 1.02%, 0.85%, 0.60% and 0.71% decrease in MCFi respectively. When comparing the DIs of participants with the highest and lowest deciles of MCFi, those with the lowest MCFi had significantly higher DIs throughout life except in early-life (Figure 4C).
E/e’ ratio
A unit increase in DI in older age increased the log odds of having elevated LV filling pressure (E/e’ >13) by 5.3 (p=0.007) in fully adjusted models. The addition of one new health deficit from early-life to older age (DI4–1) increased the multiplicative odds of having elevated LV filling pressure by 2.1, by 75.4 times for DI4–2 and 78.5 times for DI4–3. Participants with the highest decile of E/e’ had significantly higher DIs throughout life when compared to those in the lowest decile (Figure 4D).
LVEDVi
After adjustment for sex and SEP there was no association between any of the DIs and LVEDVi (Supplementary Table S4).
Sensitivity analysis
After removing cardiovascular health deficits from all the DIs, associations of early-life and/or young adulthood health deficit burden with the 4 key echocardiographic parameters, remained significant for LVmassi and MCFi (Supplementary Table S5). After additional adjustment for childhood SEP, β coefficients were slightly attenuated but DIs mostly retained significant associations with echocardiographic parameters (Supplementary Table S6). There was no significant difference (p>0.10) between the echocardiographic parameters of participants included in the DI analysis compared to those excluded due to >20% missing data (Supplementary Table S7). There was no significant difference in cardiac values (p>0.50) between those with complete or some missing echocardiographic parameters (Supplementary Table S8). Associations between DIs and echocardiographic parameters persisted in fully adjusted models after multiple imputation (Supplementary Table S9). Associations between PCA wDIs and LVmassi, and MCFi persisted both in young adulthood and older age, but were lost for EF and E/e’ (Supplementary Table S10). The cluster analysis revealed no clear separation of data in terms of health deficit patterns in any of the age-intervals. (Supplementary Figures S4-S7).
DISCUSSION
Life-course data from the MRC NSHD British birth cohort indicates that health deficits accumulated during childhood, young adulthood, middle age and older age, and the accumulation of new health deficits between these periods, associate with an adverse cardiac phenotype in older age–one that is characterized by hypertrophy, elevated filling pressure and reduced systolic function.
A health deficit index objectively captures the vulnerability of an individual by counting health deficits. Although theoretically a DI can have any value from 0 to 1, values higher than 0.5 were rare amongst NSHD participants (<1%), as a high burden of health deficits limits survival14. The relatively large proportion of individuals at 0-16 and 19-44 years with low health deficit burden (Figure 2) can be explained by the broad time periods over which health deficits were scored (~20 years) as well as the existence of broadly defined health deficits encompassing multiple clinical entities (Supplementary Table S3).
The burden of health deficits in all 4 time periods of the life-course, including early-life and young adulthood, was associated with higher LVmassi and lower MCFi in later-life (60-64 years). A possible explanation may be that health deficits insidiously strain the heart which eventually manifests as maladaptive myocardial hypertrophy. Even after removing all cardiovascular-related health deficits from the DIs, the associations remained significant, particularly in early-life, confirming the biological relationship linking any kind of global index of health deficits with cardiac function parameters in older age.
We found that NSHD participants with decreasing DIs across the life-course, had fitter hearts (higher EF) when compared to the rest of the cohort. Recognizing the potential long-term cardiovascular sequelae of health deficit accrual in these early formative years opens a window of opportunity for managing health-deficits earlier rather than later. This means firmer implementation of public health strategies aimed at disease prevention, the integrated treatment of conditions as soon as they arise, and provision of more holistic and integrated models of healthcare to ensure that optimal care is provided for all co-existing conditions. Importantly, there is a need to ensure that the priorities and preferences of patients are integrated in all of these strategies to maximize what may well become their life-course engagement with the healthcare system.
Whole-of-life burden of health deficits (DIsum and DImean) was associated with lower EF and MCFi, and higher LVmassi in later-life, more so for DImean. The latter is likely due to the over-inflation of the DI brought about by summation (in DIsum) which absorbs all health deficits including reversible ones. Whole-of-life DIs were not significantly associated with elevated LV filling pressures, but the degree of health deficit accumulation between DI60_64 and any of the preceding time periods (DI4–1, DI4–2 and DI4–3) were. In a similar way, these step changes in health deficit accumulation but also those occurring earlier in life (DI2-1, DI3-1) were associated with lower LV systolic function at 60-64. These data suggest that regardless of one’s baseline health deficit burden, any rapid accumulation of multiple health deficits increases one’s risk of cardiac dysfunction in the senior years. It reinforces the point that even apparently mild de novo clinical issues should be taken seriously, to prevent them from having a long-term destabilizing effect on the myocardium. This is especially important given that with modern medicine, people are living longer, and previously deadly diseases are now being turned into chronic manageable problems contributing to multi-morbidity. It is already well-established that multi-morbidity erodes our well-being and quality of life, predisposing to further health deficit accumulation, premature death, mental health problems, negative health behaviours and reduced self-care conduct15. What the current work adds to the multi-morbidity knowledgebase, is this notion that rapid multiple health deficit accumulation is harmful for the heart.
A strength of the study is the implicit age homogeneity of birth cohort participants enabling age-matching at all subsequent assessments for DIs. Participants were exposed to similar risk factors and had access to similar diagnostic and treatment facilities over time. However, such secular effects may also have led to the underestimation of early-life DIs through underdiagnoses. Conversely, the burden of certain medical conditions may have been over-represented in early-life DIs compared to a modern-day population because treatments exist now that may not have been available then.
A limitation inherent to quantitative health deficit approaches3 is the equal weighting assigned to heterogenous health deficits that are unlikely to be straining the organism in a commensurate way. However, we performed a PCA and showed that the associations persisted for both LVmassi, and MCFi, even when the deficit index was formulated as a weighted average of health deficits. The associations were lost for EF and E/e’ because the weighted index diluted the contributions of cardiovascular deficits upon which EF and E/E’ strongly rely, as previously revealed by the sensitivity analysis (Supplementary Table S5). In addition, we did not account for the possibility of some health deficits spilling over (through causation or modulation) into other conditions. Another limitation is the potential of random zeros due to missing data, but we have shown that the prevalence of this was negligibly low in the analysed sample. The timing of echocardiographic assessment in the 7th decade of life excludes participants that already passed away who might have had the highest DIs and most adverse cardiac phenotype. Moreover, selective follow-up could also be a potential confounder as changes in LV morphology could be a protective mechanism against mortality. There are recognized precision limitations in the GE Vivid I system used to capture all echocardiographic data in this study. Whilst the current data indicate a significant association between health deficit burden and echocardiographic parameters at 60-64 years, the evidence is insufficient to claim a causal line. Despite the association between high E/e’ and elevated LV filling pressure, there are instances where it is unreliable16. Although the prediction of normal and abnormal LV filling pressure is most reliable when the E/e’ is <8 or >15, cut-off of >13 was chosen to avoid sensitivity loss17.
CONCLUSION
The accumulation of health deficits during the life-course associates with a maladaptive cardiac phenotype in older age characterized by hypertrophy and poorer function. It could be that health deficits accumulated during childhood, young adulthood, middle age and older age, insidiously strain the myocardium potentially paving the way to future cardiac dysfunction in susceptible individuals.
Data Availability
NSHD data is available from: https://www.nshd.mrc.ac.uk/data.
SOURCES OF FUNDING
This study was funded by the UK Medical Research Council (program codes MC_UU_12019/1; MC_UU_12019/4; MC_UU_12019/5). G.C. is supported by British Heart Foundation (MyoFit46 Special Programme Grant SP/20/2/34841), the National Institute for Health Research Rare Diseases Translational Research Collaboration (NIHR RD-TRC) and by the NIHR UCL Hospitals Biomedical Research Center. J.C.M. is directly and indirectly supported by the UCL Hospitals NIHR BRC and Biomedical Research Unit at Barts Hospital respectively.
DISCLOSURES
The views expressed in this article are those of the authors who declare that they have no conflict of interest.
ACKNOWLEDGEMENTS
The authors would like to thank all the NSHD members for their participation and continuous engagement with follow-up and all NSHD scientific and data collection teams. The authors are also grateful to Andrew Wong and Adam Moore at the MRC Unit for Lifelong Health and Ageing at UCL for their assistance with NSHD data access.
REFERENCES
