Skip to main content
  • Research article
  • Open access
  • Published:

Changes in the proportion of clinical clusters contribute to the phenotypic evolution of Behçet’s disease in Japan



We hypothesized that Behçet’s disease (BD) consists of several clinical subtypes with different severity, resulting in heterogeneity of the disease. Here, we conducted a study to identify clinical clusters of BD.


A total of 657 patients registered in the Yokohama City University (YCU) regional BD registry between 1990 and 2018, as well as 6754 patients who were initially registered in the Japanese Ministry of Health, Labour and Welfare (MHLW) database between 2003 and 2014, were investigated. The YCU registry data regarding the clinical manifestations of BD, human leukocyte antigen (HLA) status, treatments, and hospitalizations were analyzed first, followed by similar analyses of the MHLW for validation. A hierarchical cluster analysis was independently performed in both patient groups.


A hierarchical cluster analysis determined five independent clinical clusters in the YCU cohort. Individual counterparts of the YCU clusters were confirmed in the MHLW registry. Recent phenotypical evolutions of BD in Japan, such as increased gastrointestinal (GI) involvement, reduced complete type according to the Japan Criteria, and reduced HLA-B51 positivity were associated with chronologically changing proportions of the clinical clusters.


In this study, we identified independent clinical clusters among BD patients in Japan and found that the proportion of each cluster varied over time. We propose five independent clusters namely “mucocutaneous”, “mucocutaneous with arthritis”, “neuro”, “GI”, and “eye.”


Behçet’s disease (BD) is a rare systemic inflammatory disease that causes periodic ocular and mucocutaneous inflammation [1, 2]. Environmental and genetic factors such as human leukocyte antigen (HLA)-B51 have been implicated in its pathogenesis [3,4,5]. The 2018 update EULAR recommendation emphasized individualized treatment for BD [6]: intensive treatment is particularly necessary for ocular, vascular, neurological, and gastrointestinal (GI) symptoms, which are associated with poor outcomes [6, 7]. Therefore, identifying predictive factors for these unfavorable BD subtypes is essential.

One of the primary obstacles facing personalized medicine is the heterogeneity of BD patients [2]. Serious organ involvement such as ocular and vascular lesions in young men [8, 9], genital ulcers in young women [9], and GI symptoms in pediatric patients [9, 10] have been reported. Geographic and genetic disease variations are other features of BD, including distinct clinical manifestations in countries with prevalent disease [1, 2, 11]. GI involvement is more frequent in Far East-Asian countries, including Japan [9, 12,13,14], whereas vascular involvement is more common in Mediterranean and Middle Eastern countries [15, 16].

The clinical features of BD have evolved in Japan despite the preservation of a homogenous genetic background because of low immigration rates. Increased GI manifestations, decreased eye involvement, and fewer HLA-B51-positive patients over the past 3 decades have been reported [13, 17, 18]. The same trend was observed in Korea, where GI involvement is also prevalent [14]. Interestingly, an explosive increase was also observed in inflammatory bowel diseases in East Asian countries [19, 20]. It is plausible that the phenotypic changes are caused by preferential expansion of a particular subset of BD in Japan; this idea is supported by previous meta-analyses that revealed a strong association between HLA-B51 and both eye involvement and male sex, whereas its association with GI involvement was weak [21, 22].

Because of the diverse manifestations of this disease, we aimed to stratify patients with varying clinical manifestations and different prognoses. For example, the genetic structures of anti-citrullinated protein antibody-positive and -negative patients with rheumatoid arthritis are different; an HLA-shared epitope is more common in patients of the former group, who also experience progressive bone destruction [23,24,25]. Previous studies have shown that patients with BD are categorized into several distinct phenotypic clusters associated with varying prognoses [26,27,28,29]. Clustering analyses have delineated several possible BD subsets, including vascular disease, while another subset is an acne-arthritis-enthesitis cluster in Turkey [26]. Interestingly, the proposed acne-arthritis-enthesitis cluster showed familial aggregation, suggesting a possible contribution of subset-specific genetic factors [30].

Prognostic factors for BD are necessary, but the disease heterogeneity is precluding the development of precision medicine. In this study, clustering analysis was performed on patients with BD using data from nationwide and regional registries in Japan to stratify complex phenotype of BD.



Data from patients with BD treated in Japan were analyzed retrospectively. All patients met the 1987 revised diagnostic criteria for BD as designated by the Behçet’s Disease Research Committee (BDRC) of Japan [31]. According to Japanese BD diagnostic criteria, oral ulcers, skin lesions, eye inflammation, and genital ulcers are “major symptoms;” arthritis, epididymitis, gastrointestinal symptoms, vascular lesions, and neurological manifestations are “minor symptoms.” Patients in one of the following categories meet criteria for BD: (1) all four major symptoms, (2) three major symptoms, (3) two major and two minor symptoms, (4) typical recurrent ocular inflammation and one or more major symptoms, or (5) typical recurrent ocular inflammation and two minor symptoms meet criteria for BD diagnosis. “GI symptoms” were noted only for patients who were found to have deep ulcers in the ileocecum area, whereas patients with “neurological manifestations” included those with parenchymal neuro-BD.

