Skip to main content

Using hierarchical similarity to examine the genetics of Behçet’s disease



Behçet’s disease (BD) is a multisystem inflammatory disease that affects patients along the historic silk road. Thus far, the pathogenesis of the disease has proved elusive due to the complex genetic interactions of the disease. In this paper, we seek to clarify the genetic factors of the disease while also uncovering other diseases of interest that present with a similar genotype as BD.


To do this, we employ a convergent functional genomics approach by leveraging the hierarchical similarity tool available in Geneweaver. Through our analysis, we were able to ascertain 7 BD consensus genes and 16 autoimmune diseases with genetic overlap with BD. The results of our study will inform further research into the pathogenesis of Behçet’s disease.


Behçet’s disease (BD) is a multi-system inflammatory disease, demonstrating lesions of the mouth, the genitalia, and the eye [1, 2]. The disease is thought to be triggered in genetically predisposed individuals by a currently unknown combination of viral infection and environmental factors [3].

The genetics of BD are of particular interest. The HLA-B51 gene is the most closely associated risk factor for BD. However, the sizable presence of HLA-B51 negative BD patients indicates the role of at least one other genetic loci in BD pathogenesis [4, 5]. The lack of a clear mechanism of BD pathogenesis and the genetic variation present in BD patients makes it difficult to separate true genetic factors of the disease from background genetic noise.

The autoimmune nature of BD and the discovery of isolated shared genetic factors between other autoimmune diseases and BD is intriguing. For instance, ERAP1 has been shown to be a susceptibility gene in BD, psoriasis, and ankylosing spondylitis [6,7,8]. The molecular and genetic overlap between autoimmune diseases presents a significant area for research, as gaining a deeper understanding of the overlap can have important implications for diagnosis, treatments, and future research [9].

Addition research posits an intriguing approach towards uncovering this overlap. In addiction research, single genes rarely impact clinical phenotype, and hundreds of variants are needed to fully explain underlying genetics of the diseases. Convergent functional genomics (CFG) is leveraged to isolate important signals against this background. The premise of CFG is straight forward: the more lines of evidence for a gene, the higher it is prioritized as a gene of interest [9, 10]. This approach allows data from a variety of genomic sources to be synthesized as evidence to ascertain the impact of genes on phenotype [10].

To implement the CFG approach, we leverage the computational power of the Geneweaver Hierarchical Similarity (HiSIM) graph [11]. This tool creates a graph that is a hierarchical network of multiway geneset intersections. This enables users to find genes connected to all populated subsets of an input set of gene lists (genesets) [12]. Nodes at the top of the resulting graph contain a few strongly associated genes while nodes at the bottom contain weaker associated genes. In this work, we use the CFG approach to first find candidate genes of interest for BD and then evaluate other autoimmune diseases for genetic overlap with BD.

Main text


Data collection

Genes associated with BD were collected as genesets on Geneweaver [12]. Each of the sixteen genesets collected contained a record of BD associated genes from a single source. Eleven of the genesets originated from GWAS studies and were collected by searching the GWAS Catalog and publicly available curated genesets from Geneweaver [12, 13]. The combined GWAS data came from a global population that included samples from Iranian, Japanese, Turkish, Korean, Spanish, Western European, Middle Eastern, and Han Chinese populations. [14,15,16,17, 17,18,19,20,21,22,23]. One of the GWAS studies collected data regarding BD and a special type of BD that effects the GI tract, Intestinal BD (IBD) [17], and was split into two genesets for purposes of this study. The remaining BD genesets came from the NCBI gene2mesh tool [24], the Online Mendelian Inheritance in Man (OMIM) database tag “autoinflammatory, familial, Behçet’s-like” [25], Malacards [26], a BD review paper that consolidated genes of interest proposed in the literature, and a BD gene expression profile paper [27, 28]. In total, the data present in the 16 genesets, including the data consolidated by the review article geneset, came from 73 unique data sources and contained 319 genes (Additional file 1: Table S1). The union of all 16 genesets was subsequently collected using the Boolean Algebra tool on Geneweaver and stored as another geneset [12].

