Analysis of families in the multiple autoimmune disease genetics consortium (MADGC) collection: the PTPN22 620W allele associates with multiple autoimmune phenotypes

Autoimmune disorders constitute a diverse group of phenotypes with overlapping features and a tendency toward familial aggregation. It is likely that common underlying genes are involved in these disorders. Until very recently, no specific alleles--aside from a few common human leukocyte antigen class II genes--had been identified that clearly associate with multiple different autoimmune diseases. In this study, we describe a unique collection of 265 multiplex families assembled by the Multiple Autoimmune Disease Genetics Consortium (MADGC). At least two of nine "core" autoimmune diseases are present in each of these families. These core diseases include rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), type 1 diabetes (T1D), multiple sclerosis (MS), autoimmune thyroid disease (Hashimoto thyroiditis or Graves disease), juvenile RA, inflammatory bowel disease (Crohn disease or ulcerative colitis), psoriasis, and primary Sjogren syndrome. We report that a recently described functional single-nucleotide polymorphism (rs2476601, encoding R620W) in the intracellular tyrosine phosphatase (PTPN22) confers risk of four separate autoimmune phenotypes in these families: T1D, RA, SLE, and Hashimoto thyroiditis. MS did not show association with the PTPN22 risk allele. These findings suggest a common underlying etiologic pathway for some, but not all, autoimmune disorders, and they suggest that MS may have a pathogenesis that is distinct from RA, SLE, and T1D. DNA and clinical data for the MADGC families are available to the scientific community; these data will provide a valuable resource for the dissection of the complex genetic factors that underlie the various autoimmune phenotypes.     

Selective IgA deficiency in autoimmune diseases

Selective immunoglobulin A deficiency (IgAD) is the most common primary immunodeficiency in Caucasians. It has previously been suggested to be associated with a variety of concomitant autoimmune diseases. In this review, we present data on the prevalence of IgAD in patients with Graves disease (GD), systemic lupus erythematosus (SLE), type 1 diabetes (T1D), celiac disease (CD), myasthenia gravis (MG) and rheumatoid arthritis (RA) on the basis of both our own recent large-scale screening results and literature data. Genetic factors are important for the development of both IgAD and various autoimmune disorders, including GD, SLE, T1D, CD, MG and RA, and a strong association with the major histocompatibility complex (MHC) region has been reported. In addition, non-MHC genes, such as interferon-induced helicase 1 (IFIH1) and c-type lectin domain family 16, member A (CLEC16A), are also associated with the development of IgAD and some of the above diseases. This indicates a possible common genetic background. In this review, we present suggestive evidence for a shared genetic predisposition between these disorders.     

A comprehensive analysis of shared loci between systemic lupus erythematosus (SLE) and sixteen autoimmune diseases reveals limited genetic overlap

In spite of the well-known clustering of multiple autoimmune disorders in families, analyses of specific shared genes and polymorphisms between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) have been limited. Therefore, we comprehensively tested autoimmune variants for association with SLE, aiming to identify pleiotropic genetic associations between these diseases. We compiled a list of 446 non–Major Histocompatibility Complex (MHC) variants identified in genome-wide association studies (GWAS) of populations of European ancestry across 17 ADs. We then tested these variants in our combined Caucasian SLE cohorts of 1,500 cases and 5,706 controls. We tested a subset of these polymorphisms in an independent Caucasian replication cohort of 2,085 SLE cases and 2,854 controls, allowing the computation of a meta-analysis between all cohorts. We have uncovered novel shared SLE loci that passed multiple comparisons adjustment, including the VTCN1 (rs12046117, P = 2.02×10⁻⁰⁶) region. We observed that the loci shared among the most ADs include IL23R, OLIG3/TNFAIP3, and IL2RA. Given the lack of a universal autoimmune risk locus outside of the MHC and variable specificities for different diseases, our data suggests partial pleiotropy among ADs. Hierarchical clustering of ADs suggested that the most genetically related ADs appear to be type 1 diabetes with rheumatoid arthritis and Crohn''s disease with ulcerative colitis. These findings support a relatively distinct genetic susceptibility for SLE. For many of the shared GWAS autoimmune loci, we found no evidence for association with SLE, including IL23R. Also, several established SLE loci are apparently not associated with other ADs, including the ITGAM-ITGAX and TNFSF4 regions. This study represents the most comprehensive evaluation of shared autoimmune loci to date, supports a relatively distinct non–MHC genetic susceptibility for SLE, provides further evidence for previously and newly identified shared genes in SLE, and highlights the value of studies of potentially pleiotropic genes in autoimmune diseases.     

