Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Patterns of Multimorbidity

View ORCID ProfileKien Wei Siah, Chi Heem Wong, Jerry Gupta, View ORCID ProfileAndrew W. Lo
doi: https://doi.org/10.1101/2021.05.13.21256888
Kien Wei Siah
1Laboratory for Financial Engineering, Sloan School of Management, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America
2Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Kien Wei Siah
Chi Heem Wong
1Laboratory for Financial Engineering, Sloan School of Management, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America
2Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jerry Gupta
3Digital Catalyst, Swiss Re, Cambridge, Massachusetts, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrew W. Lo
1Laboratory for Financial Engineering, Sloan School of Management, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America
2Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America
4Sante Fe Institute, Santa Fe, New Mexico, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Andrew W. Lo
  • For correspondence: alo-admin@mit.edu
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Supplementary material
  • Data/Code
  • Preview PDF
Loading

Abstract

Background With multimorbidity becoming the norm rather than the exception, the management of multiple chronic diseases is a major challenge facing healthcare systems worldwide.

Methods Using a large, nationally representative database of electronic medical records from the United Kingdom spanning the years 2005 to 2016 and consisting over 4.5 million patients, we apply statistical methods and network analysis to identify comorbid pairs and triads of diseases and identify clusters of chronic conditions across different demographic groups. Unlike many previous studies, which generally adopt cross-sectional designs, we examine temporal changes in the patterns of multimorbidity. In addition, we perform survival analysis to examine the impact of multimorbidity on mortality.

Results The proportion of the population with multimorbidity has increased by approximately 2.5 percentage points over the last decade, with more than 17% having at least two chronic morbidities. We find that the prevalence and the severity of multimorbidity increase progressively with age. Stratifying by socioeconomic status, we find that people living in more deprived areas are more likely to be multimorbid compared to those living in more affluent areas at all ages. The same trend holds consistently for all years in our data. In addition to a number of strongly associated comorbid pairs (e.g., cardiac-vascular and cardiac-metabolic disorders), we identify three principal clusters: a respiratory cluster, a cardiovascular cluster, and a mixed cardiovascular-renal-metabolic cluster. These are supported by established pathophysiological mechanisms and shared risk factors, and are largely consistent with existing studies in the medical literature.

Conclusions In this paper, we use data-driven methods to characterize multimorbidity patterns in different demographic groups and their evolution over the past decade. Our findings contribute to the better understanding of the epidemiology of multimorbidity that is needed to develop more effective primary care for multimorbid patients.

1 Introduction

Multimorbidity, defined as the coexistence of two or more chronic medical conditions in an individual patient1, is a growing public health concern for healthcare systems worldwide. It has been found to be associated with adverse health outcomes, including a higher risk of mortality, a lower quality of life, increased utilization of healthcare, and correspondingly higher healthcare costs.2–14 It is most prevalent in the elderly population, as organs gradually lose full function with the aging process.8,15–17 With an increasing life expectancy and an aging population, the number of people with multiple health conditions is set to rise, as is public expenditure on long-term medical care. Unfortunately, current healthcare systems are largely designed to treat single diseases, resulting in the need to use multiple services to manage multimorbidity.11,18–20 Due to poor coordination and integration in medical care, causing a lack of continuity in treatment, disorders not designated as the primary condition are often undertreated.21

In order to align primary care more closely to the needs of patients suffering from multiple health conditions, a better understanding of the epidemiology of multimorbidity in the general population is necessary. However, previous studies on multimorbidity have mostly been limited in focus to elderly patients, due to its high prevalence in that population.2,19,20,22– 25 Since this data is typically collected through self-administered questionnaires and research interviews, it may be subject to self-reporting bias. Analyses based on small sample sizes from selected populations are known not to generalize well. Furthermore, many studies employ only a narrow range of methods to study multimorbidity patterns.

In this paper, we aim to characterize multimorbidity patterns not only in older patients, but also across groups with different demographic and socioeconomic statuses, using a large, nationally representative primary care electronic medical records database. We apply standard statistical methods to identify common comorbid pairs and triads of diseases, and use network analysis algorithms to identify clusters of chronic conditions. Unlike many previous studies, which adopted a cross-sectional design (i.e., analyzing data at a single point in time), we examine temporal changes in the patterns of multimorbidity across a decade of patient data. In addition, we analyze the impact of multimorbidity on mortality using survival models. In the discussion, we compare our findings with related works in literature.

2 Methods

a Data

We use anonymized electronic medical records from The Health Improvement Network (THIN)26 database for our analysis. The database contains longitudinal patient data collected at primary care clinics throughout the UK, covering approximately 6% of the UK population. We extract demographic information (e.g., date of birth, sex, geographical location, and socioeconomic group), baseline vitals (e.g., smoking and alcohol status), and medical history (e.g., medical condition and date of diagnosis) from patient records between 2005 and 2016. We categorize the subjects into seven mutually exclusive age groups based on Medical Subject Headings (MeSH) definitions (see Appendix A).27

Diagnoses are recorded in the THIN database using Read Codes, a coded thesaurus of clinical terms used by the UK National Health Service since 1985.28 There is no standard method for the selection and definition of morbidities in the literature. After consulting with medical officers and Life & Health (L&H) actuaries at Swiss Re, we identify chronic conditions in the records, that is, diseases that are either permanent, caused by nonreversible pathological alterations, or require long periods of rehabilitation and care,19,29 and map them to a list of 46 higher level morbidities. Furthermore, we classify the morbidities into 14 System Organ Classes (SOCs) as defined in the Medical Dictionary for Regulatory Activities (MedDRA) dictionary. (See Appendix A for list of morbidities and classifications.) As in similar studies, we define multimorbidity as the presence of at least 2 of the 46 morbidities in a patient.

b Statistical analysis

