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

Replication study confirms the association between UBAC2 and Behçet's disease in two independent Chinese sets of patients and controls

Abstract

Introduction

The purpose of this study was to replicate genetic factors associated with the susceptibility to Behçet's disease (BD). We conducted a two-stage candidate genes association and functional study, involving 477 BD patients and 1,334 normal controls of Chinese Han descent.

Methods

The genotyping of five candidate genes/loci, including LOC100129342, KIAA1529, CPVL, UBASH3B and UBAC2, were performed using TaqMan single nucleotide polymorphism (SNP) assays. Real-time PCR and luciferase reporter assay were performed to test the function of the identified promoter polymorphism. The main outcome measures were genotype frequencies and expression levels in BD patients.

Results

The first-stage study results showed that UBAC2 (rs9513584, Pc = 0.018, OR = 1.4), but not LOC100129342, KIAA1529, CPVL, UBASH3B was associated with the susceptibility to BD in Chinese Han. The fine-mapping association study of UBAC2 identified six risk SNPs for BD in the Chinese cohort; three of them were verified in validation study (rs3825427, first-stage Pc = 2.2 × 10-3, second-stage Pc = 9.3 × 10-3, combined Pc = 6.9 × 10-6; rs9517668, first-stage Pc = 1.7 × 10-3, second-stage Pc = 0.03, combined Pc = 3.3 × 10-4; rs9517701, first-stage Pc = 5.1 × 10-3, second-stage Pc = 9.0 × 10-3, combined Pc = 2.9 × 10-5; respectively). Functional analysis showed that the risk T allele of the promoter polymorphism rs3825427 had a significantly lower promoter activity than the non-risk G allele (P = 0.002) and a decreased expression of UBAC2 transcript variant 1 in peripheral blood mononuclear cells (PBMCs) and skin of normal controls carrying the risk T allele than that in individuals with the G allele (P = 0.045, P = 0.025; respectively). The mRNA expression of UBAC2 transcript variant 1 was significantly decreased in PBMCs and skin of BD patients as compared with controls (P = 0.025; P = 0.047, respectively). The mRNA expression of UBAC2 transcript variant 2 was significantly increased in skin of BD patients as compared with controls (P = 0.004).

Conclusions

This study replicates a predisposition gene to BD, UBAC2, and suggests that UBAC2 may be involved in the development of BD through its transcriptional modulation.

Introduction

Behçet's disease (BD) is generally considered as a refractory multisystem disorder characterized by recurrent oral ulceration, genital ulceration, recurrent uveitis and multiple skin lesions [1, 2]. BD presents a significant geographical distribution, mainly seen in the countries along the ancient silk route [3]. It is one of the most common and severe sight-threatening uveitis entities in China [3]. Although the etiology of this disease remains unclear, it is currently thought that a genetic predisposition coupled with a triggering event seems to lead to its development. Although HLA-B51 has been found to be strongly associated with BD [47], it is estimated that HLA-B51 only accounts for 20% of the risk genes for this disease [8]. A whole genome linkage study has implicated several potential non-HLA contributory loci/genes for BD [9]. Candidate gene association studies have identified a number of susceptible genes for BD, including IL23R [10], CTLA4 [11], eNOS [12], SUMO4 [13, 14], NOD2 [15], IL1A [16], manganese superoxide dismutase (SOD) [17] and so on. These studies tremendously enhanced our understanding of the genetic etiology of BD. However, these identified risk genes only contribute to about 30% of the predisposition for BD [18].

Fei and coworkers [19] conducted a genome-wide association study (GWAS) on 152 Turkish Behçet's disease patients using the DNA pooling approach and reported the association of BD with five novel loci/genes, including LOC100129342, KIAA1529, CPVL, UBASH3B and ubiquitin-associated domain containing 2 (UBAC2), among which the genetic association between UBAC2 and BD was lately confirmed and functional significance was tested by the same study group [20]. As far as we know, the association result needs to be confirmed by further replication studies, particularly in other ethnic groups. We, therefore, conducted a candidate association to investigate the association between the identified risk genes/loci and BD in the non-Turkish ethnic population. Our result revealed that UBAC2 was associated with BD in Han Chinese. Fine mapping and validation studies further confirmed UBAC2 as susceptibility gene for Behçet's disease. Further functional study showed that the risk T allele of rs3825427 in UBAC2 was associated with the decreased promoter activity and its mRNA expression. These genetic associations and functional studies suggest that UBAC2 is a risk gene for Behçet's disease.

