Next Article in Journal
COVID-19: A Cross-Sectional Study of Healthcare Students’ Perceptions of Life during the Pandemic in the United States and Brazil
Next Article in Special Issue
Academic Promotion of Physicians in Medical Schools: A Special Focus on Primary Health Care in Taiwan
Previous Article in Journal
Missing In-Center Hemodialysis Sessions among Patients with End Stage Renal Disease in Banda Aceh, Indonesia
Previous Article in Special Issue
Family Physicians Working at Hospitals: A 20-Year Nationwide Trend Analysis in Taiwan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combined Multimorbidity and Polypharmacy Patterns in the Elderly: A Cross-Sectional Study in Primary Health Care

by
Grant Stafford
1,2,
Noemí Villén
3,4,
Albert Roso-Llorach
2,4,5,
Amelia Troncoso-Mariño
3,6,
Mònica Monteagudo
2,5 and
Concepción Violán
5,7,*
1
Programa de Máster en Salud Pública, Universitat Pompeu Fabra, 08003 Barcelona, Spain
2
Unitat Transversal de Recerca (UTR), Fundació Institut Universitari per a la Recerca a l’Atenció Primària de Salut Jordi Gol i Gurina (IDIAPJGol), 08007 Barcelona, Spain
3
Àrea del Medicament i Servei de Farmàcia, Atenció Primària Barcelona Ciutat, Institut Català de la Salut (ICS), 08015 Barcelona, Spain
4
Programa de Doctorat en Metodologia de la Recerca Biomèdica i Salut Pública, Universitat Autònoma de Barcelona, Bellaterra (Cerdanyola del Vallès), 08193 Barcelona, Spain
5
Universitat Autònoma de Barcelona, Bellaterra (Cerdanyola del Vallès), 08193 Barcelona, Spain
6
Department of Clinical Sciences, University of Barcelona and IDIBELL, L’Hospitalet de Llobregat, 08907 Barcelona, Spain
7
Unitat de Suport a la Recerca Metropolitana Nord, Fundació Institut Universitaria per a la Recerca a l’Atenció Primària de Salut Jordi Gol i Gurina (IDIAPJGol), 08303 Mataró, Spain
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(17), 9216; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18179216
Submission received: 22 June 2021 / Revised: 25 August 2021 / Accepted: 26 August 2021 / Published: 1 September 2021
(This article belongs to the Special Issue Primary Healthcare)

Abstract

:
(1) Background: The acquisition of multiple chronic diseases, known as multimorbidity, is common in the elderly population, and it is often treated with the simultaneous consumption of several prescription drugs, known as polypharmacy. These two concepts are inherently related and cause an undue burden on the individual. The aim of this study was to identify combined multimorbidity and polypharmacy patterns for the elderly population in Catalonia. (2) Methods: A cross-sectional study using electronic health records from 2012 was conducted. A mapping process was performed linking chronic disease categories to the drug categories indicated for their treatment. A soft clustering technique was then carried out on the final mapped categories. (3) Results: 916,619 individuals were included, with 93.1% meeting the authors’ criteria for multimorbidity and 49.9% for polypharmacy. A seven-cluster solution was identified: one non-specific (Cluster 1) and six specific, corresponding to diabetes (Cluster 2), neurological and musculoskeletal, female dominant (Clusters 3 and 4) and cardiovascular, cerebrovascular and renal diseases (Clusters 5 and 6), and multi-system diseases (Cluster 7). (4) Conclusions: This study utilized a mapping process combined with a soft clustering technique to determine combined patterns of multimorbidity and polypharmacy in the elderly population, identifying overrepresentation in six of the seven clusters with chronic disease and chronic disease-drug categories. These results could be applied to clinical practice guidelines in order to better attend to patient needs. This study can serve as the foundation for future longitudinal regarding relationships between multimorbidity and polypharmacy.

1. Introduction

The global life expectancy at birth has increased from 52.6 years in 1960 to 72.6 years in 2018 [1]. While it is certain that people are living longer on average, this does not necessarily mean they are living healthier lives, as an increase in life expectancy anticipates an increase in morbidity [2,3]. As individuals age, the body changes and experiences a state of physical decline, resulting in weaker defenses and easier acquisition of chronic illnesses in the later years of life [2,4]. The diagnosis of two or more chronic diseases in the same individual is referred to as multimorbidity [5].
Multimorbid individuals tend to be prescribed a high number of medications in order to combat their diagnosed chronic illnesses. Consumption of prescribed drugs holds a higher prevalence and relevance in older adults, and complications could include potentially inappropriate prescribing [6]. While a homogeneous operational definitional is lacking throughout the field, literature supports the definition of polypharmacy as the consumption of five or more drugs daily in the same individual [7]. Polypharmacy is considered a critical public health problem that is related to drug-drug and drug-disease interactions, adverse drug events [8,9,10,11,12], falls, hospital admissions and mortality [13,14]. Polypharmacy has been on the rise over the past several decades [11] and is highly associated to multimorbidity [8].
In a world with an ageing population, the burdens of multimorbidity and polypharmacy have undue individual and system-wide impacts on health. While there exists a growing amount of literature regarding multimorbidity and polypharmacy, the vast majority of studies analyze polypharmacy as descriptive drugs in multimorbidity patterns, focus almost exclusively on one topic or the other without meaningfully connecting the two, or examine the disease rather than the individual as the unit of analysis [15,16,17]. Furthermore, medication is considered a proxy variable to disease [15,18,19], and, for this reason, jointly analyzing multimorbidity and polypharmacy can produce an overestimation error due to the fact that people with prevalent diseases such as diabetes or cardiovascular diseases are treated with many medications for both clinical conditions and risk factors and, for this reason, are overestimated. To avoid this, drug groups can be analyzed according to their associated disease, thereby preventing prevalent diseases from being overestimated. This type of approach would permit a better understanding of the patient groups and, at the same time, facilitate strategies aimed at prevention, diagnosis, and treatment because it includes diseases with or without drug treatment.
As far as we understand, little research has been completed regarding methods that simultaneously analyze the combined patterns of multimorbidity and polypharmacy at an individual level. Machine-learning soft clustering models are a robust tool capable of performing such an analysis. Cluster analysis involves assigning individuals to a certain cluster so that the items (i.e., units of analysis—diseases and drugs) are as similar as possible, while individuals in different clusters are as least similar as possible. Cluster identification is based on similarity measures, and their choice is reliant upon the data and/or the reason for analysis [20]. Hard clustering forces each individual to belong to only one cluster, while soft clustering (also called fuzzy clustering) grants varying degrees of membership, thus allowing for the individual to pertain to multiple clusters [20]. The aim of this study was to determine combined patterns of multimorbidity and polypharmacy in the Catalan population 65–99 years of age through a machine-learning soft clustering technique that incorporates the research team’s mapping of chronic disease and drug associations.

2. Materials and Methods

2.1. Setting, Design, and Population

Catalonia, an autonomous community of Spain, is a Mediterranean region with 7,515,398 reported inhabitants for the year 2012 [21]. Universal health coverage is established for residents in Spain by the National Health Service and is implemented in a decentralized fashion through each of the seventeen autonomous communities [22]. In Catalonia, the Catalan Health Institute (CHI) manages over 283 primary care centers, offering health services to over six million residents [23].
A cross-sectional study was performed on the baseline year (2012) of a longitudinal study (2012–2016) using electronic health records (EHRs) in Catalonia. Inclusion criteria for the cross-sectional study population allowed for individuals 65–99 years of age on 31 December 2011, who survived until 31 December 2012, and had at least one visit to a CHI-managed primary care center during the longitudinal study period (2012–2016). No new entries were permitted in the study, and dropouts were due to either death or transfer to another primary care center outside of CHI governance. A total of 916,619 eligible individuals were included at the baseline year (Figure 1).

2.2. Data Source

The Information System for Research in Primary Care (SIDIAP) contains EHRs from the primary care centers managed by the CHI [24]. The SIDIAP database, in addition to clinical information, contains demographic, laboratory, and invoiced drug information, with every datapoint linked to the individual via an anonymous and unique personal identifier.

2.3. Variables

The SIDIAP database was the single source of information for all variables. All variables were analyzed exclusively within the study period (2012).

2.3.1. Chronic Diseases and Multimorbidity

All diseases in the SIDIAP database were coded according to the International Classification of Diseases, Version 10 (ICD-10). An operational definition for multimorbidity was based on the 60 chronic disease categories determined by Calderón-Larrañaga et al. in the Swedish National study of Aging and Care in Kungsholmen (SNAC-K) [25]. Each chronic disease category was included as an individual binary variable, and multimorbidity was defined via a dichotomous variable as the presence of two or more diagnoses from the 60 chronic disease categories. However, only chronic disease categories with ≥2% prevalence in the study population were included for final analysis, thus leaving 47 SNAC-K chronic disease categories in total (Appendix A Table A1).

2.3.2. Drugs and Classification