Integrated Enrichment Analysis of Variants and Pathways in Genome-Wide Association Studies Indicates Central Role for IL-2 Signaling Genes in Type 1 Diabetes,and Cytokine Signaling Genes in Crohn's Disease

Pathway analyses of genome-wide association studies aggregate information over sets of related genes, such as genes in common pathways, to identify gene sets that are enriched for variants associated with disease. We develop a model-based approach to pathway analysis, and apply this approach to data from the Wellcome Trust Case Control Consortium (WTCCC) studies. Our method offers several benefits over existing approaches. First, our method not only interrogates pathways for enrichment of disease associations, but also estimates the level of enrichment, which yields a coherent way to promote variants in enriched pathways, enhancing discovery of genes underlying disease. Second, our approach allows for multiple enriched pathways, a feature that leads to novel findings in two diseases where the major histocompatibility complex (MHC) is a major determinant of disease susceptibility. Third, by modeling disease as the combined effect of multiple markers, our method automatically accounts for linkage disequilibrium among variants. Interrogation of pathways from eight pathway databases yields strong support for enriched pathways, indicating links between Crohn''s disease (CD) and cytokine-driven networks that modulate immune responses; between rheumatoid arthritis (RA) and “Measles” pathway genes involved in immune responses triggered by measles infection; and between type 1 diabetes (T1D) and IL2-mediated signaling genes. Prioritizing variants in these enriched pathways yields many additional putative disease associations compared to analyses without enrichment. For CD and RA, 7 of 8 additional non-MHC associations are corroborated by other studies, providing validation for our approach. For T1D, prioritization of IL-2 signaling genes yields strong evidence for 7 additional non-MHC candidate disease loci, as well as suggestive evidence for several more. Of the 7 strongest associations, 4 are validated by other studies, and 3 (near IL-2 signaling genes RAF1, MAPK14, and FYN) constitute novel putative T1D loci for further study.     

Integrative Pathway-Based Approach for Genome-Wide Association Studies: Identification of New Pathways for Rheumatoid Arthritis and Type 1 Diabetes

Genome-wide association studies (GWAS) led to the identification of numerous novel loci for a number of complex diseases. Pathway-based approaches using genotypic data provide tangible leads which cannot be identified by single marker approaches as implemented in GWAS. The available pathway analysis approaches mainly differ in the employed databases and in the applied statistics for determining the significance of the associated disease markers.So far, pathway-based approaches using GWAS data failed to consider the overlapping of genes among different pathways or the influence of protein–interactions. We performed a multistage integrative pathway (MIP) analysis on three common diseases - Crohn''s disease (CD), rheumatoid arthritis (RA) and type 1 diabetes (T1D) - incorporating genotypic, pathway, protein- and domain-interaction data to identify novel associations between these diseases and pathways. Additionally, we assessed the sensitivity of our method by studying the influence of the most significant SNPs on the pathway analysis by removing those and comparing the corresponding pathway analysis results. Apart from confirming many previously published associations between pathways and RA, CD and T1D, our MIP approach was able to identify three new associations between disease phenotypes and pathways. This includes a relation between the influenza-A pathway and RA, as well as a relation between T1D and the phagosome and toxoplasmosis pathways. These results provide new leads to understand the molecular underpinnings of these diseases.The developed software herein used is available at http://www.cogsys.cs.uni-tuebingen.de/software/GWASPathwayIdentifier/index.htm.     

Identification of Potential Pleiotropic Genes for Immune and Skeletal Diseases Using Multivariate MetaCCA Analysis