Materials and methods

Ethics statement

Written informed consent was obtained from each participant, and this study was approved by the Clinical Research Ethics Committee of the First Affiliated Hospital of Chongqing Medical University and Zhongshan Ophthalmic Center, Sun Yat-sen University (Permit Number: 2009-201004).

Subjects

One hundred forty-seven BD patients and 951 normal controls enrolled in the first-stage were recruited from April 2005 to March 2008 in Southern China at Zhongshan Ophthalmic Center, Sun Yat-sen University. The replication cohort is comprised of 330 Han Chinese BD patients and 383 Han Chinese normal controls, which were recruited from April 2008 to March 2010 in Western China at the First Affiliated Hospital of Chongqing Medical University (Table 1). The diagnosis of BD was based on the criteria of the International Study Group for BD [2]. Twelve patients with typical nongranulomatous uveitis, multiform skin lesions or genital ulcers in association with arthritis or other clinical features, were diagnosed as having BD according to the criteria of the Behçet's Disease Research Committee of Japan and were also included in this study (Table 2) [21]. The control populations consisted of unrelated healthy individuals from the same geographical regions as the BD patients and were age, sex and ethnically matched with the patients.

Table 1 Summary characteristic of Behçet's patients and normal control subjects in this study
Table 2 The clinical characteristics of Behcet patients

SNPs selection

We chose the enrolled 25 SNPs based on the following principle: we focused on the UBAC2 gene and ran Haploview software (Daly Lab at the Broad Institute, Cambridge, MA, USA) using Chinese Han Beijing data to screen out tagSNPs including rs3825427, rs1927726, rs9517668, rs9554581, rs4636771, rs7332161, rs9517701, rs7325747, rs912129, rs4772190 and also included the reported BD associated SNP rs9513584 [19]. In order to enhance the accuracy rating of genotyping, we also included 16 SNPs with high heterozygosity (allele frequency > 0.1).

The three SNPs, including rs3825427, rs9517668 and rs9517701, were enrolled in the replication study because of their smaller Bonferroni corrected P-value. The call rate of SNP rs9517699 did not meet the quality standard in the replication study (call rate < 0.80), so SNP rs9517699 was excluded in the replication analysis (Additional file 1).

SNP genotyping

Genomic DNA was extracted from the peripheral blood of patients and controls using the QIAamp DNA Blood Mini Kit (QIAGEN Inc., Hilden, Germany) according to the manufacturer's instructions. The Genotyping was determined by TaqMan® SNP Genotyping Assays (Applied Biosystems, Foster City, CA, USA) on the Applied Biosystems 7500 Real-Time PCR System according to the manufacturer's instructions. All SNPs tested in this study had a genotyping success rate > 98% and accuracy > 99% as judged by random re-sequencing of 20% of samples in all subjects.

Real-time quantitative PCR analysis

Anticoagulated blood samples were obtained using vacuum tubes with EDTA. PBMCs were prepared from venous blood of BD patients (n = 8) and normal controls (n = 52) by Ficoll-Hypaque density-gradient centrifugation. Skin was obtained from the skin lesions of patients (n = 4) and the scrotum of normal controls (n = 22). Total RNA was isolated from skin and PBMCs of the patients or controls by using the QIAamp® RNA Blood Mini kit (QIAGEN Inc., Hilden, Germany) or Qiagen RNeasy Fibrous Tissue Mini Kit (QIAGEN Inc., Hilden, Germany) with treatment of Dnase I, according to the manufacturer's instructions and reversed into cDNA according to the Superscript protocol (SuperScript III First-Strand Synthesize System, Invitrogen, Carlsbad, CA, USA). To compare the mRNA expression of the two protein-coding transcript variants of UBAC2, Real-time Quantitative PCR was performed using the Applied Biosystems 7500 System (Applied Biosystems) with the following primers (UBAC2F:5'CCG GCT CCA GTG GGC TCT ACA3'/UBAC2R:5'GGG CGA GCA GGA GGG AGA GG3') for the UBAC2 transcript variant 1 to generate a 81 bp product and (UBAC2F: 5'TAGGAAGTCGTGGCGAGGGAGC3'/UBAC2R: 5'GCCTTGTCTGCTGAC CACCGCT3') for the UBAC2 transcript variant 2 to generate a 84 bp product. The β-actin (β-actinF: 5'GGA TGC AGA AGG AGA TCA CTG3'/β-actinR:5'CGA TCC ACA CGG AGT ACT TG3') was chosen as the internal reference gene to normalize UBAC2 expression. RT-PCR conditions were one cycle of 95°C for 10 minutes, followed by 40 cycles in which each cycle included 95°C for 15 sec, 60°C for 1 minute, and 95°C for 15 sec, 60°C for 1 minute, 95°C for 15 sec, 60°C for 15 sec.

