Skip to main content
Advertisement
  • Loading metrics

Estimating and characterizing the burden of multimorbidity in the community: A comprehensive multistep analysis of two large nationwide representative surveys in France

Abstract

Background

Given the increasing burden of chronic conditions, multimorbidity is now a priority for healthcare and public health systems worldwide. Appropriate methodological approaches for assessing the phenomenon have not yet been established, resulting in inconsistent and incomplete descriptions. We aimed to estimate and characterize the burden of multimorbidity in the adult population in France in terms of number and type of conditions, type of underlying mechanisms, and analysis of the joint effects for identifying combinations with the most deleterious interaction effects on health status.

Methods and findings

We used a multistep approach to analyze cross-sectional and longitudinal data from 2 large nationwide representative surveys: 2010/2014 waves of the Health, Health Care, and Insurance Survey (ESPS 2010–2014) and Disability Healthcare Household Survey 2008 (HSM 2008), that collected similar data on 61 chronic or recurrent conditions. Adults aged ≥25 years in either ESPS 2010 (14,875) or HSM 2008 (23,348) were considered (participation rates were 65% and 62%, respectively). Longitudinal analyses included 7,438 participants of ESPS 2010 with follow-up for mortality (97%) of whom 3,798 were reinterviewed in 2014 (52%). Mortality, activity limitation, self-reported health, difficulties in activities/instrumental activities of daily living, and Medical Outcomes Study Short-Form 12-Item Health Survey were the health status measures. Multiple regression models were used to estimate the impact of chronic or recurrent conditions and multimorbid associations (dyads, triads, and tetrads) on health status. Etiological pathways explaining associations were investigated, and joint effects and interactions between conditions on health status measures were evaluated using both additive and multiplicative scales.

Forty-eight chronic or recurrent conditions had an independent impact on mortality, activity limitations, or perceived heath. Multimorbidity prevalence varied between 30% (1-year time frame) and 39% (lifetime frame), and more markedly according to sex (higher in women), age (with greatest increases in middle-aged), and socioeconomic status (higher in less educated and low-income individuals and manual workers). We identified various multimorbid combinations, mostly involving vasculometabolic and musculoskeletal conditions and mental disorders, which could be explained by direct causation, shared or associated risk factors, or less frequently, confounding or chance. Combinations with the highest health impacts included diseases with complications but also associations of conditions affecting systems involved in locomotion and sensorial functions (impact on activity limitations), and associations including mental disorders (impact on perceived health). The interaction effects of the associated conditions varied on a continuum from subadditive and additive (associations involving cardiometabolic conditions, low back pain, osteoporosis, injury sequelae, depression, and anxiety) to multiplicative and supermultiplicative (associations involving obesity, chronic obstructive pulmonary disease, migraine, and certain osteoarticular pathologies). Study limitations included self-reported information on chronic conditions and the insufficient power of some analyses.

Conclusions

Multimorbidity assessments should move beyond simply counting conditions and take into account the variable impacts on health status, etiological pathways, and joint effects of associated conditions. In particular, the multimorbid combinations with substantial health impacts or shared risk factors deserve closer attention. Our findings also suggest that multimorbidity assessment and management may be beneficial already in midlife and probably earlier in disadvantaged groups.

Author summary

Why was this study done?

  • The growing burden of chronic conditions has made multimorbidity a priority for healthcare and public health systems worldwide.
  • Over the past 2 decades, multimorbidity has been increasingly investigated in various settings and populations, but the large variations in estimates can partly be attributed to the heterogeneity of settings, designs, and analytical methods.
  • The impact and joint effects of the most frequently associated chronic or recurrent conditions on health have not yet been quantified, and the underlying mechanisms explaining these associations have been seldom explored.

What did the researchers do and find?

  • A comprehensive and multistep analysis of 2 large nationwide representative surveys allowed the estimation and characterization of the burden of multimorbidity in the general population in France.
  • This study clarified several issues regarding the assessment and measurement of multimorbidity, notably in terms of the number and type of conditions to be considered, the type of underlying mechanisms to explain the associations, and the analysis of the joint effects of associated conditions to identify those with most deleterious interaction effects on health.

What do these findings mean?

  • Multimorbidity assessments should move beyond simply counting chronic conditions and include a wide spectrum of chronic conditions. They should especially consider the variable impacts and joint effects of associated conditions.
  • Multimorbid combinations with large health impacts or multiplicative interactions should also be appropriately prioritized for monitoring and management. As multimorbidity does not only affect older people, its assessment and management should begin in midlife and probably earlier in disadvantaged groups.

Introduction

The gains in life expectancy observed worldwide have increased the risk of living with chronic conditions [1]. Over the past 2 decades, multimorbidity, defined as the co-occurrence of at least 2 chronic or recurrent conditions [2], has been increasingly investigated and estimated in various populations, settings, and geographical locations [36]. The broad variations observed in the available estimates of the burden of multimorbidity can be attributed, at least partly, to heterogeneous study methodologies in terms of the number and type of conditions (e.g., diseases, risk factors, symptoms, specific conditions) as well as the thresholds used (e.g., two or three) [7,8]. In particular, the identification of conditions for inclusion in multimorbidity assessments has largely relied on the prevalence, or more commonly, the availability of conditions in existing datasets [8]. The type and meaning of multimorbid associations have subsequently received attention [9], although the etiological pathways have seldom been explored [10]. Significantly, the quantified impact of multimorbid associations on relevant health indicators such as mortality or limitations in daily activities has scarcely been estimated [11,12]. Furthermore, the metrics of the joint effects of the associated conditions (i.e., additive or multiplicative scale), which is a critical issue for measuring multimorbidity, has been largely ignored.

To address this gap, we aimed to estimate and characterize the burden of multimorbidity in the French general population using a comprehensive and multistep analytical approach in order to (1) identify the relevant chronic or recurrent conditions for inclusion in multimorbidity assessments based on their impact on health status; (2) determine the frequency of the co-occurrence of these conditions and their underlying etiological pathways; and (3) evaluate the impact of the number of conditions and their associations as well as the joint effects of associated conditions on health. This study further aimed to contribute to our knowledge of multimorbidity in France where general population studies are rare [6].

Materials and methods

Survey designs and study populations