Invoiced drugs recorded in the SIDIAP database included drugs dispensed in pharmacies. Drugs received in hospital and/or dispensed by a hospital pharmacy and all other drugs not subsidized through the national health system were excluded from this study. Drugs were coded according to the Anatomical Therapeutic Chemical (ATC) Classification System, which categorizes drugs through various levels of specificity into groups (hereafter referred to as “drug categories”) according to the targeted organ/system and their chemical, pharmacological, and therapeutic properties [26]. Drug categories with ≥1% prevalence in the study population were included for final analysis. Chronic use for invoiced drugs was determined for individuals who were invoiced three or more packages of the same drug categories during the study period. While not meeting the prevalence or chronic use criteria, the drug category Other Drugs Affecting Bone Structure and Mineralization (ATC 4th level code “M05BX”) was added to the study due to its twice-a-year treatment regimen for chronic diseases. Each drug category was included as an individual binary variable, and polypharmacy was defined via a dichotomous variable as chronic use in the same individual for five or more different drug categories (ATC 4th level) from the eighty-nine drug categories outlined in the following section.

2.3.3. Grouping of Drugs and Mapping to Chronic Diseases

The research team identified 89 different drug categories (ATC 4th level) (Appendix A Table A2) associated to the 60 chronic disease categories mentioned prior via a thorough revision of several databases that are well-known international clinical guidelines for each disease [27,28,29]. A mapping was done of all 60 SNAC-K chronic disease categories [25] and 89 drug categories. Drug categories were mapped to the SNAC-K chronic disease categories for which they are prescribed to treat. Chronic disease-drug categories were then created in the form of dichotomous variables for individuals who were diagnosed with a SNAC-K chronic disease category and, depending on disease management criteria, invoiced 0, 1, or 1+ of the mapped drug categories. (For example, the Allergy disease category was mapped with seven drug categories. An individual would qualify in this category if diagnosed with a disease pertaining to the Allergy disease category and invoiced at least one of the seven drug categories). A total of seven of the final 47 SNAC-K chronic diseases categories require non-pharmacological treatments or treatment with drugs excluded from this study and, therefore, could not be mapped. The remaining 40 SNAC-K chronic disease categories were mapped to drug categories, resulting in 29 chronic disease-drug categories containing ≥ 2% prevalence in the study population (see Appendix A Table A3 for mapping example; see Supplementary Materials Table S1 for complete mapping process). The seven chronic disease categories requiring non-pharmacological treatment or treatment with drugs excluded from this study and the 29 chronic disease-drug categories, all with prevalence ≥ 2%, were included to determine combined patterns of multimorbidity and polypharmacy (Figure 2).

2.3.4. Other Variables

Pertinent demographic data analyzed in the study includes age (measured in years), sex (female or male), and socioeconomic status (measured by the MEDEA [Mortality in Spanish Small Areas and Socioeconomic and Environmental Inequalities] Index via quintiles from least deprived to most deprived for urban areas, while rural areas were sorted into an independent category [30]).

2.4. Statistical Analysis

Descriptive statistics were applied to summarize preliminary findings.
Due to the high dimensionality of the SIDIAP database, dimension-reduction techniques were exercised through PCA Mix, an application of principal component analysis (PCA) for numeric original variables and multiple correspondence analysis (MCA) for binary variables. This method reduced the size of the database while maintaining the complexity of the original data. The Karlis–Saporta–Spinaki rule was applied in order to select the appropriate number of dimensions to preserve [31].
Using the reduced database, the combined multimorbidity and polypharmacy patterns were determined through a fuzzy c-means (FCM) clustering algorithm [32], incorporating the twenty-nine chronic disease-drug categories and the seven chronic disease categories with non-pharmacological treatment or treatment with drugs excluded from this study, all of which satisfied the prevalence threshold of ≥2% in the study population. To obtain a range with the ideal number of clusters, validation indices [33] were calculated (Supplementary Materials Figure S1). The outcome of the FCM clustering algorithm was a determined number of models with different numbers of clusters in each model, as calculated by the validation indices. Each model contained varying degrees of disease-medication association, and the final model of clusters was determined by the research team according to clinical relevance.
The clusters were described in two parts: (1) observed/expected ratios (O/E ratios) were calculated by dividing the prevalence of a chronic disease or chronic disease-drug category in a specific cluster by the prevalence of the same chronic disease or chronic disease-drug category in the entire study population; (2) exclusivity was determined by dividing the number of individuals with a chronic disease or chronic disease-drug category in a specific cluster by the amount of all individuals with the same chronic disease or chronic disease-drug category in the entire study population. A threshold of two for the O/E ratio was set in order for a disease/medication to be considered a relevant part of a cluster [34,35]. An exclusivity threshold of 30% was a secondary, but not determining, factor when evaluating the chronic disease or chronic disease-drug categories association with a cluster.
All analyses were performed in R version 4.0.3 and Stata version 15. Specifically, R was used to run the PCA mix and FCM clustering algorithm; Stata was used for data management.

3. Results

Of the 916,619 eligible individuals 65 years and over (women: 57.8%; mean age: 75.4; standard deviation; 7.4), 853,085 (93.1%) satisfied the criteria for multimorbidity, and 457,576 (49.9%) for polypharmacy (Figure 3). The most frequent chronic disease categories in the population were hypertension (71%), dyslipidemia (50.9%), osteoarthritis and other degenerative joint diseases (32.8%), obesity (28.7%), and diabetes (25.1%) [Appendix A Table A1], with a median of six chronic diseases (interquartile range [IQR] 4.0–8.0) per person. The most prevalent drug categories included proton pump inhibitors (44.3%), HMG CoA reductase inhibitors (38.2%), anilides (28.4%), platelet aggregation inhibitors, excluding heparin (35.6%), and benzodiazepine derivatives (20.9%) [Appendix A Table A2], with a median of 5 drugs (IQR 2.0–8.0) per person.
The authors identified a seven-cluster solution for combined patterns of multimorbidity and polypharmacy. Cluster 2 to Cluster 7 contained over 99% multimorbidity in each cluster and reported higher overrepresentation values (O/E ratio > 2) for at least one of the chronic disease or chronic disease-drug categories. Characteristics of the study participants in each cluster are detailed in Appendix A Table A4. Principal results and the most frequent chronic disease and chronic disease-drug categories by cluster (Table 1) are highlighted below:
Cluster 1 (non-specific) included a substantial number of individuals that do not present any overrepresented chronic disease or chronic disease-drug category (O/E ratios are below two and exclusivity values are below 30%). Cluster 1 was also the cluster with the lowest average age (74.20 years, SD 7.47) and the lowest percentage of individuals with multimorbidity (81.74%) and polypharmacy (15.38%).
Diabetes (Cluster 2): The only category that surpassed the O/E ratio threshold was chronic disease-drug group for “diabetes” (O/E ratio 2.15), with exclusivity of 41.93%.
Neurological and musculoskeletal, female dominant (Cluster 3): 69.05% of the individuals in this cluster were female, with exclusivity for the chronic disease-drug groups for “peripheral neuropathy”, “dorsopathies”, and “other musculoskeletal and joint diseases” all surpassing 30%. Similar to Cluster 3 is behavioral, neurological, and musculoskeletal, female dominant (Cluster 4), which also has a high proportion of females (75.08%) and contains many of the chronic disease-drug groups from the neurological and musculoskeletal categories. However, the differentiating factor is the emphasis on behavioral chronic disease-drug groups, with several groups pertaining to this category as well.
Cardio-cerebrovascular and renal (Cluster 5): The vast majority (six out of eight) of the significant groups in this cluster are cardiovascular related, with exclusivity of chronic disease-drug groups for “peripheral vascular disease” and “ischemic heart disease” both exceeding 30%. This cluster contained the highest proportion of polypharmacy (83.65%).
Cardiovascular, renal, inflammatory, and respiratory (Cluster 6): Cardiovascular categories occupy the top three spots, with chronic disease-drug categories for “atrial fibrillation” (O/E ratio 5.99) and “heart failure” (O/E ratio 5.70), and chronic disease category “bradycardias and conduction disorders” (O/E ratio 4.14), all exceeding an exclusivity of 30%. This cluster contained the highest proportions for individuals diagnosed with 10 or more chronic diseases (30.82%), individuals prescribed 10 or more chronic drugs (27.71%), and individuals with 10 or more primary care visits (75.38%). Polypharmacy was nearly that of Cluster 5 (83.54%).
The multisystem pattern identified in Cluster 7 represents the smallest group of individuals from the study population (5.41%) and contains several overrepresented diseases corresponding to multiple systems, in which individuals with digestive disease dominate the cluster.
Socioeconomic status (measured by MEDEA index) between categories appears to remain relatively stable in Cluster 2 to Cluster 7, with most deprived fluctuating between 13.63–15.79% and least deprived between 14.90–16.57%. Those in rural health settings fluctuate between 19.34–23.83% in Cluster 2 to Cluster 7.