Luciferase Reporter Assay

The whole genome synthesized promoter sequences of UBAC2 carrying the T allele or the G allele of SNP rs3825427 were cloned into pGL3-basic vector (Promega, Madison, WI, USA). This vector was then transiently transfected into HEK293 cells. Reporter plasmid was transfected to cells using Lipofectamine reagent (Life Technologies, Grand Island, NY, USA). Transfection efficiency was standardized by cotransfecting with pRL-SV40 (Promega). Luciferase activity was determined after 24 hours incubation using a Luciferase Assay system (Beyotime, Jiangsu, China). For each plasmid construct, three independent transfection experiments were performed, and each was done in triplicate.

Statistical analysis

Hardy-Weinberg equilibrium (HWE) was tested using the chi-square test and no SNPs showed significant deviation from HWE (P > 0.05). Odds ratios (OR) and 95% confidence intervals (95% CI) were calculated by using using SPSS version 17.0 (Chicago, IL, USA) to estimate disease risk. Linkage disequilibrium (LD) was examined by using Haploview (version 4.2, Daly Lab at the Broad Institute, Cambridge, MA, USA). To account for multiple testing, the Bonferroni correction was applied. Review Manager 4.2 (Cochrane Collaboration, Oxford, UK) was used to perform meta-analysis. The conditional logistic regression analysis (SAS, 9.13, Cary, NC, USA) was performed to localize the effect and identify the number of independent effects.

Results

The clinical findings of the BD patients included in our study are shown in Table 2. No statistical difference in the distribution of age and gender was observed between BD patients and controls (P > 0.05).

The identification of risk gene for BD by candidate genes association study in Southern China (first-stage study)

Five genes/loci, that is, LOC100129342, KIAA1529, CPVL, UBASH3B and UBAC2, have been identified as risk factors for BD in Turkish. These five candidate genes/loci were not validated by other ethnic populations. We conducted a candidate gene association study to confirm the association result using 147 BD patients and 951 controls. Our results showed that the rs9513584 polymorphism inUBAC2 was associated with the susceptibility to BD (Pc = 0.018, OR = 1.4, meta-analysis OR = 1.5) (Table 3). There was no association of LOC100129342, KIAA1529, CPVL and UBASH3B with BD in the Chinese cohorts (Table 3).

Table 3 The meta- analysis result combined our data and Turkish data

Fine-mapping study and validation study of UBAC2 association in Chinese cohort

The variant in UBAC2 showed a consistent association in the Turkish and Chinese cohort, therefore, UBAC2 was considered as the most important candidate among the five loci. We, therefore, fine-mapped the UBAC2 using 25 SNPs and found that 6 SNPs, that is, rs3825427, rs1927726, rs9517668, rs9554581, rs9517701 and rs9517699, were associated with BD (P ≤0.05/25 = 0.002) (Table 4, Figure 1 and Additional file 1). We next checked the independence of the evidence for multiple associations within the UBAC2 locus by carrying out a conditional logistic regression analysis of the 25 SNPs in UBAC2. After control for the genetic effect of rs3825427, none of the 24 SNPs remained significantly associated with Behçet's disease after correction for the number of SNPs tested by the conditional analysis (the smallest P = 0.0024 > 0.05/25 = 0.002). Therefore, the multiple associations within the UBAC2 locus are not independent.

Table 4 Summary of the association of UBAC2 SNPs with Behcet's disease in Han Chinese population
Figure 1
figure 1