Publicly available Geneweaver genesets for 27 conditions were compiled in order to test BD ’s relation to other autoimmune diseases. For each of the 27 conditions, Homo Sapien genesets were found using the GeneWeaver search feature using keywords. The top 25 genesets per condition were included in the study, with priority given to manually curated genesets (Tier III) with a focused scope of study. Genesets associated with primary data (Tier I) were prioritized next, followed by derivative (Tier II) data sets. For example, a tier 1 geneset documenting a psoriasis GWAS study was given preference over a tier 2 study documenting a psoriasis and Crohn’s disease GWAS study. Disease genesets were created by using the Geneweaver Boolean Algebra tool on publicly available curated genesets (Additional file 1: Table S2) [12]. These 27 genesets were the product of 309 unique and individual genesets spanning over 4000 genes created using the Boolean Algebra tool on Geneweaver [11]

BD HiSIM run

To find BD consensus genes, the Geneweaver HiSIM graph was run with the parameters reported in Additional file 1: Table S3 and with the 16 collected BD genesets as input.

Autoimmune HiSIM run

To determine the genetic overlap between BD and other autoimmune diseases, the HiSIM graph tool was run with the parameters documented in Additional file 1: Table S3 on 27 autoimmune disease genesets constructed using the Boolean Algebra tool and the BD Union geneset. The data from the Autoimmune HiSIM run was then filtered to find relevant connections between BD and other autoimmune diseases.

Jaccard geneset analysis

A Jaccard geneset analysis was conducted on autoimmune diseases identified with a high BD genetic similarity to determine the statistical significance of the overlap. The Jaccard geneset analysis was run using the Geneweaver Jaccard Similarity tool [12].

Neighbor joining analysis

In order to find the top five most genetically similar diseases to BD, a neighbor joining tree was created using the ape R package [29]. The distance between any two diseases was defined using the normalized Jaccard formula seen in Formula 1.

$$\begin{aligned} 1 - \frac{A \cap B}{min(size(A), size(B))}. \end{aligned}$$

Formula 1 took into account the differences in geneset size between two genesets when determining their overlap. By taking the size of the intersection between the two sets and dividing it by the largest possible overlap that could occur between the two sets, the formula allowed sets of different sizes to be compared. The distance between all possible pair combinations of the identified diseases was calculated and stored in a distance matrix. The distance matrix was used as input to the ape neighbor joining program [29, 30].


BD HiSIM run

Running the HiSIM graph on the BD collected genesets resulted in the graph shown in Fig. 1. HLA-B was identified in 7 genesets making it the most common gene among all tested genes. It was followed by IL-10 in 6 genesets, IL23R in 5 genesets, and HLA-A, STAT4, MICA, and ERAP1 in 4 genesets (Additional file 1: Table S4).

Fig. 1

HiSIM graph between 17 BD gensets

Autoimmune disease run

After collecting autoimmune genesets on Geneweaver, the HiSIM graph was run on 27 unique autoimmune diseases and BD. Figure 2 displays the genetic overlap between the autoimmune diseases and BD as determined by the HiSIM graph. Out of the 27 conditions tested, 16 were found to have some genetic overlap with BD and were subsequently used to preform the Jaccard Analysis.

Fig. 2

Graph displaying identified HiSIM nodes with the largest gene overlap with BD. 16 conditions were found to have some overlap with BD out of the 27 tested

Jaccard geneset analysis

Using the 16 identified conditions from the Autoimmune HiSIM run, a Jaccard Geneset Analysis was run. The results from the Jaccard geneset analysis found no statistically significant overlap between any of the tested pairs of conditions (Additional file 1: Figure S1).

Neighbor joining analysis

The genetic overlap between BD and the 16 identified autoimmune diseases was then normalized using Formula 1 and used as input for the Neighbor Joining Tree (Fig. 3). The tree indicates that BD is closest to Sarcoidosis, Uveitis, Sjögren’s Syndrome, Hemolytic Anemia, and Myositis.

Fig. 3

Results of the neighbor joining tree


The results of the BD HiSIM graph run are consistent with previous BD research. All top scoring genes have previously been identified as associated with BD [6, 27, 31]. The emergence of HLA-B as the top scoring gene is unsurprising given that HLA-B51 is the most closely associated risk factor of BD [4, 31, 32].

The results of the BD autoimmune disease run identified 16 diseases with a BD genetic overlap. However, these associations were unable to achieve statistical significance based on the results of our Jaccard Similarity run (Additional file 2: Figure S1). To account for this, we proposed a modified Jaccard similarity formula (Formula 1) and constructed a neighbor joining tree that allowed us to quantify the genetic similarity of each of the 16 identified diseases to BD.

The results of our CFG approach have the benefit of identifying autoimmune diseases with a shared genetic architecture to BD. The inclusion of certain diseases in the 16 identified diseases was expected: Uveitis is symptom of BD [32], Sjögren’s Syndrome often presents with but is not associated with BD [33,34,35], and Crohn’s disease presents with similar symptoms to IBD [36].

The inclusion of psoriasis in our list was noteworthy. Psoriatic arthritis has been documented to sometimes be confused with BD articular involvement [37]. Recently, Hahn et al. found that psoriasis patients were twice as likely to be diagnosed with BD compared to controls [38].

Encouragingly, there have been documented case reports of patients presenting simultaneously with both BD and one of conditions identified in this study [39,40,41,42,43,44,45,46,47,48,49,50,51,52,53]. The only exception to this was Alopecia Areata which to our knowledge has not been reported as presenting simultaneously with BD. These findings, coupled with our genetic analysis, might provide evidence for an underlying genetic mechanism related to the pathogenesis of these diseases. It will be left up to future work to fully clarify these relationships.


Behçet’s disease is a complex multi-system inflammatory disease whose exact pathogenesis continues to elude researchers [2, 3]. Genetically, the identification of HLA-B51 as a major, but not sole, susceptibility gene has led to the hunt for other genetic factors of the disease [4]. The resulting identification of multiple genes of interest, however, does not explicitly establish the contributions of each of these genes towards the overall presentation of the disease [6, 7, 27, 32, 54,55,56,57]. Interestingly, research into some of these genes has uncovered roles of these genes in the progression of other autoimmune diseases [6,7,8, 58]. Therefore, researching autoimmune diseases in relation to other autoimmune diseases might be needed to fully understand human disease.

In this study, we employed a functional convergent genomics approach to discover (1) the genetic factors and (2) related autoimmune conditions of Behçet’s disease. The power of this approach lies in its ability to synthesize information from multiple genomics data sources [9] and in its recognition of shared genetic factors of autoimmune diseases [6, 7, 56]. The Geneweaver HiSIM graph’s ability to synthesize this information presents a valuable opportunity for further discovery in this line of research [11]. Furthermore, our results using this approach confirmed existing BD genetics research and identified 16 autoimmune diseases that share an underlying genetic relationship to BD. Almost all of these associations have documented clinical findings linking them with BD, further providing evidence towards our results. It will be left up to further research to fully uncover the complex genetic interactions underlying these diseases and the shared genetic mechanisms between them.


A limitation to our study was that our selection criteria excluded non-autoimmune conditions, such as oral cavity cancer, that might mimic BD phenotype presentations. Another limitation to this study was that the Jaccard Similarity test used suffers from dataset size limitations [59].

Availability of data and materials

The genesets used for this study were exported from Geneweaver and included as Additional files 3 and 4.



Behçet’s disease


Convergent functional genomics


  1. 1.

    Behçet H, Matteson EL. On relapsing, aphthous ulcers of the mouth, eye and genitalia caused by a virus. 1937. Clin Exp Rheumatol. 2010;28(4 Suppl 60):S2-5.

    PubMed  Google Scholar 

  2. 2.

    Davatchi F. Behcet’s disease. Int J Rheumat Dis. 2014;17(4):355–7.

    Article  Google Scholar 

  3. 3.

    Bulur I, Onder M. Behçet disease: new aspects. Clin Dermatol. 2017;35(5):421–34.

    Article  Google Scholar 

  4. 4.

    Verity DH, Marr JE, Ohno S, Wallace GR, Stanford MR. Behçet’s disease, the Silk Road and HLA-B51: historical and geographical perspectives: Behçet’s disease, the Silk Road and HLA-B51. Tissue Antigens. 1999;54(3):213–20.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Muruganandam M, Rolle NA, Sibbitt WL Jr, Cook GB, Emil NS, Fangtham M, et al. Characteristics of Behcet’s disease in the American Southwest. Seminars in Arthritis and Rheumatism. 2019;49(2):296–302. Available from:

  6. 6.

    Kirino Y, Nakajima H. Clinical and genetic aspects of Behçet’s disease in Japan. Int Med. 2019;58(9):1199–1207. Available from:

  7. 7.

    Conde-Jaldón M, Montes-Cano MA, García-Lozano JR, Ortiz-Fernández L, Ortego-Centeno N, González-León R, et al. Epistatic Interaction of ERAP1 and HLA-B in Behçet disease: a replication study in the Spanish population. PLoS ONE. 2014;9(7):e102100.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Tran TM, Colbert RA. Endoplasmic reticulum aminopeptidase 1 and rheumatic disease: functional variation. Curr Opin Rheumatol. 2015;27(4):357–363. Available from:

  9. 9.

    Spanagel R. Convergent functional genomics in addiction research—a translational approach to study candidate genes and gene networks. Silico Pharmacol. 2013;1:18.

    Article  Google Scholar 

  10. 10.

    Niculescu AB. Convergent functional genomics of psychiatric disorders. Am J Med Genet Part B Neuropsychiatr Genet. 2013;162(7):587–94.

    CAS  Article  Google Scholar 

  11. 11.

    Baker E, Bubier JA, Reynolds T, Langston MA, Chesler EJ. GeneWeaver: data driven alignment of cross-species genomics in biology and disease. Nucleic Acids Res. 2016;44(D1):D555–9.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Baker EJ, Jay JJ, Bubier JA, Langston MA, Chesler EJ. GeneWeaver: a web-based system for integrative functional genomics. Nucleic Acids Res. 2012;40(Database issue):D1067-1076.

    CAS  Article  Google Scholar 

  13. 13.

    MacArthur J, Bowler E, Cerezo M, Gil L, Hall P, Hastings E, et al. The new NHGRI-EBI Catalog of published genome-wide association studies (GWAS Catalog). Nucleic Acids Res. 2017;45(D1):D896–901.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Lee YJ, Horie Y, Wallace GR, Choi YS, Park JA, Choi JY, et al. Genome-wide association study identifies GIMAP as a novel susceptibility locus for Behçet’s disease. Ann Rheumat Dis. 2013;72(9):1510–6.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Fei Y, Webb R, Cobb BL, Direskeneli H, Saruhan-Direskeneli G, Sawalha AH. Identification of novel genetic susceptibility loci for Behçet’s disease using a genome-wide association study. Arthritis Res Ther. 2009;11(3):R66.

    Article  Google Scholar 

  16. 16.

    Ortiz-Fernández L, Carmona FD, Montes-Cano MA, García-Lozano JR, Conde-Jaldón M, Ortego-Centeno N, et al. Genetic analysis with the immunochip platform in Behçet disease identification of residues associated in the HLA class I region and new susceptibility Loci. PloS One. 2016;11(8):e0161305.

    Article  Google Scholar 

  17. 17.

    Kim SW, Jung YS, Ahn JB, Shin ES, Jang HW, Lee HJ, et al. Identification of genetic susceptibility loci for intestinal Behçet’s disease. Sci Rep. 2017;7:39850.

    CAS  Article  Google Scholar 

  18. 18.

    Mizuki N, Meguro A, Ota M, Ohno S, Shiota T, Kawagoe T, et al. Genome-wide association studies identify IL23R-IL12RB2 and IL10 as Behçet’s disease susceptibility loci. Nat Genet. 2010;42(8):703–6.

    CAS  Article  Google Scholar 

  19. 19.

    Kirino Y, Bertsias G, Ishigatsubo Y, Mizuki N, Tugal-Tutkun I, Seyahi E, et al. Genome-wide association analysis identifies new susceptibility loci for Behçet’s disease and epistasis between HLA-B*51 and ERAP1. Nat Genet. 2013;45(2):202–7.

    CAS  Article  Google Scholar 

  20. 20.

    Xavier JM, Shahram F, Sousa I, Davatchi F, Matos M, Abdollahi BS, et al. FUT2: filling the gap between genes and environment in Behçet’s disease? Ann Rheumat Dis. 2015;74(3):618–24.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Hou S, Yang Z, Du L, Jiang Z, Shu Q, Chen Y, et al. Identification of a susceptibility locus in STAT4 for Behçet’s disease in Han Chinese in a genome-wide association study. Arthritis Rheumat. 2012;64(12):4104–13.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Kappen JH, Medina-Gomez C, van Hagen PM, Stolk L, Estrada K, Rivadeneira F, et al. Genome-wide association study in an admixed case series reveals IL12A as a new candidate in Behçet disease. PloS One. 2015;10(3):e0119085.

    Article  Google Scholar 

  23. 23.

    Remmers EF, Cosan F, Kirino Y, Ombrello MJ, Abaci N, Satorius C, et al. Genome-wide association study identifies variants in the MHC class I, IL10, and IL23R-IL12RB2 regions associated with Behçet’s disease. Nat Genet. 2010;42(8):698–702.

    CAS  Article  Google Scholar 

  24. 24.

    Athey BD, Cavalcoli JD, Jagadish HV, Omenn GS, Mirel B, Kretzler M, et al. The NIH National Center for Integrative Biomedical Informatics (NCIBI). J Am Med Inform Assoc. 2012;19(2):166–70.

    Article  PubMed  Google Scholar 

  25. 25.

    Amberger JS, Bocchini CA, Scott AF, Hamosh A. leveraging knowledge across phenotype–gene relationships. Nucleic Acids Res. 2019;47(D1):D1038–D1043. Available from:

  26. 26.

    Rappaport N, Nativ N, Stelzer G, Twik M, Guan-Golan Y, Iny Stein T, et al. MalaCards: an integrated compendium for diseases and their annotation. Database. 2013;2013. Available from:

  27. 27.

    Deng Y, Zhu W, Zhou X. Immune Regulatory Genes Are Major Genetic Factors to Behcet Disease: Systematic Review. The Open Rheumatology Journal. 2018;12(1):70–85. Available from:

  28. 28.

    Puccetti A, Fiore PF, Pelosi A, Tinazzi E, Patuzzo G, Argentino G, et al. Gene expression profiling in Behcet’s disease indicates an autoimmune component in the pathogenesis of the disease and opens new avenues for targeted therapy. J Immunol Res. 2018;2018:1–18. Available from:

  29. 29.

    Paradis E, Schliep K. Ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics. 2019;35:526–8.

    CAS  Article  Google Scholar 

  30. 30.

    R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2013. Available from:

  31. 31.

    Mizuki N, Inoko H, Ohno S. Molecular genetics (HLA) of Behçet’s disease. Yonsei Med J. 1997;38(6):333–49.

    CAS  Article  Google Scholar 

  32. 32.

    Salmaninejad A, Gowhari A, Hosseini S, Aslani S, Yousefi M, Bahrami T, et al. Genetics and immunodysfunction underlying Behçet’s disease and immunomodulant treatment approaches. J Immunotoxicol. 2017;14(1):137–51.

    CAS  Article  Google Scholar 

  33. 33.

    Nair J, Singh T. Sjogren’s syndrome: review of the aetiology, pathophysiology a potential therapeutic interventions. J Clin Exp Dentistry. 2017;0–0. Available from:

  34. 34.

    Günaydin I, Ustündağ C, Kaner G, Pazarli H, Yurdakul S, Hamuryudan V, et al. The prevalence of Sjögren’s syndrome in Behçet’s syndrome. J Rheumatol. 1994;21(9):1662–4.

    PubMed  Google Scholar 

  35. 35.

    Ju FH, Xu TZ, Hong HH, Mao H, Wang M, Wang Z. Behcet disease combined with Sjogren syndrome: a unique case report and literature review. Medicine. 2018;97(12):e0138.

    Article  Google Scholar 

  36. 36.

    Nordstrom E, Fischer M. The great masquerader: Behcet’s disease. Case Reports. 2014;2014(apr 19 1):bcr2013202919–bcr2013202919.

    Article  Google Scholar 

  37. 37.

    Kokturk A. Clinical and pathological manifestations with differential diagnosis in Behçet’s disease. Pathol Res Int. 2012;2012:1–9. Available from:

  38. 38.

    Hahn HJ, Kwak SG, Kim DK, Kim JY. Association of Behçet disease with psoriasis and psoriatic arthritis. Sci Rep. 2021;11(1):2531. Available from:

  39. 39.

    Louro M, Vaio T, Crespo J, Santos R, Carvalho A. Pulmonary sarcoidosis in Behçet’s disease treated with adalimumab. Eur J Case Rep Intern Med. 2017;4(4):000576.

    PubMed  PubMed Central  Google Scholar 

  40. 40.

    Watanabe H, Yashiro M, Asano T, Sato S, Takahashi A, Katakura K, et al. A case of Behçet’s disease and systemic sclerosis developing after an earthquake disaster. Fukushima J Med Sci. 2015;61(1):86–90. Available from:

  41. 41.

    Ichimura Y, Asano Y, Takahashi T, Toyama T, Taniguchi T, Tamaki Z, et al. Development of systemic sclerosis in patients with Behçet’s disease: remission of Behçet’s disease in parallel with the progression of skin sclerosis. J Dermatol. 2014;41(12):1113–4.

    Article  PubMed  Google Scholar 

  42. 42.

    Ogose T, Tamaki W, Shinahara K, Kaneko M, Wakata Y, Naruse K. A case of recurrent myositis as the main manifestation of Behçet disease. Pediatr Int. 2010;52(2):e101–4.

    Article  PubMed  Google Scholar 

  43. 43.

    Dursun D, Akova Y, Yücel E. Myositis and scleritis associated with Behcet’s disease: an atypical presentation. Ocular Immunol Inflam. 2004;12(4):329–32.

    Article  Google Scholar 

  44. 44.

    Sarui H. Necrotising myositis in Behcet’s disease: characteristic features on magnetic resonance imaging and a review of the literature. Annal Rheumat Dis. 2002;61(8):751–2.

    CAS  Article  Google Scholar 

  45. 45.

    Carlos Arteaga R, Otto J Hernandez F, Renato Puppi M, Olga Judith Hernández F. Polymyositis in Adamantiades-Behçet’s Disease. Int J Neurol Neurother. 2020;7(1). Available from:

  46. 46.

    Mignemi G, Facchini C, Raimondo D, Montanari G, Ferrini G, Seracchioli R. A Case report of nasal endometriosis in a patient affected by Behcet’s disease. J Minimally Invasive Gynecol. 2012;19(4):514–516. Available from:

  47. 47.

    Joshi H, Shahriar I, Sharma P, Sagi SV, Oyibo SO. A rare coexistence of Behcet’s disease and Graves’ thyrotoxicosis in a young man: a case report. Oxford Med Case Rep. 2020;2020(1):132.

    Article  Google Scholar 

  48. 48.

    Ertaş R, Özyurt K, Avcı A, Ketenci Ertas S, Atasoy M. Case report: Behçet’s disease accompanied with vitiligo. F1000Research. 2017;6:310. Available from:

  49. 49.

    Guney E, Akcali G, Akcay BI, Unlu C, Erdogan G, Bozkurt TK, et al. Vitiligo in a patient treated with interferon alpha-2a for Behçet’s Disease. Case Rep Med. 2012;2012:1–3. Available from:

  50. 50.

    Chyuma Y, Utsunomiya A, Saito T, Hanada S, Nishimata H, Arima T. Hemolytic anemia complicated with Behçet’s disease and myelodysplastic syndrome] [Rinsho Ketsueki]. Japan J Clin Hematol. 1992;33(3):333–7.

    CAS  Google Scholar 

  51. 51.

    Chung SY, Oh KK, Chang HS. Sonographic findings of tuberculous thyroiditis in a patient with Behçet’s syndrome: tuberculous thyroiditis. J Clin Ultrasound. 2002;30(3):184–8.

    Article  PubMed  Google Scholar 

  52. 52.

    Carvalho FP, Maia B, Marinho A, Veloso A, Cunha I, Leite I, et al. Subacute thyroiditis in association with psoriasis and behcet’s disease. Endocrine Abstracts. 2018;Available from:

  53. 53.

    Iwadate H. A case of primary biliary cirrhosis complicated by Behçet’s disease and palmoplantar pustulosis. World J Gastroenterol. 2006;12(13):2136. Available from:

  54. 54.

    Piga M, Mathieu A. Genetic susceptibility to Behcet’s disease: role of genes belonging to the MHC region. Rheumatology. 2011;50(2):299–310.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Leccese P, Alpsoy E. Behçet’s disease: an overview of etiopathogenesis. Front Immunol. 2019;10:1067.

    CAS  Article  Google Scholar 

  56. 56.

    Takeuchi M, Ombrello MJ, Kirino Y, Erer B, Tugal-Tutkun I, Seyahi E, et al. A single endoplasmic reticulum aminopeptidase-1 protein allotype is a strong risk factor for Behçet’s disease in HLA-B*51 carriers. Annal Rheumat Dis. 2016;75(12):2208–11.

    CAS  Article  Google Scholar 

  57. 57.

    Sousa I, Shahram F, Francisco D, Davatchi F, Abdollahi BS, Ghaderibarmi F, et al. Brief report: association of CCR1, KLRC4, IL12A-AS1, STAT4, and ERAP1 With Behçet’s disease in Iranians. Arthritis Rheumatol. 2015;67(10):2742–8.

    Article  PubMed  Google Scholar 

  58. 58.

    Gibson AW, Edberg JC, Wu J, Kimberly RP. The role of IL-10 in autoimmune pathology. Landes Biosci. 2013. Publication Title: Madame Curie Bioscience Database [Internet]. Available from:

  59. 59.

    Salvatore S, Dagestad Rand K, Grytten I, Ferkingstad E, Domanska D, Holden L, et al. Beware the Jaccard: the choice of similarity measure is important and non-trivial in genomic colocalisation analysis. Briefings in Bioinformatics. 2020;21(5):1523–1530. Available from:

Download references


We thank our colleagues from Baylor University who provided insight and expertise that greatly assisted the research, although they may not agree with all of the interpretations/conclusions of this paper.


Not applicable.

Author information




All authors provided substantial contributions to the conception and design of the study. SS carried out the experiments and wrote the manuscript. EB supervised the project and revised the manuscript critically for important intellectual content. Both authors read and approved the final manuscript.

Corresponding author

Correspondence to Erich J. Baker.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Table S1 contains a list of the 27 autoimmune diseases tested,the number of genesets collected for each condition, the total number of genes collected for condition, and theassociated publications for the genesets. Table S2 contains the parameters used to run the HiSIM graph. Table S3contains the top results from the BD HiSIM run. Table S4 contains the sources used to create the review articlegeneset used for the BD HiSIM run.

Additional file 2: Figure S1

shows the results from the Jaccard Similarity run.

Additional file 3:

The genesets included in this text (.txt) file reference studies looking into the genetics of BD. These genesets wereused to conduct the BD HiSIM Run. This text file was created using the export genesets functionality inGeneweaver. All genesets included are publicly available and searchable on Geneweaver.

Additional file 4:

The genesets included in this text (.txt) file were created using the Boolean Algebra tool on publicly availablegenesets on Geneweaver. These genesets were used to conduct the BD Auotimmune Disease Run. This text file wascreated using the export genesets functionality in Geneweaver. All genesets included are publicly available andsearchable on Geneweaver.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Shenoi, S.J., Baker, E.J. Using hierarchical similarity to examine the genetics of Behçet’s disease. BMC Res Notes 14, 353 (2021).

Download citation


  • Behçet syndrome
  • Genomics
  • Autoimmune disease