Dense genotyping of immune-related loci implicates host responses to microbial exposure in Behçet's disease susceptibility.

We analyzed 1,900 Turkish Behçet's disease cases and 1,779 controls genotyped with the Immunochip. The most significantly associated SNP was rs1050502, a tag SNP for HLA-B*51. In the Turkish discovery set, we identified three new risk loci, IL1A-IL1B, IRF8, and CEBPB-PTPN1, with genome-wide significance (P < 5 × 10-8) by direct genotyping and ADO-EGR2 by imputation. We replicated the ADO-EGR2, IRF8, and CEBPB-PTPN1 loci by genotyping 969 Iranian cases and 826 controls. Imputed data in 608 Japanese cases and 737 controls further replicated ADO-EGR2 and IRF8, and meta-analysis additionally identified RIPK2 and LACC1. The disease-associated allele of rs4402765, the lead marker at IL1A-IL1B, was associated with both decreased IL-1α and increased IL-1ß production. ABO non-secretor genotypes for two ancestry-specific FUT2 SNPs showed strong disease association (P = 5.89 × 10-15). Our findings extend the list of susceptibility genes shared with Crohn's disease and leprosy and implicate mucosal factors and the innate immune response to microbial exposure in Behçet's disease susceptibility.

IL10 low-frequency variants in Behçet's disease patients.

AIM: To explain the missing heritability after the genome-wide association studies era, sequencing studies allow the identification of low-frequency variants with a stronger effect on disease risk. Common variants in the interleukin 10 gene (IL10) have been consistently associated with Behçet's disease (BD) and the goal of this study is to investigate the role of low-frequency IL10 variants in BD susceptibility. METHODS: To identify IL10 low-frequency variants, a discovery group of 50 Portuguese BD patients were Sanger-sequenced in a 7.7 kb genomic region encompassing the complete IL10 gene, 0.9 kb upstream and 2 kb downstream, and two conserved regions in the putative promoter. To assess if the novel variants are BD- and/or Portuguese-specific, they were assayed in an additional group of BD patients (26 Portuguese and 964 Iranian) and controls (104 Portuguese and 823 Iranian). RESULTS: Rare IL10 coding variants were not detected in BD patients, but we identified 28 known single nucleotide polymorphisms with minor allele frequencies ranging from 0.010 to 0.390, and five novel non-coding variants in five heterozygous cases. ss836185595, located in the IL10 3' untranslated region, was also detected in one Iranian control individual and therefore is not specific to BD. The remaining novel IL10 variants (ss836185596 and ss836185602 in intron 3, ss836185598 and ss836185604 in the putative promoter region) were not found in the replication dataset. CONCLUSION: This study highlights the importance of screening the whole gene and regulatory regions when searching for novel variants associated with complex diseases, and the need to develop bioinformatics tools to predict the functional impact of non-coding variants and statistical tests which incorporate these predictions.

Characterization of the major histocompatibility complex locus association with Behçet's disease in Iran.

INTRODUCTION: The aim of this study was to characterize the association of human leukocyte antigen (HLA) B alleles and major histocompatibility complex (MHC) single nucleotide polymorphisms (SNPs) with Behçet's disease (BD) in an Iranian dataset. METHODS: The association of three SNPs in the MHC region previously identified as the most associated in high-density genotyping studies was tested in a case-control study on 973 BD patients and 825 controls from Iran, and the association of HLA-B alleles was tested in a subset of 681 patients and 414 controls. RESULTS: We found that HLA-B*51 (P = 4.11 × 10(-41), OR [95% CI] = 4.63[3.66-5.85]) and B*15 confer risk for BD (P = 2.83 × 10(-2), OR [95% CI] = 1.75[1.08-2.84]) in Iranian, and in B*51 negative individuals, only the B*15 allele is significantly associated with BD (P = 2.51 × 10(-3), OR [95% CI] = 2.40[1.37-4.20]). rs76546355, formerly known as rs116799036, located between HLA-B and MICA (MHC class I polypeptide-related sequence A), demonstrated the same level of association with BD as HLA-B*51 (P adj = 1.78 × 10(-46), OR [95% CI] = 5.46[4.21-7.09], and P adj = 8.34 × 10(-48), OR [95% CI] = 5.44[4.20-7.05], respectively) in the HLA-B allelotyped subset, while rs2848713 was less associated (P adj = 7.14 × 10(-35), OR [95% CI] = 3.73[2.97-4.69]) and rs9260997 was not associated (P adj = 1.00 × 10(-1)). Additionally, we found that B*51 genotype-phenotype correlations do not survive Bonferroni correction, while carriers of the rs76546355 risk allele predominate in BD cases with genital ulcers, positive pathergy test and positive BD family history (2.31 × 10(-4) ≤ P ≤ 1.59 × 10(-3)). CONCLUSIONS: We found that the HLA-B*51 allele and the rs76546355/rs116799036 MHC SNP are independent genetic risk factors for BD in Iranian, and that positivity for the rs76546355/rs116799036 risk allele, but not for B*51, does correlate with specific demographic characteristics or clinical manifestations in BD patients.