Genetic association analysis in UBAC2 gene region. Overview of SNPs across the UBAC2 gene region in the Chinese cohort. Linkage disequilibrium (r2) to the most significantly SNP (rs9517668, red diamond) is color-coded (red: r2 > 0.8; orange: r2 = 0.5-0.8; yellow: r2 = 0.2-0.5; white: r2 < 0.2). Recombination rates across each region in HapMap CHB are shown in light blue (right y axis). The chromosomal locations and relative positions of genes according to hg18 are shown (x axis).

To further validate the fine-mapping result, we performed a replication study on the selected three SNPs, including rs3825427, rs9517668 and rs9517701, using the West-Southern China cohort. The result showed that all of the three selected SNPs were associated with BD (rs3825427, First-stage Pc = 2.2 × 10-3, OR = 1.5, Second-stage Pc = 9.3 × 10-3, OR = 1.4, combined Pc = 6.9 × 10-6, OR = 1.5; rs9517668, First-stage Pc = 1.7 × 10-3, OR = 1.6, Second-stage Pc = 0.03, OR = 1.3, combined Pc = 3.3 × 10-4, OR = 1.4; rs9517701, First-stage Pc = 5.1 × 10-3, OR = 1.5, Second-stage Pc = 9.0 × 10-3, OR = 1.4, combined Pc = 2.9 × 10-5, OR = 1.4; respectively) (Table 4).

Effects of promoter polymorphism of rs3825427 on UBAC2 transcriptional level

Since the rs3825427 polymorphism is located in the promoter region of UBAC2 and bioinformatics analysis revealed that this G/T SNP changed the binding site for transcriptional factors, we performed a series of functional analyses to examine whether rs3825427 polymorphism affected the expression of UBAC2. We tested two protein-coding transcript variants of UBAC2 in mRNA level with the presence of different genotype (TT/GT/GG) in rs3825427. The mRNA level of UBAC2 transcript variant 1 was found decreased in PBMCs and the skin of normal individuals with the TT genotype compared with those with the GG genotype (Bonferroni corrected P = 0.045, P = 0.025; respectively) (Figure 2). No significant difference was observed in the expression of UBAC2 transcript variant 2 between the various genotypes (Figure 2).

Figure 2
figure 2

Comparison of UBAC2 variants transcriptional level among different genotypes of rs3825427. Five GG genotype, four GT genotype, three TT genotype PBMCs samples and four GG genotype, three GT genotype and five TT genotype skin samples were used to examine the expression of UBAC2 transcript variant 1 and 2 and each sample were assayed three times. The mean ± SD is given for each genotype from three experiments. Statistical significance was taken when Bonferroni corrected two-tailed P < 0.05 using SPSS17.0.

We subsequently performed a Dual Luciferase Reporter Gene Assay to evaluate whether the promoter sequence carrying different alleles have different promoter activity. The promoter sequences of UBAC2 carrying T allele or G allele of SNP rs3825427 were synthesized by Sangon Biotech (Shanghai, P.R. China). Sequencing was performed to validate the result of synthesis and the result showed that only one variant was observed in the location of the rs3825427 polymorphism. The luciferase reporter expression was found to be decreased in the T allele as compared to cells carrying the G allele (P = 0.002) (Figure 3).

Figure 3
figure 3

Effects of the rs3825427 genotype in UBAC2 on luciferase activity in cultured HEK 293 cells. pGL3 luciferase reporter recombinant plasmids containing an UBAC2 promoter sequence with the risk allele T or wild-type G allele at SNP rs3825427 were transfected into HEK293 cells. Renilla luciferase plasmid pTK-SV40 was cotransfected with each construct as an internal control for normalization. The mean ± SD is given for each construct from three experiments. Statistical significance was taken when two-tailed P < 0.05 using SPSS17.0 independent sample t test.

Down-regulated expression of UBAC2 in PBMCs and skin of BD patients

A further study was performed to test whether the mRNA expression of UBAC2 was altered in BD. The results showed that the expression of UBAC2 was significantly decreased in the PBMCs and skin from BD patients as compared with that observed in normal controls (P = 0.025, P = 0.047; respectively) (Figure 4). The mRNA expression of UBAC2 transcript variant 2 was significantly increased in skin of BD patients as compared with controls (P = 0.004) (Figure 4).