4. Discussion

4.1. Key Results

In this article, patterns for multimorbidity and polypharmacy were analyzed in a joint manner by means of creating a variable relating each drug to one or more diseases for which they are indicated. Such patterns were obtained via specified criteria that allowed for the identification of multimorbid individuals that also qualified as polymedicated (those who were invoiced more than five distinct drugs in one year). This study ultimately identified distinct patterns with singular characteristics, with some more profound in women than in men. This is a key point in order to carry out a stricter follow up of the individuals who have a higher risk of presenting secondary effects caused by drugs or drug interaction.
Cluster 1 (non-specific) did not overrepresent any disease and Cluster 7 (multisystemic) overrepresented diseases from many systems are less specific than Clusters 2 to 6. These cluster types have been identified in other studies [35,36]. In our study, we were able to identify overrepresented diseases that do not have any prescribed drug in primary healthcare, such as neoplastic disease, deafness and hearing loss, and bradycardias and conduction diseases. These data suggest that diseases and drugs should be considered in combination to more efficiently identify multimorbidity patterns [37].
It was observed that the chronic disease-drug categories “hypertension” and “dyslipidemia” are distributed throughout the seven obtained clusters. Individuals with these common pathologies are therefore not grouped together in any one specific cluster; rather, they are spread throughout all seven. The fuzzy c-means method is a technique that classifies individuals based on cluster probability membership, and, given that there were many individuals with these chronic diseases and associated drugs in the study, the results are homogenously distributed throughout the study population [37].
The patterns from Cluster 2 to Cluster 6 group individuals who share similar chronic disease-drug categories. Cluster 2 is classified as diabetes because 41.93% of individuals included in the cluster are diabetic. Although not reaching the O/E threshold of two, associations are still observed between diabetes and obesity, cataracts, and neoplasms, all diseases that are frequently presented in diabetic individuals [38,39,40].
Cluster 3 (neurological and musculoskeletal) and Cluster 4 (behavioral, neurological, and musculoskeletal) include patterns predominantly in women. The clustering model identified these two groups of individuals that present similar health problems with singular characteristics. The patterns within behavioral, neurological, and musculoskeletal show an association with autoimmune diseases that are also more frequent in women [41].
Cluster 5 (cardio-cerebrovascular and renal) and Cluster 6 (cardiovascular, renal, inflammatory, and respiratory) include chronic disease-drug patterns that primarily treat cardiac, cardiovascular, cerebrovascular, and renal pathologies. This association, which has been observed in a previous study that indicated a pattern of increased mortality [16], concerns closely related pathologies that also share risk factors and even treatments. This study allows for the grouping of individuals with similar chronic disease and chronic disease-drug categories into two differentiated and specific clusters, given that Cluster 6 (cardiovascular, renal, inflammatory, and respiratory) includes a chronic disease-drug category for inflammatory pathologies as well as a chronic disease-drug group for respiratory diseases. Chronic inflammation as a mechanism in atherosclerosis carries a higher risk for cardiovascular and cerebrovascular events [42]. Autoimmune diseases and arthrosis have also been linked to an increase in cardiovascular disease, although this could be due to adverse effects of the treatments for these diseases, such as corticoids [43,44].
A small group of individuals who presented polypharmacy did not satisfy the condition for multimorbidity. This could be due to the individuals having received multiple treatments during the study period without them coinciding. On the other hand, there is also the possibility that some pathologies could be treated with five or more drugs simultaneously for certain patients, with neuropathy being a prime example.
The inclusion of the same study variable for chronic diseases and drugs indicated for their respective diseases’ treatment is an advantage that allows for a better identification of the individuals and, thereby, allows to better orient the clinical management of these groups of individuals.
The methodology used for mapping drugs and chronic diseases permitted the identification of a new chronic diseases-drug category. This new variable applied to fuzzy clustering methods is less susceptible to outliers in the data, choice of distance measure, and the inclusion of inappropriate or irrelevant variables [35,45,46].
This mapping process between chronic diseases and drugs, together with the application of the fuzzy method, will be very useful for multimorbidity studies where drugs are used as a proxy variable.
Soft clustering methods offer a new methodological approach towards understanding the relationships between specific diseases or drugs in individuals. This is an essential step in improving the care of patients and the quality of health systems. Analyzing multimorbidity patterns based on drugs permits the identification of patient subgroups with different clinical approaches and attention. Our analysis focuses on groups of patients with specific diseases and drugs as opposed to other studies centered in diseases.
This methodology can be applied in a large variety of studies using electronic health records to study polypharmacy and multimorbidity.

4.2. Comparison to the Literature

Studies performed with clustering methods concerning multimorbidity and polypharmacy are novel but becoming more frequent [47]; however, these studies generally approach multimorbidity with a polypharmacy perspective or vice versa. To the authors’ knowledge, this study is the first in conveying a combined approach to multimorbidity and polypharmacy by combining both topics into one inclusive variable to determine joint outcomes. Although, as previously mentioned, some of the patterns obtained (neurological and musculoskeletal, female dominant; behavioral, neurological, and musculoskeletal, female dominant; and cardio-cerebrovascular and renal) coincide with the literature [17,37,48].

4.3. Strengths and Limitations

The sample size, with nearly one million individuals, is a clear advantage of the study. Considering individuals as the primary variable for analysis rather than diseases [34], the study employs a large, high-quality database composed of primary health care records representative of the Catalan population aged ≥ 65 years [24], an age group that is more susceptible to present health problems related to drugs. The inclusion of all drugs that are indicated for different diseases with a minimum number of packages facilitates the identification of those individuals who should be more thoroughly monitored. Regarding the method, soft clustering offers a methodologic focus in understanding the relations between specific diseases and individuals. The analysis of multimorbidity and polypharmacy patterns can identify subgroups of patients with different associated diseases and drugs. The extensive mapping process of drug categories to their respective SNAC-K chronic disease categories provides a further level of detail than most other clustering techniques within the literature. This can ultimately serve as a framework for future studies that wish to map diseases and drugs under certain conditions to determine combined patterns.
Some limitations of this study should be considered. Individuals who met initial selection criteria but sought care outside of CHI governance were ineligible, possibly introducing selection bias for individuals who chose to seek care in a private healthcare facility. However, this is a small group of the population, and the results of this study can be applied to the general population. Drugs were recorded via invoices, and hospital drugs were not included, which could have influenced the calculation for polypharmacy. Concerning the clusters, the results obtained in this study are similar to those obtained in other studies without prior mapping; nevertheless, we consider that this methodology more precisely determines which pathologies and drugs are overrepresented in each cluster, thus more adequately defining patient profiles. Finally, while the fuzzy c-means clustering technique used in this study is an unsupervised, exploratory method, the authors believe that, combined with a vigorous internal validation system, this technique produces robust results and minimizes potential pitfalls.

5. Conclusions

This study utilized a mapping process combined with a soft clustering technique to determine combined patterns of multimorbidity and polypharmacy in the elderly Catalan population in 2012, identifying overrepresentation in 6 of the 7 clusters with chronic disease and chronic disease-drug categories. Cluster 2 to Cluster 6 provided recognizable patterns, predominantly in diabetes; neurological and musculoskeletal, female dominant and behavioral, neurological, and musculoskeletal, female dominant; and cardio-cerebrovascular and renal and cardiovascular, renal, inflammatory, and respiratory. These patterns further highlight the differences between sexes, specifically within neurological and musculoskeletal, female dominant and behavioral, neurological, and musculoskeletal, female dominant.
The combined patterns of multimorbidity and polypharmacy identified in this study will contribute key information to the evaluation of multimorbid and polymedicated individuals, facilitating the identification of these subgroups of individuals that require specific attention. The obtained results could be applied to clinical practice guidelines differentiating distinct population groups, such as multimorbid individuals with and without associated diseases and/or polypharmacy. Due to the relationship between multimorbidity and polypharmacy over long periods of time, this analysis could serve as the base to deepen this relationship in further longitudinal studies. The patterns obtained in this research will allow for an in-depth study of the prescription of multiple medications in elderly people in relation to medication-related problems.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijerph18179216/s1, Figure S1: validation indices, Table S1: mapping of SNAC-K chronic disease categories and their associated ATC drugs categories.

Author Contributions

Conceptualization, A.R.-L. and C.V.; methodology, G.S., N.V., A.R.-L., A.T.-M., and C.V.; software, G.S. and A.R.-L.; validation, N.V. and A.T.-M.; formal analysis, G.S. and A.R.-L.; investigation, G.S., N.V., A.R.-L., A.T.-M., and C.V.; resources, A.R.-L. and C.V.; data curation, G.S. and A.R.-L.; writing—original draft preparation, G.S., N.V., and C.V.; writing—review and editing, G.S., N.V., A.R.-L., A.T.-M., M.M., and C.V.; visualization, G.S. and A.R.-L.; supervision, A.R.-L. and C.V.; project administration, C.V.; funding acquisition, C.V. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a research grant from the Instituto de Salud Carlos III, under the auspices of the Ministry of Science, Innovation and Universities, awarded in the 2016 call under the Health Strategy Action 2013–2016, within the National Research Program oriented to Societal Challenges, within the Technical, Scientific and Innovation Research National Plan 2013–2016 (Grant Number PI16/00639), co-funded with the EU European Regional Development Fund and the Department of Health of the Catalan Government, in the 2017 call that awarded subsidies for the Strategic Plan for Research in Health (Pla Estratègic de Recerca i Innovació en Salut, PERIS) 2016–2020, modality research oriented to primary care (Grant Number SLT002/16/00058) and by the Catalan Government (Grant Number AGAUR 2017 SGR 578). This project has been partially funded by the Carlos III Health Institute (Ministry of Economy and Competitiveness, Spain) through the Network for Prevention and Health Promotion in Primary Care (redIAPP, RD16/0007/0001) and by European Union European Regional Development Fund (ERDF) funds (European Regional Development Fund).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Clinical Research Ethics Committee, Fundació Institut Universitari per a la Recerca a l’Atenció Primària de Salut Jordi Gol i Gurina (IDIAPJGol) (Protocol No: P17/080). All data were anonymized in agreement with national and international law.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets are not available, since researchers signed an agreement with the Information System for the Development of Research in Primary Care (SIDIAP) concerning confidentiality and security of the dataset, which forbids providing data to third parties. The SIDIAP is subject to periodic audits.

Acknowledgments