FUT2: filling the gap between genes and environment in Behçet's disease?

OBJECTIVES: To identify new susceptibility loci for Behçet's disease (BD), we performed a genome-wide association study (GWAS) using DNA pooling. METHODS: Two replicate pools of 292 Iranian BD cases and of 294 age- and sex-matched controls were allelotyped in quadruplicate on the Affymetrix Genome-Wide Human SNP Array 6.0. Of the 51 top markers, 47 were technically validated through individually genotyping. Replication of validated single nucleotide polymorphisms (SNPs) was performed in an independent Iranian dataset (684 cases and 532 controls). RESULTS: In addition to the well-established HLA-B locus, rs7528842 in a gene desert on chromosome 1p21.2, and rs632111 at the 3'UTR of FUT2 were associated in both the discovery and replication datasets (individually and in combination). However, only the FUT2 SNP was associated in a previous GWAS for BD in Turkish people. Fine-mapping of FUT2 in the full Iranian dataset showed additional associations in five coding SNPs (2.97E-06

Gene expression profiling and association studies implicate the neuregulin signaling pathway in Behçet's disease susceptibility.

Behçet's disease (BD) is a complex disease with genetic and environmental risk factors implicated in its etiology; however, its pathophysiology is poorly understood. To decipher BD's genetic underpinnings, we combined gene expression profiling with pathway analysis and association studies. We compared the gene expression profiles in peripheral blood mononuclear cells (PBMCs) of 15 patients and 14 matched controls using Affymetrix microarrays and found that the neuregulin signaling pathway was over-represented among the differentially expressed genes. The Epiregulin (EREG), Amphiregulin (AREG), and Neuregulin-1 (NRG1) genes of this pathway stand out as they are also among the top differentially expressed genes. Twelve haplotype tagging SNPs at the EREG-AREG locus and 15 SNPs in NRG1 found associated in at least one published BD genome-wide association study were tested for association with BD in a dataset of 976 Iranian patients and 839 controls. We found a novel association with BD for the rs6845297 SNP located downstream of EREG, and replicated three associations at NRG1 (rs4489285, rs383632, and rs1462891). Multifactor dimensionality reduction analysis indicated the existence of epistatic interactions between EREG and NRG1 variants. EREG-AREG and NRG1, which are members of the epidermal growth factor (EGF) family, seem to modulate BD susceptibility through main effects and gene-gene interactions. These association findings support a role for the EGF/ErbB signaling pathway in BD pathogenesis that warrants further investigation and highlight the importance of combining genetic and genomic approaches to dissect the genetic architecture of complex diseases.

Association study of IL10 and IL23R-IL12RB2 in Iranian patients with Behçet's disease.

OBJECTIVE: Independent replication of the findings from genome-wide association studies (GWAS) remains the gold standard for results validation. Our aim was to test the association of Behçet's disease (BD) with the interleukin-10 gene (IL10) and the IL-23 receptor-IL-12 receptor ß2 (IL23R-IL12RB2) locus, each of which has been previously identified as a risk factor for BD in 2 different GWAS. METHODS: Six haplotype-tagging single-nucleotide polymorphisms (SNPs) in IL10 and 42 in IL23R-IL12RB2 were genotyped in 973 Iranian patients with BD and 637 non-BD controls. Population stratification was assessed using a panel of 86 ancestry-informative markers. RESULTS: Subtle evidence of population stratification was found in our data set. In IL10, rs1518111 was nominally associated with BD before and after adjustment for population stratification (odds ratio [OR] for T allele 1.20, 95% confidence interval [95% CI] 1.02-1.40, unadjusted P [P(unadj) ] = 2.53 × 10(-2) ; adjusted P [P(adj) ] = 1.43 × 10(-2) ), and rs1554286 demonstrated a trend toward association (P(unadj) = 6.14 × 10(-2) ; P(adj) = 3.21 × 10(-2) ). Six SNPs in IL23R-IL12RB2 were found to be associated with BD after Bonferroni correction for multiple testing, the most significant of which were rs17375018 (OR for G allele 1.51, 95% CI 1.27-1.78, P(unadj) = 1.93 × 10(-6) ), rs7517847 (OR for T allele 1.48, 95% CI 1.26-1.74, P(unadj) = 1.23 × 10(-6) ), and rs924080 (OR for T allele 1.29, 95% CI 1.20-1.39, P = 1.78 × 10(-5) ). SNPs rs10489629, rs1343151, and rs1495965 were also significantly associated with BD in all tests performed. Results of meta-analyses of our data combined with data from other populations further confirmed the role of rs1518111, rs17375018, rs7517847, and rs924080 in the risk of BD, but no epistatic interactions between IL10 and IL23R-IL12RB2 were detected. Results of imputation analysis highlighted the importance of IL23R regulatory regions in the susceptibility to BD. CONCLUSION: These findings independently confirm, extend, and refine the association of BD with IL10 and IL23R-IL12RB2. These associations warrant further validation and investigation in patients with BD, as they may have implications for the development of novel therapies (e.g., immunosuppressive therapy targeted at IL-23p19).