Background An Immunochip study recently identified the association of a number of new genetic loci with Behcet’s disease (BD).
Objective To confirm the association between new genetic loci reported in an Immunochip study and BD in a Han Chinese population.
Methods A two-stage association study was carried out in 1238 patients with BD and 1458 healthy controls. Twenty-two candidate single nucleotide polymorphisms (SNPs) were selected for genotyping by iPLEXGold genotyping or TaqMan SNP assays and a meta-analysis was performed for significantly associated markers.
Results The results showed that four SNPs (LACC1/rs9316059, CEBPB-PTPN1/rs913678, ADO-EGR2/rs224127 and RIPK2/rs10094579) were associated with BD in an allelic association test (rs9316059 T allele: pc=4.95×10−8, OR=0.687; rs913678 C allele: pc=3.01×10−4, OR=1.297; rs224127 A allele: pc=3.77×10−4, OR=1.274; rs10094579 A allele: pc=6.93×10−4, OR=1.302). For four SNPs tested by meta-analysis, the association with BD was strengthened and all exceeded genome-wide significance (rs9316059: p=2.96×10−16; rs913678: p=2.09×10−16; rs224127: p=5.28×10−13; rs10094579: p=9.21×10−11).
Conclusions Our findings confirmed the association of four loci (LACC1, CEBPB-PTPN1, ADO-EGR2 and RIPK2) in Chinese Han patients with BD.
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Statistics from Altmetric.com
BD is a chronic systemic vasculitis that mainly presents with recurrent uveitis, oral ulcers, genital ulcers and multiple skin lesions.1 BD is more common among countries along the ‘silk route’ from the Mediterranean, Middle East, China and Japan, but is rare in the USA and Europe.2 Although the aetiology and pathogenesis of BD remain unclear, it is currently thought that both genetic and environmental factors contribute to disease occurrence and development. In addition to HLA-B*51 which has been shown to have the strongest association with BD,3–7 a series of genome-wide association studies in different populations have identified a number of non-human leucocyte antigen susceptibility loci for BD, including IL23R-IL12RB2, IL10, STAT4, CCR1-CCR3, KLRC4, ERAP1, TNFAIP3, IL12A and FUT2.3–11 These findings have increased our understanding of immunogenetic factors involved in the disease. However, these identified genetic risk loci do not fully explain the genetic aetiology of BD, and other genetic factors remain to be identified.
Recently, Takeuchi et al 12 conducted a genetic association study using the Immunochip genotyping array in a Turkish cohort (1900 patients with BD and 1779 controls) and further confirmed the associations of BD with HLA-B*51, IL23R-IL12RB2, IL10, CCR1, KLRC4, ERAP1, IL12A and FUT2. More significantly, the same study identified six new BD risk loci (IL1A-IL1B, IRF8, CEBPB-PTPN1, ADO-EGR2, RIPK2 and LACC1) with genome-wide significance (p<5×10−8) and a number of new loci with a suggestive disease association (p<5×10−5). Takeuchi et al also performed a replication study12 and confirmed the association for some, but not all, loci in Iranian and Japanese cohorts. These findings indicated that the genetic background of BD may differ among different ethnic groups. Since the novel loci mentioned above have not yet been tested in other populations, we decided to perform a replication study to assess whether the findings of the Takeuchi study could be confirmed in Chinese Han patients with BD.
Materials and methods
A total of 1238 patients with BD and 1458 healthy controls were included in the present study. All the patients were Han Chinese and recruited from the Department of Ophthalmology of the First Affiliated Hospital of Chongqing Medical University (Chongqing, China) from June 2008 to December 2016. BD was strictly diagnosed based on the criteria of the International Study Group for BD13 and all patients had uveitis. All control subjects were matched for age, sex, ethnicity (Han Chinese) and geographic origin with patients with BD. The present study was conducted in accordance with the tenets of the Declaration of Helsinki. Informed consent was obtained from all participating individuals.
We selected 27 candidate SNPs in 20 susceptibility loci from the BD Immunochip association study.12 Criteria used were as follows: we selected the lead SNPs in the potential susceptibility loci with a p value less than 5×10−5. In order to increase the reliability of our results, the other susceptibility SNPs that were identified by meta-analysis and imputation were also included in this study.
Some SNPs were excluded from this study for the following reasons. (1) The SNPs in the loci that have been reported previously by our team (IL10, IL23R-IL12RB2, CCR1, ERAP1, KLRC4, IRF8, FOXP1).14–18 (2) The SNPs that were not polymorphic in Han Chinese. (3) The SNPs that were in high linkage disequilibrium (LD) with the lead SNPs in Chinese Han (r2>0.8, Han Chinese in Beijing, HCB) (table 1).
Genotyping of the selected 22 SNPs was performed in 478 BD cases and 662 controls drawn from a Chinese Han population in a first-stage study, and another independent cohort including 760 BD cases and 796 controls was examined in a second-stage confirmation study. Genotyping of 21 SNPs was performed using the MassARRAY platform (Sequenom, California, USA) and iPLEX Gold Assay. The SNPs rs1047781was genotyped by TaqMan SNP Genotyping Assay (Applied Biosystems, Foster City, California, USA) and the probe fluorescence signal was detected using the 7500 Real-Time PCR System (Applied Biosystems, USA). All SNPs tested had a high call rate (≥95% in all individual) and conformed to Hardy-Weinberg equilibrium (HWE) in the normal controls (p for HWE ≥0.05).
The χ2 test was applied for the evaluation of the HWE. Genotype and allele frequencies were compared between patients with BD and normal controls by the χ2 test or Fisher’s exact test using SPSS V.17.0. The Bonferroni correction method was applied for the correction of p values for multiple comparisons. HWE was tested by the SHEsis website. Meta-analysis in multiple populations was performed using STATA software V.12.0. The p value of heterogeneity and I2 were calculated to evaluate heterogeneity between populations. Phet<0.05 and I2 >0.5 were considered to be significant. Statistical power was estimated from effect size in the original Turkish data sets from the ‘Takeuchi’ study,12 allele frequency and sample size in the Chinese Han population. Statistical power analysis was performed using PS Power and Sample Size Calculations software (V.3.1.2; Department of Biostatistics, Vanderbilt University, Nashville, Tennessee, USA).
Clinical characteristics of patients with BDs
All 1238 patients with BD had uveitis, of which 26% patients had hypopyon. The most frequent type of uveitis was panuveitis (93.2%), followed by posterior uveitis (4.9%) and anterior uveitis (1.9%). Oral ulcers (94.4%) were the most frequent extraocular manifestation, followed by skin lesions (75.2%) and genital ulcers (55.4%). The distribution of age and gender and clinical features of the enrolled patients with BD and the healthy controls in this study are shown in table 2.
Association test of examined SNPs in the first phase
In the first-stage study, 22 SNPs were genotyped in 478 patients with BD and 662 normal controls. Significant higher frequencies of the CEBPB -PTPN1/rs913678 C allele (pc=1.01×10−2, OR=1.382), RIPK2/rs10094579 A allele (pc=3.16×10−2, OR=1.368) and ADO-EGR2/rs224127 A allele (pc=3.36×10−2, OR=1.318) were observed in patients with BD. In addition, the frequency of the T allele for LACC1/rs9316059 was significantly lower in patients with BD (pc=1.53×10−5; OR=0.613) (table 3). However, there was no association between the remaining SNPs and BD in this Chinese Han cohort (online supplementary table 1).
Association test of examined SNPs in the second phase and combined study
To further confirm the outcome of the first-stage study, we enrolled a separate set of 760 patients with BD and 796 healthy individuals for a second-stage test. We only tested the SNPs that showed a significant association in the first phase. The frequency of the T allele for LACC1/rs9316059 in patients with BD was confirmed to be significantly lower (pc=1.34×10−3, OR=0.745). In addition, the result again demonstrated significantly higher frequencies of the CEBPB-PTPN1/rs913678 C allele (pc=4.43×10−2, OR=1.222), RIPK2/rs10094579 A allele (pc=2.74×10−2, OR=1.252) and ADO-EGR2/rs224127 A allele (pc=1.88×10−2, OR=1.232) in patients with BD (table 3).
Combination of the data from the first-stage and second-stage study showed that four SNPs (rs913678, rs9316059, rs10094579 and rs224127) were significantly associated with BD (rs913678 C allele: p=1.37×10−5, pc=3.01× 10−4, OR=1.297; rs9316059 T allele: p=2.25×10−9, pc=4.95×10−8, OR=0.687; rs10094579 A allele: p=3.15×10−5, pc=6.93×10−4, OR=1.302; rs224127 A allele: p=1.71×10−5, pc=3.77×10−4, OR=1.274) (table 3).
Stratified analysis for rs913678, rs9316059, rs10094579 and rs224127
We also analysed whether rs913678, rs9316059, rs10094579 and rs224127 showed an association with the main clinical features of BD. The results did not show any significant association between the four tested SNPs and groups of patients with BD divided according to their clinical features (online supplementary table 2).
To further investigate the risk conferred by the SNPs (rs913678, rs9316059, rs10094579 and rs224127) associated with BD, we performed a meta-analysis of the genetic polymorphisms for which data were available from the Takeuchi association study and our study data sets. The results showed that the disease association of tested SNPs were reinforced after meta-analysis and all exceeded genome-wide significance (rs913678: p=2.09×10−16; rs9316059: p=2.96×10−16; rs10094579: p=9.21×10−11; rs224127: p=5.28×10−13) (table 4).
In the present study, we performed a replication study in a Han Chinese BD cohort for 22 candidate SNPs identified with an association p value <5× 10−5 with BD in a recent Immunochip study.12 The results showed that four SNPs (rs913678 in CEBPB - PTPN1, rs9316059 in LACC1, rs10094579 in RIPK2, rs224127 in ADO-EGR2) contribute to the genetic susceptibility to BD in a Chinese Han population.
LACC1/rs9316059, the most significantly associated SNP with BD in our study, displayed genome-wide significant association (table 3). The Immunochip study was performed in Turkish patients, and the findings concerning LACC1 rs9316059 was also confirmed in a Japanese cohort.12 In the same study, the association with SNP rs2121033 in LACC1 was identified by meta-analysis in three populations (Turkish, Iranian and Japanese). These two SNPs are in strong LD with each other (r2=0.891, in HCB). This study in combination with ours indicates that the protective LACC1 locus is a commonly associated gene for BD in all the populations tested including Chinese Han, Turkish, Iranian and Japanese. In addition, we identified the susceptibility SNP rs10094579 in RIPK2 in our cohort (table 3). To our knowledge, this is the first report showing that rs10094579 in RIPK2 confers risk to BD. This SNP only showed a suggestive disease association with a Turkish BD cohort but was not confirmed both in an Iranian and Japanese BD cohort.12 Moreover, the association we found for the other two SNPs (CEBPB-PTPN1/rs913678 and ADO-EGR2/rs224127) are in agreement with data in an Iranian population for CEBPB-PTPN1/rs913678 and a Japanese population for ADO-EGR2/rs224127 (table 3). In a meta-analysis of populations (table 4 and figure 1), we show that ADO-EGR2/rs224127, CEBPB-PTPN1/rs913678 and PIPK2/s10094579 all exceeded genome-wide significance. Based on these findings, we propose that LACC1, CEBPB–PTPN1, RIPK2 and ADO-EGR2 constitute BD susceptibility genes in Chinese Han, together with other established loci such as IL10, IL23R–IL12RB2, CCR1, IRF8, KLRC4, STAT4, ERAP1, TNFAIP3, TNFSF4, UBAC2, IL-37, IL-18RAP, GAS6, PROS1, CD6, CD11c, ATG5, TRAF5, TRAF3IP2, JAK1, MIF, PDGFRL, CD40, CIITA, NOD1, NOS3, REL and TLR2.14–34
Of note, we also did not find evidence supporting the disease association with two other reported loci (IL1A-IL1B and FUT2), which both showed a genome-wide association (p<5×10−8) in the original ‘Immunochip’ report.12 IL1A-IL1B/rs3783550 that was identified in Turkish patients could also not be confirmed in Iranian as well as Japanese patients with BD. We could also not confirm the association with the FUT2/rs1047781 (T) allele, which is an ancestry-specific FUT2 non-secretor mutation (in Japanese and Han Chinese) with a significant association with BD in Japanese.12 The fact that the two SNPs mentioned above show a lack of association in Chinese Han is probably due to the different genetic background between Chinese Han and Japanese and Turkish populations, since our sample size was large enough to find a possible existing association (Power>0.8) (online supplementary table 1). Populations between different continents show between 16% and 19% genetic differences. Even within a continent, populations may differ genetically, whereby Japanese and Chinese, differ by 6.78%.35 Similarly, we did not confirm an association with the other two SNPs (rs1509966 and rs7075773) in ADO-EGR2 and BD susceptibility. Although they had a relatively high statistic power of 0.641 and 0.738, respectively, we could not exclude a false-negative disease association of rs1509966 and rs7075773 in our study. In addition, the difference of the associations for three SNPs in ADO-EGR2 among the nationalities may be partly explained by the population genetic heterogeneity. Likewise, RIPK2/rs2230801, another SNP identified in Turkish patients and confirmed in Japanese patients, did not show a significant association with BD in our study. Given that rs2230801is a rare variant in China (minor allele frequency<0.05), the statistical power to confirm this finding was low (0.23) (online supplementary table 1).
A recent study reported that the LACC1/rs3764147 (p.Ile254Val) is in high LD with rs9316059 (r2=0.892, HCB) and leads to impaired protein function.36 Furthermore, Lacc1−/− mice produce decreased IL-1β in response to lipopolysaccharide treatment, consistent with a role for IL-1β in BD pathogenesis.36 The minor allele of the other three SNPs increase the risk for BD. The CEBPB-PTPN1/rs913678 C allele is associated with decreased gene expression, and Cebpb−/− mice show increased susceptibility to pathogens.37 38 The RIPK2 kinase transduces signalling downstream of the intracellular peptidoglycan sensors NOD1 and NOD2 to promote a productive inflammatory response.39 However, excessive NOD2 signalling has been associated with numerous diseases, including inflammatory bowel disease (IBD), sarcoidosis and inflammatory arthritis.40–42 Interestingly, ADO-EGR2 has also been identified as a risk for Vogt-Koyanagi-Harada (VKH) syndrome by a previous genome-wide association study of our team.43 ADO and EGR2 were all expressed in the iris, whereas EGR2 was also expressed in ciliary body and choroid.43 VKH syndrome and BD are two of the most common types of uveitis in Chinese Han, and the fact that they share common susceptibility loci suggests that ADO-EGR2 may be a common genetic locus for uveitis, which may provide a theoretical basis for prevention and treatment of other type of uveitis.
The four novel BD susceptibility loci that we could confirm in Chinese patients have also been reported to be associated with other immune disorders. ADO-EGR2, LACC1 and CEBPB-PTPN1 are shared by BD and IBD.41 44–46 RIPK2, ADO-EGR2 and LACC1 have been shown to be associated with leprosy.47–50 These observations suggest that these diseases, whether being autoinflammatory (BD and IBD) or infectious (leprosy), may share molecular pathways although their exact role (protective or susceptibility) may differ markedly. The C allele of rs913678 in CEBPB-PTPN1 confers risk of BD but was protective for IBD and ulcerative colitis (UC).41 46 The minor allele of LACC1/rs9316059 confers protection for BD but is in high LD with a common coding variant, rs3764147 (r2=0.892, HCB), which increases risk for IBD, Crohn’s disease (CD) and leprosy.41 49 A similar discrepancy is also seen in RIPK2, where the minor allele of rs10094579 conferred risk of BD but is in high LD with rs7015630 (r2=0.976, HCB) that was protective for CD and leprosy.41 49 These discordant observations suggests that these genes which are involved in various signalling pathways may play opposite roles in BD as compared with IBD and leprosy. Further functional investigations may help to elucidate the molecular mechanisms underlying BD development and increase our understanding on the impact of these loci in the pathogenesis of autoinflammatory and infectious diseases.
Our study has several limitations. We only replicated the lead SNPs from the loci identified by a previous Immunochip study12 and we cannot exclude that other suggestive SNPs may show an association with BD in Chinese Han. In addition, there were four genes that did not show informative results. These included LONRF2, RBM6 and C5orf56 since the only SNP included was not polymorphic, and IL12A, where the two SNPs included, had one that was not polymorphic and data analysis of the other SNP IL12A/rs2647935 showed a statistical power that was too low to obtain a meaningful conclusion. Further, fine mapping for these gene regions is needed in Chinese Han to definitely show whether this is not a false-negative association. It should also be noted that all the patients with BD in the present study suffered from uveitis (100%), whereas the Turkish, Iranian and Japanese patients from the ‘Immunochip’ association study showed a lower uveitis incidence (39.4%, 56.3%, 86.9%, respectively),12 which indicates that there may be a selection bias towards ocular BD in our study. One should also be aware of the fact that most of our patients with BD are male (80.9%), which is in agreement with the previous reports from countries along the ancient Silk Road,51–53 whereas studies from Europe or the USA often show an almost equal gender distribution.54–56
In conclusion, our study not only confirms the association of LACC1/rs9316059, CEBPB-PTPN1/rs913678 and ADO-EGR2/rs224127 with BD but also identifies a novel RIPK2/rs10094579 polymorphism that affects BD susceptibility in Chinese Han. Our findings are an addition to the growing body of data from different ethnic populations, thereby gradually revealing the genetic risk landscape of BD. Further investigations on how these gene polymorphisms exactly affect BD are needed and may provide future targets for its treatment.
The authors thank the patients and controls for their participation in this study.
PW, LD and SH contributed equally.
Contributors PW, LD and SH are joint first authors and anlaysed the data. PW and PY designed the study. PW, LY and JH collected the data. PY made clinical diagnoses. GS, LY, JH, JD, QC and CZ collected the samples. LY, JH, JD, QC and GY extracted the blood DNA. PW drafted the manuscript. PY, SH, GS and AK helped revise the manuscript. All authors have read and approved the final manuscript.
Funding This work was supported by Major Research Development Program of China (2016YFC0904000), Natural Science Foundation Major International (Regional) Joint Research Project (81720108009), Natural Science Foundation Major International (Regional) Joint Research Project (81320108009), National Natural Science Foundation Project (81522013), Chongqing Outstanding Youth Grant (cstc2014jcyjjq10005), Chongqing Key Laboratory of Ophthalmology (CSTC2008CA5003), Chongqing Science and Technology Platform and Base Construction Program (cstc2014pt-sy10002) and Research fund for Traditional Chinese Medicine of Chongqing Health and FamilyPlanning Commission (ZY201401013). Thanks to all donors enrolled in the present study.
Competing interests None declared.
Patient consent Parental/guardian consent obtained.
Ethics approval This study was approved by the ethics committee of the First Affiliated Hospital of Chongqing Medical University, Chongqing, China (Permit Number 2009–201008).
Provenance and peer review Not commissioned; externally peer reviewed.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.