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

© 2015 Xavier et al.; licensee BioMed Central. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.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.This work was supported by the Portuguese Fundação para a Ciência e a Tecnologia (grants PTDC/SAU-GMG/098937/2008, PTDC/IIM-GES/5015/2012 and CMUP-ERI/TPE/0028/2013, fellowships SFRH/BD/43895/2008 to JMX, SFRH/BPD/35737/2007 to PA, SFRH/BPD/70008/2010 to IS, a Ciência and an Investigator-FCT contract to SAO), and the Research Committee of the Tehran University of Medical Sciences (grant 132/714).info:eu-repo/semantics/publishedVersio

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

© Springer-Verlag Berlin Heidelberg 2013Behç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.This research was supported by the Research Committee of the Tehran University of Medical Sciences (grant 132/714), the Portuguese Fundação para a Ciência e a Tecnologia (grant PTDC/SAU-GMG/098937/2008, doctoral fellowship SFRH/BD/43895/2008 to JMX, and a Ciência contract to SAO), and the Portuguese Instituto do Emprego e Formação Profissional (fellowship to JMX, TK, BVF).info:eu-repo/semantics/publishedVersio

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

© 2014 Asia Pacific League of Associations for Rheumatology and Wiley Publishing Asia Pty Ltd.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.This work was supported by the Portuguese Fundação para a Ciência e a Tecnologia (grant PTDC/SAU-GMG/098937/2008, fellowships SFRH/BD/43895/2008 to JMX, SFRH/BPD/35737/2007 to PA, SFRH/BPD/70008/2010 to IS, a Ciência and an Investigator-FCT contract to SAO) and the Research Committee of the Tehran University of Medical Sciences (grant 132/714).info:eu-repo/semantics/publishedVersio

Dense genotyping of immune-related loci implicates host responses to microbial exposure in Behcet's disease susceptibility

We analyzed 1,900 Turkish Behcet'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-IL18, IRF8, and CEBPB-PTPN1, with genomewide significance (P < 5 x 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 R1PK2 and LACC1. The disease associated allele of rs4402765, the lead marker at 11.1A-ILIB, was associated with both decreased IL-1 alpha and increased IL-1 beta production. ABO non-secretor genotypes for two ancestry specific FUT2 SNPs showed strong disease association (P = 5.89 x 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 Behcet's disease susceptibility

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 single nucleotide polymorphism (SNP) was rs1050502, a tag SNP for HLA-B*51. In the Turkish discovery set, we identified three novel loci, IL1A-IL1B, IRF8, and CEBPB-PTPN1, with genome-wide significance (P<5×10(−8)) by direct genotyping, and ADO-EGR2 by imputation. ADO-EGR2, IRF8, and CEBPB-PTPN1 replicated by genotyping 969 Iranian cases and 826 controls. Imputed data in 608 Japanese cases and 737 controls replicated ADO-EGR2 and IRF8 and meta-analysis additionally identified RIPK2 and LACC1. The disease-associated allele of rs4402765, the lead marker of the IL1A-IL1B locus, was associated with both decreased interleukin-1α and increased interleukin-1β production. ABO non-secretor genotypes of two ancestry-specific FUT2 SNPs showed strong disease association (P=5.89×10(−15)). Our findings extend shared susceptibility genes with Crohn’s disease and leprosy, and implicate mucosal factors and the innate immune response to microbial exposure in Behçet’s disease susceptibility