BackgroundImmune and skeletal systems physiologically and pathologically interact with each other. Immune and skeletal diseases may share potential pleiotropic genetics factors, but the shared specific genes are largely unknown.ObjectiveThis study aimed to investigate the overlapping genetic factors between multiple diseases (including rheumatoid arthritis (RA), psoriasis, osteoporosis, osteoarthritis, sarcopenia, and fracture).MethodsThe canonical correlation analysis (metaCCA) approach was used to identify the shared genes for six diseases by integrating genome-wide association study (GWAS)-derived summary statistics. The versatile Gene-based Association Study (VEGAS2) method was further applied to refine and validate the putative pleiotropic genes identified by metaCCA.ResultsAbout 157 (p<8.19E-6), 319 (p<3.90E-6), and 77 (p<9.72E-6) potential pleiotropic genes were identified shared by two immune diseases, four skeletal diseases, and all of the six diseases, respectively. The top three significant putative pleiotropic genes shared by both immune and skeletal diseases, including HLA-B, TSBP1, and TSBP1-AS1 (p<E-300), were located in the major histocompatibility complex (MHC) region. Nineteen of 77 putative pleiotropic genes identified by metaCCA analysis were associated with at least one disease in the VEGAS2 analysis. Specifically, the majority (18) of these 19 putative validated pleiotropic genes were associated with RA.ConclusionThe metaCCA method identified some pleiotropic genes shared by the immune and skeletal diseases. These findings help to improve our understanding of the shared genetic mechanisms and signaling pathways underlying immune and skeletal diseases.     

Pathway Analysis of GWAS Provides New Insights into Genetic Susceptibility to 3 Inflammatory Diseases

Although the introduction of genome-wide association studies (GWAS) have greatly increased the number of genes associated with common diseases, only a small proportion of the predicted genetic contribution has so far been elucidated. Studying the cumulative variation of polymorphisms in multiple genes acting in functional pathways may provide a complementary approach to the more common single SNP association approach in understanding genetic determinants of common disease. We developed a novel pathway-based method to assess the combined contribution of multiple genetic variants acting within canonical biological pathways and applied it to data from 14,000 UK individuals with 7 common diseases. We tested inflammatory pathways for association with Crohn''s disease (CD), rheumatoid arthritis (RA) and type 1 diabetes (T1D) with 4 non-inflammatory diseases as controls. Using a variable selection algorithm, we identified variants responsible for the pathway association and evaluated their use for disease prediction using a 10 fold cross-validation framework in order to calculate out-of-sample area under the Receiver Operating Curve (AUC). The generalisability of these predictive models was tested on an independent birth cohort from Northern Finland. Multiple canonical inflammatory pathways showed highly significant associations (p 10⁻³–10⁻²⁰) with CD, T1D and RA. Variable selection identified on average a set of 205 SNPs (149 genes) for T1D, 350 SNPs (189 genes) for RA and 493 SNPs (277 genes) for CD. The pattern of polymorphisms at these SNPS were found to be highly predictive of T1D (91% AUC) and RA (85% AUC), and weakly predictive of CD (60% AUC). The predictive ability of the T1D model (without any parameter refitting) had good predictive ability (79% AUC) in the Finnish cohort. Our analysis suggests that genetic contribution to common inflammatory diseases operates through multiple genes interacting in functional pathways.     

Comorbidities of Psoriasis - Exploring the Links by Network Approach

Increasing epidemiological studies in patients with psoriasis report the frequent occurrence of one or more associated disorders. Psoriasis is associated with multiple comorbidities including autoimmune disease, neurological disorders, cardiometabolic diseases and inflammatory-bowel disease. An integrated system biology approach is utilized to decipher the molecular alliance of psoriasis with its comorbidities. An unbiased integrative network medicine methodology is adopted for the investigation of diseasome, biological process and pathways of five most common psoriasis associated comorbidities. A significant overlap was observed between genes acting in similar direction in psoriasis and its comorbidities proving the mandatory occurrence of either one of its comorbidities. The biological processes involved in inflammatory response and cell signaling formed a common basis between psoriasis and its associated comorbidities. The pathway analysis revealed the presence of few common pathways such as angiogenesis and few uncommon pathways which includes CCKR signaling map and gonadotrophin-realising hormone receptor pathway overlapping in all the comorbidities. The work shed light on few common genes and pathways that were previously overlooked. These fruitful targets may serve as a starting point for diagnosis and/or treatment of psoriasis comorbidities. The current research provides an evidence for the existence of shared component hypothesis between psoriasis and its comorbidities.     

IL2RA Genetic Heterogeneity in Multiple Sclerosis and Type 1 Diabetes Susceptibility and Soluble Interleukin-2 Receptor Production