Figure 4
figure 4

Real-time RT-PCR analysis of UBAC2 transcript variant 1 and 2 mRNA levels between Behçet's patients and normal controls. UBAC2 transcript variant 1 and 2 mRNA levels derived from PBMCs of 23 normal controls and 6 patients or from skin of 9 normal controls and 4 patients. Each RT-PCR was assayed in triplicate. Significance was examined by using SPSS's independent sample t test.

Discussion

In this study we performed a candidate gene analysis combined with a fine mapping study and the result showed a link of several SNPs in the UBAC2 gene with the susceptibility to Behçet's disease in Chinese Han patients. Functional analysis revealed that the risk-associated T allele of rs3825427 significantly down-regulated the expression of UBAC2 mRNA. A significantly decreased expression of UBAC2 was observed in PBMCs and skin of BD patients compared to normal controls. Our study confirms the association between a single SNP rs9513584 of the UBAC2 gene and BD, which was recently reported in a group of Turkish BD patients. Of interest was the observation that the G allele of rs9513584 was associated with the absence of eye disease in the study in Turkish patients. We could not address this observation since all our patients had eye disease.

The UBAC2 gene, also known as PHGDHL1, is located at 13q32.3. The UBAC2 gene encodes an ubiquitination related structural domain. The study reported by Fei et al. [19] and their subsequent replication showed the association of UBAC2 gene with Behçet's disease and revealed the possible contribution of UBAC2 in the pathogenesis of BD [20]. Our present study successfully replicated the association between UBAC2 and BD in a Chinese cohort. Two other GWAS did not identify the UBAC2 gene as the risk gene for Behçet's disease [22, 23]. The inconsistent result may be partly explained by the population heterogeneity and the disparate SNP array (Affymetrix SNP array, Santa Clara, CA, USA vs Illumina SNP array, San Diego, CA, USA). In a previous study, we also identified a polymorphism in a ubiquitin-related gene, SUMO4, that was associated with the susceptibility to Behçet's disease in a Chinese Han population [13]. Another ubiquitination related gene, UBE2QL1, was also found to be associated with BD in a Turkish cohort [23]. Previous studies suggest that ubiquitination reactions are involved in the regulation of receptor tyrosine kinase signaling and may play important roles in the TNF-α, IL-1β, and TCR-mediated NF-κB activation pathway [24, 25]. NF-κB has been demonstrated to play a crucial role in the pathogenesis of BD through regulating the apoptosis-related factors and increasing the resistance of T cells to apoptosis [26]. All these findings indicate that the ubiquitination-related pathway may have a protective effect in the development of Behçet's disease.

SNP rs3825427 in the UBAC2 gene encodes a promoter polymorphism. Our study showed a close association of this SNP with the susceptibility to BD. To determine whether the promoter polymorphism may modulate the expression of this gene, dual luciferase reporter gene assays and real-time PCR analysis were performed and the results showed that the T allele was associated with decreased promoter activity and may down-regulate the expression of UBAC2. These findings suggest that the rs3825427 polymorphism may lead to the down-regulation of UBAC2 expression. More interestingly, UBAC2 transcription was found to be down-regulated in PBMCs and skin of BD patients as compared to controls. Our results are consistent with that observed in mice deficient in the homologous UBA domain as reported by Carpino et al. [27] They found that UBA knockout mice showed increased cytokine production and also exhibited an increased incidence and severity of experimental autoimmune encephalomyelitis (EAE) as compared with wild-type mice [27].

Several possible limitations of the present study merit particular consideration. First, the patients enrolled in our study were recruited from the Eye Department. Because BD usually affects multiple systems, the patients enrolled in this study might, therefore, represent a separate disease population. Second, the patients enrolled in the replication study were recruited only from Chinese Han individuals and the sample size is relatively small. The results presented here need to be confirmed in other ethnic populations and using larger samples. Finally, although molecular biology study provides functional evidence for the rs3825427 polymorphism, this result is generally suggestive and does not exactly explain how the genetic variant translates into physiologic processes and then affects disease susceptibility. Therefore, the association results presented here should be investigated further using more functional experiments.

Conclusion

