Comorbidity patterns have become a major source of information to explore shared mechanisms of pathogenesis between disorders. In hypothesis-free exploration of comorbid conditions, disease-disease networks are usually identified by pairwise methods. However, interpretation of the results is hindered by several confounders. In particular a very large number of pairwise associations can arise indirectly through other comorbidity associations and they increase exponentially with the increasing breadth of the investigated diseases. To investigate and filter this effect, we computed and compared pairwise approaches with a systems-based method, which constructs a sparse Bayesian direct multimorbidity map (BDMM) by systematically eliminating disease-mediated comorbidity relations. Additionally, focusing on depression-related parts of the BDMM, we evaluated correspondence with results from logistic regression, text-mining and molecular-level measures for comorbidities such as genetic overlap and the interactome-based association score. We used a subset of the UK Biobank Resource, a cross-sectional dataset including 247 diseases and 117,392 participants who filled out a detailed questionnaire about mental health. The sparse comorbidity map confirmed that depressed patients frequently suffer from both psychiatric and somatic comorbid disorders. Notably, anxiety and obesity show strong and direct relationships with depression. The BDMM identified further directly co-morbid somatic disorders, e. g. irritable bowel syndrome, fibromyalgia, or migraine. Using the subnetwork of depression and metabolic disorders for functional analysis, the interactome-based system-level score showed the best agreement with the sparse disease network. This indicates that these epidemiologically strong disease-disease relations have improved correspondence with expected molecular-level mechanisms. The substantially fewer number of comorbidity relations in the BDMM compared to pairwise methods implies that biologically meaningful comorbid relations may be less frequent than earlier pairwise methods suggested. The computed interactive comprehensive multimorbidity views over the diseasome are available on the web at Co=MorNet: bioinformatics. mit. bme. hu/UKBNetworks.
It has long been recognised that medical disorders frequently co-occur in the same individual [1] but the significance of comorbidity in revealing shared mechanisms of pathogenesis and outcome is a more recent realisation [2–4]. For a given disease or for a focused disease group, the exploration of comorbidities is largely hypothesis driven, together with the cautious selection and management of potential confounders [4–7]. The availablility of large health data sets with full multimorbidity information provides an unprecedented opportunity to understand the overall network of dependencies underpinning complex multimorbidities. These multivariate dependencies in turn become new targets for drug development and other therapies for multimorbid conditions, particularly relevant in aging societies [8–10]. However, the dissection of comorbidity relations is hindered by myriads of confounding factors [11–13]. Following the characterization from Bagley et al. [12], epidemiological co-occurrences can arise through different routes: 1) shared genetic background, 2) disease interactions (a disorder directly causes another), 3) common environmental cause and 4) different biases (diagnosis artifacts, selection biases). Earlier diseasome-wide works focused on the exploration of shared genetic background (1) behind comorbidities [2,3, 12,13] and the underlying molecular networks [3,14–18]. These works relied on pairwise comorbid relations partly controlled for potential confounding factors such as age (for controlling with disease onset see e. g. [2], for incidence-based control see e. g. [12]). However, these approaches do not address the issue of apparent comorbidity mediated by intervening associations with other diseases; a problem of indirect relations that is already attracting attention in other areas of network science [19]. In this paper we demonstrate that probabilistic graphical models (PGMs) in the Bayesian statistical framework provide a principled, unified solution for filtering such disease-mediated indirect relations, for correcting for potential external confounders and for coping with limitations and uncertainty of the data. Specifically, we construct sparse multimorbidity maps by applying PGMs for all diseases, i. e. for the whole diseasome. To our knowledge, this method has not been applied for the diseasome so far, despite the unique ability of PGMs to represent maximally sparse models, demonstrated on, for example genomic datasets [19–25]. Our diseasome-wide evaluations show that this approach efficiently scores and discriminates direct and disease-mediated indirect comorbidity relations and has resulted in a loss of more than 80% of comorbidity relations from prevailing pairwise methods. We made a more detailed investigation of BDMMs in the subset of psychiatric and metabolic disorders of the diseaseome. We focused on depression, which is a common psychiatric disorder with a complex neurobiological and psychosocial background [26,27], with approximately 10% prevalence worldwide, and according to forecasts depression will be the largest contributor to the disease burden in the middle- and high-income countries by 2030 according to the World Health Organization [28–30]. Many epidemiologic studies have reported high comorbidity between mental illnesses [5,31], which was partly explained by shared heritability between psychiatric disorders [32,33]. Less is known about the complex biopsychosocial mechanisms which underlie associations between somatic and common psychiatric disorders: depression frequently co-occurs with a wide range of somatic disorders, for example with migraine [34], with other disorders causing chronic pain [35,36], and with cardio-metabolic syndromes [37]. It has been also demonstrated that patients with depression have increased number of diagnosed disorders compared to non-depressed patients [5,38] and depression worsens the treatment outcome of the comorbid conditions [39] and is an independent predictor of increased mortality rate [40,41]. Therefore besides the exploration of further comorbidities of depression, it is equally important to discriminate its comorbidities as direct and indirect comorbidities. Such discrimination could reveal more specific pathophysiological subgroups of this heterogeneous condition and thus transform the power of genetic and epidemiological studies to advance precision medicine in psychiatry and metabolic disorders. We also evaluated the correspondence of BDMMs with molecular-level measures and relations, such as genetic overlap and the interactome-based association score for a depression-related subset of diseases. Focusing on depression we identified a direct multimorbid neighbourhood and confirmed that direct comorbidities correspond to direct relationships in the molecular interactome.
To explore the direct and indirect status of comorbidities in the diseasome and to investigate the effect of filtering disease-mediated, indirect comorbidity relations by BDMMs, we computed multiple pairwise measures. The prevailing non-systems based approaches usually use a pairwise measure such as odds ratio (OR), Pearson’s correlation coefficient (Φ) or logistic regression to determine the epidemiologic relationship of disorders [2,3]. We computed these most often used statistical descriptors of comorbidity for all the investigated diseases (see S2 Dataset for all computed pairwise measures for all possible disease pairs). We also computed BDMMs and cross-compared with the prevailing pairwise approaches [3,11,12,18]. The Fig 3 illustrates the sparsity of the direct (systems-based) map compared to a non-systems based network over the diseasome. The BDMM approach resulted in 320 direct connections (BDMM edge posterior > 0. 95) whereas applying the χ2 independence test 1714 disorder-disorder relations have significant p-values with a threshold of 0. 05 after Bonferroni correction. We also investigated the sufficiency of the sample size on different approaches. Fig 4 shows characteristics of the most significant results in each approaches. We transformed the different scores to the [0,1] interval to make them comparable (see Materials and methods). Descriptors and statistics are available on the web (Co=MorNet: bioinformatics. mit. bme. hu/UKBNetworks), also as an interactive tool to visualise networks of direct and mediated connections of selected diseases. The examination of the filtering capacity of the BDMM approach confirms that the BDMM edge posteriors strongly differentiate the direct connections from the mediated ones, e. g. there are only a few disease pairs which have a BDMM edge posterior between 0. 05 and 0. 95 (Fig 5). For the definition of BDMM edge and BDMM (structural) association relations see Table S1 in S1 Appendix. On Fig 5, we mark the direct comorbid connections with high BDMM edge posteriors (> 0. 95) in red, showing that all of them have significant Bonferroni-corrected χ2 p-values. Additionally, focusing on a subset of relations with BDMM edge posterior less than 0. 95, we examined the BDMM association posteriors versus the parametric association (see Fig 6). It confirms that almost all such strongly significant parametric association has high BDMM association posterior. Note, that relations with high BDMM edge posterior have high BDMM association posteriors as well. Furthermore, Fig 6 also illustrates that BDMM association posterior indicates many more structurally distant -presumably parametrically weaker- associations which cannot be inferred by parametric association tests such as Bonferroni-corrected χ2 p-values. To evaluate the genetic relevance of direct versus indirect comorbidity relations, we extended the epidemiological-level analysis with molecular-level approaches using a genetic overlap [3] and the interactome-based separation score [16]. This analysis included depression, metabolic syndromes and hypertension, see Fig 7, S2 Fig and Table S3 in S1 Appendix. The genetic overlap measures for comorbidities were computed based on manually curated databases, using the DisGeNet [44] and NCBI PheGenI [45] (for details, see S1 Appendix, for a related earlier work, see [2]). The interactome-based connectivity/separation scores for comorbidities were computed by the supplied method and data from Menche et al. [16] (see S1 Appendix). Beside of the investigation of direct and indirect status of comorbidities and the cross-comparison of pairwise and BDMM measures, we performed a more detailed medical evaluation of BDMMs on psychiatric and metabolic disorders, especially on depression and its comorbidities.
First of all, results indicate that the UK Biobank dataset is sufficiently large for the construction of BDMMs for this variable set, as the BDMM edge posteriors are peaked at 0 and 1 (see upper part in Fig 5), indeed, with the thresholds 0. 05 and 0. 95 we can efficiently separate the statistically significant comorbid relations to direct (co=morbid) and indirect relations. Notably, BDMMs eliminated more than 80% of comorbidity relations as indirect ones (320 direct connections from 1714 candidate relations). Interestingly, a recent work using a priori disease categories for restricting comorbidity relations, controlling for confounding with incidence characteristics and using two independent data sets similarly reported nearly 90% elimination ratio [12]. Note that all co=morbid connections were also confirmed by standard statistical methods (the right-lower quadrant is empty in Fig 5), which implies the technical condition that the distribution of the BDMM posteriors is stable [48–50]. In Fig 6 we further evaluated the connections with significant Bonferroni-corrected χ2 p-values but below the 0. 95 edge posterior threshold. This shows that most of these connections have high BDMM structural association posterior, which suggest that BDMM indeed filtered mediated and confounded relations. There are 31 connections which have a significant pairwise association score but no structural association (BDDM structural association posterior < 0. 1). These disorder pairs rarely occur together in patients (mean co-occurrence: 8. 45, with standard deviation: 4. 9, and quantiles: 5,7, 10. 5 for the 25%, 50% and 75% respectively). In case of the χ2 test we applied Yates’ continuity correction but for such weak connections even with the large UK Biobank dataset the BDMM approach was not able to catch that weak dependency structure. Our results demonstrated that the interactome-based score provided similar maps as BDMM co=morbid diseases, in sharp contrast to the associative genetic overlap scores which followed the pattern of the pairwise disease relative risk (see Fig 7, S2 Fig and Table S3 in S1 Appendix). Note that the interactome-based score and co=morbidity are analogous as both use a systems-based approach on different levels (molecular and epidemiologic level respectively). For detailed description of the molecular level methods and results see S1 Appendix and [2,16]. The high comorbidity between mood disorders and anxiety or stress related disorders is well known, and twin studies suggested that these comorbidities originated mainly from shared genetic risk factors [52,53]. Our results showed another expected aspect, namely that this relationship is independent of the order of the onset of these disorders. This observation is in line with a longitudinal study which showed that generalised anxiety disorder (GAD) and major depressive disorder (MDD) are strongly comorbid with an equal probability of GAD or MDD occurring first or simultaneously suggesting they might not be distinct disorders [54]. Although overlapping genetic risk factors for anxiety and depression have not yet been identified, common genetic vulnerability has been found for other comorbid psychiatric disorders [32,33,55]. Our BDMM further indicate strong and stable co=morbidity between depression, anxiety, stress, postnatal depression and nervous breakdown, pointing toward interactome-level overlaps; this reinforces the need to find potential common biological mechanisms [56]. Epidemiologic studies repeatedly report high comorbidity between depression and metabolic disorders [57], depression and diabetes [58], depression and cardiovascular disorders [37], depression and hypertension [59,60], and depression and obesity [61,62]. However, there have been several contradictory results, and this suggests a more complex relationship. Indeed, recent GWAS results found no shared genetic risk between these disorders and depression [55]. Nevertheless, based on the UK Biobank cohort data, obesity was co=morbid with depression while cardiovascular disorders, hypertension, high cholesterol and diabetes, including type 2 diabetes were only indirectly related to depression. When we excluded occurrences after the onset of depression the direct relationship between obesity and depression remained as expected but entirely new links with high posterior probability emerged suggesting a strong relationship between the consequences of metabolic syndromes and depression. Studies of the genetic relationship between obesity and depression suggest that atypical depression, characterised by increase in appetite and weight, is associated with genetic risk factors and polygenic risk scores of increased BMI and triglycerides, while typical depression, with decreased appetite and weight, show more similarities with other psychiatric disorders [63,64]. Thus, in line with our results comorbid obesity and metabolic disorders may identify a specific subtype of depression with a distinct biological background. Caution is required in inferring shared biology for co=morbid disorders given our current lack of knowledge about relevant GxE interactions. Although metabolic disorders are only co=morbid with depression when they precede it, a variety of non-biological associations of mediators may be at play. We cannot currently exclude the possibility that lifestyle factors, such as diet, physical activity and stress, or medication used to treat hypertension, hypercholesterolaemia and obesity may contribute to the later development of depression [65]. As a specific example, a previous study demonstrated that current psychological distress amplified the effect of genetic risk of high BMI [66]. Patients with increased genetic risk to become overweight showed worse physical outcome (higher BMI), and quite probably more comorbid psychological symptoms, when life stress was present. Furthermore, it has been reported that statins, drugs with cholesterol-lowering effect, have antidepressant effect in patients with comorbid depression and coronary artery disease while the same drugs can have pro-depressive effect or no effect on depression when comorbidities and depression subtypes were not taken into account [67]. Migraine [34], IBS [68], FM [69], and CFS [70] are highly comorbid with depression based on epidemiologic studies. It is therefore puzzling that they involve different etiological mechanisms. In addition, their symptoms often overlap making it difficult to apply diagnostic categories. We found that these disorders were not relevant when they occurred before depression but were highly co=morbid in the full analysis. The probable explanation is that in general, these disorders are related to consequences of depression and only specific subtypes of these disorders can be expected to have causal relations, e. g. shared biological background with depression. For example, a genetic risk score analysis demonstrated that migraine with comorbid depression was more genetically related to depression than to pure migraine, which suggests that migraine might develop as a consequence of different polygenic backgrounds [71]. Similarly, a large general population cohort study confirmed that FM, CFS and IBS increase the odds of depression and anxiety but that most patients who suffer from FM, CFS and IBS have no mood or anxiety disorder [72]. One of the main limitations is that all disorders were self-reported, although trained nurses evaluated and corrected all entries during face-to-face interviews. The second one is that the applied treatments or medications were not included in the analysis which could highlight comorbidities due to the side effect of treatments. We will address this problem in follow-up studies. Note that we only used a subset of the UK Biobank dataset selecting those participants who filled out the Mental Health Questionnaire and provided online dietary information, which may introduce confounding through selection bias. However, limiting our study to this subpopulation enabled us to test different definitions of depression and will allow us to connect this comorbidity network to relevant environmental risk factors. The use of large-scale health data sets, such as the UK Biobank dataset hold the promise of complementing and guiding the molecular level research of complex diseases. Adopting an intermediary approach between statistical association analysis and causal discovery we investigated the use of Bayesian networks in the Bayesian model averaging framework to explore direct probabilistic relations with respect to a given set of variables, i. e. to eliminate confounders and mediatory effects by a systems-based approach. We demonstrated the applicability of BDMMs, especially their principled capability of discriminating direct and indirect comorbidities. In summary, PGMs offer maximally sparse dependency models and utilize the omic nature of the epidemiologic data jointly modelling all the morbidities; while the Bayesian approach through posteriors provides an explicit representation for the uncertainties in a dataset. Thus the Bayesian direct morbidity maps provide sparse, systems-based, omic-wide perspectives. From a clinical perspective, our results also highlight that the direct and indirect subtypes of comorbidities support a finer biological interpretation, namely an interactome-based detailed interpretation using molecular mechanisms corresponding to direct relations, whereas genetic overlap using associative gene sets may only reflect indirect comorbidities. In addition, re-running the analysis by including only instances of disorders which preceded depression, we delineated comorbid disorders of depression with more refined causal roles that could specify subgroups of depressed patients with more homogenous background. The investigation of Bayesian direct morbidity maps also demonstrated, that even large-scale datasets such as UK Biobank, are still limited for non-ambiguous identification of complex dependency patterns such as multimorbidities [73]. However, the applied Bayesian statistical framework offers an automated, normative solution for the multiple hypothesis testing problem and the application of probabilistic graphical models in the Bayesian framework supports the versatile post-processing of the results and their efficient communication and sharing. The results of our research highlight the advantages of Bayesian systems-based modelling, which could be vital to cope with the growing heterogeneity of new health data sets containing full personal genetic information with high dimensional data about lifestyle, environmental factors and sequential decisions on drug therapies [8,74].
In the present study we used the UK Biobank (http: //www. ukbiobank. ac. uk/ and [75] cohort where subjects’ chronic illness history together with onset age were ascertained by trained nurses during face-to-face interviews and were processed by experts resulting in 525 different disease categories. The investigated subset in this study consisted of 117,392 participants (female: 64,320; male: 53,072) who provided the extended Mental Health Questionnaire (http: //biobank. ctsu. ox. ac. uk/crystal/label. cgi? id=100060) and the online diet questionnaire data together with the extensive baseline dataset. We used the UK Biobank original disease categories with at least 1‰ prevalence in the selected subset, which resulted in n = 247 diseases including depression (n = 6040). In addition, we coded obesity in cases where BMI were equal or greater than 30kg/m2, for further analyses. For statistical analysis, sex was included into the data set, and age was binned into 3 equal frequency categories with thresholds 60 and 68 years. Then we applied the different pairwise measures and logistic regression together with Bayesian systems-based modelling to compare the models computed on these datasets (see below and in S1 Appendix). To investigate the effect of disease onset, self-reported disease onset data was used to filter the dataset. 6,040 patients affected by depression, provided onset data whose comorbid illnesses were eliminated if they occurred after the onset of depression. After removal of diseases with prevalence less than 1‰ the dataset contained 241 diseases. We extended the dataset with sex, age and BMI-based obesity. The data were analysed using same statistical methods as with the non-filtered dataset. To test the stability of comorbid relationships with depression we also used an alternative depression definition instead of self-reported depressive disorder. Depression and its severity was defined by the Mental Health Questionnaire data [46], for definition see S1 Appendix. These alternative depression categories were analysed with Bayesian systems-based modelling. We applied text-mining and conventional statistical methods to explore comorbid relations, see S1 Appendix. For these computations we used in-house written R scripts together with the statistical programs included in the stats package of R [76]. To overcome the limitations of these conventional methods, we applied a Bayesian network Markov Chain Monte Carlo (BN-MCMC) method to explore the overall system of dependencies-independencies, visualized as an undirected graph with weighted edges [20,21,23–25,42]. The weighted edges correspond the a posteriori probabilities (Pr) of direct, nonmediated “co=morbidity” relations, the weights are in the [0,1] interval and the higher values show stronger relationship. The systems-based approach using Bayesian networks prunes the indirect, mediated connections between morbidities, thus resulting in a sparse co=morbidity map compared to pairwise association networks (for detailed description of the method and for further types of dependency relations, see S1 Appendix and Table S1 in S1 Appendix).
Depression is one of the most common of psychiatric disorders and its causation is correspondingly multifactorial, complex and heterogeneous. It occurs in combination with a number of physical illnesses far more commonly than expected by chance. Such comorbidities may be important clues pointing to shared environmental and genetic risk factors and could identify different causal types of depression. However, a method is still needed to weed out statistically significant pairings that nevertheless arise through indirect routes involving comorbidities between other diseases. We examined the pairwise associations among 247 diseases of 117,392 participants recorded in the UK Biobank database. We found that the great majority of disease associations were indirect consequences of a sparse network of ‘direct’ comorbidities (‘sparse diseaseome’) constructed using probabilistic graphical models (PGMs) within the Bayesian statistical framework. In a depression-related subset of illnesses, we found that several pairwise associations of depression were indirect and due to their comorbidities with obesity which had a strong direct connection with depression. Furthermore, the direct comorbidities in a depression-related subset of disorders, but not the pairwise associations, strongly mapped onto an underlying molecular network (‘interactome’) suggesting that this approach significantly improved correspondence with molecular reality.
body weight
medicine and health sciences
pathology and laboratory medicine
migraine
headaches
physiological parameters
signs and symptoms
obesity
mood disorders
genetic epidemiology
epidemiology
hypertension
blood pressure
mental health and psychiatry
metabolic disorders
clinical genetics
diagnostic medicine
physiology
genetics
biology and life sciences
depression
vascular medicine
genetics of disease
Comorbidities in the diseasome are more apparent than real: What Bayesian filtering reveals about the comorbidities of depression