Multiple sclerosis (MS) and type 1 diabetes (T1D) are organ-specific autoimmune disorders with significant heritability, part of which is conferred by shared alleles. For decades, the Human Leukocyte Antigen (HLA) complex was the only known susceptibility locus for both T1D and MS, but loci outside the HLA complex harboring risk alleles have been discovered and fully replicated. A genome-wide association scan for MS risk genes and candidate gene association studies have previously described the IL2RA gene region as a shared autoimmune locus. In order to investigate whether autoimmunity risk at IL2RA was due to distinct or shared alleles, we performed a genetic association study of three IL2RA variants in a DNA collection of up to 9,407 healthy controls, 2,420 MS, and 6,425 T1D subjects as well as 1,303 MS parent/child trios. Here, we report “allelic heterogeneity” at the IL2RA region between MS and T1D. We observe an allele associated with susceptibility to one disease and risk to the other, an allele that confers susceptibility to both diseases, and an allele that may only confer susceptibility to T1D. In addition, we tested the levels of soluble interleukin-2 receptor (sIL-2RA) in the serum from up to 69 healthy control subjects, 285 MS, and 1,317 T1D subjects. We demonstrate that multiple variants independently correlate with sIL-2RA levels.     

The Crosstalk of Pathways Involved in Immune Response Maybe the Shared Molecular Basis of Rheumatoid Arthritis and Type 2 Diabetes

Rheumatoid arthritis (RA) and Type 2 diabetes (T2D) are both systemic diseases linked with altered immune response, moderate mortality when present together. The treatment for both RA and T2D are not satisfied, partly because of the linkage between them has not yet been appreciated. A comprehensive study for the potential associations between the two disorders is needed. In this study, we used RNA sequencing to explore the differently expressed genes (DEGs) in peripheral blood mononuclear cells (PBMC) of 10 RA and 10 T2D patients comparing with 10 healthy volunteers (control). We used bioinformatics analysis and the Ingenuity Pathways Analysis (IPA) to predict the commonalities on signaling pathways and molecular networks between those two diseases. 212 DEGs in RA and 114 DEGs in T2D patients were identified compared with healthy controls, respectively. 32 DEGs were shared between the two comparisons. The top 10 shared pathways interacted in cross-talking networks, regulated by 5 shared predicted upstream regulators, leading to the activated immune response were explored, which was considered as partly of the association mechanism of this two disorders. These discoveries would be considered as new understanding on the associations between RA and T2D, and provide novel treatment or prevention strategy.     

Autoimmune Disease Classification by Inverse Association with SNP Alleles

With multiple genome-wide association studies (GWAS) performed across autoimmune diseases, there is a great opportunity to study the homogeneity of genetic architectures across autoimmune disease. Previous approaches have been limited in the scope of their analysis and have failed to properly incorporate the direction of allele-specific disease associations for SNPs. In this work, we refine the notion of a genetic variation profile for a given disease to capture strength of association with multiple SNPs in an allele-specific fashion. We apply this method to compare genetic variation profiles of six autoimmune diseases: multiple sclerosis (MS), ankylosing spondylitis (AS), autoimmune thyroid disease (ATD), rheumatoid arthritis (RA), Crohn''s disease (CD), and type 1 diabetes (T1D), as well as five non-autoimmune diseases. We quantify pair-wise relationships between these diseases and find two broad clusters of autoimmune disease where SNPs that make an individual susceptible to one class of autoimmune disease also protect from diseases in the other autoimmune class. We find that RA and AS form one such class, and MS and ATD another. We identify specific SNPs and genes with opposite risk profiles for these two classes. We furthermore explore individual SNPs that play an important role in defining similarities and differences between disease pairs. We present a novel, systematic, cross-platform approach to identify allele-specific relationships between disease pairs based on genetic variation as well as the individual SNPs which drive the relationships. While recognizing similarities between diseases might lead to identifying novel treatment options, detecting differences between diseases previously thought to be similar may point to key novel disease-specific genes and pathways.     

Knowledge-based analysis of genetic associations of rheumatoid arthritis to inform studies searching for pleiotropic genes: a literature review and network analysis

IntroductionPleiotropy describes the genetic effect of a single gene on multiple phenotypic traits. Gene variants directly affect the normal processes of a series of physiological and biochemical reactions, and therefore cause a variety of diseases traits to be changed accordingly. Moreover, a shared genetic susceptibility mechanism may exist between different diseases. Therefore, shared genes, with pleiotropic effects, are important to understand the sharing pathogenesis and hence the mechanisms underlying comorbidity.MethodsIn this study, we proposed combining genome-wide association studies (GWAS) and public knowledge databases to search for potential pleiotropic genes associated with rheumatoid arthritis (RA) and eight other related diseases. Here, a GWAS-based network analysis is used to recognize risk genes significantly associated with RA. These RA risk genes are re-extracted as potential pleiotropic genes if they have been proved to be susceptible genes for at least one of eight other diseases in the OMIM or PubMed databases.ResultsIn total, we extracted 116 potential functional pleiotropic genes for RA and eight other diseases, including five hub pleiotropic genes, BTNL2, HLA-DRA, NOTCH4, TNXB, and C6orf10, where BTNL2, NOTCH4, and C6orf10 are novel pleiotropic genes identified by our analysis.ConclusionsThis study demonstrates that pleiotropy is a common property of genes associated with disease traits. Our results ascertained the shared genetic risk profiles that predisposed individuals to RA and other diseases, which could have implications for identification of molecular targets for drug development, and classification of diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0715-1) contains supplementary material, which is available to authorized users.     

Recent findings on genes associated with inflammatory disease

Inflammatory diseases encompass a variety of medical conditions. In this chapter, autoimmune diseases and allergic disorders will be our focus. The autoimmune diseases include organ-specific autoimmunities, such as type I diabetes mellitus and autoimmune thyroiditis (AITD), and organ non-specific disorders such as systemic lupus erythematosus (SLE). All of them seem to share aspects of aberrant immunologic tolerance toward self-antigens. Asthma and atopic diathesis are among the allergies. Crohn disease and SLE are relatively rare with a prevalence of 10-50 per 100,000, and rheumatoid arthritis (RA), psoriasis, AITD and asthma are commoner with a prevalence of 500 per 100,000 or much higher. The difference among ethnic groups is not prominent for rheumatoid arthritis, psoriasis, AITD or asthma, but Crohn disease and SLE affect some ethnic populations more than others. Although all of these disorders have some environmental component, asthma and atopy seem most affected by environmental factors, as is suggested by the significant increase in their incidence over the last several decades with changes in various environmental factors, especially in developed countries. Over the last 10 years, multiple linkage studies revealed many disease-linked loci throughout the genome with various consistencies. As implicated by some pathophysiological studies of inflammatory immune system related disorders, certain loci are involved in multiple disorders. In the following sections, reports on the identification of disease-associated genes or markers will be summarized for individual diseases (cytotoxic T lymphocyte-associated 4 (CTLA4), CARD15, DLG5, SLC22A4/A5, programmed cell death 1 (PDCD1), RUNX1, SLC9A3R1/NAT9, PADI4, ADAM33, DPP10, PHF11 and GPRA), followed by a discussion of the genes that have been implicated in multiple disorders.     

Defining the role of the MHC in autoimmunity: a review and pooled analysis

The major histocompatibility complex (MHC) is one of the most extensively studied regions in the human genome because of the association of variants at this locus with autoimmune, infectious, and inflammatory diseases. However, identification of causal variants within the MHC for the majority of these diseases has remained difficult due to the great variability and extensive linkage disequilibrium (LD) that exists among alleles throughout this locus, coupled with inadequate study design whereby only a limited subset of about 20 from a total of approximately 250 genes have been studied in small cohorts of predominantly European origin. We have performed a review and pooled analysis of the past 30 years of research on the role of the MHC in six genetically complex disease traits – multiple sclerosis (MS), type 1 diabetes (T1D), systemic lupus erythematosus (SLE), ulcerative colitis (UC), Crohn's disease (CD), and rheumatoid arthritis (RA) – in order to consolidate and evaluate the current literature regarding MHC genetics in these common autoimmune and inflammatory diseases. We corroborate established MHC disease associations and identify predisposing variants that previously have not been appreciated. Furthermore, we find a number of interesting commonalities and differences across diseases that implicate both general and disease-specific pathogenetic mechanisms in autoimmunity.     

Pervasive sharing of genetic effects in autoimmune disease