We examine the distribution of multimorbidity in relation to age and socioeconomic status, as done in Barnett et al.11 However, we use the Index of Multiple Deprivation (IMD) as a proxy for socioeconomic status. The IMD is a widely used measure of relative deprivation or poverty of wards and districts in the UK. It is computed using census data as a weighted index of deprivation in seven domains, including income, employment, education, health, crime, barriers to housing and services, and living environment.30 (IMD data was available only for a subset of the patients. See Appendix B for the sample sizes used in this analysis.) We note that the same approach, defining socioeconomic status by the area of residence, has been used in previous studies.11,31

For each age group, we also compute the observed prevalence for all individual, pairs, and triplets of morbidities. By the assumptions of probability theory, we expect diseases that are independent to co-occur at a rate close to the product of the observed prevalence of each individual constituent disease (i.e., the expected prevalence). Therefore, by comparing the ratio of the observed prevalence versus the expected prevalence (i.e., the lift), we can identify pairs and triads of diseases that occur together more frequently than expected by chance, possibly driven by an underlying pathophysiological mechanism. As a second metric, we estimate the odds ratio using logistic regression models to determine the association between each pair of diseases, both without adjustment and adjusted by age, sex, and all other diseases.

Next, we construct multimorbidity networks to study the natural clustering of diseases in the dataset. We consider diseases as nodes with sizes proportional to their observed prevalence. For each pair of diseases, we connect their nodes with an undirected edge weighted by the estimated lift, a measure of the strength of the association between the comorbid pair. This creates a dense network where each node is linked to almost every other node. This density, however, makes visualization and inference difficult. As a pre-processing step for subsequent analysis, we extract the main graph structure by removing edges from the adjacency matrix that are peripheral and relatively unimportant. We prune the edges between nodes that have joint prevalence below the 90th percentile, and keep only the edges that have a lift above 2.0, i.e., those edges between pairs that co-occur two times more frequently than expected by chance. Similar thresholds have been used in related studies.2,32–35

We compute measures of centrality to identify the most important vertices in the multimorbidity network. In particular, for each node, we compute the degree centrality, which is defined as the number of links incident on a node, a direct measure of the connectivity of a node. In this context, a disease with high degree centrality is important because it often co-occurs with a large number of pathologies. We also estimate the eigenvector centrality, a measure of the transitive influence of nodes. To calculate the eigencentrality, each node is assigned a score that is proportional to the sum of the scores of all of its neighbors. Nodes with high eigencentrality either have many connections, or are connected to important neighbors. In addition, we compute the graph clustering coefficient (also known as the transitivity) as a quantitative measure of the network’s tendency to aggregate in smaller subgroups. To identify any clusters embedded in the multimorbidity networks, we apply a community detection algorithm based on modularity maximization36– 38 to partition nodes into groups that have dense intra-group connections and sparse inter-group connections.

To gain insight into temporal disease associations, we construct directed multimorbidity networks. We extract from each patient’s medical history a sequence of diseases ordered by the time of diagnosis. Using these trajectories, we can derive the probability of any given disease conditional on some prior diagnosis, i.e., Prob(Disease B given Disease A). We use these probabilities as weights of the directed edges in the network. As before, we prune the network based on node prevalence and edge weights. Since these connections are directed, we can compute the in-degree and out-degree centralities, defined as the number of edges directed to the node, and the number of edges directed from the node to others, respectively. A node with a high in-degree centrality is often diagnosed following other diseases; a node with a high out-degree centrality often leads to subsequent diagnoses in other diseases. These metrics are useful for understanding disease progression, and any causal or contributory relationships between diseases.

Finally, we examine the impact of multimorbidity on mortality by performing survival analysis on the dataset. We use the five-year overall survival as the primary outcome variable, and consider in our models a range of features, including demographic group, baseline vitals, baseline medical history, the severity of multimorbidity as quantified by the number of co-occurring chronic conditions, and the presence of any of the top ten most prevalent pairs and triplets of morbidities as observed in the aged and elderly age groups. We exclude those subjects aged 65 or less from this part of the analysis, as younger age groups have five-year overall mortality rates close to zero.

We explore three standard methods used in survival modeling, the Cox proportional hazards model39, the regularized Cox model, and the accelerated failure time model, and additionally, we apply a nonlinear and non-parametric neural network survival model.40 For model estimation and validation, we randomly split the original dataset into two disjoint sets, a training set that comprises 70% of the data, and a testing set that comprises the remaining 30%. We use the training set to estimate our models, and keep the testing set as an out-of-sample dataset for performance validation. We use the concordance index (C-index) as the metric for model performance. This metric is commonly used in survival analysis to evaluate its predictive power.41 It is a measure of the concordance between orderings of observed survival times and the predicted times or risks. (A C-index of 0.5 corresponds to a random model, while a value of 1.0 corresponds to a perfect model.) We use cross-validation to tune the hyperparameters of the models.

In addition to discriminative power, we assess the calibration of our models by comparing the actual and the predicted survival probabilities at 36, 48, and 60 months of overall survival. For each time cutoff, we divide the test set into quintiles based on the predicted risk scores. We then compute the average predicted score and the true survival probability observed in each of the quintiles. Lastly, we create calibration plots by plotting the observed probabilities against the predicted probabilities. In the ideal case, the points should lie as close as possible to the diagonal line, which represents perfect calibration.

3 Results

We summarize the demographic statistics of the study population in Table 1. On average, the dataset consists of approximately 4.6 million patients each year, with an even mix of both sexes in all years. Most of the patient records were collected in England, which makes up the largest part of the population of the United Kingdom. However, the distribution in geographical location has evolved over the years, shifting towards other regions in the country. Over 60% of the patients are in the Adult (19 to 45 years old) and Middle-Aged (45 to 65 years old) age groups, as defined by the MeSH classification (see Appendix A). Approximately 15% are over 65 years old.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 1. Demographics of the dataset between 2005 and 2016.

