University of Groningen Contribution of chronic disease to the burden of disability

Background: Population ageing is expected to lead to strong increases in the number of persons with one or more disabilities, which may result in substantial declines in the quality of life. To reduce the burden of disability and to prevent concomitant declines in the quality of life, one of the first steps is to establish which diseases contribute most to the burden. Therefore, this paper aims to determine the contribution of specific diseases to the prevalence of disability and to years lived with disability, and to assess whether large contributions are due to a high disease prevalence or a high disabling impact. Methodology/Principal Findings: Data from the Dutch POLS-survey (Permanent Onderzoek Leefsituatie, 2001–2007) were analyzed. Using additive regression and accounting for co-morbidity, the disabling impact of selected chronic diseases was calculated, and the prevalence and years lived with ADL and mobility disabilities were partitioned into contributions of specific disease. Musculoskeletal and cardiovascular disease contributed most to the burden of disability, but chronic nonspecific lung disease (males) and diabetes (females) also contributed much. Within the musculoskeletal and cardiovascular disease groups, back pain, peripheral vascular disease and stroke contributed particularly by their high disabling impact. Arthritis and heart disease were less disabling but contributed substantially because of their high prevalence. The disabling impact of diseases was particularly high among persons older than 80. Conclusions/Significance: To reduce the burden of disability, the extent diseases such as back pain, peripheral vascular disease and stroke lead to disability should be reduced, particularly among the oldest old. But also moderately disabling diseases with a high prevalence, such as arthritis and heart disease, should be targeted. Citation: Klijs B, Nusselder WJ, Looman CW, Mackenbach JP (2011) Contribution of Chronic Disease to the Burden of Disability. PLoS ONE 6(9): e25325. doi:10.1371/journal.pone.0025325 Editor: Jos H. Verbeek, Finnish Institute of Occupational Health, Finland Received February 18, 2011; Accepted September 1, 2011; Published September 22, 2011 Copyright: 2011 Klijs et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: This study was part of the project "Living longer in good health," which was financially supported by Netspar, an independent network for research, education and knowledge exchange in the field of pensions, aging and retirement. This work was supported by the European Public Health Programme (2006109) and the Erasmus MC, University Medical Centre. The funders had no role in study design, data collection, and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. * E-mail: b.klijs@erasmusmc.nl


Introduction
Population ageing is expected to lead to a sharp increase in the occurrence of disability. Disability is associated with an increased need for social services, e.g. healthcare, and a loss in quality of life [1,2,3]. To enable the future health care system to cope with increasing demands, and to avoid strong decrements in the quality of life, it is crucial to develop strategies that effectively lead to reductions in the burden of disability. One of the first steps crucial in developing these strategies is to identify which diseases contribute most to the total burden of disability, but also to clarify whether large contributions are related with a high prevalence of disease or with a high disabling impact, i.e. a high extent the disease leads to disability.
To date, only a limited number of studies have assessed the contribution of specific diseases to the burden of disability [4,5,6,7,8,9,10,11,12,13,14]. Unfortunately, most of these studies were based on relatively old data [4,5,6,7,8,9]. As both the prevalence of chronic diseases and their disabling impacts change over time, the results of these studies may be outdated [13,14].
Furthermore, only two studies explicitly addressed the role of disease prevalence and disabling impact in contributions to disability [5,6]. The most comprehensive study assessing burden of disease has probably been WHOs Global Burden of Disease study (GBD) [12]. In this study, however, disease burden is quantified on the basis of panel valuations of health states rather than actual presence of physical or mental disabilities [12].
Most of the previous studies were based on cause elimination techniques, which provide outcome measures that reflect the reduction in the prevalence of disability if the disease would no longer be present [4,7,8,11]. Drawbacks of this method are that the results may be inconsistent in a situation of co-morbidity, as they depend on the ordering of the elimination, and that contributions of all diseases do not add up to the total disability prevalence. Recently, Nusselder et al. developed a methodology, based on an additive regression technique, that enables exact partitioning of the burden of disability into additive contributions of disease in the presence of co-morbidity [5,15].
The aim of this study was to investigate the current contribution of specific chronic diseases to the total burden of disability in mobility and activities of daily living among the elderly, both in terms of disability prevalence and years lived with disability. This study is among the first to highlight the role of both the prevalence and the disabling impact of specific diseases to determine their contribution to the burden of disability.

Study population
The study population consisted of subjects from seven successive years (2001)(2002)(2003)(2004)(2005)(2006)(2007) of the POLS health and labor survey, which is being conducted by Statistics Netherlands. The survey does not include the institutionalized population. To account for selective non-response and to ensure representativeness for the Dutch noninstitutionalized population, weights were used that were attached to the data.
Information on disabilities in mobility and activities of daily living was collected through face-to-face interviews and information on the presence of chronic disease through written questionnaires. From 2001From -2007,766 subjects were approached and the response was 62%. For our analyses selected elderly subjects who were 55 years and older (n = 17,404) were selected. 22% of these subjects could not be included in the study population because they lacked disease information. Table 1 provides further detailed information on numbers of persons in the study population and Figure 1 on the difference between the source and study population.

Disability
Subjects were asked if they were able to 'walk up and down the stairs', 'walk outside', 'enter/leave the house', 'sit down/get up from a chair', 'move around on the same floor', 'get in/out of bed',  'eat/drink', 'get dressed/undressed', 'wash face/hands' and 'wash completely' and could answer with 'without difficulty', 'with minor difficulty', 'with major difficulty' and 'only with help'. Someone was considered disabled if he or she opted for one of the latter two answers at least once.

Definition of disease groups
Information on the presence of a range of diseases was collected in the questionnaire. From the original questions 'cardiovascular disease' (CVD) was compiled, containing 'stroke', 'heart disease' and 'peripheral vascular disease (PVD)', and 'musculoskeletal disease' was compiled, containing 'back pain', 'arthritis' and 'disorder neck/arm'. Furthermore, 'diabetes mellitus (DM)', 'cancer', 'chronic non-specific lung disease (CNSLD)' and 'other' were distinguished. Appendix S1 shows the original questions and the diseases that were distinguished. Skin cancer was not included because it is not associated with disability. PVD did not include vascular disease of the upper extremity. In subjects who had suffered from back pain or disorders of the neck/arm but indicated that they currently did not suffer anymore, the condition was regarded as no longer present. If a disease was defined on the basis of multiple questions and information on any of the questions was missing while none of the questions indicated the presence of a disease, the disease information was considered missing.