We thank the Catalan Institute of Health and the SIDIAP, which provided the database for the study.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Prevalence of the 60 chronic diseases included in the study in individuals aged 65–99 years (n = 916,619, Catalonia, 2012). In the last column, list of diseases included by prevalence cut off of 2%.
Table A1. Prevalence of the 60 chronic diseases included in the study in individuals aged 65–99 years (n = 916,619, Catalonia, 2012). In the last column, list of diseases included by prevalence cut off of 2%.
RankChronic ConditionsFrequencyPercentage (%)Cut Off of 2%
1Hypertension650,89971.0
2Dyslipidemia466,58550.9
3Osteoarthritis and other degenerative joint diseases300,80332.8
4Obesity262,88828.7
5Diabetes230,46025.1
6Anemia167,57718.3
7Cataract and other lens diseases156,62217.1
8Chronic kidney diseases153,75616.8
9Prostate diseases153,63516.8
10Osteoporosis151,84716.6
11Depression and mood diseases148,75116.2
12Solid neoplasms137,04515.0
13Colitis and related diseases131,51214.4
14Venous and lymphatic diseases126,99713.9
15Other musculoskeletal and joint diseases124,76513.6
16Dorsopathies124,60313.6
17Neurotic, stress-related and somatoform diseases123,39513.5
18COPD, emphysema, chronic bronchitis109,60312.0
19Ischemic heart disease95,43410.4
20Deafness, hearing impairment90,2619.9
21Sleep disorders88,7399.7
22Thyroid diseases88,4459.7
23Other genitourinary diseases85,4689.3
24Cerebrovascular disease80,2648.8
25Atrial fibrillation80,2478.8
26Esophagus, stomach and duodenum diseases80,0438.7
27Heart failure74,0778.1
28Other eye diseases68,9397.5
29Glaucoma66,1627.2
30Inflammatory arthropathies62,4506.8
31Dementia59,2136.5
32Cardiac valve diseases52,1005.7
33Peripheral neuropathy49,1275.4
34Other psychiatric and behavioral diseases46,8415.1
35Asthma43,6634.8
36Allergy40,3944.4
37Autoimmune diseases39,3504.3
38Ear, nose, throat diseases38,7524.2
39Peripheral vascular disease30,6743.4
40Other neurological diseases28,5413.1
41Chronic pancreas, biliary tract and gallbladder diseases27,3213.0
42Migraine and facial pain syndromes25,9992.8
43Bradycardias and conduction diseases25,4762.8
44Chronic liver diseases22,6332.5
45Other digestive diseases22,0222.4
46Parkinson and parkinsonism20,8332.3
47Other metabolic diseases 18,9972.1
48Other cardiovascular diseases16,8331.8
49Other skin diseases15,3631.7
50Chronic ulcer of the skin13,8691.5
51Blood and blood forming organ diseases 13,5751.5
52Other respiratory diseases99741.1
53Epilepsy89811.0
54Hematological neoplasms81740.9
55Chronic infectious diseases66470.7
56Inflammatory bowel diseases55490.6
57Schizophrenia and delusional diseases47920.5
58Blindness, visual impairment 47720.5
59Multiple sclerosis 5760.1
60Chromosomal abnormalities770.0
Abbreviations: COPD: chronic obstructive pulmonary disease. The black shading in the lower part of the table IS necessary, as it indicates a cut off of 2%, as defined in the first line.
Table A2. Prevalence of the 89 medication categories diseases included in the study in individuals aged 65–99 years (n = 916,619, Catalonia, 2012).
Table A2. Prevalence of the 89 medication categories diseases included in the study in individuals aged 65–99 years (n = 916,619, Catalonia, 2012).
RankATC-5 CodeMedication CategoriesFrequencyPercentage (%)
1A02BCProton pump inhibitors405,94244.29
2C10AAHMG CoA reductase inhibitors349,67638.15
3N02BEAnilides260,01828.37
4B01ACPlatelet aggregation inhibitors excl. heparin234,30625.56
5N05BABenzodiazepine derivatives191,87020.93
6C09AAACE inhibitors, plain182,90619.95
7A10BABiguanides119,95513.09
8C08CADihydropyridine derivatives116,71312.73
9N06ABSelective serotonin reuptake inhibitors116,49712.71
10C07ABBeta blocking agents, selective102,44211.18
11C03CASulfonamides, plain101,93311.12
12C09CAAngiotensin II antagonists, plain99,74210.88
13A12AXCalcium, combinations with vitamin D and/or other drugs93,58110.21
14C03AAThiazides, plain86,4149.43
15C09DAAngiotensin II antagonists and diuretics80,1158.74
16M01AEPropionic acid derivatives79,3938.66
17C09BAACE inhibitors and diuretics74,4948.13
18G04CAAlpha-adrenoreceptor antagonists71,0287.75
19B01AAVitamin K antagonists62,3316.80
20M05BABisphosphonates62,3136.80
21R03BBAnticholinergics53,1595.80
22N05CDBenzodiazepine derivatives52,8595.77
23A10BBSulfonylureas50,8125.54
24H03AAThyroid hormones50,6995.53
25R03AKAdrenergics in combination with corticosteroids or other drugs, excl. anticholinergics50,6445.53
26N02AXOther opioids48,3885.28
27M04AAPreparations inhibiting uric acid production47,6295.20
28N06AXOther antidepressants45,3574.95
29R03ACSelective beta-2-adrenoreceptor agonists43,9724.80
30S01EDBeta blocking agents42,5134.64
31C01DAOrganic nitrates41,6744.55
32S01EEProstaglandin analogues40,6434.43
33N03AXOther antiepileptics37,8734.13
34B03AAIron bivalent, oral preparations35,133.83
35M01AXOther anti-inflammatory and antirheumatic agents, non-steroids34,4903.76
36M01ABAcetic acid derivatives and related substances30,8483.37
37N02BBPyrazolones30,6003.34
38C02CAAlpha-adrenoreceptor antagonists26,3332.87
39N07CAAntivertigo preparations25,2902.76
40C01EBOther cardiac preparations24,2722.65
41G04BDDrugs for urinary frequency and incontinence23,6462.58
42D01ACImidazole and triazole derivatives23,4592.56
43A03FAPropulsives22,4612.45
44A10AEInsulins and analogues for injection, long-acting21,6782.36
45H02ABGlucocorticoids21,4932.34
46C01AADigitalis glycosides21,1902.31
47D07ACCorticosteroids, potent (Group III)20,5992.25
48G04CBTestosterone-5-alpha reductase inhibitors19,3792.11
49N06DAAnticholinesterases19,3662.11
50C08DBBenzothiazepine derivatives19,2112.10
51R03BAGlucocorticoids19,1432.09
52R06AXOther antihistamines for systemic use18,6092.03
53C07AGAlpha and beta blocking agents18,1981.99
54A02BAH2-receptor antagonists18,1581.98
55C03BASulfonamides, plain17,9981.96
56C10ABFibrates17,8741.95
57N05CFBenzodiazepine related drugs17,9141.95
58G03CANatural and semisynthetic estrogens, plain16,5131.80
59C03DAAldosterone antagonists16,3671.79
60N06DXOther anti-dementia drugs16,2921.78
61A10BDCombinations of oral blood glucose lowering drugs15,3971.68
62N06AANon-selective monoamine reuptake inhibitors15,3641.68
63N05AHDiazepines, oxazepines, thiazepines and oxepines14,2491.55
64R01ADCorticosteroids13,8801.51
65G04CXOther drugs used in benign prostatic hypertrophy13,6841.49
66J01MAFluoroquinolones13,6101.48
67C03EALow-ceiling diuretics and potassium-sparing agents13,1791.44
68A10BXOther blood glucose lowering drugs, excl. insulins12,9331.41
69A12AACalcium12,8091.40
70C09DBAngiotensin II antagonists and calcium channel blockers12,6981.39
71A12BAPotassium12,6251.38
72J01CRCombinations of penicillins, incl. beta-lactamase inhibitors12,6581.38
73C04ADPurine derivatives12,5891.37
74C10AXOther lipid modifying agents12,5791.37
75M01AHCoxibs12,5371.37
76N04BADopa and dopa derivatives12,2801.34
77D01AEOther antifungals for topical use11,5141.26
78A11CCVitamin D and analogues11,2621.23
79A10ADInsulins and analogues for injection, intermediate- or long-acting combined with fast-acting11,0561.21
80A10ACInsulins and analogues for injection, intermediate-acting10,9251.19
81N02ABPhenylpiperidine derivatives10,8021.18
82B03BAVitamin B12 (cyanocobalamin and analogues)10,6131.16
83C01BDAntiarrhythmics, class III10,5111.15
84N05AXOther antipsychotics10,4041.14
85C07AABeta blocking agents, non-selective10,2311.12
86S01ECCarbonic anhydrase inhibitors97371.06
87B03ABIron trivalent, oral preparations95041.04
88N03AEBenzodiazepine derivatives92111.00
89M05BXOther drugs affecting bone structure and mineralization76520.83
Abbreviations: HMGCoA-reductase: 3-Hidroxi-3-metil-glutaril-CoA reductase; ACE inhibitors: angiotensin converting enzyme inhibitors.
Table A3. Mapping of SNAC-K chronic disease categories and their associated ATC drugs categories (first two categories).
Table A3. Mapping of SNAC-K chronic disease categories and their associated ATC drugs categories (first two categories).
ALLERGY
ICD-10 codes and labels included in chronic disease categoryATC-5 codes and labels for drugs associated to chronic disease category. Includes none or any of the following:
J301Allergic rhinitis due to pollenR01ADCorticosteroids (nasal use)
J302Other seasonal allergic rhinitisR03ACSelective beta-2-adrenoreceptor agonists
J303Other allergic rhinitisR03AKAdrenergics in combination with corticosteroids or other drugs, excl. anticholinergics
J304Allergic rhinitis, unspecifiedR03BAGlucocorticoids (inhalation)
J450Predominantly allergic asthmaR03BBAnticholinergics
K522Allergic and dietetic gastroenteritis and colitisR06AXOther antihistamines for systemic use
L20Atopic dermatitisD07ACCorticosteroids, potent (Group III) (topical use)
L23Allergic contact dermatitis
L500Allergic urticaria
Z516Desensitization to allergens
ANEMIA
ICD-10 codes and labels included in chronic disease categoryATC-5 codes and labels for drugs associated to chronic disease category. Includes none or any of the following:
D50Iron deficiency anaemiaB03AAIron bivalent, oral preparations
D51Vitamin B12 deficiency anaemiaB03ABIron trivalent, oral preparations
D52Folate deficiency anaemiaB03BAVitamin B12 (cyanocobalamin and analogues)
D53Other nutritional anaemiasH02ABGlucocorticoids (systemic use, plain)
D55Anaemia due to enzyme disorders
D56Thalassaemia
D57Sickle-cell disorders
D58Other hereditary haemolytic anaemias
D59Acquired haemolytic anaemia
D60Acquired pure red cell aplasia [erythroblastopenia]
D61Other aplastic anaemias
D63Anaemia in chronic diseases classified elsewhere
D64Other anaemias
Only chronic diseases and their associated ATC drug categories resulting in chronic disease or chronic disease-drug categories ≥ 2% prevalence are included in this table. For the sake of brevity, only the first two mapped categories are shown. The complete mapping table can be found in Supplementary Materials Table S1. Remaining categories with <2% prevalence can be sent to readers upon request.
Table A4. Variables characterizing each cluster for chronic disease or chronic disease-drug category threshold of 2% (n = 916,619).
Table A4. Variables characterizing each cluster for chronic disease or chronic disease-drug category threshold of 2% (n = 916,619).
1. Non-Specific2. Diabetes3. Neurological and Musculoskeletal, Female Dominant4. Behavioral, Neurological, and Musculoskeletal, Female Dominant5. Cardio- Cerebrovascular and Renal6. Cardiovascular, Renal, Inflammatory, and Respiratory7. MultisystemStudy Population (All)
Number of people, n344,958178,457102,75090,28780,85569,72049,592916,619
Multimorbidity, n (%)281,952 (81.74)178,412 (99.97)102,629 (99.88)90,245 (99.95)80,819 (99.96)69,680 (99.94)49,347 (99.51)853,084 (93.07)
Polypharmacy, n (%)53,070 (15.38)105,323 (59.02)69,503 (67.64)68,276 (75.62)67,635 (83.65)58,245 (83.54)35,523 (71.63)457,575 (49.92)
Women, n (%)187,691 (54.41)98,251 (55.06)70,946 (69.05)67,785 (75.08)39,649 (49.04)35,130 (50.39)29,679 (59.85)529,131 (57.73)
Men, n (%)157,267 (45.59)80,206 (44.94)31,803 (30.95)22,502 (24.92)41,207 (50.96)34,590 (49.61)19,913 (40.15)387,488 (42.27)
Age, mean (SD)74.28 (7.47)74.61 (6.86)75.82 (7.07)75.43 (7.06)77.06 (7.34)78.47 (7.30)78.27 (7.52)75.41 (7.39)
Age (categories), n (%)
[65,70)119,921 (34.76)51,741 (28.99)24,050 (23.41)23,121 (25.61)15,649 (19.35)9796 (14.05)7901 (15.93)252,178 (27.51)
[70,80)140,724 (40.79)83,142 (46.59)47,314 (46.05)41,082 (45.5)34,708 (42.93)28,197 (40.44)19,419 (39.16)394,586 (43.05)
[80,90)71,317 (20.67)39,345 (22.05)27,970 (27.22)23,398 (25.92)26,574 (32.87)27,204 (39.02)18,936 (38.18)234,744 (25.61)
[90,99]12,997 (3.77)4229 (2.37)3416 (3.32)2686 (2.97)3924 (4.85)4524 (6.49)3336 (6.73)35,111 (3.83)
MEDEA * index, n (%)
R71,007 (21.98)35,761 (21.22)18,739 (19.54)16,239 (19.34)15,426 (21.01)14,514 (23.28)10,563 (23.83)182,249 (19.88)
U163,526 (19.67)25,114 (14.90)14,580 (15.20)13,131 (15.64)11,684 (15.91)9413 (15.10)7344 (16.57)144,792 (15.80)
U253,337 (16.51)26,606 (15.79)15,239 (15.89)13,307 (15.85)11,546 (15.72)9563 (15.34)6832 (15.41)136,430 (14.88)
U351,420 (15.92)28,004 (16.62)15,988 (16.67)13,976 (16.64)11,928 (16.24)9839 (15.78)7068 (15.94)138,223 (15.08)
U446,598 (14.43)28,059 (16.65)16,264 (16.96)14,066 (16.75)11,795 (16.06)9726 (15.60)6480 (14.62)132,988 (14.51)
U537,128 (11.49)24,998 (14.83)15,109 (15.75)13,264 (15.79)11,050 (15.05)9290 (14.90)6043 (13.63)116,882 (12.75)
Number of chronic diseases, median [IQR]3 [2–5]6 [4–7]7 [5–9]8 [6–10]8 [6–10]8 [6–10]7 [6–10]6 [4–8]
Number of chronic diseases (categories), n (%)
025,380 (7.36)0 (0)0 (0)0 (0)0 (0)0 (0)0 (0)25,380 (2.77)
137,626 (10.91)45 (0.03)121 (0.12)42 (0.05)36 (0.04)40 (0.06)244 (0.49)38,154 (4.16)
[2,5)176,529 (51.17)50,912 (28.53)12,934 (12.59)7896 (8.75)7990 (9.88)6133 (8.8)6441 (12.99)268,836 (29.33)
[5,10)99,930 (28.97)114,840 (64.35)67,586 (65.78)58,154 (64.41)51,082 (63.18)42,058 (60.32)30,059 (60.61)463,709 (50.59)
≥105493 (1.59)12,660 (7.09)22,109 (21.52)24,195 (26.8)21,747 (26.9)21,489 (30.82)12,847 (25.91)120,540 (13.15)
Number of medications, median [IQR]2 [0–3]5 [3–7]6 [4–8]6 [5–9]7 [5–9]7 [5–10]6 [4–9]4 [2–7]
Number of medications, (categories), n (%)
0113,702 (32.96)536 (0.3)1414 (1.38)132 (0.15)24 (0.03)180 (0.26)1801 (3.63)117,789 (12.85)
150,694 (14.7)5232 (2.93)3057 (2.97)1282 (1.42)581 (0.72)602 (0.86)1364 (2.75)62,812 (6.85)
[2,5)127,492 (36.96)67,366 (37.75)28,777 (28.01)20,597 (22.81)12,615 (15.6)10,692 (15.34)10,904 (21.99)278,442 (30.38)
[5,10)50,510 (14.64)90,914 (50.94)55,102 (53.63)52,401 (58.04)47,710 (59.01)38,927 (55.83)25,494 (51.41)361,058 (39.39)
≥102560 (0.74)14,410 (8.07)14,401 (14.02)15,875 (17.58)19,924 (24.64)19,318 (27.71)10,029 (20.22)96,518 (10.53)
Number of visits 2012, median [IQR]6 [2–10]10 [6–16]12 [7–18]12 [7–19]14 [8–23]18 [10–29]13 [7–23]9 [5–16]
Number of visits 2012 (categories), n (%)
044,954 (13.03)825 (0.46)637 (0.62)317 (0.35)242 (0.3)204 (0.29)766 (1.54)47,945 (5.23)
125,389 (7.36)3155 (1.77)1463 (1.42)1213 (1.34)1026 (1.27)683 (0.98)954 (1.92)33,884 (3.7)
[2,5)75,895 (22)20,685 (11.59)9889 (9.62)8168 (9.05)6708 (8.3)4336 (6.22)4759 (9.6)130,439 (14.23)
[5,10)101,438 (29.41)56,518 (31.67)27,564 (26.83)22,805 (25.26)18,077 (22.36)11,938 (17.12)11,008 (22.2)249,349 (27.2)
≥1097,283 (28.2)97,274 (54.51)63,197 (61.51)57,783 (64)54,802 (67.78)52,558 (75.38)32,105 (64.74)455,002 (49.64)
For the sake of simplicity, all numbers in the table were rounded to its closest natural number. * MEDEA index starts with U1 (urban setting, least deprived) and ends with U5 (urban setting, most deprived). Individuals in rural settings are classified in the variable R. MEDEA index n = 851,564.