The proportion of the population with multimorbidity has increased by approximately 2.5 percentage points over the last decade, with more than 17% of all patients having at least two chronic morbidities in 2016. We find that the prevalence and the severity of multimorbidity increase progressively with age (see Figure 1). By age 60, approximately half the population has been diagnosed with at least one chronic condition, after which we observe a steep rise in multimorbidity, with close to 1 in 3 patients suffering from at least two morbidities by age 70. Stratifying the prevalence of multimorbidity by IMD in Figure 2, we find that people living in more deprived areas are more likely to be multimorbid compared to those living in more affluent areas at all ages. The same trend holds consistently for all years in our data.

Figure 1.
  • Download figure
  • Open in new tab
Figure 1.

Number of co-occurring chronic conditions by age in 2016. See Appendix C for other years.

Figure 2.
  • Download figure
  • Open in new tab
Figure 2.

Prevalence of multimorbidity by age and IMD in 2016. See Appendix C for other years. Higher score corresponds to greater deprivation.

We characterize the epidemiology of individual diseases by plotting heat maps of disease prevalence in different age groups separately. We find that asthma and respiratory conditions have high prevalence across all age groups, with the former occurring especially frequently in the Adolescent age group (13 to 19 years old). We observe the onset of metabolic and cardiovascular diseases in the Middle-Aged and older age groups, in particular, of diabetes and hypertension. Not surprisingly, diseases such as dementia, kidney diseases, and stroke occur most frequently in the oldest patients (65 years and above). We observe an increasing trend in prevalence for some diseases. For example, the prevalence of diabetes in the Aged age group (65 to 80 years old) has increased by almost 35% in the past decade (see Figure 3). In contrast, the prevalence of diseases such as angina has been falling over the study period.

Figure 3.
  • Download figure
  • Open in new tab
Figure 3.

Single disease prevalence in the Aged subgroup between 2005 and 2016. See Figures 7 and 8 for disorder to index mapping and MedDRA SOC abbreviations. See Appendix C for other age groups.

In Table 2, we summarize the lift and odds ratio of the top ten most frequently co-occurring pairs of diseases in the Aged age group in 2016. (See Appendix C for other age groups and years.) In all age groups, asthma occurs in combination with other respiratory-related diseases approximately twice more often than expected by chance (i.e., the lift is greater than 2.0). Additionally, the estimated odds ratios, both unadjusted and adjusted, indicate that patients with asthma are at least twice as likely to suffer from other respiratory conditions at the same time, and vice versa.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 2. Lift and odds ratio of the top 10 most prevalent multimorbidity pairs in the Aged subgroup in 2016. See Appendix C for other age groups and years.

Hypertension is most associated with a second condition in the older age groups, although most pairs do not necessarily occur more frequently than by chance. The combination of hypertension and diabetes stands out with a relatively high lift and an odds ratio that is greater than 2.0, suggesting that this disease pair may be comorbid. Angina and coronary artery disease (CAD) also demonstrate a strong association in the Aged and Elderly age groups with unusually high lift and odds ratio.

To better visualize the data, we plot the lift of all combinations of disease pairs in heat maps, stratified by MedDRA system organ classes. (See Figure 4 and Appendix C for other age groups and years.) The co-occurrence of cardiac-cardiac and cardiac-respiratory disorders is a major risk across all age groups. We observe significant coupling between cardiac and hepatobiliary disorders in the Adolescent and Child (2 to 13 years old) age groups. On the other hand, combinations of cardiac-vascular and cardiac-metabolic disorders are the most dominant in the Middle-Aged and older age groups. We observe the same general patterns across time.

Figure 4.
  • Download figure
  • Open in new tab
Figure 4.

Heat map of lift of multimorbidity pairs in the Aged subgroup in 2016. See Figures 7 and 8 for disorder to index mapping and MedDRA SOC abbreviations. See Appendix C for other age groups and years.

As shown in Figure 1, the proportion of patients with three or more co-occurring disorders is small in the younger age groups. For patients aged 45 years and older, triplets involving angina, CAD, hypertension, diabetes and myocardial infarction (MI) occur most frequently with high lift, suggesting strong correlations between these diseases (see Table 3).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 3. Lift of the top 10 most prevalent multimorbidity triplets in the Aged subgroup in 2016. See Appendix C for other age groups and years.

In Figures 5 and 6, we plot the undirected and directed multimorbidity networks observed in the Aged age group in 2016. (See Appendix C for other age groups and years.) Instead of a force-directed layout, we place the nodes in fixed positions around a circle to allow easy visualization of temporal changes in connections and clusters when comparing plots from different years. The edge thickness is proportional to the lift between each disease pair. Apart from single-node clusters, the communities detected using modularity maximization are given different colors.

Figure 5.
  • Download figure
  • Open in new tab
Figure 5.

Undirected multimorbidity network in the Aged subgroup in 2016. Edge thickness is proportional to the lift between each disease pair. Intra-group edges and inter-group edges are represented by solid lines and dashed lines, respectively. Only communities with more than one node are colored. See Figures 7 and 8 for mapping of disorder to index and MedDRA SOC abbreviations. See Appendix C for other age groups and years.

Figure 6.
  • Download figure
  • Open in new tab
Figure 6.

Directed multimorbidity network in the Aged subgroup in 2016. Edge thickness is proportional to the lift between each disease pair. Intra-group edges and inter-group edges are represented by solid lines and dashed lines, respectively. Only communities with more than one node are colored. See Figures 7 and 8 for mapping of disorder to index and MedDRA SOC abbreviations. See Appendix C for other age groups and years.

Figure 7.
  • Download figure
  • Open in new tab
Figure 7.

Mapping between index and chronic conditions.

Figure 8.
  • Download figure
  • Open in new tab
Figure 8.

Abbreviations for MedDRA SOCs used in figures.

In Tables 4 and 5, we identify clusters that remain relatively stable throughout the years in undirected and directed multimorbidity networks, respectively. We find between one and four clusters for each age group. The number of diseases in each cluster ranges between two and twelve. In general, the communities found in Adolescent and younger patients can vary greatly from year to year compared to older age groups, where the clusters evolve very little over time. This is expected, given that only a small proportion of the former cohort suffers from more than two co-occurring disorders (see Figure 1), so the results are sensitive to small changes in prevalence each year.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 4. Clusters identified in undirected multimorbidity networks in different age groups.
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 5. Clusters identified in directed multimorbidity networks in different age groups.

A respiratory cluster of asthma, chronic obstructive pulmonary disease (COPD), and respiratory-related diseases appears to be present in all age groups in both undirected and directed graphs. Similarly, a vascular-metabolic-hepatobiliary-renal cluster that is characterized by hypertension, diabetes, liver diseases, and kidney diseases, with the occasional appearance of cardiac disorders, is also present in almost all cohorts. As observed in previous analyses, we also find several clusters dominated by cardiovascular disorders such as angina, CAD, myocardial infarction (MI), atrial fibrillation, cardiac arrhythmia, heart failure, heart valve disorder (HVD), stroke, peripheral artery disease (PAD), and transient ischemic attack (TIA).

In Tables 6 and 7, we summarize the top five diseases for each centrality measure. (See Appendix C for the full set of results.) The degree centrality and eigencentrality for hypertension, diabetes, CAD, and angina are the highest when all age groups are aggregated in undirected multimorbidity networks. In the Adolescent and younger age groups, kidney disease shows both high degree centrality and eigencentrality. Other important nodes include respiratory-related diseases and HVD, which have high degree centrality and high eigencentrality, respectively. For the Adult and Middle-Aged age groups, hypertension and diabetes are the most central nodes with respect to both measures. In the Aged and Elderly age groups, we find that cardiac disorders make up all of the top five most connected nodes.

View this table:
  • View inline
  • View popup
Table 6. Centrality measures for top five diseases in undirected multimorbidity networks with mean computed over time.
View this table:
  • View inline
  • View popup
Table 7. Centrality measures for top five diseases in directed multimorbidity networks with mean computed over time. We exclude eigencentralities that are close to zero.

We observe similar results in directed networks. In general, hypertension, diabetes, and respiratory-related diseases demonstrate high in-degree centrality and eigencentrality, while cardiac disorders show high out-degree centrality. In the Middle-Aged and younger age groups, asthma emerges as a new central node with high in-degree centrality, while the top five diseases for the Aged and Elderly age groups remain dominated by cardiovascular diseases.

We summarize the dataset used for survival analysis in Table 8. The sample consists of approximately 390,000 patients in the Aged and Elderly age groups for each year between 2010 and 2012. More than 50% of the patients are multimorbid. In terms of predicting five-year overall survival, we find the performance of the linear and nonlinear survival models explored to be very similar (see Appendix D). We focus on the Cox model here due to its ease of interpretability. The model achieves a promising C-index of 0.81 (95% CI 0.80–0.81) on out-of-sample data in 2012. In addition, its calibration curves lay close to the ideal diagonal, indicating that the model is well calibrated, i.e., the model does not systematically overestimate or underestimate survival rates in any of the quintiles. (See Appendix C for plots.)

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 8. Summary statistics of dataset used for survival analysis.

In Table 9, we extract the top ten coefficients in the Cox model to identify specific risk factors. To correct for multiple testing, we perform the Benjamini-Hochberg adjustment with a 5% false discovery rate for identifying significant factors. Apart from cancers, we find the presence of multimorbidity to be a strong adverse risk factor, i.e., the higher the number of co-occurring chronic conditions, the greater the mortality risk. For example, the hazard ratio of having four or more chronic conditions is 2.44 (95% CI 2.22–2.69). We also find a high IMD, corresponding to a lower socioeconomic status, to be significantly associated with increased risk, although this factor is not in the top ten coefficients.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 9. Top ten risk factors of the Cox Proportional Hazards model for 2012. Asterisks denote factors that are significant after the Benjamini-Hochberg adjustment.

4 Discussion

With multimorbidity becoming the norm rather than the exception2,12,17,22,42,43, the management of multiple chronic diseases in older adults is a major challenge facing healthcare systems worldwide. It is clear that a better understanding of the epidemiology of multimorbidity is required to develop more effective preventive interventions and better primary medical care for multimorbid patients. In this paper, we use data-driven methods to characterize multimorbidity patterns in different demographic groups and their evolution over the past decade, using a large, representative electronic medical records database consisting of over 4.5 million patients.

Consistent with other studies, we find that the prevalence and severity of multimorbidity increase substantially with age. In addition, we observe social inequalities in multimorbidity, with patients in socioeconomically deprived areas more likely to be multimorbid.11,12,31,42,44– 46 Our findings also support the role of hypertension as an important risk factor in older adults, as reported in the literature.2,33,47 Hypertension is one of the most prevalent and most central chronic conditions in our dataset, and one that serves as an important bridge between many diseases in our networks. Other trends identified in our analysis, such as the falling prevalence of angina48–51 and the growing prevalence of diabetes52, are also well documented in previously published population studies.

In our pairwise analysis, we find strong association between multiple pairs of chronic conditions, including between asthma and respiratory-related diseases53,54 in the Adolescent age group, between hypertension and diabetes25,55–57 and between CAD and angina58 among older patients, and between cardiovascular and respiratory disorders in all age groups.2 Triplets involving cardiovascular and metabolic disorders, such as CAD, hypertension, and diabetes, also occurred more frequently than expected by chance.2,22,25,59

Our network analysis further identified several meaningful communities that are common across all demographics, including a respiratory cluster (e.g., asthma and COPD)60, a cardiovascular cluster19,61, and a mixed cardiovascular-renal-metabolic cluster32,62–64, all of which are supported by either established pathophysiological mechanisms or shared risk factors. For example, it is well known that cardiovascular diseases are one of the most common complications of diabetes. While we do not find any particular multimorbidity pattern to have a significant effect on mortality, our models do indeed verify the substantial burden of multimorbidity (as quantified by the number of co-occurring chronic conditions) on overall survival in older patients.7,12,23,65,66

However, we must emphasize that our results do not necessarily imply any causal link between diseases identified to be in the same cluster. The association might be attributable to shared risk factors (e.g., smoking) or other adverse events, and any temporal relationships to be inferred from the multimorbidity directed networks might be administrative in nature (e.g., incomplete medical records that are rectified in subsequent visits) or biased by delayed diagnosis.

In general, the lack of an accepted standard for defining multimorbidity makes it difficult for any meaningful comparison of results across different studies.67,687 Moreover, because results can be highly dependent on the study population, the disease ontology used, and the number of chronic conditions considered, it is not uncommon for studies to report seemingly conflicting findings. In this paper, we consider a wide range of demographic groups and a total of 46 morbidities, which is more than most similar studies,11 and well above the minimum of 11 to 12 as recommended by systematic reviews in this field of research.68,69 In addition, our findings are largely consistent with existing studies in the medical literature.

Lastly, we note that cancer appears to be under represented in the THIN database. This is because many cancer patients are treated separately in cancer centers under the care of specialized clinical teams. Unfortunately, data on such patients rarely make their way back to the primary care clinics where the THIN data is collected, leading to a gap in this area.

5 Conclusions

Current healthcare systems are largely centered on single-disease approaches to treatment, resulting in the fragmentation of care and a lack of continuity in the management of multiple diseases. Even most clinical trials exclude multimorbid patients. Because multimorbidity is more common in disadvantaged groups, the current structure exacerbates health inequalities in society. In this paper, we apply statistical methods and network analysis to characterize multimorbidity associations in the general UK population using a large electronic medical records database spanning the years 2005 to 2016. Our findings contribute to a better understanding of multimorbidity that may be useful for the early detection and prevention of comorbidities, for example, recommending that hypertension patients reduce sugar intake as a preventive measure for diabetes.

There is a pressing need for a universal framework that standardizes the way that multimorbidity is assessed (e.g., the minimum number of diseases and the choice of chronic conditions to include) in order to facilitate comparisons between studies and populations. With the “Omics” revolution, the combination of phenotypic, genomic, and epigenomic data has the potential to provide deeper insights into the underlying pathophysiological associations between comorbid diseases. Unfortunately, the availability of such linked datasets remains very limited. Further research is also needed to better understand the impact of multimorbidity on different health outcomes, such as the quality of life and healthcare costs, in order to align the primary healthcare system more closely to the needs to multimorbid patients.

Data Availability

The data supporting the current study have not been deposited in a public repository due to its proprietary nature. The data are available from The Health Improvement Network (THIN) at https://www.the-health-improvement-network.com/en/. Restrictions apply to the availability of the data, which were used under license for this study.

Declarations

Ethics approval and consent to participate

This study was approved by the Scientific Review Committee at IQVIA.

Consent for publication

Not applicable

Availability of data and materials