To minimize patient heterogeneity and increase comparability with other studies, we also applied the International Study Group (ISG) and International Team for the Revision of the International Criteria for Behçet’s Disease (ITR-ICBD) criteria [32, 33].

Hospital-based regional BD registry

The Yokohama City University (YCU) registry includes BD patients who were treated at any of the seven hospitals within the Kanagawa district in mid-Japan from July 1991 to December 2018. All the patients were Japanese except eight patients (one Filipino, one Syrian, one Korean, two Chinese, three dual heritage of Pakistani, Finnish, and French). Retrospectively, we obtained the following data: date of birth, sex, date of BD onset, follow-up period, individual symptoms of BD, HLA-B51 positivity, episode of hospitalization, blindness, death caused by BD, and medications. Blindness is defined as visual acuity of less than 3/60, or a corresponding visual field loss to less than 10°, even with just one eye.

Japan national BD registry

Since 1972, the Japan MHLW has been granted public subsidies for patients’ medical expenses for treatment of intractable diseases, including BD. To be eligible for this assistance, patients must meet the BRDC criteria and provide documentation of their date of birth, gender, date of onset, the presence of BD symptoms included in the Japanese BD diagnostic criteria (as filed by the attending physician), questions about laboratory data, admission episodes, and medical treatment. We extracted data recorded in the MHLW national registry between 2003 and 2014; to adjust for the observation period, only newly submitted cases were included. We selected newly registered patients in order to adjust for the observation period, given that manifestations can accumulate over one’s clinical course and change the BD-associated phenotype [17]. Moreover, the effects of treatment on the phenotypes were marginal in this group because new registrations were linked to financial support regarding medical costs. Some of the patients were diagnosed before 2003 although they were initially registered after that year.

Statistical analysis

R version 3.6.2 ( was used for our analyses. We performed Ward’s hierarchical cluster analysis with dissimilarity measures of the square Euclidean distance for both registry data. The validity of the number of clusters was assessed by several methods based on inter- and intra-cluster “scatter” indicators [34,35,36]. In these methods, the indices CH, indicator diffH, and indicator KL/minimization indicator H, which maximizes the number of clusters, were selected as arguments. The Kruskal-Wallis and chi-squared were performed to compare groups. A P value < 0.05 was considered statistically significant.



Among 707 patients with BD in the YCU registry and 8282 in the MHLW registry, 657 and 6754 patients were included after excluding those having missing data (Supplementary Figure 1). Table 1 shows the clinical characteristics of the two patient groups. These characteristics were similar to those previously reported in Japan, albeit with a lower incidence of eye involvement in the MHLW registry [1].

Table 1 Characteristic of patients with Behçet’s disease who were enrolled in the study

Discovery of distinct clusters in the hospital-based BD registry

We first performed cluster analysis using 8 variables (4 major and 4 minor symptoms as described in the Japanese criteria except for epididymitis). A dendrogram revealed 5 major clusters assigned as clusters 1 to 5 (Supplementary Figure 2a). The validity of the number of clusters was supported by statistical analysis (Supplementary Table 1). We then compared the clinical features of each cluster (Table 2). In brief, all patients in cluster 1 showed mucocutaneous symptoms without vital organ involvement except in the eyes. Cluster 2 was characterized by the highest proportion of GI involvement (75.7%), and the lowest of eye involvement (34.3%). In cluster 3, mucocutaneous symptoms without vital organ involvement were main features like cluster 1, but arthritis was absent. In cluster 4, all patients had eye involvement with low frequency of other mucocutaneous manifestations. Cluster 5 was associated with high rates of neurological involvement (98.5%).

Table 2 Characteristics of the cluster with Behçet’s Disease patients in Yokohama City University registry

We compared other features including demography, HLA-B51 positivity, fulfillment of the ISG and the ITR-ICBD criteria, therapy, and hospitalization events among the 5 clusters. More female patients were found in clusters 1 and 3, whereas male was predominant in cluster 4. Frequency of HLA-B51-positive individuals was comparable among clusters, except cluster 2 (GI cluster) that showed significantly lower positivity (33.0%) than any other clusters. Less than quarter of patients did not meet the ISG criteria in clusters 2, 4, and 5, whereas all patients met the ITR-ICBD criteria except 12 patients (8.6%) in cluster 2. However, background characteristics of individual clusters were unchanged in patients meeting the ISG or ITR-ICBD criteria (Supplementary Tables 2 and 3). Besides clinical manifestations, the data showed therapeutic approaches were different among the clusters. The use of glucocorticoids, particularly high dosage, and immunosuppressants were more common in cluster 2 and cluster 5 than the others, due to the gastrointestinal and neurological involvement, respectively. Introduction of TNF inhibitors was the most common and the earliest in cluster 4 because of the eye involvement (periods from the diagnosis to administration; 2.96 ± 4.46 years). The cluster was associated with low incidence of mucocutaneous symptoms, and high unsatisfactory rate of ISG criteria (22.3%). On the other hand, patients in clusters 2 and 5 needed frequent hospitalization compared with the others, suggesting that they had more severe physical condition. Blindness is more in clusters 4 and 5, but most cases with blindness were diagnosed before 1999 (Supplementary Table 4). We thought that it is due to the advance of therapy, such as TNF inhibitors. The number of death caused by BD was very small (3 cases).

Validation in the Japanese national BD registry

To validate the findings, the same clustering analysis was conducted using the MHLW dataset (Supplementary Figure 2b). Although approximately 200 patients in the YCU registry overlapped with the MHLW database, these patients were not excluded because of complete anonymization. A hierarchical cluster analysis of eight BD clinical manifestations identified seven different clusters that were statistically valid (Supplementary Table 1) and designated as clusters A to G (Table 3). The clustering patterns between the YCU Registry and the MHLW database corresponded to clusters 1 and E, clusters 2 and A, clusters 3 and G, clusters 4, D and F, and clusters 5 and B, respectively, but there was some discrepancy. The exception was cluster C, which had the lowest frequency of oral ulcers, suggesting that the MHLW database may contain suspicious BD cases. Otherwise, the features of each cluster in the MHLW database were also shared with counterparts of the clusters in the YCU registry. For example, cluster A in the MHLW database was characterized by the frequent GI and uncommon eye involvement, low HLA-B51 positivity, low fulfillment of the ISG and ITR-ICBD criteria, and high frequency of hospitalization. These characteristics were also found in the counterpart cluster, cluster 2 in the YCU registry.

Table 3 Characteristics of the cluster with Behçet’s disease patients in MHLW registry

Thus, our data showed that clusters in the hospital-based BD registry were reproduced in the nationwide registry in which patients with a short disease duration were enriched, indicating that the clustering patterns represent those in BD patients in Japan.

Chronological changes of clinical clusters

We previously reported that Japanese BD patient phenotypes have been evolving over the past 30 years, with ocular involvement becoming less common while the incidence of intestinal BD continues to rise [9, 13]. Similar evolution of clinical phenotypes has also been reported in Korea [14]. In Japan, the phenotypical changes were also associated with decreased HLA-B51 positivity in the patients [9, 13]. Therefore, we hypothesized that the recent phenotypic evolution of BD is caused by alternating proportion of individual disease clusters in Japan.

As in a previous study [13], we compared the proportion of clusters in the three generation groups according to the year of onset in the YCU registry. We found a continuous increase in clusters 2 and 4, but a gradual decrease in the other clusters (Fig. 1). Particularly, expansion of cluster 2 characterized by GI involvement and low HLA-B51 positivity is concordant with recent epidemiological changes in Japan [13]. Another increment was found in cluster 4, which was characterized by the most common and earliest use of TNF inhibitors and a lower incidence of mucocutaneous symptoms (Table 2), particularly after 2010 (Fig. 1). Proportions of the other clusters had been relatively shrinking.

Fig. 1
figure 1

Evolution of clusters in the Yokohama City University registry. The Y-axis represents the percentage of each cluster to the total number of specified periods


Because BD involves variable disease manifestations, many investigators refer to it as “Behçet’s syndrome” [6, 28]. Our data strongly support the notion that BD patients can be categorized into several clinical subtypes based on clinical presentations. A hierarchical clustering analysis of the hospital-based, YCU registry identified five clinical clusters, which consist of “mucocutaneous,” “mucocutaneous with arthritis,” “neuro,” “GI,” and “eye without other involvements” subtypes. The similar clustering pattern was also reproduced in the national registry, the MHLW database. Moreover, the present data showed that recent phenotypical evolution of BD was associated with altering the proportion of individual clinical clusters in the YCU registry. The chronological proportional changes of clinical subtypes were considered to cause recent epidemiological changes of BD patients in Japan, which are characterized by an increase in GI involvement and decreases complete type according to the Japan Criteria and HLA-B51 positivity.

Both genetic and environmental factors have been implicated in the formation of clinical clusters. However, due to the relative stability of the genetic background in Japan, recent changes in the clinical picture of BD patients over time may be dependent solely on environmental factors. Chronological analysis of the YCU registry showed two clinical clusters which had been recently expanded: one is cluster 2, which is a cluster with the highest proportion of GI patients and has expanded over the past 30 years in the YCU registry. Counterpart of the GI cluster was considered to correspond to cluster A in the MHLW database. The GI involvement is much more frequent in patients in East Asian countries including Japan and Korea than any regions, but it is unclear what is responsible for the geographical concentration. In addition to similar genetic backgrounds, the preferential use of the Japanese BD diagnostic criteria in these 2 countries [31], which incorporate the intestinal symptoms, is a possible reason. Interestingly, a westernized lifestyle is considered to contribute to recent rapid increase in inflammatory bowel diseases such as Crohn disease and ulcerative colitis in both countries [20]. Likewise, environmental changes might be involved in recent expansion of the GI cluster in BD patients in the East Asia.

In the YCU, another expansion was found in cluster 4, which was characterized by the highest incidence of, eye involvement, low incidence of mucocutaneous symptoms, and most common and earliest use of TNF inhibitors. The total number of patients with ocular lesions has decreased (Supplementary Table 4), but cluster 4 has increased. We consider that it is because TNF inhibitors for ocular involvement have a decreasing effect on mucocutaneous symptoms. This cluster was considered to correspond to clusters D and F in the MHLW database. Unlike the GI cluster, the cluster expansion was prominent after 2010. This is not likely coincident, because the use of TNF inhibitors was rapidly increased after official approval of infliximab and adalimumab for the BD related uveitis at 2007 and 2016, respectively, in Japan. Notably, the YCU has a Behçet’s Disease Medical Research Center with experienced rheumatologists and ophthalmologists, prompting this trend because ophthalmology department of YCU had many severe BD uveitis patients and they were administered TNF inhibitors. As disease modifying effects have been shown with early introduction of azathioprine [37], it is plausible that early and intensive treatments including biologics for the ocular involvement suppress extraocular manifestations, resulting in low incidence of mucocutaneous symptoms and low satisfactory rate of the ISG criteria in the cluster. However, besides the treatment related disease modifying effects, other factors also contribute to the formation of this cluster, because substantial number of patients was categorized into the cluster in the pre-biologics era. On the other hand, clusters 1, 3, and 5, identified in the YCU registry, are steadily shrinking in Japan (Fig. 1). These clusters include patients with the ocular and neurological involvement. Besides, a recent study also found that ocular disease is becoming milder in Japan [38], though the reason remains uncertain.

A study from Japan has shown that neurological involvement is more likely to develop in male smokers [39], suggesting that the recent shrinkage in cluster 5 is relate to fewer male smokers in Japan, though positive or negative effects of smoking on BD have been controversial [39,40,41,42]. Alternatively, improvement of oral hygiene [43] and altered microbiota [44] can also contribute to lowering the rate of serious BD phenotypes. It is important to determine environmental factors that are alterable unlike genetic backgrounds, especially when considering prophylactic strategies.

Most of the findings were consistent between the YCU and MHLW databases, but there were some discrepancies. For example, cluster 2 of the YCU had a high frequency of complications of GI and vascular involvements [2], but cluster A of the MHLW database did not show a significant association between the two. We suspect that different periods of observation are one of the reasons for this discrepancy (YCU, 13.7 ± 11.8 years; MHLW, 4.1 ± 7.6 years). In addition, while data were collected cross-sectionally in the MHLW, data were accumulated over time in the YCU registry. As shown in previous studies, symptoms of BD accumulate over time, and therefore, additional symptoms are likely to appear after the end of the observation period in the MHLW database [13, 17]. In fact, the incidence of GI and vascular lesions differed, especially vascular lesions, which were only 2.2% in the MHLW database, compared to 8.4% in the YCU registry (Table 1). We suspect that there is some relationship between GI and vascular lesions, but we are unable to establish that point in this study because of the different periods of observation in two cohorts. Therefore, it is necessary to validate in a population of patients with similar disease duration.

There is accumulating evidence that clinical presentations of BD are different among ethnic groups, suggesting that each ethnic has a different clustering pattern [28]. The acne-arthritis-enthesitis cluster in Turkey, which likely belongs on the spectrum of seronegative spondyloarthritis [28], was not apparent in our population, though the clinical manifestations were not sufficiently assessed for comparison. This is partly due to the lower prevalence of enthesitis, spondyloarthritis, and HLA-B27 positivity in Japan than in Europe [45]. Both mucocutaneous manifestations and arthritis were enriched in cluster 1 of the YCU registry and cluster E of the MHLW database; these clusters are possible counterparts to the acne-arthritis-enthesitis cluster in Turkey. Further studies with precise phenotypic data are necessary to assess whether the acne-arthritis-enthesitis cluster exists in Japan; moreover, international comparative studies on BD cluster patterns might shed light on ethnic and geographic variations in the disease.

To establish precision medicine, it is important to focus more on genetic contributions rather than environmental factors. We previously reported multiple BD susceptibility loci through genome-wide association studies (GWAS) such as HLA-A26, IL23R, and others [46,47,48,49]. As shown in Tables 2 and 3, HLA-B51 positivity was significantly less prevalent in cluster 2 (YCU) and cluster A (MHLW) than in other clusters, suggesting that the GI cluster is associated with a distinct genetic background [50]. Identification of the diverse genetic characteristics of the intestinal variant will facilitate its management as a separate disease entity [51]. Moreover, it was possible that autoinflammatory disease (AID) associated with single gene mutation had been diagnosed as BD. For example, “haploinsufficiency of A20” which is caused by TNFAIP3 mutations has a phenotype similar to that of patients in the GI cluster of BD [52].

Several studies have attempted to subcategorize BD based on clinical manifestations [26, 27, 53]; however, a patient’s natural disease course [13, 17] as well as treatment-related disease modification may prevent the emergence of other phenotypes [54]. Thus, because phenotype-based analysis alone is not necessarily sufficient, the inclusion of genetic analysis would be more instructive.

Several limitations should be noted. First, our results are not necessarily be generalizable to other countries because most of the patients were Japanese who had relatively unique and homogenous genetic backgrounds. Second, selection bias was unavoidable because of the study design; however, we minimized this bias by evaluating two distinct patient populations with unique profiles. Finally, there are only a few reports of standardized index for severity of BD [55, 56]. Therefore, we would like to investigate this point by establishing a large-scale BD patient’s registry in Japan.

Based on our current findings and our clinical experience, we propose four distinct clinical clusters namely “mucocutaneous,” which consists of either with or without arthritis, “neuro,” “GI,” and “eye”(Fig. 2). To establish this concept, it is important to further accumulate future data by analyzing genetic and biological samples. A multidisciplinary approach is useful for the early detection of clinical clusters. Ideally, treatment should be optimized for individual clinical subtypes, as the prognosis varies between clusters.

Fig. 2
figure 2

Proposing Behçet’s disease clusters based on the present and previous studies


We identified independent clinical clusters among BD patients in Japan and found that the proportion of each cluster varied over time.

Availability of data and materials

The datasets used and analyzed during the current study are available from the corresponding author on a reasonable request.



Autoinflammatory disease


Behçet’s disease


Behçet’s Disease Research Committee




Genome-wide association studies


Haploinsufficiency of A20


Human leukocyte antigen


International Study Group


International Team for the Revision of the International Criteria for Behçet’s Disease


Ministry of Health, Labour and Welfare


Tumor necrosis factor


Yokohama City University


  1. Sakane T, Takeno M, Suzuki N, Inaba G. Behcet’s disease. N Engl J Med. 1999;341(17):1284–91.

    Article  CAS  PubMed  Google Scholar 

  2. Yazici H, Seyahi E, Hatemi G, Yazici Y. Behcet syndrome: a contemporary view. Nat Rev Rheumatol. 2018;14(2):107–19.

    Article  CAS  PubMed  Google Scholar 

  3. de Menthon M, Lavalley MP, Maldini C, Guillevin L, Mahr A. HLA-B51/B5 and the risk of Behcet’s disease: a systematic review and meta-analysis of case-control genetic association studies. Arthritis Rheum. 2009;61(10):1287–96.

    Article  PubMed  CAS  Google Scholar 

  4. Ono S, Aoki K, Sugiura S, Nakayama E, Itakura K. Letter: HL-A5 and Behcet’s disease. Lancet. 1973;2(7842):1383–4.

    Article  Google Scholar 

  5. Ombrello MJ, Kirino Y, de Bakker PI, Gul A, Kastner DL, Remmers EF. Behcet disease-associated MHC class I residues implicate antigen binding and regulation of cell-mediated cytotoxicity. Proc Natl Acad Sci U S A. 2014;111(24):8867–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Hatemi G, Christensen R, Bang D, Bodaghi B, Celik AF, Fortune F, Gaudric J, Gul A, Kotter I, Leccese P, et al. 2018 update of the EULAR recommendations for the management of Behcet’s syndrome. Ann Rheum Dis. 2018;77(6):808–18.

    PubMed  Google Scholar 

  7. Desbois AC, Vallet H, Domont F, Comarmond C, Cacoub P, Saadoun D. Management of severe complications in Behcet’s disease with TNF inhibitors. Expert Opin Biol Ther. 2017;17(7):853–9.

    Article  CAS  PubMed  Google Scholar 

  8. Yazici H, Tuzun Y, Pazarli H, Yurdakul S, Ozyazgan Y, Ozdogan H, Serdaroglu S, Ersanli M, Ulku BY, Muftuoglu AU. Influence of age of onset and patient’s sex on the prevalence and severity of manifestations of Behcet’s syndrome. Ann Rheum Dis. 1984;43(6):783–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Ishido T, Horita N, Takeuchi M, Kawagoe T, Shibuya E, Yamane T, Hayashi T, Meguro A, Ishido M, Minegishi K, et al. Clinical manifestations of Behcet’s disease depending on sex and age: results from Japanese nationwide registration. Rheumatology (Oxford). 2017;56(11):1918–27.

    Article  Google Scholar 

  10. Kone-Paut I, Shahram F, Darce-Bello M, Cantarini L, Cimaz R, Gattorno M, Anton J, Hofer M, Chkirate B, Bouayed K, et al. Consensus classification criteria for paediatric Behcet’s disease from a prospective observational cohort: PEDBD. Ann Rheum Dis. 2016;75(6):958–64.

    Article  PubMed  Google Scholar 

  11. Ishigatsubo Y, Takeno M. Overview. In: Ishigatsubo Y, editor. Behçet’s disease: from genetics to therapies. Tokyo: Springer Japan; 2015. p. 1–20.

    Google Scholar 

  12. Ideguchi H, Suda A, Takeno M, Miyagi R, Ueda A, Ohno S, Ishigatsubo Y. Gastrointestinal manifestations of Behcet’s disease in Japan: a study of 43 patients. Rheumatol Int. 2014;34(6):851–6.

    Article  PubMed  Google Scholar 

  13. Kirino Y, Ideguchi H, Takeno M, Suda A, Higashitani K, Kunishita Y, Takase-Minegishi K, Tamura M, Watanabe T, Asami Y, et al. Continuous evolution of clinical phenotype in 578 Japanese patients with Behcet’s disease: a retrospective observational study. Arthritis Res Ther. 2016;18(1):217.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Kim DY, Choi MJ, Cho S, Kim DW, Bang D. Changing clinical expression of Behcet disease in Korea during three decades (1983–2012): chronological analysis of 3674 hospital-based patients. Br J Dermatol. 2014;170(2):458–61.

    Article  CAS  PubMed  Google Scholar 

  15. Sarica-Kucukoglu R, Akdag-Kose A, Kayabal IM, Yazganoglu KD, Disci R, Erzengin D, Azizlerli G. Vascular involvement in Behcet’s disease: a retrospective analysis of 2319 cases. Int J Dermatol. 2006;45(8):919–21.

    Article  PubMed  Google Scholar 

  16. Davatchi F, Shahram F, Chams-Davatchi C, Shams H, Abdolahi BS, Nadji A, Faezi T, Akhlaghi M, Ghodsi Z, Karimi N, et al. Behcet’s disease in Iran: analysis of 7641 cases. Mod Rheumatol. 2019;29(6):1023–30.

    Article  PubMed  Google Scholar 

  17. Ideguchi H, Suda A, Takeno M, Ueda A, Ohno S, Ishigatsubo Y. Behcet disease: evolution of clinical manifestations. Medicine (Baltimore). 2011;90(2):125–3.

  18. Mizuki Y, Horita N, Horie Y, Takeuchi M, Ishido T, Mizuki R, Kawagoe T, Shibuya E, Yuda K, Ishido M, et al. The influence of HLA-B51 on clinical manifestations among Japanese patients with Behcet’s disease: a nationwide survey. Mod Rheumatol. 2020;30(4):708–14.

    Article  CAS  PubMed  Google Scholar 

  19. Kuwahara E, Murakami Y, Nakamura T, Inoue N, Nagahori M, Matsui T, Watanabe M, Suzuki Y, Nishiwaki Y. Factors associated with exacerbation of newly diagnosed mild ulcerative colitis based on a nationwide registry in Japan. J Gastroenterol. 2017;52(2):185–93.

    Article  PubMed  Google Scholar 

  20. Ng WK, Wong SH, Ng SC. Changing epidemiological trends of inflammatory bowel disease in Asia. Intest Res. 2016;14(2):111–9.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Maldini C, Lavalley MP, Cheminant M, de Menthon M, Mahr A. Relationships of HLA-B51 or B5 genotype with Behcet’s disease clinical characteristics: systematic review and meta-analyses of observational studies. Rheumatology (Oxford). 2012;51(5):887–900.

    Article  Google Scholar 

  22. Suwa A, Horita N, Ishido T, Takeuchi M, Kawagoe T, Shibuya E, Yamane T, Hayashi T, Meguro A, Ishido M, et al. The ocular involvement did not accompany with the genital ulcer or the gastrointestinal symptoms at the early stage of Behcet’s disease. Mod Rheumatol. 2019;29(2):357–62.

    Article  PubMed  Google Scholar 

  23. Kroot EJ, de Jong BA, van Leeuwen MA, Swinkels H, van den Hoogen FH, van't Hof M, van de Putte LB, van Rijswijk MH, van Venrooij WJ, van Riel PL. The prognostic value of anti-cyclic citrullinated peptide antibody in patients with recent-onset rheumatoid arthritis. Arthritis Rheum. 2000;43(8):1831–5.

    Article  CAS  PubMed  Google Scholar 

  24. Ding B, Padyukov L, Lundstrom E, Seielstad M, Plenge RM, Oksenberg JR, Gregersen PK, Alfredsson L, Klareskog L. Different patterns of associations with anti-citrullinated protein antibody-positive and anti-citrullinated protein antibody-negative rheumatoid arthritis in the extended major histocompatibility complex region. Arthritis Rheum. 2009;60(1):30–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Holoshitz J. The rheumatoid arthritis HLA-DRB1 shared epitope. Curr Opin Rheumatol. 2010;22(3):293–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Hatemi G, Fresko I, Tascilar K, Yazici H. Increased enthesopathy among Behcet’s syndrome patients with acne and arthritis: an ultrasonography study. Arthritis Rheum. 2008;58(5):1539–45.

    Article  PubMed  Google Scholar 

  27. Jun Hayashi CT, Sugano Y, Nonaka N, Enomoto M, Ishii O, Nakamura T, Matsuda T, Sugimmori H. Study on clinical manifestations of Behcet’s disease by cluster analysis. Clin Rheumatol Relat Res. 2005;17(2):99–105.

    Google Scholar 

  28. Yazici H, Ugurlu S, Seyahi E. Behcet syndrome: is it one condition? Clin Rev Allergy Immunol. 2012;43(3):275–80.

    Article  CAS  PubMed  Google Scholar 

  29. Seyahi E. Phenotypes in Behcet’s syndrome. Intern Emerg Med. 2019;14(5):677–89.

    Article  PubMed  Google Scholar 

  30. Karaca M, Hatemi G, Sut N, Yazici H. The papulopustular lesion/arthritis cluster of Behcet’s syndrome also clusters in families. Rheumatology (Oxford). 2012;51(6):1053–60.

    Article  Google Scholar 

  31. Mizushima Y, Inaba G, Mimura Y. Guide for the diagnosis of Behçet’s disease. Report of Behçet’s Disease Research Committee, Japan; 1987. p. 8–17.

    Google Scholar 

  32. Criteria for diagnosis of Behcet’s disease. International study group for Behcet’s disease. Lancet. 1990;335(8697):1078–80.

    Google Scholar 

  33. The International Criteria for Behcet's Disease. (ICBD): a collaborative study of 27 countries on the sensitivity and specificity of the new criteria. J Eur Acad Dermatol Venereol. 2014;28(3):338–47.

    Article  Google Scholar 

  34. Caliński T, J Harabasz J. A dendrite method for cluster analysis. Communications in Statistics. 1974;3(1):1–27.

  35. Hartigan JA. Clustering analysis. Wiley; 1975.

  36. Krzanowski WJ, Lai YT. A criterion for determining the number of clusters in a data set. Biometrics. 1985;44(1):23–34.

    Article  Google Scholar 

  37. Hamuryudan V, Ozyazgan Y, Hizli N, Mat C, Yurdakul S, Tuzun Y, Senocak M, Yazici H. Azathioprine in Behcet’s syndrome: effects on long-term prognosis. Arthritis Rheum. 1997;40(4):769–74.

    Article  CAS  PubMed  Google Scholar 

  38. Nakahara H, Kaburaki T, Tanaka R, Yoshida A, Takamoto M, Kawata M, Fujino Y, Kawashima H, Aihara M. Comparisons of clinical features in Japanese patients with Behce’s uveitis treated in the 1990s and the 2000s. Ocul Immunol Inflamm. 2020;28(2):262–9.

    Article  CAS  PubMed  Google Scholar 

  39. Aramaki K, Kikuchi H, Hirohata S. HLA-B51 and cigarette smoking as risk factors for chronic progressive neurological manifestations in Behcet’s disease. Mod Rheumatol. 2007;17(1):81–2.

    Article  PubMed  Google Scholar 

  40. Malek Mahdavi A, Khabbazi A, Yaaghoobian B, Ghojazadeh M, Agamohammadi R, Kheyrollahiyan A, Rashtchizadeh N. Cigarette smoking and risk of Behcet’s disease: a propensity score matching analysis. Mod Rheumatol. 2019;29(4):633–9.

    Article  PubMed  Google Scholar 

  41. Lee YB, Lee JH, Lee SY, Yu DS, Han KD, Park YG. Association between smoking and Behcet’s disease: a nationwide population-based study in Korea. J Eur Acad Dermatol Venereol. 2019;33(11):2114–22.

    Article  CAS  PubMed  Google Scholar 

  42. Rizvi SW, McGrath H Jr. The therapeutic effect of cigarette smoking on oral/genital aphthosis and other manifestations of Behcet’s disease. Clin Exp Rheumatol. 2001;19(5 Suppl 24):S77–8.

    CAS  PubMed  Google Scholar 

  43. Direskeneli H, Mumcu G. A possible decline in the incidence and severity of Behcet’s disease: implications for an infectious etiology and oral health. Clin Exp Rheumatol. 2010;28(4 Suppl 60):S86–90.

    PubMed  Google Scholar 

  44. Consolandi C, Turroni S, Emmi G, Severgnini M, Fiori J, Peano C, Biagi E, Grassi A, Rampelli S, Silvestri E, et al. Behcet’s syndrome patients exhibit specific microbiome signature. Autoimmun Rev. 2015;14(4):269–76.

    Article  PubMed  Google Scholar 

  45. Kishimoto M, Yoshida K, Ichikawa N, Inoue H, Kaneko Y, Kawasaki T, Matsui K, Morita M, Suda M, Tada K, et al. Clinical characteristics of patients with spondyloarthritis in Japan in comparison with other regions of the world. J Rheumatol. 2019;46(8):896–903.

    Article  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Mizuki N, Meguro A, Ota M, Ohno S, Shiota T, Kawagoe T, Ito N, Kera J, Okada E, Yatsu K, et al. Genome-wide association studies identify IL23R-IL12RB2 and IL10 as Behcet’s disease susceptibility loci. Nat Genet. 2010;42(8):703–6.

  49. Gul A. Genome-wide association studies in Behcet’s disease: expectations and promises. Clin Exp Rheumatol. 2011;29(4 Suppl 67):S3–5.

    PubMed  Google Scholar 

  50. Kim SW, Jung YS, Ahn JB, Shin ES, Jang HW, Lee HJ, Il Kim T, Kim DY, Bang D, Kim WH, et al. Identification of genetic susceptibility loci for intestinal Behcet’s disease. Sci Rep. 2017;7:39850.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Kirino Y, Nakajima H. Clinical and genetic aspects of Behcet’s disease in Japan. Intern Med. 2019;58(9):1199–207.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Tsuchida N, Kirino Y, Soejima Y, Onodera M, Arai K, Tamura E, Ishikawa T, Kawai T, Uchiyama T, Nomura S, et al. Haploinsufficiency of A20 caused by a novel nonsense variant or entire deletion of TNFAIP3 is clinically distinct from Behcet’s disease. Arthritis Res Ther. 2019;21(1):137.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Arida A, Vaiopoulos G, Markomichelakis N, Kaklamanis P, Sfikakis PP. Are clusters of patients with distinct clinical expression present in Behcet’s disease? Clin Exp Rheumatol. 2009;27(2 Suppl 53):S48–51.

  54. Yazici H, Pazarli H, Barnes CG, Tuzun Y, Ozyazgan Y, Silman A, Serdaroglu S, Oguz V, Yurdakul S, Lovatt GE, et al. A controlled trial of azathioprine in Behcet’s syndrome. N Engl J Med. 1990;322(5):281–5.

    Article  CAS  PubMed  Google Scholar 

  55. Krause I, Rosen Y, Kaplan I, Milo G, Guedj D, Molad Y, Weinberger A. Recurrent aphthous stomatitis in Behçet’s disease: clinical features and correlation with systemic disease expression and severity. J Oral Pathol Med. 1999;28(5):193–6.

  56. Krause I, Mader R, Sulkes J, Paul M, Uziel Y, Adawi M, Weinberger A. Behçet’s disease in Israel: the influence of ethnic origin on disease expression and severity. J Rheumatol. 2001;28(5):1033–6.

Download references


Not applicable.


The study was supported by the Japan Society for the Promotion of Science (Y.K.19H03700, and M.T. 17 K09990) and partly supported by the Health and Labour Sciences Research Grants (Research on Intractable Diseases) from the Ministry of Health, Labour and Welfare of Japan. The funding source had no role in the design of the study, analysis, interpretation of the data, or the decision to approve publication of the final manuscript. We would like to thank Editage for English language editing.

Author information

Authors and Affiliations



Conception and design: Y. Kirino, M. Takeno. Analysis and interpretation of the data: Y. Kirino, M. Takeno, Y. Soejima, and M. Kurosawa. Drafting of the article: Y. Kirino, M. Takeno, Y. Soejima, and M. Kurosawa. Critical revision for important intellectual content: Y. Kirino, M. Takeno, Y. Soejima, M. Kurosawa, M. Takeuchi, R. Yoshimi, Y. Sugiyama, S. Ohno, Y. Asami, A. Sekiguchi, T. Igarashi, S. Nagaoka, Y. Ishigatsubo, H. Nakajima, and N. Mizuki. Final approval of the article: Y. Kirino, M. Takeno, Y. Soejima, M. Kurosawa, M. Takeuchi, R. Yoshimi, Y. Sugiyama, S. Ohno, Y. Asami, A. Sekiguchi, T. Igarashi, S. Nagaoka, Y. Ishigatsubo, H. Nakajima, and N. Mizuki. Obtaining of funding: Y. Kirino, M. Takeno, and N. Mizuki. Collection and assembly of data: Y. Kirino, M. Takeno, Y. Soejima, and M. Kurosawa. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Yohei Kirino.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee of YCU (A191100002). This survey is based on the application of the opt-out strategy in accordance with Japan MHLW Epidemiological Research Ethics Guidelines.

Consent for publication

Not applicable.

Competing interests

Y. Kirino: Amgen (speaking and consultation fees); M. Takeno: Amgen (advisor for the clinical trial of Apremilast); and S. Nagaoka: Chugai, Bristol-Myers Squibb, Gilead (consultation fees).

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.

Supplementary file.

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Soejima, Y., Kirino, Y., Takeno, M. et al. Changes in the proportion of clinical clusters contribute to the phenotypic evolution of Behçet’s disease in Japan. Arthritis Res Ther 23, 49 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: