A single nucleotide polymorphism in the IL1RL1 gene is associated with Behcet’s disease in a Chinese Han population

Xin-Shu Liu, Zi-Yan Wu, Si Chen, Chan Zhao, Fei Gao, Ming-Hang Pei, Shan-Shan Jia,Yong-Zhe Li, Pei-Zeng Yang, Mei-Fen Zhang

1Department of Ophthalmology, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing 100730, China

2Department of Rheumatology and Clinical Immunology,Peking Union Medical College Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College,Key Laboratory of Rheumatology and Clinical Immunology,Ministry of Education, Beijing 100730, China

3Department of Clinical Laboratory, Beijing Anzhen Hospital,Capital Medical University, Beijing 100029, China

4Department of Clinical Laboratory, Peking Union Medical College Hospital, Peking Union Medical College and Chinese Academy of Medical Sciences, Beijing 100730, China

5The First Affiliated Hospital of Chongqing Medical University, Chongqing Key Laboratory of Ophthalmology and Chongqing Eye Institute, Chongqing 400016, China

Abstract

● KEYWORDS: Behcet’s disease; single nucleotide polymorphism; Chinese Han population; IL33; IL1RL1

INTRODUCTION

B ehcet’s disease (BD) is defined as a kind of chronic recurrent systemic vasculitis; its most common manifestations are aphthous ulceration, skin lesions, genital ulcers and ocular inflammation[1]. BD distributes worldwide and is especially prevalent in Mediterranean countries,the Middle East and Southeast Asia[2]. The etiology of BD remains poorly understood. Considerable evidence indicates that the immunopathogenesis of BD is critical for clarifying the initiation and progression of this disease. Cytokines involved in Th1 and Th17, such as IFN-γ, IL-12, and IL-17,were exhibited according to remissions and exacerbations of inflammation in BD patients[3]. Currently, it is supposed that the effect of environmental risk factors on genetically susceptible individuals may be a trigger of this pathological process.Polymorphisms in theIL10,IL12,IL23andIL37genes have been discovered to be associated with the disease[4-5].IL-33 belongs to the IL1 cytokine family. After binding with

ST2, the complex activates the downstream NF-κB and MAPK pathways[6]. Recent studies indicate the degree of participation of the IL-33/ST2 axis in various diseases, especially in immune and inflammatory disorders, such as rheumatoid arthritis(RA)[7], giant cell arteritis (GCA)[8], Grave’s disease (GD)[9],inflammatory bowel disease (IBD)[10]and systemic sclerosis(SSc)[11]. Correspondingly, genetic polymorphism studies have identified IL33 and/or IL1RL1 loci as susceptibility genes in RA[12], GCA[13], autoimmune thyroid diseases (AITD)[14]and IBD[15]. Several studies have displayed a significant change in IL-33 and/or ST2 levels in the peripheral circulation and/or at inflammatory sites[16], but few studies have been conducted to explore the genetic predisposition of IL33/IL1RL1 to BD.

In this study, we hypothesized that IL33/IL1RL1 gene polymorphisms may be associated with genetic susceptibility to BD in the Chinese Han population.

SUBJECTS AND METHODS

Ethical ApprovalThis study was approved by the Institutional Review Board of the Peking Union Medical College Hospital and the First Affiliated Hospital of Chongqing Medical University and adhered to the tenets of the Declaration of Helsinki. All participants signed written informed consent forms.

SubjectsA total of 783 BD patients and 701 ethnically matched healthy controls who visited Peking Union Medical College Hospital and the First Affiliated Hospital of Chongqing Medical University between October 2011 and October 2015,were consecutively recruited in this case-control study. All subjects were Han nationality Chinese and were not related to one another. Patients who fulfilled the criteria for the diagnosis of BD[17]were enrolled as cases, while those concomitant with other autoimmune or inflammatory diseases, such as systemic lupus erythematosus (SLE) and RA, were excluded. Healthy controls without any autoimmune or inflammatory disorders were included during their physical examination.

Selection of Single Nucleotide PolymorphismsFour single nucleotide polymorphisms (SNPs; rs11792633, rs7025417,rs10975519, and rs1048274) in the IL33 gene, which had previously shown associations with BD[18]or other autoimmune diseases[13,19], and four tag SNPs (rs2310220, rs12712142,rs13424006, and rs3821204) in the IL1RL1 gene were selected for subsequent analyses.

Tag SNPs in the IL1RL1 gene (chr2:102294394-102334929)were identified by Haploview 4.2 software from HapMap CHB data (HapMap Data Rel 27 PhaseII+III, Feb09), with a pairwise linkage disequilibrium (LD) ofr2≥0.8 and minor allele frequency (MAF) values ≥0.1.

GenotypingGenomic DNA of the participants was extracted from EDTA peripheral venous blood by using a DNA isolation kit (Tiangen, Beijing, China) following the manufacturer’s instructions. The SNPs were genotyped by the Sequenom MassARRAY system (San Diego, CA, USA) according to standard procedures. Specifically, first of all, designed the primers for multiplex polymerase chain reaction (PCR) and locus-specific single-base extension with the MassARRAY Assay Design 4.0 software; second, carried out the PCRs, the products were used for locus-specific single-base extension reactions; third, desalted and transferred the final products to a 384-element Spectro CHIP array for allele detection, which was performed by matrix-assisted laser desorption ionizationtime-of-flight mass spectrometry (MALDI-TOF MS). Finally,the resultant mass spectrometry data was analyzed by using MassARRAY Typer 4.0 software.

Statistical AnalysesStatistical analyses were mainly accomplished by PLINKv1.07 software (http://pngu.mgh.harvard.edu/purcell/plink/)[20]. Hardy-Weinberg equilibrium(HWE) in the control population was assessed by the Chisquare (χ2) test for each SNP. Any SNP with significant deviation from HWE was excluded from subsequent analyses.The basic analyses of allele frequencies and genotype distributions were performed by theχ2test. For additional genotype analyses under additive, dominant and recessive model, the Logistic regression test was used. Statistical power was calculated by a freely available power and sample size calculation program (PSv.3.1.2, http://biostat.mc.vanderbilt.edu/wiki/Main/PowerSampleSize)[21]. Haplotype analyses were performed by Haploviewv4.2 software (http://www.broadinstitute.org/haploview)[22].Pvalues less than 0.05 were considered statistically significant.Pvalues for multiple comparisons were corrected by the Bonferroni method(Pc=P×n, n was the number of tested SNPs).

RESULTS

Clinical Features of the ParticipantsThe baseline demographic and clinical features of the participants were displayed in Table 1. A total of 783 BD patients (63.6% male)and 701 (55.5% male) ethnically matched healthy controls were recruited in the present study. The mean ages of the patients and controls were 35.8±11.2 and 38.4±10.2 years old,respectively.

Association Analyses of the Single Nucleotide PolymorphismsIL1RL1rs3821204 was excluded from further analyses because of deviation from HWE in the control group (P＜0.05). The remaining seven SNPs were all in HWE.The average genotyping rate of the seven SNPs was over 98%. The sample size provided a statistical power of 80.1%(α=0.05) for detecting the association between rs12712142 and BD based on the odds ratio (OR) and MAF value of the present study.

The allele frequencies and genotype distributions of SNPs in theIL33andIL1RL1gene regions in BD patients and controls are presented in Table 2. The frequency ofIL1RL1rs12712142variant allele A was significantly lower in BD patients than that in controls (29.4%vs34.1%, OR=0.8, 95%CI: 0.69-0.94,Pc=0.039; Table 2). The genotype distribution of rs12712142 was also significantly different between patients and controls(Pc=0.043; Table 2). Further Logistic regression analyses under additive, dominant and recessive model are displayed in Table 3, and rs12712142 was associated with BD based on the additive and dominant model (OR=0.8, 95%CI: 0.69-0.94,Pc=0.040 and OR=0.72, 95%CI: 0.58-0.88,Pc=0.011; Table 3).However, no significant difference was observed in allele frequencies or genotype distributions in theIL33SNPs(rs11792633, rs7025417, rs10975519 and rs1048274) and otherIL1RL1SNPs (rs2310220 and rs13424006) between BD patients and controls (allPc＞0.05; Table 2), and further Logistic regression analyses based on additive, dominant and recessive model revealed no significant associations of the above SNPs with BD patients either (allPc＞0.05; Table 3).

Table 1 Demographic and clinical data of BD patients and controls

Haplotype Analyses of the IL33 Single Nucleotide PolymorphismsThe haplotype distributions of SNPs inIL33were analyzed by Haploview software. Pairwise LD was observed for rs11792633, rs10975519, and rs1048274 (r2＞0.8)in our data, which is shown in Figure 1. However, none of the distributions of the three haplotypes (TTA, CCG, TCG)formed by the above SNPs indicated any significant difference between BD patients and controls (allPc＞0.05; Table 4).

DISCUSSION

Figure 1 LD analyses of the SNPs in the IL33 gene region The LD plots were generated by Haploview software v4.2 using our data.The numbers (divided by 100) in the small squares represent r2 value and range from 0 to 1. The three SNPs (rs11792633, rs10975519, and rs1048274) in IL33 reside in one LD block.

To the best of our knowledge, this is the first hospitalbased case-control study conducted in China describing the relationship betweenIL33/IL1RL1gene polymorphisms and BD. In our study, the results demonstrated that theIL1RL1rs12712142 polymorphism was associated with BD in a Chinese Han population for the first time. In our cohort, the frequency of variant allele A ofIL1RL1rs12712142 in BD patients was significantly lower than that in healthy controls.Accordingly, the basic genotype distribution and analyses under additive and dominant models also showed significant differences between BD patients and controls. This result suggested that variant allele A ofIL1RL1rs12712142 seemed to be protective against BD.IL-33, encoded by theIL33gene, is expressed constitutively in endothelial and epithelial cells nucleus[23]. ST2, encoded by theIL1RL1gene, was identified as a receptor of IL-33[24]. ST2 has two isoforms: ST2 (a membrane-bound form) is expressed in immune cells such as dendritic cells, Th1 cells and Th2 cells,binds to IL-33 and then activates the NF-κB pathway, while sST2 (a soluble form) acts as a decoy receptor[24-25]. Emerging evidence suggests that the IL-33/ST2 pathway plays a vital role in autoimmune and inflammatory diseases. Increased levels of IL-33 were observed in the synovial fluid of RA patients[7].Animal models showed that at the onset of collagen-induced arthritis, disease severity could be attenuated by administration of anti-ST2 antibody; at the meanwhile, joint destruction could be reduced and a marked decrease in IFN-γ production was observed[26]. Two of the previous studies showed elevated serum IL-33 and sST2 in BD patients compared to those in healthy controls; moreover, serum IL-33 levels were higher in active BD patients than those in inactive BD patients[27-28].SNPs (rs1342326, rs7044343, and rs11792633) in theIL33gene region were also shown to be associated with BD[18,29].

In our cohort, the frequency of variant allele A ofIL1RL1rs12712142 in BD patients was significantly lower than that in healthy controls. We assumed that the protective role of variant allele A ofIL1RL1rs12712142 was mediated by decreasing level of ST2 and/or increasing level of sST2. A recent study supports our assumption that genetic polymorphisms may affect the expression of ST2. The study firstly showed a lower frequency of theIL18R1rs12987977 variant allele G in BD patients, then revealed a downregulation ofIL1RL1in carriersof the protective homozygous rs12987977/GG genotype compared with the TT genotype in functional experiments[30].Moreover, Hoet al[31]identified that genetic factors determined 45% sST2 production variation and that genetic variation inIL1RL1could lead to increased level of sST2. Higher decoy sST2 may be protective against BD, and increased level of serum sST2 in BD patients[28]may be the result of a compensatory protective reaction of the human body.

Table 2 Allele frequencies and genotype distributions of the IL33 and IL1RL1 gene markers in BD patients and controls

Table 3 Analyses of the seven SNPs based on additive, dominant, and recessive genetic models

Table 4 Haplotype analyses of IL33 SNPs between BD patients and controls

None of the SNPs in theIL33region showed a significant association with BD in our study at the allelic, genotypic or haplotypic levels. Although the variant allele T of rs11792633,which has been reported to be protective for BD in a Turkish population[18], was also lower in BD patients than in the controls from our data, the difference was not statistically significant. This inconsistency may be attributed to the genetic heterogeneity of different ethnic groups.

Despite the relatively large sample size in the current study,we only assessed four SNPs inIL33that have been previously reported to be associated with BD or other autoimmune diseases and four tag SNPs inIL1RL1that may not completely represent the whole genetic region. We did not examine the function ofIL1RL1rs12712142in vivoorin vitro, either. More genetic and mechanistic studies are warranted to determine the role ofIL33/IL1RL1in BD pathogenesis.

In conclusion, this study clearly demonstrates the association of theIL1RL1polymorphism with BD in a Chinese Han population. The A variant inIL1RL1rs12712142 is correlated with a decreased risk of BD, which may suggest a protective role of theIL1RL1gene in the pathogenesis of BD.

ACKNOWLEDGEMENTS

Foundation:Supported by the National Natural Science Foundation of China (No.81770917).

Conflicts of Interest: Liu XS,None;Wu ZY,None;Chen S,None;Zhao C,None;Gao F,None;Pei MH,None;Jia SS,None;Li YZ,None;Yang PZ,None;Zhang MF,None.