Statistical methods
Disability prevalence by cause was estimated from individual information on the presence or absence of disability, the presence or absence of the selected disease groups, gender and age. A method based on a multivariate additive regression model was used, which is described in more detail elsewhere [5,6,16]. This method takes into account that persons who do not report a disease may be disabled (this risk is referred to as ''background'') and that persons can have more than one disease (co-morbidity). Disability in persons without a reported disease is entirely attributed to background. Disability in persons with at least one disease is attributed partly to background and partly to the disease(s). We assume that causes of disability (diseases and background risk) act as independently competing causes. Assuming independence, the hazard of someone having disease A and B (and zero background risk) equals hazard A + hazard B. In a situation of no competing risk (one disease and zero background risk) the hazard can be easily converted to a probability of being disabled (1-exp(-disease hazard)). When more causes are competing, the probability of being disabled from the specific disease is lower than in the situation of no competition. In the footnote of table 2, a calculation example is given of how was dealt with co-morbidity.
To estimate cause-specific disability prevalence from crosssectional data the following assumptions were made. First, the distribution of disability by cause is explained entirely by diseases that are (still) present at the time of the survey and the risk of disability in absence of any reported diseases. Second, this distribution is proportional to the distribution of the risk of becoming disabled in the time-period preceding the survey. Thirdly, causes of disability (diseases and background) act as independently competing causes.
The regression model is specified as follows: whereŷ y is the estimated probability that the person has disability, e is the base of the natural logarithm and g the linear predictor. The latter is defined as the sum of the background rate by age (a a ) and the cause-specific rates of disability (b d , labeled as ''disabling impact'') for the disease groups (d) that are present in the respondent (given by the dummy variables X d ). Background was handled as a cause that is prevalent for everyone; the rate of background is age dependent (5-year age groups). The disabling impact b d may also vary by age. As the full age-interaction term would require n (number of age classes) times m (number of diseases) different parameters, the rank of the interaction was reduced to one, which means that the age-specific disabling impact of each disease, b da , is estimated as the product of an age patternc a which is equal for each disease, and a disease effect d d , which varies by disease, but not by age (Reduced Rank Regression) [17,18]. While a one rank solution restricts the age pattern to be the same for all diseases, also second rank solutions were fitted (b da~ca1 :d d1 zc a2 :d d2 ) and scaled deviances were compared to test for differences in age patterns for different diseases. Adding a one rank interaction improved the fit of the model (log-likelihood ratio test with P-value of 0.05), indicating that the disabling impact varies by age. Because adding a second rank did not further improve the fit of the model, the same age pattern for all diseases was used. Models were fitted using a quasi-Newton method that was programmed for the statistical package R version 2.7.1. [19]. All analyses were done seperately for men and women, as the loglikelihood ratio test indicated that both the rates of background and the disease-specific rates given sex-specific rates of background differed signifcantly by sex. The significance of the differences in the disabling impact between males and females was assessed by assuming normal distribution of the parameters with the mean of the original ones as the standard error.

Calculation of number of disabled by cause across subgroups
Disability prevalence by cause depends on the prevalence of the disease (X d ) and the disabling impact of the disease (b d or b da ). Analogous to using the proportional distribution of mortality rates to obtain probabilities of death in the presence of competing causes (so called ''crude probabilities''), the attribution of disease d is (b d X d =g):y and of background is (a a =g):y [5,6]. Applying these formulas gives for every individual the probability of being disabled caused by background or disease (if present). Adding the cause-specific probabilities of an individual gives the probability of being disabled for that individual. Adding the cause-specific probabilities of all persons in the dataset, or in a specific age group, gives the total number of disabled by cause in the population, or in that age group. Dividing the number of disabled persons by cause by the total number of persons gives the proportion of disability by cause.
Years lived with disability at age 55 by cause were obtained using the Sullivan method [20]. The Sullivan method uses the prevalence of disability in each age group to divide the number of person-years into years lived with and without disability [20]. Instead of using the age and sex specific prevalence of disability, we used age, sex and cause specific prevalence of disability, yielding the years with disability by each cause. Adding the years with disability by each cause yields the life expectancy with disability, as the sum of the cause-specific disability add to the total prevalence. A life table for the Dutch population 2001-2007, available from the EHEMU database, was used [21].
Confidence intervals around the estimates of disabling impacts, prevalences of disability by cause and contributions to LED at age 55 were obtained with bootstrapping based on 1000 replicas [22].
For the disabling impacts and prevalences by cause we used nonparametric bootstrapping. For the bootstrap of life expectancy, which is used in combination with prevalence by cause to calculate LED, we used parameterized bootstrap, assuming Poisson distribution of the numbers of deaths.
The software for additive regression and for calculation of disability prevalence by cause is available from the authors on request.

Results
In males, the prevalence of disability increased from 4% at ages 55-59 to 20% at ages older than 80. In females, this was from 6% to 37% (Figure 1). Also persons without any disease reported disability.
The contribution of diseases to the prevalence of disability depends both on the prevalence and disabling impact of disease. Among males, arthritis, heart disease and other, and at younger ages also back pain, had the highest prevalences (Table 2). Among females these were arthritis, other, disorder neck/arm and back pain, and at older ages also DM and heart disease.
Among males, PVD had the highest disabling impact (Table 2). Furthermore, stroke, CNSLD and back pain showed high disability risks. Among females, back pain had the highest disabling impact, but stroke and PVD also showed high impacts. In both sexes, cancer did not lead to much disability, as was the case for DM and heart disease in males. The disabling impacts of DM, heart disease, arthritis and disorder neck/arm were significantly higher among females than among males (p,0.05). Although the disabling impact of cancer was low, it was significantly higher among males. The disabling impact of the diseases increased with age (P,0.05).
In males, the most important contributors to the prevalence of disability were musculoskeletal disease and CVD (Table 3, Figure 2). Musculoskeletal disease accounted for 40% of the prevalence of disability below age 65, which was about twice the contribution of CVD. At older ages the two conditions contributed equally. Most of the disability attributed to CVD was caused by Abbreviations: CNSLD = chronic non-specific lung disease; DM = diabetes mellitus; PVD = peripheral vascular disease (upper extremity excluded). 5 year age groups were used to calculate the disabling impact of diseases and the prevalence of disability by cause (table 3). To summarize and as is shown in the table, disabling impacts were also calculated using three age-aggregated groups (55-64.9, 65-79.9, . = 80). 'Background', representing presence of disability irrespective of disease presence, is preferably modelled according to 5 year age groups and could therefore not be shown in the table.
The disabling impact represents the rate of disability from a specific cause given that the disease is present. Adding these specific disability rates for the diseases present and the background rate of disability (by age and gender) gives the total disability rate for a specific exposure group. The proportion of the cause-specific rate in total rate is used to divide the probability of disability in this group by cause. For example, for males aged 75-79 with PVD and arthritis, adding the background rate (0.03), the rate for PVD (0.17), and the rate for arthritis (0.07) yields a total disability rate of 0.27 and a total probability of disability of 0 stroke and PVD. In males younger than 65, most disability attributed to musculoskeletal disease was caused by back pain. At older ages, the largest part was caused by arthritis. CNSLD contributed less than musculoskeletal disease and CVD, but was still responsible for 10-15%. In females, musculoskeletal disease was by far the most important contributor and accounted for 40-50% of the disability burden. Disability attributed to musculoskeletal disease was mostly caused by arthritis, but back pain was also an important cause. The second important cause was CVD, contributing more than 15% in females aged 80 and older. Stroke, heart disease and PVD all three contributed importantly to the disability attributed to CVD. DM was also important and contributed 5-8%. Diseases contributions to years lived with disability were similar to contributions to the prevalence of disability (Table 3, Figure 3). Of the 3.01 years with disability in males, 0.85 were contributed by CVD and 0.86 by musculoskeletal disease. CNSLD contributed 0.40 years. In females, musculoskeletal diseases were responsible for 2.93 of the 6.64 life years with disability. The contribution of arthritis alone (1.75 years) was about 0.8 years more than the contribution of all CVDs together (0.97 years). DM contributed 0.44 years.

Discussion
This study is among the first to investigate contributions of various chronic diseases to the prevalence of disability and is the first to present years lived with disability by cause. Musculoskeletal disease is the main contributor and CVD, particularly important among males, is a second. CNSLD is the third contributor among males and DM among females. Within the group of musculoskeletal disease, arthritis-disorder neck arm contributes mostly by a high prevalence and back pain by a high disabling impact, although the prevalence of this condition is also high. Within the group of CVD, heart disease contributes mostly by its high prevalence and PVD and stroke by its high disabling impact. The disabling impact of all diseases increases with age.

Evaluation of data and methods/limitations
Selection bias may limit the external validity of the results. Possible selection bias caused by non-response (38%) was minimized by using individual weights to adjust for selection effects by age, gender, marital status, urbanization grade, province, employment, health-and smoking status [23].
Twenty two percent of all subjects aged 55 and older lacked disease information and were excluded from the analysis, which led to a slight underestimation of disability in our study (Figure 1). The higher prevalence of disability among non-responders may suggest a slight underestimation of the prevalence and/or disabling impact of some diseases, and hence, that our results should be regarded as conservative.
Our results may not be generalizable to the institutionalized population, which has a higher prevalence of disability and also the relative contribution of specific diseases to the total burden of disability inside institutions may differ from our estimates [24,25]. However, in the Netherlands only a minor part of elderly people lives in an institution (i.e. 90% of those aged 80-85 still live at home), hence, bias due to excluding the institutionalized population is probably small and negligible at younger ages [25].
Due to differences in the extent diseases have remained undiagnosed in the population and due to differences in the reference periods used in the questionnaire, some variation may exist in the extent the prevalences derived from the POLS survey reflect true prevalences. Previous literature showed that self-report of most chronic conditions is fairly accurate, except for arthritis, which may be underestimated as well as overestimated [26]. If the prevalence of arthritis or another disease in our study was underestimated, subjects with less severe forms of the disease would be most likely to be uncounted and, hence, the disabling impact would be overestimated. This implies that although there may be some bias in the estimates of disease prevalence, this is nullified by a bias in opposite direction in the disabling impact. Hence, the bias in the estimates of contributions of specific diseases to the prevalence of disability and life expectancy with disability is likely to have remained small. Self-report of ADL disability has been found to correlate well with performance based measures and, hence, is expected not to have affected the results substantially [27]. Contributions of specific diseases to the prevalence of disability were estimated on the basis of diseases prevalence and disabling impact in the study sample from the POLS health and labor survey, the Netherlands, 2001-2007. The disabling impact represents the rate of disability from a specific cause given that the disease is present. Adding specific disability rates for the diseases present and the background rate of disability (by age and gender) gives the total disability rate for a specific exposure group. The contributions of specific diseases presented in the figure add up to the total prevalence of disability. doi:10.1371/journal.pone.0025325.g002 Using a less stringent definition of disability, defined as one or more items answered with ''at least minor difficulty'' resulted in a lower percentage of disability explained by the diseases included. Our substantive conclusions regarding which diseases contributed most remained unaffected. Due to a lack of power, using a more stringent cut-level could not be evaluated.
The cross-sectional nature our methods and data did not allow identifying cases in which disability was present prior to the onset of the disease. In these cases the disability might be falsely attributed to the disease.
It was decided not to exclude conditions such as 'dizziness with falling' and 'involuntary loss of urine', which have an intermediate position between diseases and disability, but to add them in a separate category. Adding this group mainly reduced the contribution of background, but did not affect the contribution of the other diseases substantially. In both sexes only 'dizziness with falling' and 'involuntary loss of urine' had a high disabling impact and therefore accounted for most of the contribution of other (data available on request).
Mental health conditions are a strong predictor of disease burden and disability onset [12,28]. Unfortunately, mental conditions were not included in the checklist of diseases in the POLS survey. To obtain an impression of the extent mental health issues contribute to the burden of disability, in an additional analysis, we included a score of 60 or lower on the RAND Mental Health Inventory (MHI-5) to the diseases, representing (light or severe) mental health problems [29]. In this analysis, 8-13% of the prevalence of disability was attributed to mental health conditions, at the expense of contributions of most other disease groups and background ( Figure  S1). Together with findings from previous studies, these results suggest that a substantial part of the burden of disability may be associated with mental health problems [10,12]. However, as the MHI-5 is particularly useful for assessing mental health status of the general population and may lack validity to diagnose individual cases, the results need to be interpreted with caution. We decided to provide the results including mental health as supplementary material instead of presenting them in the main analysis.

Comparison with previous studies
In most previous studies musculoskeletal and cardiovascular disease were reported as the most important contributors to disability [4,5,6,7,8,9,11]. Our study confirmed these results also for the oldest old (80+). Additionally, it was shown that the high prevalence of disability at this age is caused by a high prevalence of diseases, but even more by the high extent these diseases lead to disability at older age. This age dependence was not studied in most earlier studies. Furthermore, the current study was the first to identify PVD as an important contributor to the burden of disability, which is associated with much disability mainly by a high disabling impact.
Compared with WHOs Global Burden of Disease study, our method was substantially different [12]. Most importantly, WHOs approach uses disability weights to quantify burden of disease that are based on panel valuations and do not refer to physical (or mental) disabilities only, but represent a broader spectrum of health loss. Additionally, in this approach, each case of disease is assigned a disability weight, irrespective of individual factors such as age and sex which may affect disabling consequences of diseases, and, hence, the estimated burden of disease. In the Global Burden of Disease study, unipolar depressive disorders, alcohol use disorders, hearing loss (adult onset) and Alzheimer and other dementias were associated with most years lost due to disability in high income countries [12]. Differences in methodology and diseases included hampers further comparison with our results.

Interpretation of results
The high prevalence of disability at older ages reflects both an increase in the presence of disabling diseases with increasing age and an increase of the disabling impact. The increasing contributions of background suggest that at older ages also frailty or age-related diseases not included in our study may have contributed [30]. This could for instance be dementia, of which the prevalence by age shows considerable resemblance with the pattern of disability attributed to background in our study [31]. Also conditions that caused persisting disability but are no longer 'present', e.g. falls, might have contributed [32].
Gender differences in the disabling impact of heart disease, DM, cancer, arthritis and disorder neck/arm were found. In addition to differences by gender in self-reporting behavior, differences in physiological, psychosocial and environmental factors may explain these variations. For heart disease, a higher disabling impact among females may be related to a more pronounced decline in lethality among females than males during the last decades, or to gender differences in the nature of heart disease causing better survival among disabled females than among disabled males [33,34]. For DM, a more negative interference of the disease with protective mechanisms in the vascular wall causing thrombogenesis, and a negative influence of female gender on the effect of some cardiovascular risk factors are potential causes of the greater risk for vascular complications in females than in males [35,36,37]. Via cardiovascular complications, these mechanisms may also be responsible for the larger disabling impact among females. Due to differences in hormonal factors, coping styles and anxiety, females are more sensitive to pain than males [38,39,40]. These differences therefore may also explain the differences in the disabling impact of arthritis and disorder neck/arm.

Implications; conclusion
This study clearly shows that diseases contribute to the burden of disability by high disabling impacts, e.g. stroke in males, by moderate disabling impacts but high prevalences, e.g. arthritis, or by both high prevalences and disabling impacts, e.g. back pain in females. Diseases that only have a high disabling impact, e.g. stroke in females, or only a high prevalence, e.g. heart disease in males, do not necessarily contribute much to the burden of disability. The current results showed that the largest contributors are musculoskeletal disorders, particularly arthritis, and CVD. For policy makers this means that the largest reductions in the burden of disability can be obtained by interventions that prevent the primary cause of disability, i.e. that prevent disease onset. Further reductions can be achieved by diminishing disabling impacts. Evidence is accumulating that effective interventions to reduce the extent diseases cause disability, such as home visit and exercise programs, are increasingly available [41,42]. As frail elderly are particularly vulnerable to disability, this group should receive priority [43,44]. The coming decades, the population of oldest old will increase massively, i.e. in 2050 about a quarter of the population of 50 years and older is expected to be older than 80 [45]. The large burden of disability at this age shows that reductions are urgently needed.

Supporting Information
Appendix S1 (DOC) Figure S1 Prevalence of disability by cause, including ill mental health. Abbreviations: CNSLD = chronic non-specific lung disease; CVD = cardiovascular disease; DM = diabetes mellitus; PVD = peripheral vascular disease (upper extremity excluded); MHI-5 = RAND mental health inventory. Contributions of specific diseases to the prevalence of disability were estimated on the basis of diseases prevalence and disabling impact in the study sample from the POLS health and labor survey, the Netherlands, 2001-2007. The disabling impact represents the rate of disability from a specific cause given that the disease is present. Adding specific disability rates for the diseases present and the background rate of disability (by age and gender) gives the total disability rate for a specific exposure group. The contributions of specific diseases presented in the figure add up to the total prevalence of disability. The total prevalence for females aged. = 80 is higher than in the original analysis, which may be related with exclusion of subjects who had information missing on items for MHI-5. (TIF)