We used data from 2 large nationwide representative surveys conducted in France between 2008 and 2014, which used similar interview and sampling methods as well as similar lists of chronic conditions and health indicators: the Health, Health Care, and Insurance Survey from 2010 and 2014 (Enquête Santé et Protection Sociale, ESPS) and the Disability Healthcare Household Survey from 2008 (Enquête Handicap–Santé Ménages, HSM). Despite several differences, notably in terms of the time frame used to investigate morbidities that prevented cross-validation or replication in the strictest sense, there was sufficient commonality between the studies to look for convergent evidence based on their parallel analyses.

ESPS is a longitudinal health survey representative of individuals living in households in France (95% of the total population). It collects information about their health status through telephone and face-to-face interviews conducted by specially trained interviewers as well as self-administered questionnaires [13]. In 2010, the participation rate was 65%, resulting in 14,875 participants aged ≥25 years. Half of participants (N = 7,727, one per household) were scheduled to be followed-up to be reinterviewed regarding health status measures in 2014. Linkage with vital statistics, for these participants only, also allowed for the assessment of mortality up to and including 2014.

HSM is a purely cross-sectional two-stage survey conducted in 2008 with a focus on health, disability, and dependency. The participation rates for the first stage to screen individuals with disabilities and the second stage were 80% and 77%, respectively, leading to 23,348 participants aged ≥25 years residing in France being evaluated in face-to-face interviews and self-administered questionnaires in the second stage [14].

Both HSM 2008 and ESPS 2010–2014 received institutional review board approval, and participants provided written informed consent.

Chronic and recurrent conditions

Both surveys used similar interviews and questionnaires, although there were minor variations in the lists and labels of self-reported conditions. A total of 61 chronic or recurrent conditions were recorded, which is at the top end of the range (25 to 75) recommended by Holzer and colleagues for multimorbidy assessement [15]. This list included conditions typically found in previous multimorbidity studies [15] (Table 1). The ESPS survey considered conditions occurring within the past 1-year period, whereas the HSM survey considered lifetime occurrence. Both time frames are relevant, as several conditions have lifelong implications (e.g., cancer, vascular diseases), while others have more temporary effects (e.g., depression, asthma) [5].

thumbnail
Table 1. Frequency and associations of conditions with activity limitations, perceived health, or mortality in the ESPS and HSM surveys.

An X indicates that the condition was independently associated with the health status measure. All estimates are weighted to represent the French general population estimates. The 48 conditions selected for further analysis are highlighted in bold.

https://doi.org/10.1371/journal.pmed.1003584.t001

Health status measures

In addition to the mortality recorded in the ESPS survey, several measurements relating to functioning and perceived health were taken into account using valid and reliable instruments. Both surveys included the European Union Global Activity Limitation Indicator (GALI; extent to which the participant is limited because of a health problem in their daily activities: severely limited/limited but not severely/not limited at all) and Self-Reported Health (SRH; how the participant evaluates his/her general health: very good/good/fair/bad/very bad) [16]. Difficulties in activities of daily living (ADL, N = 7) [17], instrumental activities of daily living (IADL, N = 12) [18], and the Medical Outcomes Study Short-Form 12-Item Health Survey (SF-12) questionnaire were only included in the HSM survey [19]. Change scores were calculated for functioning and perceived health measures in ESPS participants when repeated measurements were available (GALI, SRH). Fig 1 presents the conceptual model assumed in this study. In line with a large body of research in the field of disability and perceived health and quality of life [2022], the pathological processes were considered to cause disability and nonoptimal health; age, sex, and socioeconomic status were considered to be the determinants of morbid conditions, disability, and nonoptimal health. Under this model, it is possible to evaluate and quantify the impacts of morbid conditions on health status.

thumbnail
Fig 1. Schematic diagram of the conceptual model assumed in this study.

https://doi.org/10.1371/journal.pmed.1003584.g001

Statistical analysis

We implemented a sequential multistep analytical approach to (1) identify the relevant conditions for inclusion in the multimorbidity analysis; (2) determine the frequency and type of associations of conditions identified as relevant; and (3) estimate the impact of these multimorbid associations on health status.

Identifying the relevant conditions associated with impacts on health status.

We aimed to select conditions that were consistently associated with independent impacts on health status across surveys and/or within surveys across indicators (functioning, perceived health, mortality) and across thresholds in the case of categorical indicators (see below). Multiple regression models were used to assess the independent effects on health status measures of each recorded condition adjusted for age and sex: binary logistic regression for mortality, new limitation in GALI, and health deterioration in SRH; polytomous logistic regression for categorical indicators in GALI (two thresholds: any and only severe limitation), SRH (two thresholds: “fair health” and “good health”), and ADL-IADL (two thresholds: any and limitation in ≥3 ADLs or ≥2 IADLs); and linear regression for SF-12 summary scores. As recommended for over a decade [2325], lasso methods were used to select the conditions to be retained in the final models. Interactions with age and sex were tested for each condition included in the final models. We also compared the explanatory power of the raw count of included conditions relative to the inclusion of binary indicators for all conditions; we used restricted cubic splines to model nonlinear relationships between the count and the health status measures using 5 knots: 2, 3, 4, 6, and 9, respectively, located at the 35 or 50th, 60 or 70th, 88th, 95th, and 99th percentiles of the number of morbidities according to the study.

Determining the frequency and type of associated conditions.

We described the frequency of dyads, triads, and tetrads of the associated conditions retained in the previous step and investigated the etiological pathways that most plausibly explain the most frequent associations (≥0.25%) using the typology by Valderas and colleagues [10] who retained 5 pathways: direct causation, shared risk factors, associated risk factors, confounding by another condition, and only chance (Fig 2). This categorization was practically implemented in a standardized way, through the following steps:

  1. The age- and sex-adjusted odds ratio (OR) for the dyad and the fully adjusted OR, including all conditions independently associated with both components of the dyad, were computed using binary logistic models (with the first component of the dyad as the dependent variable). If the 95% confidence interval (CI) of the age- and sex-adjusted OR included 1, we considered that the dyad could be “explained by chance.” If the 95% CI of the age- and sex-adjusted OR did not include 1, while the 95% CI of the fully adjusted OR included 1, we considered that the dyad could be explained by confounding by another condition. These analyses were conducted in both ESPS and HSM surveys.
  2. When the 95% CI of the fully adjusted OR for the dyad did not include 1 (independent association) in both surveys, a thorough literature search was performed by one of us (JC) on MEDLINE and Google Scholar search engines using the following terms: “association,” “cause,” “determinant,” “complication,” “risk factor,” or “epidemiology” along with each component of the dyad. Systematic reviews as well as primary studies were considered. If a reasonable consensus was found in the literature regarding a causal link between the conditions, we considered that the dyad could be explained by a causal relationship. We subsequently specified whether this relation was bidirectional or due to a complication. Otherwise, if there was a reasonable consensus regarding the shared risk factors, we considered that the dyad could be explained by “shared risk factors” and thus mentioned them. Otherwise, if there was a reasonable consensus regarding the association of risk factors, we considered that the dyad could be explained by “associated risk factors” and thus mentioned them. Several exceptions were made to this categorization to indicate “strongly related pathological processes” such as hip and knee osteoarthritis, ischemic heart disease and peripheral arterial disease (atherosclerosis), or anxiety and depression.
thumbnail
Fig 2. Typology of etiological associations between comorbid conditions.

https://doi.org/10.1371/journal.pmed.1003584.g002

Estimating the impact and joint effects of multimorbid associations.

We described the health impacts of dyads, triads, and tetrads of associated conditions on activities, perceived health, and mortality. Following VanderWeele [26], the joint effects of multimorbid associations were conceptualized as the sum of the effects of each condition and their interaction, with the latter considered on a continuum ranging from subadditive (negative additive interaction) to supermultiplicative (positive multiplicative interaction). The placement on the continuum depended on the relative magnitude of the probability of negative health status in the group in terms of the associated conditions: Associations at the (super)multiplicative end are associated with worse health status than those at the (sub)additive end (Fig 3).

thumbnail
Fig 3. Interaction continuum and 5 categories of interaction on additive and multiplicative scales considered in this study.

RoOR, ratio of odds ratios; RERI, relative excess risk due to interaction.

https://doi.org/10.1371/journal.pmed.1003584.g003

Joint effects and interactions of associated conditions on functioning (severely limited in GALI; limited in ≥3 ADLs or ≥2 IADLs; new limitation) and perceived health (bad or very bad; health deterioration) were thus evaluated on both additive and multiplicative scales by estimating the relative excess risk due to interaction (RERI) [27] and the ratio of odds ratios (RoOR) in the logistic models, including the associations of conditions, age, sex, and all the conditions independently associated with the predicted health indicator as identified in the first step. Due to a limited number of events, it was not possible to study the joints effects on mortality. To ensure that the RoOR and RERI estimates (based on the ORs) were valid, joint-effect analyses were only conducted for health indicators with a prevalence ≤10% (the ORs approximate the relative risks in this context). CIs of RERI were computed according to Zou [28].

All the analyses (descriptive and analytic) were performed separately for the 2 surveys using SAS, version 9.2 software (SAS Institute, Cary, NC, USA) and R software, version 3.6.1 (lasso methods). Appropriate weights were used to provide valid estimates for the French population, while taking into account the unequal probabilities of selection resulting from sample design, nonresponse, and noncoverage in the 2 surveys [13,14]. Due to the data collection methodology used (mostly personal interviews in the 2 surveys), there were very few missing data (<0.005) either for conditions, standard sociodemographic variables, or health status measures. Due to the complex sampling designs, however, classical standard goodness of fit and predictive ability statistics including the Hosmer–Lemeshow statistic and c-index were not available. The Akaike information criterion (AIC) was used for model comparison.

The main analyses were planned before being conducted. Further analyses were only performed to respond to the reviewers’ comments, notably regarding the characteristics of ESPS participants in the longitudinal analyses and the risk of multimorbidity associated with sociodemographic variables, with no influence on the main findings. No data-driven changes were made to the analyses.

This study is reported in line with the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines (S1 STROBE Checklist).

Results

All 14,875 ESPS and 23,348 HSM participants aged ≥25 were included in the cross-sectional analyses. Nearly half of the ESPS sample (N = 7,438, 97% of scheduled participants) was followed up for mortality until 2014, when 3,798 (52%) agreed to participate again in the survey for longitudinal analyses. A description of the 4 analyzed samples is given in S1 Table. There were no important differences between them.

Chronic or recurrent conditions with a negative impact on health status

Overall, 48 out of 61 recorded chronic or recurrent conditions consistently affected activity limitations, perceived health, or mortality (Table 1 and S2 Table). The remaining conditions were not associated with any of the health status measures, with the exception of nutritional anemias, which had an inconsistent impact on GALI. There was no interaction of sex with any assessed risk. An interaction of age was observed for low back pain, depression, and injury sequelae: Younger participants with these conditions had a higher risk of limitations or bad health than older ones. S3 Table details the impact of the studied chronic conditions on functioning and perceived health within the 2 reference time frames (1-year versus lifetime). Although generally similar, the proportion of participants reporting more severe impacts tended to be higher within the lifetime frame (a difference not explained by the age of the affected participants).

The raw count of these 48 conditions was associated with all health status measures (except mortality), although it was less informative than the separate inclusion of all the conditions (AIC 5% to 10% higher than those obtained with the final models). Moreover, the changes in slope (negative) observed for all indicators except for one at 3 and 6 conditions in spline regressions suggest that the health impacts rapidly stabilize, probably due to ceiling effects, when using this count (see below).

Prevalence of, and factors associated with multimorbidity

Overall, multimorbidity was observed in 30.4% of participants (95% CI: 29.6% to 31.3%) using the 1-year time frame of the ESPS survey and 39.0% (95% CI: 38.0% to 40.0%) using the lifetime frame of the HSM survey (Table 2). Regardless of the time frame, the prevalence of multimorbidity increased almost linearly with age and was higher in women than in men: Using the 1-year time frame, it ranged from 7% among 25- to 34-year-old men to 67% among women over 85 years. The greatest increases in prevalence were observed between the age groups of 45 to 54 and 55 to 64 years, and 55 to 64 and 65 to 74 years. Prevalence was also higher in less educated and low-income participants and manual workers, particularly in women and middle-aged participants. Multiple logistic regression analyses confirmed the independent and “dose-effect” associations of age, sex, education level, occupation, and household income with multimorbidity (S4 Table).

thumbnail
Table 2. Weighted frequency of multimorbidity, defined as having at least 2 of the 48 selected conditions, analyzed according to age, sex, and 3 socioeconomic status indicators (education level, occupation, and household income).

Conditions present during the last 12 months (ESPS Survey) and lifetime presence (HSM Survey).

https://doi.org/10.1371/journal.pmed.1003584.t002

The vast majority of participants reporting multiple conditions had up to 4 conditions (74% and 78% in the ESPS and HSM surveys, respectively) (S5 Table). The presence of 6 or more concomitant conditions was uncommon before 75 years but was observed in one-quarter of women after 85 years.

Regardless of participant age and survey time frame, about three-quarters of conditions were more often observed in multimorbid combinations than alone, while 24 and 4 conditions had a median of 2 and 3 associated conditions, respectively (S6 Table).

Associations of conditions and type of multimorbid combinations

The most frequent dyads (N = 223), triads and tetrads (N = 76) of chronic conditions along with the plausible etiological pathways explaining their associations are presented in Table 3 and S7 and S8 Tables. The most frequent combinations include hypertension, low back pain, obesity, osteoarthritis (knee, hip, or other peripheral joints), migraine, diabetes, anxiety, depression, and ear ailments. Fig 4 summarizes the two-by-two independent associations (dyads) between the 48 selected chronic conditions by visualizing 3 main aggregates of associations that are consistent across time frames: (i) vasculometabolic, including diseases of the circulatory system, diabetes, and obesity; (ii) diseases of the musculoskeletal system; and (iii) mental disorders. However, more scattered associations also exist, namely, anxiety and depression with various somatic conditions, especially musculoskeletal disorders, or the latter with vasculometabolic conditions.

thumbnail
Fig 4. Independent associations between the 48 selected chronic conditions.

Associations were identified in logistic models built with forward and backward stepwise procedures to identify the conditions independently associated with each condition while adjusting for age and sex. Upper right panel, above the diagonal: ESPS survey and 1-year time frame; lower left panel, below the diagonal: HSM survey and lifetime frame. Obesity is categorized according to the standard BMI criteria (obese: BMI 30–35; morbidly obese: BMI >35) and analyzed as a 3-category variable. A light blue box indicates an OR ≤2. A star sign (*) indicates an OR between 2 and 4; 2 stars (**) indicate an OR between 4 and 10; 3 stars (***) indicate an OR ≥10. An orange box indicates a negative association (OR <1). Frames indicate conditions within the same ICD-10 chapter. BMI, body mass index; ESPS, Enquête Santé et Protection Sociale; HIV, human immunodeficiency virus; HSM, Enquête Handicap–Santé Ménages; ICD-10, 10th revision of the International Statistical Classification of Diseases and Related Health Problems; OR, odds ratio.

https://doi.org/10.1371/journal.pmed.1003584.g004

thumbnail
Table 3. Weighted frequency of dyads of conditions, strength of associations, and plausible etiological pathways explaining the associations.

All dyads with a frequency of ≥1% in both survey samples are considered. Dyads are presented in decreasing order of frequency (mean frequency based on the 2 surveys).

https://doi.org/10.1371/journal.pmed.1003584.t003

Overall, 44% of the most frequent two-by-two associations were found plausibly explicable by direct causation, either unidirectional as in the case of a complication (e.g., myocardial infarction and heart failure or hypertension and stroke) or bidirectional (e.g., osteoarthritis and obesity or depression and migraine). Furthermore, 30% can be plausibly explained by shared risk factors (e.g., chronic obstructive pulmonary disease (COPD) and ischemic heart disease), or more rarely, by associated risk factors (e.g., cardiac rhythm disorders and low back pain). Several associations (3%) may be explained by strongly related or similar pathological processes (e.g., ischemic heart disease and myocardial infarction or hip and knee osteoarthritis). Finally, for the remaining 23% of the associations, the co-occurrence of 2 conditions can, after adjustment, be attributed to chance (e.g., hypertension and osteoporosis) or confounding factors (e.g., heart failure and low back pain).

Impact and joint effects of associations between conditions (dyads and triads)

The impact on health status varied greatly across the associated conditions (see S9 and S10 Tables for dyads and triads, respectively). With the exception of cardiovascular diseases complicated by heart failure and inflammatory arthritis by osteoarthritis, the combinations with the greatest impact on activities and general health frequently included conditions affecting different bodily systems, notably circulatory or respiratory, musculoskeletal, and nervous. Obesity, diabetes, eye and ear ailments, and mental disorders (especially depression) also contributed to the associations with the most impact. These impacts were consistent across time frames, although they were slightly greater in the lifetime frame.

Interactions between the associated conditions on health status ranged from subadditive to supermultiplicative (Fig 5 and S11 Table). Dyads with multiplicative or supermultiplicative interactions included obesity, chronic obstructive pulmonary disease, migraine accompanying other conditions, inflammatory arthritis and several forms of osteoarthritis, cardiac rhythm disorders, and stroke; dyads with additive interactions included cardiometabolic conditions, musculoskeletal conditions, and injury sequelae with other conditions. Depression and anxiety had a subadditive interaction: When associated with another condition, they had an additive interaction. The patterns of effects were similar across time frames, although the interactions shifted to the additive side for the lifetime frame compared to 1-year time frame (median RERI shift of 0.68 and 0.31 for activity limitations and perceived health, respectively). The same pattern of effects was observed within triads (S12 Table) where the joint effect of the third condition with the other two was similar to when the condition was separately associated with each of the other two.

thumbnail
Fig 5.

Interactions between the associated conditions in dyads on activity limitations (severely limited in GALI, upper right panel above the diagonal) and perceived health (bad or very bad SRH, lower left panel below the diagonal) in ESPS (Panel A) and HSM surveys (Panel B). Interactions are evaluated by the RERI (see S8 Table for values) and represented using a color scale from green (low values: subadditive or additive) to red (high values: multiplicative or supermultiplicative). Only dyads with a frequency of ≥0.25% in at least 1 survey sample are considered. BMI, body mass index; ESPS, Enquête Santé et Protection Sociale; GALI: Global Activity Limitation Indicator; HSM, Enquête Handicap–Santé Ménages; RERI, relative excess risk due to interaction; SRH: Self-Reported Health indicator.

https://doi.org/10.1371/journal.pmed.1003584.g005

Discussion

The use of a comprehensive and multistep analytical approach to analyze 2 large representative general population samples with cross-sectional and longitudinal data on numerous chronic conditions and a broad range of health indicators allowed us to estimate and characterize the burden of multimorbidity in France. This study also makes novel contributions to 3 crucial aspects of multimorbidity assessments and measurements: (1) proposing a comprehensive method for estimating and characterizing the burden based on the health impacts; (2) determining the underlying mechanisms to explain the observed associations; and (3) identifying multimorbid combinations with a high impact and the most deleterious interaction effects on health.

Burden of multimorbidity in France

Our study is the first to characterize the burden of multimorbidity in France. Estimates of the prevalence of specific conditions included in this study are similar to those previously reported in this country [29], while estimates of the prevalence of multimorbidity and specific multimorbid combinations as well as the observed higher prevalence among women are in line with other studies conducted in Europe [5,8] and worldwide [7]. Our study strongly confirms that multimorbidity is not just a problem of old age [3033]: In middle age (35 to 54 years), it affects 1 in 6 men and 1 in 4 women. This study also found marked socioeconomic gradients, with participants, especially women, with lower education and socioeconomic status exhibiting higher multimorbidity rates earlier in the life, as already observed [5,34].

Estimating and characterizing the burden of multimorbidity

This study showed that a large number of chronic conditions (almost 50) significantly and independently impact mortality, activity limitations, or perceived heath in adults. With the exception of the Global Burden of Diseases (GBD) program, which included over 200 noncommunicable diseases in its 2017 updated version [1], most public health agencies worldwide monitor much fewer chronic conditions. For example, the US Centers for Disease Control monitors 10 conditions in the Behavioral Risk Factor Surveillance System [35], while Public Health England generates prevalence estimates for 12 health problems according to “health profiles,” to which modeled estimates of 7 conditions are added [36]. Until a consensus is reached about which chronic conditions and combinations of conditions should be monitored regionally or internationally, this study pleads for the assessment of a large but also reasoned list of conditions that independently impact mortality, disability, or perceived health. This requirement for an independent impact on health enables us to move beyond the debate about which type of condition (risk factors, symptoms, syndromes, diseases) should be considered in the multimorbidity assessment [37]: Any chronic or recurrent condition found in the current nosology (as reflected in the International Classification of Diseases, which includes diseases and sometimes syndromes and symptoms) with an independent impact on health should be considered as long as its frequency in the community justifies its monitoring. The magnitude of the impact may be considered among the prioritization criteria in order to limit the list.

Multimorbid potential of conditions and meanings of the combinations

This study confirmed that several conditions have a high “multimorbid potential,” i.e., a high possibility of being associated with other conditions, which is notably the case with vasculometabolic and musculoskeletal conditions, chronic obstructive pulmonary disease, and especially mental disorders. These conditions were found in many dyads and triads in line with those retrieved by Prados-Torres and colleagues [9] and Violan and colleagues [8]. This suggests that interventions aimed at reducing the incidence of conditions with a high multimorbid potential could be an effective way to alleviate the burden of multimorbidity, which should thus be prioritized. On the contrary, as already reported [38], most cancers, multiple sclerosis, and dementia occur in a more isolated manner, although the underdiagnosis or underreporting of comorbid conditions cannot be excluded, in particular for dementia. One of the original aspects of this study, however, was the analysis of more than 200 dyads in terms of plausible etiological pathways. Using the typology of Valderas and colleagues [10], this analysis gave meaning to the combinations of conditions and provides clues for their clinical management beyond the simple mechanism of complication, which only drives a part of them. In particular, associations probably driven by bidirectional causality are especially important in the clinical setting if one condition is more treatable than another as in the case of depression and low back pain. Associations that share the same risk factors also require attention, as these factors may be more easily detected and tackled in multimorbid participants. In this regard, the major (and often causal) role played by obesity, depression, and anxiety in many associations, either as a component or as a shared risk factor of the components, should be stressed, especially as these conditions are not considered or are even excluded in several multimorbidity studies [39,40].

Impact and joint effects of multimorbid associations

This study is one of the first to analytically address the impact and joint effects of a large number of multimorbid associations on activity limitations and perceived heath. These different but related aspects largely determine the potential severity of the multimorbid associations. As they directly reflect the burden experienced by the participants, the impact of multimorbid associations on health indicators is best investigated and reported in terms of absolute risks. Diseases with complications have the greatest impact on health, but associations of conditions that affect different systems also appear to have large impacts. Associations concerning systems involved in locomotion (cardiovascular, respiratory, osteoarticular), and to a lesser extent, those affecting sensorial functions strongly impact both activity limitations and perceived health, whereas associations with mental disorders mostly impact perceived health. The highly different impacts resulting from the association of 2, 3, or even 4 conditions, which form the bulk of the multimorbidity burden, make the simple counting of conditions inadequate in order to characterize the burden of multimorbidity. This inadequacy has recently been pointed out with regard to the impact on mortality [30,33,41].

The joint effects and especially the interactions of the components of associations, which are more closely related to biological mechanisms and synergy phenomena [42], probably deserve closer attention. In this study, interactions appeared to be extremely variable on a continuum from subadditive to supermultiplicative. The associations with multiplicative or supermultiplicative interactions, which included obesity, chronic obstructive pulmonary disease, migraine, and certain osteoarticular pathologies (inflammatory arthritis and osteoarthritis), were the most deleterious in terms of joint effects. By contrast, associations involving cardiometabolic conditions, low back pain, osteoporosis, injury sequelae, depression, and anxiety were less unfavorable for this effect. These results have important implications for multimorbidity measurements: The variable impacts and joint effects of associations would seem to preclude the exclusive use of a raw count of conditions to quantify multimorbidity.

Strengths and limitations of the study

The strengths of this study include the use of 2 large and nationally representative surveys, including both cross-sectional and longitudinal data on 60 chronic or recurrent conditions and 6 health status measures, and the use of a multistep, systematic, and analytical approach to estimate the impact and the joint effects of conditions on health and to assess the underlying mechanisms explaining multimorbid associations. Moreover, the convergence and consistency of results across surveys, health indicators, and time frames further supports the robustness of our results, which was of particular importance due the type I statistical error that inevitably plagues the study of numerous morbidities.

This study also has some limitations. First, despite the relatively high levels of participation and follow-up rates achieved in the 2 nationwide surveys considered here, various marginal population groups (particularly severely ill participants) may have been underrepresented, thus limiting the external validity of the study. Second, self-reported information on chronic or recurrent conditions could be subject to an information bias due to social desirability (especially mental health conditions), age, and selective recall (symptom-based conditions are reported more frequently in surveys [43]). For many of the studied conditions, frequencies (1-year or lifetime frame) were similar to estimates from previous studies conducted in France except for obesity (when compared to studies where height and weight were measured [44]), and for conditions requiring institutionalization such as dementia or schizophrenia, which were less prevalent than expected. Underestimation of weight may represent a serious problem in terms of the control of confounding factors, since obesity is involved in many multimorbid associations. In other conditions, we presume that these biases mostly had a nondifferential effect on health indicators, possibly with the exception of the impact of disease severity on recall, as more severe conditions were more likely to be recalled and associated with poorer health status (severity was not directly assessed in this study, as in most others conducted to date on multimorbidity [21,45]). However, a lack of sensitivity in the reported information is preferable to a lack of specificity in terms of the nondifferential misclassification bias observed in association studies [46]. Third, the implementation of the approach to investigating the etiological pathways underlying multimorbid associations, using the typology proposed by Valderas and colleagues [10], suffers from several limitations: The empirical analyses used to classify into “chance only” and “confounding” categories may have been vulnerable to incomplete data (unrecorded conditions), to the cross-sectional nature of part of them, and sometimes to the lack of power; and the literature searching reflects the current knowledge, limited as regards causal associations between many conditions. Formal mediation analyses would have been preferable but were far from being achievable in, and also beyond the scope of, this study. Fourth, this study refrained from using data reduction methods such as factor or cluster analyses to investigate disease clustering into “patterns.” The implementation of these methods raises various problems, many of which still need to be resolved [47]. The consideration of disease combinations in dyads/triads/tetrads, as undertaken in this study, actually addressed patterns in three-quarters of multimorbid participants and also shed light on the process of morbidity aggregation (multimorbidity first involves two conditions, then three, then four, and so on). Fifth, unmeasured or incompletely controlled biological confounding cannot be excluded. Indeed, control for confounding is necessary when considering interactions [48]. The implications of this unquantifiable confounding would be primarily an attenuation of the associations. Finally, despite the large sample size, the power was limited to detecting small effects associated with less frequent conditions or even moderate effects in longitudinal analyses. The study of mortality was clearly underpowered due to the low number of deaths observed in the ESPS survey.

Policy and research implications

The lack of methodological consensus or consistent approaches for studying multimorbidity has been a serious impediment for the development of appropriate control and prevention strategies. The proposed approach reported in this paper can guide future contributions with the aim to gain a better understanding and awareness of the different aspects of multimorbidity. Although undertaken in a single country (France), this study is replicable, as exemplified by the use of several survey datasets. It could therefore be readily implemented in many countries with existing systems for the surveillance of chronic diseases.

Multimorbidity is definitely not only a problem of old age; identification and monitoring should be implemented as early as midlife, possibly earlier in disadvantaged groups. Multimorbidity assessments should rapidly move beyond the simple counting of chronic conditions of any kind by including dozens of chronic conditions that impact health and taking into account variable impacts and joint effects of component conditions (e.g., monitoring the most frequent and impacting dyads or triads). Approaches to reduce the incidence of these combinations and minimize their impact should also be a priority for improving the health of multimorbid participants and reducing their impact on healthcare resources.

Further research is needed to evaluate the influence of the data collection method (self-reports, medical records, administrative data) and the comprehensiveness of data sources on medical conditions, which should definitely include obesity and mental disorders. Longitudinal studies of trajectories of multimorbid aggregations according to age, sex, and socioeconomic status and their impacts are also required, along with fundamental and clinical research on pathological interactions, notably between somatic and mental disorders. The clinical implications of the (statistical) concept of the interaction continuum from subadditive to supermultiplicative should be scrutinized. Without much delay, however, multimorbid combinations with large health impacts and most deleterious interactions should be appropriately monitored and managed within the healthcare system.

Supporting information

S1 STROBE Checklist. Checklist for Strengthening the Reporting of Observational Studies in Epidemiology (STROBE).

https://doi.org/10.1371/journal.pmed.1003584.s001

(DOCX)

S1 Table. Description of the studied samples (ESPS and HSM surveys).

https://doi.org/10.1371/journal.pmed.1003584.s002

(DOCX)

S2 Table. Associations of conditions with activity limitations, perceived health, or mortality in the ESPS and HSM surveys.

https://doi.org/10.1371/journal.pmed.1003584.s003

(DOCX)

S3 Table. Impact of studied chronic conditions on activity limitations and perceived health in the ESPS and HSM surveys.

https://doi.org/10.1371/journal.pmed.1003584.s004

(DOCX)

S4 Table. Risk of multimorbidity, defined as having at least 2 of the 48 selected conditions, associated with age, sex, and socioeconomic status indicators (education level, occupation, and household income) as estimated in multiple logistic regression models.

https://doi.org/10.1371/journal.pmed.1003584.s005

(DOCX)

S5 Table. Number of selected conditions and prevalence of multimorbidity according to age and sex.

https://doi.org/10.1371/journal.pmed.1003584.s006

(DOCX)

S6 Table. Risk of multimorbidity (vs monomorbidity) and number of conditions associated with each of the 48 chronic conditions retained.

https://doi.org/10.1371/journal.pmed.1003584.s007

(DOCX)

S7 Table. Weighted frequency of dyads of conditions, strength of associations, and plausible etiological pathways explaining the associations.

https://doi.org/10.1371/journal.pmed.1003584.s008

(DOCX)

S8 Table. Weighted frequency of triads and tetrads of conditions.

https://doi.org/10.1371/journal.pmed.1003584.s009

(DOCX)

S9 Table. Impact of associated conditions (dyads) on activity limitations and perceived health in the ESPS and HSM surveys.

https://doi.org/10.1371/journal.pmed.1003584.s010

(DOCX)

S10 Table. Impact of associated conditions (triads and tetrads) on activity limitations and perceived health in the ESPS and HSM surveys.

https://doi.org/10.1371/journal.pmed.1003584.s011

(DOCX)

S11 Table. Joint effects of dyads of associated conditions on activity limitations (severely limited in GALI, limited in ≥3 ADLs or ≥2 IADLs, new limitation) and perceived health (bad or very bad SRH, new health deterioration) in the ESPS and HSM surveys when the conditions were retained in the final models presented in Table 1.

https://doi.org/10.1371/journal.pmed.1003584.s012

(DOCX)

S12 Table. Joint effects of two-by-two associated conditions within triads on activity limitations and perceived health in the ESPS and HSM surveys when the conditions were retained in the final models presented in Table 1.

https://doi.org/10.1371/journal.pmed.1003584.s013

(DOCX)

Acknowledgments

Disclaimer

The present paper represents the opinions of the authors and does not necessarily reflect the position of their employers.

References

  1. 1. GBD 2017 DALYs and HALE Collaborators. Global, regional, and national disability-adjusted life-years (DALYs) for 359 diseases and injuries and healthy life expectancy (HALE) for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018;392: 1859–922. pmid:30415748
  2. 2. van den Akker M, Buntinx F, Knottnerus AJ. Comorbidity or multimorbidity: what’s in a name? A review of literature. Eur J Gen Pract. 1996;2: 65–70.
  3. 3. van den Akker M, Buntinx F, Metsemakers JF, Roos S, Knottnerus JA. Multimorbidity in general practice: prevalence, incidence, and determinants of co-occurring chronic and recurrent diseases. J Clin Epidemiol. 1998;51: 367–75. pmid:9619963
  4. 4. Schneider KM, O’Donnell BE, Dean D. Prevalence of multiple chronic conditions in the United States’ Medicare population. Health Qual Life Outcomes. 2009;7: 82. pmid:19737412
  5. 5. Barnett K, Mercer SW, Norbury M, Watt G, Wyke S, Guthrie B. Epidemiology of multimorbidity and implications for health care, research, and medical education: a cross-sectional study. Lancet. 2012;380: 37–43. pmid:22579043
  6. 6. Le Cossec C, Perrine AL, Beltzer N, Fuhrman C, Carcaillon-Bentata L. Pre-Frailty, Frailty, and Multimorbidity: Prevalences and Associated Characteristics from Two French National Surveys. J Nutr Health Aging. 2016;20: 860–9. pmid:27709236
  7. 7. Nguyen H, Manolova G, Daskalopoulou C, Vitoratou S, Prince M, Prina AM. Prevalence of multimorbidity in community settings: A systematic review and meta-analysis of observational studies. J Comorb. 2019;9: 2235042X19870934. pmid:31489279
  8. 8. Violan C, Foguet-Boreu Q, Flores-Mateo G, Salisbury C, Blom J, Freitag M, et al. Prevalence, determinants and patterns of multimorbidity in primary care: a systematic review of observational studies. PLoS ONE. 2014;9: e102149. pmid:25048354
  9. 9. Prados-Torres A, Calderón-Larrañaga A, Hancco-Saavedra J, Poblador-Plou B, van den Akker M. Multimorbidity patterns: a systematic review. J Clin Epidemiol. 2014;67: 254–66. pmid:24472295
  10. 10. Valderas JM, Starfield B, Sibbald B, Salisbury C, Roland M. Defining comorbidity: implications for understanding health and health services. Ann Fam Med. 2009;7: 357–63. pmid:19597174
  11. 11. Pererira Nunes B, Flores TR, Mielke GI, Thumé E, Facchini LA. Multimorbidity and mortality in older adults: A systematic review and meta-analysis. Arch Gerontol Geriatr. 2016;67: 130–8. pmid:27500661
  12. 12. Vetrano DL, Rizzuto D, Calderón-Larrañaga A, Onder G, Welmer AK, Bernabei R, et al. Trajectories of functional decline in older adults with neuropsychiatric and cardiovascular multimorbidity: A Swedish cohort study. PLoS Med. 2018;15: e1002503. pmid:29509768
  13. 13. Dourgnon P, Guillaume S, Rochereau T. Enquête sur la santé et la protection sociale (ESPS) 2010. Rapports de l’IRDES. 2012;553: 1–223.
  14. 14. Dos Santos S, Makdessi Y. Une approche de l’autonomie chez les adultes et les personnes âgées. Premiers résultats de l’enquête Handicap-Santé 2008. Etudes et résultats. 2010;718: 1–8.
  15. 15. Holzer BM, Siebenhuener K, Bopp M, Minder CE. Evidence-based design recommendations for prevalence studies on multimorbidity: improving comparability of estimates. Popul Health Metr. 2017;15: 9. pmid:28270157
  16. 16. Robine JM, Jagger C, Euro-REVES Group. Creating a coherent set of indicators to monitor health across Europe: the Euro-REVES 2 project. Eur J Public Health. 2003;13: 6–14. pmid:14533742
  17. 17. Katz S, Ford AB, Moskowitz RW, Jackson BA, Jaffe MW. Studies of illness in the aged. The Index of ADL: a standardized measure of biological and psychosocial function. JAMA. 1963;185: 914–9. pmid:14044222
  18. 18. Lawton M, Brody E. Assessment of older people: Self-maintenance and Instrumental Activities of Daily Living. Gerontology. 1969, 9: 179–86.
  19. 19. Gandek B, Ware JE, Aaronson NK, Apolone G, Bjorner JB, Brazier JE, et al. Cross-validation of Item Selection and Scoring for the SF-12 Health Survey in Nine Countries: Results From the IQOLA Project. International Quality of Life Assessment. J Clin Epidemiol. 1998;51: 1171–8. pmid:9817135
  20. 20. World Health Organization. International classification of functioning, disability and health. Geneva: World Health Organization; 2001.
  21. 21. Marengoni A, Angleman S, Melis R, Mangialasche F, Karp A, Garmen A, et al. Aging with multimorbidity: A systematic review of the literature. Ageing Res Rev. 2011;10: 430–9. pmid:21402176
  22. 22. Guillemin F, Leplege A, Briançon S, Spitz E, Coste J. Perceived health and adaptation in chronic disease. 1st ed. New York: Routledge; 2017.
  23. 23. Harrell FE Jr. Regression modeling strategies: with applications to linear models, logistic and ordinal regression, and survival analysis. 2nd edition. Berlin: Springer; 2015.
  24. 24. Heinze G, Wallisch C, Dunkler D. Variable selection—A review and recommendations for the practicing statistician. Biom J. 2018;60: 431–449. pmid:29292533
  25. 25. Sauerbrei W, Perperoglou A, Schmid M, Abrahamowicz M, Becher H, Binder H, et al. for TG2 of the STRATOS initiative. State of the art in selection of variables and functional forms in multivariable analysis-outstanding issues. Diagn Progn Res. 2020;4:3. pmid:32266321
  26. 26. VanderWeele TJ. The Interaction Continuum. Epidemiology. 2019;30: 648–58. pmid:31205287
  27. 27. VanderWeele TJ, Knol MJ. A Tutorial on Interaction. Epidemiol Methods. 2014;3: 33–72.
  28. 28. Zou GY. On the estimation of additive interaction by use of the four-by-two table and beyond. Am J Epidemiol. 2008;168: 212–24. pmid:18511428
  29. 29. Santé publique France. Géodes. Observatoire cartographique des indicateurs épidémiologiques. Available from: https://geodes.santepubliquefrance.fr/ [cited 2020 November 1].
  30. 30. Jani BD, Hanlon P, Nicholl BI, McQueenie R, Gallacher KI, Lee D, et al. Relationship between multimorbidity, demographic factors and mortality: findings from the UK Biobank cohort. BMC Med. 2019;17: 74. pmid:30967141
  31. 31. Pefoyo AJ, Bronskill SE, Gruneir A, Calzavara A, Thavorn K, Petrosyan Y, et al. The increasing burden and complexity of multimorbidity. BMC Public Health. 2015;15: 415. pmid:25903064
  32. 32. Hu RH, Hsiao FY, Chen LJ, Huang PT, Hsu WW. Increasing age- and gender-specific burden and complexity of multimorbidity in Taiwan, 2003–2013: a cross-sectional study based on nationwide claims data. BMJ Open. 2019;9: e028333. pmid:31182451
  33. 33. Zhu Y, Edwards D, Mant J, Payne RA, Kiddle S. Characteristics, service use and mortality of clusters of multimorbid patients in England: a population-based study. BMC Med. 2020;18: 78. pmid:32272927
  34. 34. Afshar S, Roderick PJ, Kowal P, Dimitrov BD, Hill AG. Multimorbidity and the inequalities of global ageing: a cross-sectional study of 28 countries using the World Health Surveys. BMC Public Health. 2015;15: 776. pmid:26268536
  35. 35. Centers for Disease Control and Prevention. BRFSS Prevalence and Trends data. Available from: https://www.cdc.gov/brfss/brfssprevalence/index.html [cited 2020 November 1].
  36. 36. Public Health England. Public Health Profiles. Available from: https://fingertips.phe.org.uk/ [cited 2020 November 1].
  37. 37. Willadsen TG, Bebe A, Køster-Rasmussen R, Jarbøl DE, Guassora AD, Waldorff FB, et al. The role of diseases, risk factors and symptoms in the definition of multimorbidity—a systematic review. Scand J Prim Health Care. 2016;34: 112–21. pmid:26954365
  38. 38. Tabarés-Seisdedos R, Dumont N, Baudot A, Valderas JM, Climent J, Valencia A, et al. No paradox, no progress: inverse cancer comorbidity in people with other complex diseases. Lancet Oncol. 2011;12: 604–8. pmid:21498115
  39. 39. Nowels MA, VanderWielen LM. Comorbidity indices: a call for the integration of physical and mental health. Prim Health Care Res Dev. 2018;19: 96–8. pmid:28724469
  40. 40. Sakib MN, Shooshtari S, St John P, Menec V. The prevalence of multimorbidity and associations with lifestyle factors among middle-aged Canadians: an analysis of Canadian Longitudinal Study on Aging data. BMC Public Health. 2019;19: 243. pmid:30819126
  41. 41. Chiang JI, Hanlon P, Li TC, Jani BD, Manski-Nankervis JA, Furler J, et al. Multimorbidity, mortality, and HbA1c in type 2 diabetes: A cohort study with UK and Taiwanese cohorts. PLoS Med. 2020;17: e1003094.
  42. 42. Knol MJ, VanderWeele TJ. Recommendations for presenting analyses of effect modification and interaction. Int J Epidemiol. 2012;41: 514–20. pmid:22253321
  43. 43. Violán C, Foguet-Boreu Q, Hermosilla-Pérez E, Valderas JM, Bolíbar B, Fàbregas-Escurriola M, et al. Comparison of the information provided by electronic health records data and a population health survey to estimate prevalence of selected health conditions and multimorbidity. BMC Public Health. 2013;13: 251. pmid:23517342
  44. 44. Lecoffre C, Perrine AL, Blacher J, Olié V. Cholestérol LDL chez les adultes en France métropolitaine: concentration moyenne, connaissance et traitement en 2015, évolutions depuis 2006. Bull Epidémiol Hebd. 2018;37: 710–8.
  45. 45. Fortin M, Stewart M, Poitras ME, Almirall J, Maddocks H. A systematic review of prevalence studies on multimorbidity: toward a more uniform methodology. Ann Fam Med. 2012;10: 142–51. pmid:22412006
  46. 46. Greenland S. Basic Methods for Sensitivity Analysis of Biases. Int J Epidemiol. 1996;25: 1107–16. pmid:9027513
  47. 47. Roso-Llorach A, Violán C, Foguet-Boreu Q, Rodriguez-Blanco T, Pons-Vigués M, Pujol-Ribera E et al. Comparative analysis of methods for identifying multimorbidity patterns: a study of ‘real-world’ data BMJ Open. 2018;8:e018986. pmid:29572393
  48. 48. Pearce N. Analytical implications of epidemiological concepts of interaction. Int J Epidemiol. 1989;18: 976–80. pmid:2621035