Genome-wide association (GWA) studies have identified numerous, replicable, genetic associations between common single nucleotide polymorphisms (SNPs) and risk of common autoimmune and inflammatory (immune-mediated) diseases, some of which are shared between two diseases. Along with epidemiological and clinical evidence, this suggests that some genetic risk factors may be shared across diseases-as is the case with alleles in the Major Histocompatibility Locus. In this work we evaluate the extent of this sharing for 107 immune disease-risk SNPs in seven diseases: celiac disease, Crohn's disease, multiple sclerosis, psoriasis, rheumatoid arthritis, systemic lupus erythematosus, and type 1 diabetes. We have developed a novel statistic for Cross Phenotype Meta-Analysis (CPMA) which detects association of a SNP to multiple, but not necessarily all, phenotypes. With it, we find evidence that 47/107 (44%) immune-mediated disease risk SNPs are associated to multiple-but not all-immune-mediated diseases (SNP-wise P(CPMA)<0.01). We also show that distinct groups of interacting proteins are encoded near SNPs which predispose to the same subsets of diseases; we propose these as the mechanistic basis of shared disease risk. We are thus able to leverage genetic data across diseases to construct biological hypotheses about the underlying mechanism of pathogenesis.     

A Possible Mechanism behind Autoimmune Disorders Discovered By Genome-Wide Linkage and Association Analysis in Celiac Disease

Celiac disease is a common autoimmune disorder characterized by an intestinal inflammation triggered by gluten, a storage protein found in wheat, rye and barley. Similar to other autoimmune diseases such as type 1 diabetes, psoriasis and rheumatoid arthritis, celiac disease is the result of an immune response to self-antigens leading to tissue destruction and production of autoantibodies. Common diseases like celiac disease have a complex pattern of inheritance with inputs from both environmental as well as additive and non-additive genetic factors. In the past few years, Genome Wide Association Studies (GWAS) have been successful in finding genetic risk variants behind many common diseases and traits. To complement and add to the previous findings, we performed a GWAS including 206 trios from 97 nuclear Swedish and Norwegian families affected with celiac disease. By stratifying for HLA-DQ, we identified a new genome-wide significant risk locus covering the DUSP10 gene. To further investigate the associations from the GWAS we performed pathway analyses and two-locus interaction analyses. These analyses showed an over-representation of genes involved in type 2 diabetes and identified a set of candidate mechanisms and genes of which some were selected for mRNA expression analysis using small intestinal biopsies from 98 patients. Several genes were expressed differently in the small intestinal mucosa from patients with celiac autoimmunity compared to intestinal mucosa from control patients. From top-scoring regions we identified susceptibility genes in several categories: 1) polarity and epithelial cell functionality; 2) intestinal smooth muscle; 3) growth and energy homeostasis, including proline and glutamine metabolism; and finally 4) innate and adaptive immune system. These genes and pathways, including specific functions of DUSP10, together reveal a new potential biological mechanism that could influence the genesis of celiac disease, and possibly also other chronic disorders with an inflammatory component.     

Principal Component Analysis Characterizes Shared Pathogenetics from Genome-Wide Association Studies

Genome-wide association studies (GWASs) have recently revealed many genetic associations that are shared between different diseases. We propose a method, disPCA, for genome-wide characterization of shared and distinct risk factors between and within disease classes. It flips the conventional GWAS paradigm by analyzing the diseases themselves, across GWAS datasets, to explore their “shared pathogenetics”. The method applies principal component analysis (PCA) to gene-level significance scores across all genes and across GWASs, thereby revealing shared pathogenetics between diseases in an unsupervised fashion. Importantly, it adjusts for potential sources of heterogeneity present between GWAS which can confound investigation of shared disease etiology. We applied disPCA to 31 GWASs, including autoimmune diseases, cancers, psychiatric disorders, and neurological disorders. The leading principal components separate these disease classes, as well as inflammatory bowel diseases from other autoimmune diseases. Generally, distinct diseases from the same class tend to be less separated, which is in line with their increased shared etiology. Enrichment analysis of genes contributing to leading principal components revealed pathways that are implicated in the immune system, while also pointing to pathways that have yet to be explored before in this context. Our results point to the potential of disPCA in going beyond epidemiological findings of the co-occurrence of distinct diseases, to highlighting novel genes and pathways that unsupervised learning suggest to be key players in the variability across diseases.