References

  1. The World Bank. Life Expectancy at Birth, Total (Years). Available online: https://data.worldbank.org/indicator/SP.DYN.LE00.IN (accessed on 10 March 2021).
  2. Partridge, L.; Deelen, J.; Slagboom, P.E. Facing up to the global challenges of ageing. Nature 2018, 561, 45–56. [Google Scholar] [CrossRef] [PubMed]
  3. Murray, C.J.L.; Barber, R.M.; Foreman, K.J.; Ozgoren, A.A.; Abd-Allah, F.; Abera, S.F.; Aboyans, V.; Abraham, J.P.; Abubakar, I.; Abu-Raddad, L.J.; et al. Global, regional, and national disability-adjusted life years (DALYs) for 306 diseases and injuries and healthy life expectancy (HALE) for 188 countries, 1990–2013: Quantifying the epidemiological transition. Lancet 2015, 386, 2145–2191. [Google Scholar] [CrossRef] [Green Version]
  4. Strom, B.L.; Kimmel, S.E.; Hennessy, S. (Eds.) Basic principles of clinical pharmacology relevant to pharmacoepidemiologic studies. In Pharmacoepidemiology, 5th ed.; John Wiley & Sons: Hoboken, NJ, USA, 2011; pp. 23–37. [Google Scholar]
  5. MacMahon, S.; Calverley, P.; Chaturvedi, N.; Chen, Z.; Corner, L.; Davies, M.; Ezzati, M.; Guthrie, B.; Hanson, K.; Jha, V.; et al. Multimorbidity: A Priority for Global Health Research; The Academy of Medical Sciences: London, UK, 2018; p. 6. [Google Scholar]
  6. Cahir, C.; Fahey, T.; Teeling, M.; Teljeur, C.; Feely, J.; Bennett, K. Potentially inappropriate prescribing and cost outcomes for older people: A national population study. Br. J. Clin. Pharmacol. 2010, 69, 543–552. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Masnoon, N.; Shakib, S.; Kalisch-Ellett, L.; Caughey, G.E.; Masnoon, N.; Shakib, S.; Kalisch-Ellett, L.; Caughey, G.E. What is polypharmacy? A systematic review of definitions. BMC Geriatr. 2017, 17, 1–10. [Google Scholar] [CrossRef] [Green Version]
  8. Khezrian, M.; McNeil, C.J.; Murray, A.; Myint, P.K. An overview of prevalence, determinants and health outcomes of polypharmacy. Ther. Adv. Drug Saf. 2020, 11. [Google Scholar] [CrossRef]
  9. Thorell, K.; Midlöv, P.; Fastbom, J.; Halling, A. Use of potentially inappropriate medication and polypharmacy in older adults: A repeated cross-sectional study. BMC Geriatr. 2020, 20, 1–9. [Google Scholar] [CrossRef] [Green Version]
  10. Maher, R.L.; Hanlon, J.; Hajjar, E.R. Clinical consequences of polypharmacy in elderly. Expert Opin. Drug Saf. 2013, 13, 57–65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Wastesson, J.W.; Morin, L.; Tan, E.; Johnell, K. An update on the clinical consequences of polypharmacy in older adults: A narrative review. Expert Opin. Drug Saf. 2018, 17, 1185–1196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Cullinan, S.; O’Mahony, D.; O’Sullivan, D.; Byrne, S. Use of a frailty index to identify potentially inappropriate prescribing and adverse drug reaction risks in older patients. Age Ageing 2015, 45, 115–120. [Google Scholar] [CrossRef] [Green Version]
  13. Leelakanok, N.; Holcombe, A.L.; Lund, B.C.; Gu, X.; Schweizer, M. Association between polypharmacy and death: A systematic review and meta-analysis. J. Am. Pharm. Assoc. 2017, 57, 729–738. [Google Scholar] [CrossRef]
  14. Gutiérrez-Valencia, M.; Herce, P.A.; Lacalle-Fabo, E.; Escámez, B.C.; Veloz, B.A.C.; Martínez-Velilla, N. Prevalencia de polifarmacia y factores asociados en adultos mayores en España: Datos de la Encuesta Nacional de Salud 2017. Med. Clin. 2019, 153, 141–150. [Google Scholar] [CrossRef]
  15. Forslund, T.; Carlsson, A.C.; Ljunggren, G.; Ärnlöv, J.; Wachtler, C. Patterns of multimorbidity and pharmacotherapy: A total population cross-sectional study. Fam. Pr. 2020, 38, 132–139. [Google Scholar] [CrossRef] [PubMed]
  16. Villén, N.; Guisado-Clavero, M.; Fernández-Bertolín, S.; Troncoso-Mariño, A.; Foguet-Boreu, Q.; Amado, E.; Pons-Vigués, M.; Roso-Llorach, A.; Violán, C. Multimorbidity patterns, polypharmacy and their association with liver and kidney abnormalities in people over 65 years of age: A longitudinal study. BMC Geriatr. 2020, 20, 1–13. [Google Scholar] [CrossRef]
  17. Prados-Torres, A.; Calderón-Larrañaga, A.; Hancco-Saavedra, J.; Poblador-Plou, B.; Akker, M.V.D. Multimorbidity patterns: A systematic review. J. Clin. Epidemiol. 2014, 67, 254–266. [Google Scholar] [CrossRef]
  18. McCarthy, C.; Moriarty, F.; Wallace, E.; Smith, S.M.; Barbara Clyne for the SPPiRE Study Team. The evolution of an evidence based intervention designed to improve prescribing and reduce polypharmacy in older people with multimorbidity and significant polypharmacy in primary care (SPPiRE). J. Comorb. 2020, 10, 32974211. [Google Scholar] [CrossRef]
  19. Chang, T.I.; Park, H.; Kim, D.W.; Jeon, E.K.; Rhee, C.M.; Kalantar-Zadeh, K.; Kang, E.W.; Kang, S.-W.; Han, S.H. Polypharmacy, hospitalization, and mortality risk: A nationwide cohort study. Sci. Rep. 2020, 10, 1–9. [Google Scholar] [CrossRef]
  20. Wolfram. Fuzzy Clustering. Available online: https://reference.wolfram.com/legacy/applications/fuzzylogic/Manual/12.html (accessed on 31 July 2021).
  21. Statistical Institute of Catalonia. Population on 1 January. Provinces. Available online: https://www.idescat.cat/pub/?id=aec&n=245&lang=en (accessed on 3 March 2021).
  22. Royal Decree-Law 16/2012, of April 20th, Regarding Urgent Measures to Guarantee the Sustainability of the National Health System and to Improve the Quality and Security of its Provisions. Available online: https://www.boe.es/eli/es/rdl/2012/04/20/16 (accessed on 4 March 2021).
  23. Catalan Health Institute. The CHI. Available online: http://ics.gencat.cat/es/lics/ (accessed on 4 March 2021).
  24. García-Gil, M.D.M.; Hermosilla, E.; Prieto-Alhambra, D.; Fina, F.; Rosell, M.; Ramos, R.; Rodriguez, J.; Williams, T.; van Staa, T.; Bolibar, B. Construction and validation of a scoring system for the selection of high-quality data in a Spanish population primary care database (SIDIAP). J. Innov. Health Inform. 2011, 19, 135–145. [Google Scholar] [CrossRef]
  25. Calderón-Larrañaga, A.; Vetrano, D.L.; Onder, G.; Feliu, L.A.G.; Coscollar-Santaliestra, C.; Carfí, A.; Pisciotta, M.S.; Angleman, S.; Melis, R.J.; Santoni, G.; et al. Assessing and measuring chronic multimorbidity in the older population: A proposal for its operationalization. J. Gerontol. Ser. A 2016, 72, 1417–1423. [Google Scholar] [CrossRef]
  26. World Health Organization. The Anatomical Therapeutic Chemical Classification System with Defined Daily Doses (ATC/DDD). Available online: https://www.who.int/classifications/atcddd/en/ (accessed on 9 March 2021).
  27. UpToDate. Available online: https://www.uptodate.com/contents/search (accessed on 15 March 2021).
  28. BOT PLUS. Available online: https://botplusweb.portalfarma.com/botplus.aspx (accessed on 2 March 2021).
  29. CIMA. Available online: https://cima.aemps.es/cima/publico/home.html (accessed on 3 March 2021).
  30. Borrell, C.; Cano-Serral, G.; Esnaola, S.; Nolasco, A.; Pasarín, M.I.; Ramis, R.; Saurina, C.; Escolar-Pujolar, A. Constructing a deprivation index based on census data in large Spanish Cities. Gac. Sanit. 2008, 22, 179–187. [Google Scholar]
  31. Karlis, D.; Saporta, G.; Spinakis, A. A simple rule for the selection of principal components. Commun. Stat. Theory Methods 2003, 32, 643–666. [Google Scholar] [CrossRef]
  32. Bezdek, J.C.; Ehrlich, R.; Full, W. FCM: The fuzzy c-means clustering algorithm. Comput. Geosci. 1984, 10, 191–203. [Google Scholar] [CrossRef]
  33. Zhao, Q. Cluster Validity in Clustering Methods. Ph.D. Thesis, University of Eastern Finland, Joensuu, Finland, 25 June 2012. [Google Scholar]
  34. Violán, C.; Roso-Llorach, A.; Foguet-Boreu, Q.; Guisado-Clavero, M.; Pons-Vigués, M.; Pujol-Ribera, E.; Valderas, J.M. Multimorbidity patterns with K-means nonhierarchical cluster analysis. BMC Fam. Pr. 2018, 19, 1–11. [Google Scholar] [CrossRef]
  35. Violán, C.; Foguet-Boreu, Q.; Fernández-Bertolín, S.; Guisado-Clavero, M.; Cabrera-Bean, M.; Formiga, F.; Valderas, J.M.; Roso-Llorach, A. Soft clustering using real-world data for the identification of multimorbidity patterns in an elderly population: Cross-sectional study in a Mediterranean population. BMJ Open 2019, 9, e029594. [Google Scholar] [CrossRef] [Green Version]
  36. Vetrano, D.L.; Roso-Llorach, A.; Fernández, S.; Guisado-Clavero, M.; Violán, C.; Onder, G.; Fratiglioni, L.; Calderón-Larrañaga, A.; Marengoni, A. Twelve-year clinical trajectories of multimorbidity in a population of older adults. Nat. Commun. 2020, 11, 1–9. [Google Scholar] [CrossRef]
  37. Violán, C.; Foguet-Boreu, Q.; Flores-Mateo, G.; Salisbury, C.; Blom, J.; Freitag, M.; Glynn, L.; Muth, C.; Valderas, J.M. Prevalence, determinants and patterns of multimorbidity in primary care: A systematic review of observational studies. PLoS ONE 2014, 9, e102149. [Google Scholar] [CrossRef]
  38. Chobot, A.; Górowska-Kowolik, K.; Sokołowska, M.; Jarosz-Chobot, P. Obesity and diabetes—Not only a simple link between two epidemics. Diabetes Metab. Res. Rev. 2018, 34, e3042. [Google Scholar] [CrossRef] [Green Version]
  39. Becker, C.; Schneider, C.; Aballéa, S.; Bailey, C.; Bourne, R.; Jick, S.; Meier, C. Cataract in patients with diabetes mellitus—Incidence rates in the UK and risk factors. Eye 2018, 32, 1028–1035. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Giovannucci, E.; Harlan, D.M.; Archer, M.C.; Bergenstal, R.M.; Gapstur, S.M.; Habel, L.; Pollak, M.; Regensteiner, J.G.; Yee, D. Diabetes and cancer: A consensus report. Diabetes Care 2010, 33, 1674–1685. [Google Scholar] [CrossRef] [Green Version]
  41. Angum, F.; Khan, T.; Kaler, J.; Siddiqui, L.; Hussain, A. The prevalence of autoimmune disorders in women: A narrative review. Cureus 2020, 12, e8094. [Google Scholar] [CrossRef] [PubMed]
  42. Chen, W.; Thomas, J.; Sadatsafavi, M.; FitzGerald, J.M. Risk of cardiovascular comorbidity in patients with chronic obstructive pulmonary disease: A systematic review and meta-analysis. Lancet Respir. Med. 2015, 3, 631–639. [Google Scholar] [CrossRef]
  43. Iseme, R.; McEvoy, M.; Kelly, B.; Agnew, L.; Walker, F.R.; Handley, T.; Oldmeadow, C.; Attia, J.; Boyle, M. A role for autoantibodies in atherogenesis. Cardiovasc. Res. 2017, 113, 1102–1112. [Google Scholar] [CrossRef] [PubMed]
  44. Mathieu, S.; Couderc, M.; Tournadre, A.; Soubrier, M. Cardiovascular profile in osteoarthritis: A meta-analysis of cardiovascular events and risk factors. Jt. Bone Spine 2019, 86, 679–684. [Google Scholar] [CrossRef] [PubMed]
  45. Bora, D.J.; Gupta, A.K. A Comparative study between fuzzy clustering algorithm and hard clustering algorithm. Int. J. Comput. Trends Technol. 2014, 10, 108–113. [Google Scholar] [CrossRef] [Green Version]
  46. Everitt, B.; Landau, S.; Leese, M. Cluster Analysis, 5th ed.; Balding, D.J., Cressie, N.A., Fitzmaurice, G.M., Goldstein, H., Eds.; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  47. Muth, C.; Blom, J.W.; Smith, S.; Johnell, K.; Gonzalez-Gonzalez, A.I.; Nguyen, T.S.; Brueckle, M.; Cesari, M.; Tinetti, M.E.; Valderas, J.M. Evidence supporting the best clinical management of patients with multimorbidity and polypharmacy: A systematic guideline review and expert consensus. J. Intern. Med. 2018, 285, 272–288. [Google Scholar] [CrossRef] [Green Version]
  48. Ng, S.K.; Tawiah, R.; Sawyer, M.; Scuffham, P. Patterns of multimorbid health conditions: A systematic review of analytical methods and comparison analysis. Int. J. Epidemiol. 2018, 47, 1687–1704. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Estimated population study according to selection criteria.