The data that support the findings of this study are available from The Health Improvement Network (THIN; https://www.the-health-improvement-network.com/en/) but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of THIN.

Competing interests

K.W.S. and C.H.W. declare no competing interests. J.G. is an employee of Swiss Re and declares no competing interests. A.L. reports personal investments in private biotech companies, biotech venture capital funds, and mutual funds. A.L. is a co-founder and partner of QLS Advisors, a healthcare analytics and consulting company; an advisor to BrightEdge Ventures; a director of BridgeBio Pharma, Roivant Sciences, and Annual Reviews; chairman emeritus and senior advisor to AlphaSimplex Group; and a member of the Board of Overseers at Beth Israel Deaconess Medical Center and the NIH’s National Center for Advancing Translational Sciences Advisory Council and Cures Acceleration Network Review Board. During the most recent six-year period, A.L. has received speaking/consulting fees, honoraria, or other forms of compensation from: AIG, AlphaSimplex Group, BIS, BridgeBio Pharma, Citigroup, Chicago Mercantile Exchange, Financial Times, FONDS Professionell, Harvard University, IMF, National Bank of Belgium, Q Group, Roivant Sciences, Scotia Bank, State Street Bank, University of Chicago, and Yale University.

Funding

No direct funding was received for this study; general research support was provided by the MIT Laboratory for Financial Engineering and its sponsors. The authors were personally salaried by their institutions during the period of writing (though no specific salary was set aside or given for the writing of this paper).

Author contributions

Conceptualization, K.W.S., C.H.W., J.G., and A.L.; Resources, J.G., and A.L.; Methodology, K.W.S., C.H.W., J.G., and A.L.; Software, K.W.S.; Formal Analysis, K.W.S.; Writing – Original Draft, K.W.S. and A.L.; Writing – Review & Editing, K.W.S., C.H.W., J.G., and A.L.; Supervision, A.L.; Project Administration, J.G., and A.L.

Acknowledgments

Research support from the MIT Laboratory for Financial Engineering is gratefully acknowledged. The views and opinions expressed in this article are those of the authors only, and do not necessarily represent the views and opinions of any institution or agency, any of their affiliates or employees, or any of the individuals acknowledged above.

References

  1. 1.↵
    World Health Organization. Multimorbidity Technical Series on Safer Primary Care Multimorbidity: Technical Series on Safer Primary Care. (2016).
  2. 2.↵
    Schäfer, I. et al. Reducing complexity: A visualisation of multimorbidity by combining disease clusters and triads. BMC Public Health 14, 1285 (2014).
    OpenUrlCrossRefPubMed
  3. 3.
    Kadam, U. & Croft, P. Clinical multimorbidity and physical function in older adults: a record and health status linkage study in general practice. Fam. Pract. 24, 412–419 (2007).
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.
    Laux, G., Kuehlein, T., Rosemann, T. & Szecsenyi, J. Co- and multimorbidity patterns in primary care based on episodes of care: Results from the German CONTENT project. BMC Health Serv. Res. 8, 14 (2008).
    OpenUrlCrossRefPubMed
  5. 5.
    Fung, C. H. et al. The relationship between multimorbidity and patients’ ratings of communication. J. Gen. Intern. Med. 23, 788–793 (2008).
    OpenUrlCrossRefPubMed
  6. 6.
    Schoenberg, N. E., Kim, H., Edwards, W. & Fleming, S. T. Burden of Common Multiple-Morbidity Constellations on Out-of-Pocket Medical Expenditures Among Older Adults. Gerontologist 47, 423–437 (2007).
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    Gijsen, R. et al. Causes and consequences of comorbidity: A review. J. Clin. Epidemiol. 54, 661–674 (2001).
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.↵
    Salisbury, C., Johnson, L., Purdy, S., Valderas, J. M. & Montgomery, A. A. Epidemiology and impact of multimorbidity in primary care: a retrospective cohort study. Br. J. Gen. Pract. 61, e12–e21 (2011).
    OpenUrlAbstract/FREE Full Text
  9. 9.
    Wolff, J. L., Starfield, B. & Anderson, G. Prevalence, Expenditures, and Complications of Multiple Chronic Conditions in the Elderly. Arch. Intern. Med. 162, 2269 (2002).
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.
    Fortin, M. et al. Multimorbidity and quality of life in primary care: A systematic review. Health Qual. Life Outcomes 2, 1–12 (2004).
    OpenUrlCrossRefPubMed
  11. 11.↵
    Barnett, K. et al. Epidemiology of multimorbidity and implications for health care, research, and medical education: A cross-sectional study. Lancet 380, 37–43 (2012).
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.↵
    Marengoni, A. et al. Aging with multimorbidity: A systematic review of the literature. Ageing Res. Rev. 10, 430–439 (2011).
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.
    Loza, E., Jover, J. A., Rodriguez, L. & Carmona, L. Multimorbidity: Prevalence, Effect on Quality of Life and Daily Functioning, and Variation of This Effect When one Condition Is a Rheumatic Disease. Semin. Arthritis Rheum. 38, 312–319 (2009).
    OpenUrlCrossRefPubMed
  14. 14.↵
    Crentsil, V., Ricks, M. O., Xue, Q. L. & Fried, L. P. A pharmacoepidemiologic study of community-dwelling, disabled older women: Factors associated with medication use. Am. J. Geriatr. Pharmacother. 8, 215–224 (2010).
    OpenUrlCrossRefPubMed
  15. 15.↵
    Walker, A. E. Multiple chronic diseases and quality of life: Patterns emerging from a large national sample, Australia. Chronic Illn. 3, 202–218 (2007).
    OpenUrlCrossRefPubMed
  16. 16.
    Van den Akker, M., Buntix, F., Metsemakers, J. F. M., Roos, S. & Knottnerus, J. A. Multimorbidity in general practice: Prevalence, incidence, and determinants of co-occurring chronic and recurrent diseases. J. Clin. Epidemiol. 51, 367–375 (1998).
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    Fortin, M. Prevalence of Multimorbidity Among Adults Seen in Family Practice. Ann. Fam. Med. 3, 223–228 (2005).
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    Boyd, C. M. et al. Clinical Practice Guidelines and Quality of Care for Older Patients With Multiple Comorbid Diseases. JAMA 294, 716 (2005).
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.↵
    Marengoni, A., Rizzuto, D., Wang, H.-X., Winblad, B. & Fratiglioni, L. Patterns of Chronic Multimorbidity in the Elderly Population. J. Am. Geriatr. Soc. 57, 225–230 (2009).
    OpenUrlCrossRefPubMedWeb of Science
  20. 20.↵
    Marengoni, A. et al. In-hospital death and adverse clinical events in elderly patients according to disease clustering: The REPOSI study. Rejuvenation Res. 13, 469–477 (2010).
    OpenUrlCrossRefPubMed
  21. 21.↵
    Redelmeier, D. A., Tan, S. H. & Booth, G. L. The Treatment of Unrelated Disorders in Patients with Chronic Medical Diseases. N. Engl. J. Med. 338, 1516–1520 (1998).
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    Schäfer, I. et al. Multimorbidity Patterns in the Elderly: A New Approach of Disease Clustering Identifies Complex Interrelations between Chronic Conditions. PLoS One 5, e15941 (2010).
    OpenUrlCrossRefPubMed
  23. 23.↵
    Ferrer, A., Formiga, F., Sanz, H., Almeda, J. & Padrós, G. Multimorbidity as specific disease combinations, an important predictor factor for mortality in octogenarians: the Octabaix study. Clin. Interv. Aging 12, 223–231 (2017).
    OpenUrl
  24. 24.
    Diederichs, C., Berger, K. & Bartels, D. B. The Measurement of Multiple Chronic Diseases--A Systematic Review on Existing Multimorbidity Indices. Journals Gerontol. Ser. A Biol. Sci. Med. Sci. 66A, 301–311 (2011).
    OpenUrl
  25. 25.↵
    Kirchberger, I. et al. Patterns of Multimorbidity in the Aged Population. Results from the KORA-Age Study. PLoS One 7, e30556 (2012).
    OpenUrlCrossRefPubMed
  26. 26.↵
    The Health Improvement Network. The Health Improvement Network. https://www.the-health-improvement-network.com/en/.
  27. 27.↵
    Kastner, M., Wilczynski, N. L., Walker-Dilks, C., McKibbon, K. A. & Haynes, B. Agespecific search strategies for medline. J. Med. Internet Res. 8, (2006).
  28. 28.↵
    National Health Services. Read Codes. https://digital.nhs.uk/services/terminology-and-classifications/read-codes (2020).
  29. 29.↵
    Timmreck, T. C., Cole, G. E., James, G. & Butterworth, D. D. Health education and health promotion: A look at the jungle of supportive fields, philosophies and theoretical foundations. Health Educ. 18, 23–28 (1987).
    OpenUrlPubMed
  30. 30.↵
    Ministry of Housing Communities & Local Government. English indices of deprivation. https://www.gov.uk/government/collections/english-indices-of-deprivation (2012).
  31. 31.↵
    Orueta, J. F., García-Álvarez, A., García-Goñi, M., Paolucci, F. & Nuño-Solinís, R. Prevalence and Costs of Multimorbidity by Deprivation Levels in the Basque Country: A Population Based Study Using Health Administrative Databases. PLoS One 9, e89787 (2014).
    OpenUrl
  32. 32.↵
    Aguado, A., Moratalla-Navarro, F., López-Simarro, F. & Moreno, V. MorbiNet: multimorbidity networks in adult general population. Analysis of type 2 diabetes mellitus comorbidity. Sci. Rep. 10, 2416 (2020).
    OpenUrl
  33. 33.↵
    Feldman, K. et al. Insights into Population Health Management Through Disease Diagnoses Networks. Sci. Rep. 6, 30465 (2016).
    OpenUrlCrossRefPubMed
  34. 34.
    Leva, F. & Bitonti, D. Network analysis of comorbidity patterns in heart failure patients using administrative data. Epidemiol. Biostat. Public Heal. 15, (2018).
  35. 35.↵
    Liu, J. et al. Comorbidity Analysis According to Sex and Age in Hypertension Patients in China. Int. J. Med. Sci. 13, 99–107 (2016).
    OpenUrlCrossRefPubMed
  36. 36.↵
    Brandes, U. et al. On modularity clustering. IEEE Trans. Knowl. Data Eng. 20, 172–188 (2008).
    OpenUrlCrossRef
  37. 37.
    Girvan, M. & Newman, M. E. J. Community structure in social and biological networks. Proc. Natl. Acad. Sci. 99, 7821–7826 (2002).
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    Newman, M. E. J. Fast algorithm for detecting community structure in networks. Phys. Rev. E - Stat. Physics, Plasmas, Fluids, Relat. Interdiscip. Top. 69, 5 (2004).
    OpenUrl
  39. 39.↵
    Cox, D. R. Regression Models and Life-Tables. J. R. Stat. Soc. Ser. B 34, 187–202 (1972).
    OpenUrlCrossRef
  40. 40.↵
    Faraggi, D. & Simon, R. A neural network model for survival data. Stat. Med. 14, 73–82 (1995).
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    Harrell, F. E. Evaluating the Yield of Medical Tests. JAMA J. Am. Med. Assoc. 247, 2543 (1982).
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    Violan, C. et al. Prevalence, Determinants and Patterns of Multimorbidity in Primary Care: A Systematic Review of Observational Studies. PLoS One 9, e102149 (2014).
    OpenUrlCrossRefPubMed
  43. 43.↵
    Boyd, C. M. & Fortin, M. Future of multimorbidity research: How should understanding of multimorbidity inform health system design? Public Health Rev. 32, 451–474 (2010).
    OpenUrlCrossRefWeb of Science
  44. 44.↵
    Schiøtz, M. L., Stockmarr, A., Høst, D., Glümer, C. & Frølich, A. Social disparities in the prevalence of multimorbidity - A register-based population study. BMC Public Health 17, 422 (2017).
    OpenUrlCrossRefPubMed
  45. 45.
    Dugravot, A. et al. Social inequalities in multimorbidity, frailty, disability, and transitions to mortality: a 24-year follow-up of the Whitehall II cohort study. Lancet Public Heal. 5, e42–e50 (2020).
    OpenUrl
  46. 46.↵
    Schäfer, I. et al. The influence of age, gender and socio-economic status on multimorbidity patterns in primary care. first results from the multicare cohort study. BMC Health Serv. Res. 12, 89 (2012).
    OpenUrlCrossRefPubMed
  47. 47.↵
    Chen, Y. & Xu, R. Mining Cancer-Specific Disease Comorbidities from a Large Observational Health Database. Cancer Inform. 13s1, CIN.S13893 (2014).
  48. 48.↵
    Abdalla, S. M., Yu, S. & Galea, S. Trends in Cardiovascular Disease Prevalence by Income Level in the United States. JAMA Netw. Open 3, e2018150 (2020).
    OpenUrl
  49. 49.
    Yoon, S. S. (Sarah), Dillon, C. F., Illoh, K. & Carroll, M. Trends in the Prevalence of Coronary Heart Disease in the U.S.: National Health and Nutrition Examination Survey, 2001–2012. Am. J. Prev. Med. 51, 437–445 (2016).
    OpenUrlPubMed
  50. 50.
    Benjamin, E. J. et al. Heart Disease and Stroke Statistics—2019 Update: A Report From the American Heart Association. Circulation 139, e56–e528 (2019).
    OpenUrlCrossRefPubMed
  51. 51.↵
    Lampe, F. C. et al. Is the prevalence of coronary heart disease falling in British men? Heart 86, 499–505 (2001).
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    Boyle, J. P. et al. Projection of Diabetes Burden Through 2050: Impact of changing demography and disease prevalence in the U.S. Diabetes Care 24, 1936–1940 (2001).
    OpenUrlAbstract/FREE Full Text
  53. 53.↵
    Boulet, L. P. Influence of comorbid conditions on asthma. European Respiratory Journal vol. 33 897–906 (2009).
    OpenUrlAbstract/FREE Full Text
  54. 54.↵
    Bardin, P. G., Rangaswamy, J. & Yo, S. W. Managing comorbid conditions in severe asthma. Med. J. Aust. 209, S11-S17.e3 (2018).
    OpenUrl
  55. 55.↵
    Lago, R. M., Singh, P. P. & Nesto, R. W. Diabetes and hypertension. Nature Clinical Practice Endocrinology and Metabolism vol. 3 667 (2007).
    OpenUrlCrossRefWeb of Science
  56. 56.
    Long, A. N. & Dagogo-Jack, S. Comorbidities of Diabetes and Hypertension: Mechanisms and Approach to Target Organ Protection. Journal of Clinical Hypertension vol. 13 244–251 (2011).
    OpenUrlCrossRefPubMed
  57. 57.↵
    de Boer, I. H. et al. Diabetes and Hypertension: A Position Statement by the American Diabetes Association. Diabetes Care 40, 1273–1284 (2017).
    OpenUrlFREE Full Text
  58. 58.↵
    Centers for Disease Control and Prevention. Coronary Artery Disease. https://www.cdc.gov/heartdisease/coronary_ad.htm (2019).
  59. 59.↵
    Rana, J. S., Nieuwdorp, M., Jukema, J. W. & Kastelein, J. J. P. Cardiovascular metabolic syndrome - An interplay of, obesity, inflammation, diabetes and coronary heart disease. Diabetes, Obesity and Metabolism vol. 9 218–232 (2007).
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    Maselli, D. J. & Hanania, N. A. Asthma COPD overlap: Impact of associated comorbidities. Pulmonary Pharmacology and Therapeutics vol. 52 27–31 (2018).
    OpenUrl
  61. 61.↵
    Déruaz-Luyet, A. et al. Multimorbidity and patterns of chronic conditions in a primary care population in Switzerland: A cross-sectional study. BMJ Open 7, e013664 (2017).
    OpenUrlCrossRef
  62. 62.↵
    Cherney, D. Z. I. et al. Impact of Cardio-Renal-Metabolic Comorbidities on Cardiovascular Outcomes and Mortality in Type 2 Diabetes Mellitus. Am. J. Nephrol. 51, 74–82 (2020).
    OpenUrl
  63. 63.
    Arnold, S. V. et al. Burden of cardio-renal-metabolic conditions in adults with type 2 diabetes within the Diabetes Collaborative Registry. Diabetes, Obes. Metab. 20, 2000–2003 (2018).
    OpenUrlPubMed
  64. 64.↵
    Arnold, S. V. et al. Cardiovascular Outcomes and Mortality in Type 2 Diabetes with Associated Cardio-Renal-Metabolic Comorbidities. Diabetes 67, 1582–P (2018).
    OpenUrlCrossRef
  65. 65.↵
    Lee, S. J., Lindquist, K., Segal, M. R. & Covinsky, K. E. Development and Validation of a Prognostic Index for 4-Year Mortality in Older Adults. JAMA 295, 808 (2006).
    OpenUrl
  66. 66.↵
    Walter, L. C. et al. Development and Validation of a Prognostic Index for 1-Year Mortality in Older Adults After Hospitalization. JAMA 285, 2987 (2001).
    OpenUrlCrossRefPubMedWeb of Science
  67. 67.
    Fortin, M., Hudon, C., Haggerty, J., Akker, M. Van Den & Almirall, J. Prevalence estimates of multimorbidity: A comparative study of two sources. BMC Health Serv. Res. 10, 111 (2010).
    OpenUrlCrossRefPubMed
  68. 68.↵
    Fortin, M., Stewart, M., Poitras, M.-E., Almirall, J. & Maddocks, H. A Systematic Review of Prevalence Studies on Multimorbidity: Toward a More Uniform Methodology. Ann. Fam. Med. 10, 142–151 (2012).
    OpenUrlAbstract/FREE Full Text
  69. 69.↵
    Diederichs, C., Berger, K. & Bartels, D. B. The measurement of multiple chronic diseases - A systematic review on existing multimorbidity indices. Journals Gerontol. - Ser. A Biol. Sci. Med. Sci. 66, 301–311 (2011).
    OpenUrl
Back to top
PreviousNext
Posted May 14, 2021.
Download PDF

Supplementary Material

Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Patterns of Multimorbidity
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Patterns of Multimorbidity
Kien Wei Siah, Chi Heem Wong, Jerry Gupta, Andrew W. Lo
medRxiv 2021.05.13.21256888; doi: https://doi.org/10.1101/2021.05.13.21256888
Digg logo Reddit logo Twitter logo Facebook logo Google logo LinkedIn logo Mendeley logo
Citation Tools
Patterns of Multimorbidity
Kien Wei Siah, Chi Heem Wong, Jerry Gupta, Andrew W. Lo
medRxiv 2021.05.13.21256888; doi: https://doi.org/10.1101/2021.05.13.21256888

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Public and Global Health
Subject Areas
All Articles
  • Addiction Medicine (215)
  • Allergy and Immunology (495)
  • Anesthesia (106)
  • Cardiovascular Medicine (1093)
  • Dentistry and Oral Medicine (195)
  • Dermatology (141)
  • Emergency Medicine (274)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (499)
  • Epidemiology (9757)
  • Forensic Medicine (5)
  • Gastroenterology (480)
  • Genetic and Genomic Medicine (2303)
  • Geriatric Medicine (221)
  • Health Economics (462)
  • Health Informatics (1553)
  • Health Policy (732)
  • Health Systems and Quality Improvement (602)
  • Hematology (236)
  • HIV/AIDS (501)
  • Infectious Diseases (except HIV/AIDS) (11631)
  • Intensive Care and Critical Care Medicine (616)
  • Medical Education (236)
  • Medical Ethics (67)
  • Nephrology (256)
  • Neurology (2139)
  • Nursing (134)
  • Nutrition (335)
  • Obstetrics and Gynecology (426)
  • Occupational and Environmental Health (517)
  • Oncology (1172)
  • Ophthalmology (363)
  • Orthopedics (128)
  • Otolaryngology (220)
  • Pain Medicine (145)
  • Palliative Medicine (50)
  • Pathology (309)
  • Pediatrics (694)
  • Pharmacology and Therapeutics (298)
  • Primary Care Research (265)
  • Psychiatry and Clinical Psychology (2172)
  • Public and Global Health (4645)
  • Radiology and Imaging (775)
  • Rehabilitation Medicine and Physical Therapy (455)
  • Respiratory Medicine (623)
  • Rheumatology (274)
  • Sexual and Reproductive Health (225)
  • Sports Medicine (208)
  • Surgery (250)
  • Toxicology (43)
  • Transplantation (120)
  • Urology (94)