In conclusion, our study replicated the association of UBAC2 with BD and identified a promoter SNP of UBAC2, rs3825427, to be associated with the increased risk for BD in Chinese. The functional study showed that this SNP may be involved in the development of BD through transcriptional modulation of UBAC2.

Abbreviations

BD:

Behçet's disease

EAE:

experimental autoimmune encephalomyelitis

EDTA:

ethylenediaminetetraacetic acid

GWAS:

genome-wide association study

HWE:

Hardy-Weiberg equilibrium

LD:

Linkage disequilibrium

OR:

odds ratio

PBMCs:

peripheral blood mononuclear cells

SD:

standard deviation

SNP:

single-nucleotide polymorphism

SOD:

superoxide dismutase

UBAC2:

ubiquitin-associated domain containing 2.

References

  1. Sakane T, Takeno M, Suzuki N, Inaba G: Behcet's disease. N Engl J Med. 1999, 341: 1284-1291. 10.1056/NEJM199910213411707.

    Article  CAS  Google Scholar 

  2. Criteria for diagnosis of Behcet's disease. International Study Group for Behcet's Disease. Lancet. 1990, 335: 1078-1080.

  3. Yang P, Fang W, Meng Q, Ren Y, Xing L, Kijlstra A: Clinical features of Chinese patients with Behcet's disease. Ophthalmology. 2008, 115: 312-318 e314. 10.1016/j.ophtha.2007.04.056.

    Article  Google Scholar 

  4. Cohen R, Metzger S, Nahir M, Chajek-Shaul T: Association of the MIC-A gene and HLA-B51 with Behcet's disease in Arabs and non-Ashkenazi Jews in Israel. Ann Rheum Dis. 2002, 61: 157-160. 10.1136/ard.61.2.157.

    Article  PubMed Central  CAS  Google Scholar 

  5. Sano K, Yabuki K, Imagawa Y, Shiina T, Mizuki N, Ohno S, Kulski JK, Inoko H: The absence of disease-specific polymorphisms within the HLA-B51 gene that is the susceptible locus for Behcet's disease. Tissue Antigens. 2001, 58: 77-82. 10.1034/j.1399-0039.2001.580202.x.

    Article  CAS  Google Scholar 

  6. Kera J, Mizuki N, Ota M, Katsuyama Y, Pivetti-Pezzi P, Ohno S, Inoko H: Significant associations of HLA-B*5101 and B*5108, and lack of association of class II alleles with Behcet's disease in Italian patients. Tissue Antigens. 1999, 54: 565-571. 10.1034/j.1399-0039.1999.540605.x.

    Article  CAS  Google Scholar 

  7. Mizuki N, Ota M, Yabuki K, Katsuyama Y, Ando H, Palimeris GD, Kaklamani E, Accorinti M, Pivetti-Pezzi P, Ohno S, Inoko H: Localization of the pathogenic gene of Behcet's disease by microsatellite analysis of three different populations. Invest Ophthalmol Vis Sci. 2000, 41: 3702-3708.

    CAS  Google Scholar 

  8. Yazici H, Fresko I, Yurdakul S: Behcet's syndrome: disease manifestations, management, and advances in treatment. Nat Clin Pract Rheumatol. 2007, 3: 148-155. 10.1038/ncprheum0436.

    Article  Google Scholar 

  9. Karasneh J, Gul A, Ollier WE, Silman AJ, Worthington J: Whole-genome screening for susceptibility genes in multicase families with Behcet's disease. Arthritis Rheum. 2005, 52: 1836-1842. 10.1002/art.21060.

    Article  CAS  Google Scholar 

  10. Jiang Z, Yang P, Hou S, Du L, Xie L, Zhou H, Kijlstra A: IL-23R gene confers susceptibility to Behcet's disease in a Chinese Han population. Ann Rheum Dis. 2010, 69: 1325-1328. 10.1136/ard.2009.119420.

    Article  CAS  Google Scholar 

  11. Park KS, Baek JA, Do JE, Bang D, Lee ES: CTLA4 gene polymorphisms and soluble CTLA4 protein in Behcet's disease. Tissue Antigens. 2009, 74: 222-227. 10.1111/j.1399-0039.2009.01303.x.

    Article  CAS  Google Scholar 

  12. Kim JU, Chang HK, Lee SS, Kim JW, Kim KT, Lee SW, Chung WT: Endothelial nitric oxide synthase gene polymorphisms in Behcet's disease and rheumatic diseases with vasculitis. Ann Rheum Dis. 2003, 62: 1083-1087. 10.1136/ard.62.11.1083.

    Article  PubMed Central  CAS  Google Scholar 

  13. Hou S, Yang P, Du L, Zhou H, Lin X, Liu X, Kijlstra A: SUMO4 gene polymorphisms in Chinese Han patients with Behcet's disease. Clin Immunol. 2008, 129: 170-175. 10.1016/j.clim.2008.06.006.

    Article  CAS  Google Scholar 

  14. Kamoun M, Ben Dhifallah I, Karray E, Zakraoui L, Hamzaoui K: Association of small ubiquitin-like modifier 4 (SUMO4) polymorphisms in a Tunisian population with Behcet's disease. Clin Exp Rheumatol. 2010, 28: S45-49.

    Google Scholar 

  15. Kappen JH, Wallace GR, Stolk L, Rivadeneira F, Uitterlinden AG, van Daele PL, Laman JD, Kuijpers RW, Baarsma GS, Stanford MR, Fortune F, Madanat W, van Hagen PM, van Laar JA: Low prevalence of NOD2 snps in Behcet's disease suggests protective association in Caucasians. Rheumatology (Oxford). 2009, 48: 1375-1377. 10.1093/rheumatology/kep292.

    Article  CAS  Google Scholar 

  16. Karasneh J, Hajeer AH, Barrett J, Ollier WE, Thornhill M, Gul A: Association of specific interleukin 1 gene cluster polymorphisms with increased susceptibility for Behcet's disease. Rheumatology (Oxford). 2003, 42: 860-864. 10.1093/rheumatology/keg232R.

    Article  CAS  Google Scholar 

  17. Nakao K, Isashiki Y, Sonoda S, Uchino E, Shimonagano Y, Sakamoto T: Nitric oxide synthase and superoxide dismutase gene polymorphisms in Behcet disease. Arch Ophthalmol. 2007, 125: 246-251. 10.1001/archopht.125.2.246.

    Article  CAS  Google Scholar 

  18. Wallace G: Novel genetic analysis in Behcet's disease. Arthritis Res Ther. 2009, 11: 123-10.1186/ar2757.

    Article  PubMed Central  Google Scholar 

  19. Fei Y, Webb R, Cobb BL, Direskeneli H, Saruhan-Direskeneli G, Sawalha AH: Identification of novel genetic susceptibility loci for Behcet's disease using a genome-wide association study. Arthritis Res Ther. 2009, 11: R66-10.1186/ar2695.

    Article  PubMed Central  Google Scholar 

  20. Sawalha AH, Hughes T, Nadig A, Yilmaz V, Aksu K, Keser G, Cefle A, Yazici A, Ergen A, Alarcón-Riquelme ME, Salvarani C, Casali B, Direskeneli H, Saruhan-Direskeneli G: A putative functional variant within the UBAC2 gene is associated with increased risk of Behcet's disease. Arthritis Rheum. 2011, 63: 3607-3612. 10.1002/art.30604.

    Article  PubMed Central  CAS  Google Scholar 

  21. Japan BsDRCo: The Revised Criteria and Manual of Treatments of Behcet's Disease. 2003, Tokyo: Ministry of Health and Welfare, 11-29.

    Google Scholar 

  22. Mizuki N, Meguro A, Ota M, Ohno S, Shiota T, Kawagoe T, Ito N, Kera J, Okada E, Yatsu K, Song YW, Lee EB, Kitaichi N, Namba K, Horie Y, Takeno M, Sugita S, Mochizuki M, Bahram S, Ishigatsubo Y, Inoko H: Genome-wide association studies identify IL23R-IL12RB2 and IL10 as Behcet's disease susceptibility loci. Nat Genet. 2010, 42: 703-706. 10.1038/ng.624.

    Article  CAS  Google Scholar 

  23. Remmers EF, Cosan F, Kirino Y, Ombrello MJ, Abaci N, Satorius C, Le JM, Yang B, Korman BD, Cakiris A, Aglar O, Emrence Z, Azakli H, Ustek D, Tugal-Tutkun I, Akman-Demir G, Chen W, Amos CI, Dizon MB, Kose AA, Azizlerli G, Erer B, Brand OJ, Kaklamani VG, Kaklamanis P, Ben-Chetrit E, Stanford M, Fortune F, Ghabra M, Ollier WE, 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: 698-702. 10.1038/ng.625.

    Article  PubMed Central  CAS  Google Scholar 

  24. Deng L, Wang C, Spencer E, Yang L, Braun A, You J, Slaughter C, Pickart C, Chen ZJ: Activation of the ikappab kinase complex by TRAF6 requires a dimeric ubiquitin-conjugating enzyme complex and a unique polyubiquitin chain. Cell. 2000, 103: 351-361. 10.1016/S0092-8674(00)00126-4.

    Article  CAS  Google Scholar 

  25. Zhou H, Wertz I, O'Rourke K, Ultsch M, Seshagiri S, Eby M, Xiao W, Dixit VM: Bcl10 activates the NF-kappab pathway through ubiquitination of NEMO. Nature. 2004, 427: 167-171. 10.1038/nature02273.

    Article  CAS  Google Scholar 

  26. Todaro M, Zerilli M, Triolo G, Iovino F, Patti M, Accardo-Palumbo A, di Gaudio F, Turco MC, Petrella A, de Maria R, Stassi G: NF-kappab protects Behcet's disease T cells against CD95-induced apoptosis up-regulating antiapoptotic proteins. Arthritis Rheum. 2005, 52: 2179-2191. 10.1002/art.21145.

    Article  CAS  Google Scholar 

  27. Carpino N, Turner S, Mekala D, Takahashi Y, Zang H, Geiger TL, Doherty P, Ihle JN: Regulation of ZAP-70 activation and TCR signaling by two related proteins, Sts-1 and Sts-2. Immunity. 2004, 20: 37-46. 10.1016/S1074-7613(03)00351-0.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The samples of patients and healthy controls were partially collected in ZhongShan Ophthalmic center, Sun Yat-sen University. Thanks to Ms Hongyan Zhou and Mr Xiangkun Huang for their assistance in sample collection and technical support. Thanks to all donors enrolled in the present study.

This work was supported by the Natural Science Foundation Major International (Regional) Joint Research Project (30910103912), the National Natural Science Foundation Project (81070723), the Program for the Training of a Hundred Outstanding S&T Leaders of Chongqing Municipality, Key Project of Health Bureau of Chongqing, the Project of Medical Science and Technology of Chongqing, the Chongqing Key Laboratory of Ophthalmology (CSTC, 2008CA5003), the Key Project of Natural Science Foundation (81130019), the Research Fund for the Doctoral Program of Higher Education of China (20115503110002) and the Fund for PAR-EU Scholars Program. The funders had no role in study design, data collection and analysis, the decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Peizeng Yang.

Additional information

Competing interests

None of the authors has a proprietary or financial interest in any of the products mentioned.

Authors' contributions

HSP carried out the genotyping, Dual Luciferase Reporter Gene Assay and expression studies and drafted the manuscript. SQM and JZX carried out the genotyping study. CYY, LFZ and CFL participated in the samples collection. HSP and YPZ participated in the design of the study and performed the statistical analysis. HSP and YPZ conceived of the study, and participated in its design and coordination. AK and YPZ revised the manuscript. All authors read and approved the final manuscript.

Shengping Hou, Qinmeng Shu contributed equally to this work.

Electronic supplementary material

13075_2011_3542_MOESM1_ESM.DOC

Additional file 1: The association of UBAC2 gene with Behcet's disease and linkage disequilibrium of UBAC2 gene. Table S1 presented the association analysis of three SNPs in UBAC2 gene with Behcet's disease in Han Chinese population. Figure S1 showed linkage disequilibrium plots of UBAC2 gene based on the HapMap Phase II dataset for the Han Chinese from the Beijing population by Haploview 4.2 software. (DOC 280 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Hou, S., Shu, Q., Jiang, Z. et al. Replication study confirms the association between UBAC2 and Behçet's disease in two independent Chinese sets of patients and controls. Arthritis Res Ther 14, R70 (2012). https://doi.org/10.1186/ar3789

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/ar3789

Keywords