Figure 1. Estimated population study according to selection criteria.
Ijerph 18 09216 g001
Figure 2. Number of significant chronic disease or chronic disease-drug categories after mapping.
Figure 2. Number of significant chronic disease or chronic disease-drug categories after mapping.
Ijerph 18 09216 g002
Figure 3. Multimorbid and polymedicated individuals in the study aged 65–99 years (n = 916,619, Catalonia, 2012).
Figure 3. Multimorbid and polymedicated individuals in the study aged 65–99 years (n = 916,619, Catalonia, 2012).
Ijerph 18 09216 g003
Table 1. Most frequent 15 chronic disease or chronic disease-drug categories in individuals aged 65–99 years by cluster (n = 916,619, Catalonia, 2012).
Table 1. Most frequent 15 chronic disease or chronic disease-drug categories in individuals aged 65–99 years by cluster (n = 916,619, Catalonia, 2012).
PatternDisease or Disease-Medication CategoryOO/E RatioEX
1
Non-Specific
(n = 344,958: 37.63%)
Chronic disease group for solid neoplasms12.860.8632.37
Chronic disease-drug group for prostate diseases7.820.7729.11
Chronic disease-drug group for osteoporosis7.740.7427.74
Chronic disease group for deafness and hearing loss5.910.6022.58
Chronic disease-drug group for COPD, emphysema, and chronic bronchitis4.260.5320.03
Chronic disease-drug group for esophagus, stomach, and duodenum diseases3.770.5219.54
Chronic disease-drug group for thyroid diseases2.840.5119.34
Chronic disease group for cataract and lens diseases8.470.5018.65
Chronic disease-drug group for dementia1.560.4717.61
Chronic disease-drug group for hypertension26.410.4617.48
Chronic disease group for bradycardias and conduction diseases1.120.4015.13
Chronic disease-drug group for sleep disorders2.220.3914.64
Chronic disease group for obesity11.290.3914.82
Chronic disease-drug group for dyslipidemia12.320.3814.40
Chronic disease group for chronic pancreas, biliary tract, and gallbladder diseases1.090.3713.78
2
Diabetes
(n = 178,457: 19.47%)
Chronic disease-drug group for diabetes39.522.1541.93
Chronic disease-drug group for glaucoma10.751.7834.65
Chronic disease group for obesity49.611.7333.68
Chronic disease-drug group for dyslipidemia55.121.7133.33
Chronic disease-drug group for hypertension84.371.4828.89
Chronic disease-drug group for thyroid diseases7.441.3426.17
Chronic disease-drug group for chronic kidney diseases13.921.3225.69
Chronic disease-drug group for ischemic heart disease8.521.1121.61
Chronic disease-drug group for cerebrovascular diseases6.541.0119.69
Chronic disease group for cataract and lens diseases17.151.0019.54
Chronic disease-drug group for peripheral vascular disease2.551.0019.47
Chronic disease-drug group for prostate diseases9.580.9518.46
Chronic disease group for solid neoplasms14.060.9418.31
Chronic disease group for deafness and hearing loss8.420.8516.64
Chronic disease group for chronic pancreas, biliary tract, and gallbladder diseases2.400.8015.67
3
Neurological and Musculoskeletal, Female Dominant
(n = 102,750: 11.21%)
Chronic disease-drug group for peripheral neuropathy9.733.0834.56
Chronic disease-drug group for dorsopathies24.052.8832.24
Chronic disease-drug group for other musculoskeletal and joint diseases22.422.7330.64
Chronic disease-drug group for other genitourinary diseases8.682.4427.32
Chronic disease-drug group for glaucoma14.682.4327.25
Chronic disease-drug group for osteoarthritis and other degenerative joint diseases41.872.1624.19
Chronic disease group for deafness and hearing loss19.682.0022.40
Chronic disease-drug group for neurotic, stress-related, and somatoform diseases20.381.9922.36
Chronic disease group for cataract and lens diseases33.751.9822.14
Chronic disease-drug group for depression and mood diseases23.161.8620.90
Chronic disease-drug group for osteoporosis18.911.8020.17
Chronic disease-drug group for colitis and related diseases18.221.7920.09
Chronic disease-drug group for sleep disorders9.931.7419.51
Chronic disease-drug group for other psychiatric and behavioral diseases3.181.5817.71
Chronic disease-drug group for esophagus, stomach, and duodenum diseases11.291.5517.41
4
Behavioral, Neurological, and Musculoskeletal, Female Dominant
(n = 90,287: 9.85%)
Chronic disease-drug group for other psychiatric and behavioral diseases7.423.6936.34
Chronic disease-drug group for neurotic, stress-related, and somatoform diseases37.333.6536.00
Chronic disease-drug group for peripheral neuropathy10.933.4634.09
Chronic disease-drug group for depression and mood diseases41.753.3633.10
Chronic disease-drug group for dorsopathies27.223.2532.06
Chronic disease-drug group for other musculoskeletal and joint diseases26.653.2532.00
Chronic disease-drug group for other genitourinary diseases9.812.7627.14
Chronic disease-drug group for sleep disorders15.072.6426.02
Chronic disease-drug group for osteoarthritis and other degenerative joint diseases47.212.4323.97
Chronic disease-drug group for colitis and related diseases22.462.2121.76
Chronic disease-drug group for osteoporosis20.891.9919.58
Chronic disease-drug group for esophagus, stomach, and duodenum diseases13.801.9018.70
Chronic disease-drug group for thyroid diseases8.761.5815.60
Chronic disease-drug group for autoimmune diseases3.121.4314.13
Chronic disease group for deafness and hearing loss13.771.4013.77
5
Cardio-cerebrovascular and Renal
(n = 80,855: 8.82%)
Chronic disease-drug group for peripheral vascular disease12.024.7141.57
Chronic disease-drug group for ischemic heart disease29.853.8934.32
Chronic disease-drug group for cerebrovascular diseases19.342.9926.37
Chronic disease-drug group for heart failure21.102.8324.94
Chronic disease group for bradycardias and conduction diseases7.042.5322.33
Chronic disease-drug group for atrial fibrillation15.982.3921.06
Chronic disease-drug group for other psychiatric and behavioral diseases4.632.3020.33
Chronic disease-drug group for chronic kidney diseases21.852.0718.27
Chronic disease-drug group for COPD, emphysema, chronic bronchitis15.851.9817.45
Chronic disease-drug group for colitis and related diseases19.641.9317.04
Chronic disease-drug group for anemia11.321.9317.00
Chronic disease-drug group for neurotic, stress-related, and somatoform diseases18.281.7915.79
Chronic disease-drug group for prostate diseases17.961.7815.67
Chronic disease-drug group for depression and mood diseases21.751.7515.45
Chronic disease-drug group for sleep disorders9.501.6714.69
6
Cardiovascular, Renal, Inflammatory, and Respiratory
(n = 69,720: 7.61%)
Chronic disease-drug group for atrial fibrillation40.075.9945.54
Chronic disease-drug group for heart failure42.575.7043.38
Chronic disease group for bradycardias and conduction diseases11.504.1431.48
Chronic disease-drug group for inflammatory arthropathies11.163.3425.41
Chronic disease-drug group for autoimmune diseases7.073.2524.69
Chronic disease-drug group for anemia17.613.0022.81
Chronic disease-drug group for chronic kidney diseases31.242.9622.52
Chronic disease-drug group for COPD, emphysema, chronic bronchitis20.532.5619.49
Chronic disease-drug group for ischemic heart disease17.012.2216.86
Chronic disease-drug group for peripheral vascular disease5.362.1015.97
Chronic disease group for chronic pancreas, biliary tract, and gallbladder diseases4.881.6412.45
Chronic disease-drug group for cerebrovascular diseases10.291.5912.10
Chronic disease-drug group for colitis and related diseases15.911.5611.90
Chronic disease-drug group for prostate diseases14.981.4811.27
Chronic disease-drug group for hypertension83.431.4711.16
7
Multisystem
(n = 49,592: 5.41%)
Chronic disease group for other digestive diseases23.299.7052.46
Chronic disease-drug group for dementia21.636.4835.07
Chronic disease group for chronic pancreas, biliary tract, and gallbladder diseases16.955.6930.77
Chronic disease-drug group for autoimmune diseases10.955.0327.19
Chronic disease-drug group for inflammatory arthropathies14.804.4323.97
Chronic disease-drug group for anemia19.823.3718.25
Chronic disease-drug group for atrial fibrillation14.802.2111.96
Chronic disease-drug group for heart failure16.372.1911.87
Chronic disease-drug group for colitis and related diseases19.341.9010.29
Chronic disease-drug group for chronic kidney diseases19.731.8710.12
Chronic disease group for bradycardias and conduction diseases4.491.618.73
Chronic disease-drug group for COPD, emphysema, chronic bronchitis12.501.568.44
Chronic disease-drug group for cerebrovascular diseases10.061.558.41
Chronic disease-drug group for esophagus, stomach, and duodenum diseases9.951.377.41
Chronic disease-drug group for osteoarthritis and other degenerative joint diseases25.471.317.10
Categories highlighted in gray are chronic disease categories; those in white are chronic disease-drug categories. Categories in bold reach the O/E ratio threshold of two. Abbreviations: O: disease prevalence in the cluster; O/E ratio: observed/expected ratio; Ex: exclusivity; COPD: chronic obstructive pulmonary disease.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Stafford, G.; Villén, N.; Roso-Llorach, A.; Troncoso-Mariño, A.; Monteagudo, M.; Violán, C. Combined Multimorbidity and Polypharmacy Patterns in the Elderly: A Cross-Sectional Study in Primary Health Care. Int. J. Environ. Res. Public Health 2021, 18, 9216. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18179216

AMA Style

Stafford G, Villén N, Roso-Llorach A, Troncoso-Mariño A, Monteagudo M, Violán C. Combined Multimorbidity and Polypharmacy Patterns in the Elderly: A Cross-Sectional Study in Primary Health Care. International Journal of Environmental Research and Public Health. 2021; 18(17):9216. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18179216

Chicago/Turabian Style

Stafford, Grant, Noemí Villén, Albert Roso-Llorach, Amelia Troncoso-Mariño, Mònica Monteagudo, and Concepción Violán. 2021. "Combined Multimorbidity and Polypharmacy Patterns in the Elderly: A Cross-Sectional Study in Primary Health Care" International Journal of Environmental Research and Public Health 18, no. 17: 9216. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18179216

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop