On the genetics of the Silk Route: association analysis of HLA, IL10, and IL23R-IL12RB2 regions with Behçet's disease in an Iranian population

Risk factors, clinical features and treatment of Behçet's disease uveitis

Behçet's disease is a systemic vasculitis frequently associated with intraocular inflammation. Recent findings identified independent clinical clusters in Behçet's disease, each involving distinct combinations of affected organs. Ocular Behçet's disease, mainly manifested as uveitis, is characterized as an independent cluster with a low likelihood of association with other system involvements, such as intestinal, cardiovascular, or central nervous system. A prevailing theory suggests that the pathogenesis of the disease is multifactorial, where a variety of genetic and infectious agents may interact with each other to cause the disease. Among sporadic cases, the human leukocyte antigen (HLA) genes, including HLA-B51, HLA-A26, HLA-B15, and HLA-B5701, have been found to be a key component conferring genetic susceptibility. Outside the HLA region, a set of susceptibility variants are identified, closely related to interleukin (IL)-23/IL-17 pathway, tumor necrosis factor (TNF) signaling, and pattern recognition receptor systems. Microbial infections, such as Streptococcus sanguinis, Mycobacterium tuberculosis, and Herpes simplex virus (HSV), are linked to play the triggering of disease in immunogenetically predisposed individuals. Clinically, due to the notable relapsing-remitting course of ocular Behçet's disease, the prevention of recurrent attack would be the primary treatment goal. Combination of corticosteroids and immunomodulatory drugs, such as anti-TNF agents, interferon, and conventional immunosuppressants (e.g. cyclosporine, azathioprine), have been the mainstream regimen for the disease. Future research may focus on comparing the effectiveness of immunomodulatory drugs and identifying the most suitable subgroups for a specific drug on the basis of the knowledge of the molecular heterogeneity of the disease.

Identification of novel genes in Behcet's disease using integrated bioinformatic analysis

Behcet’s disease (BD) is a chronic vascular inflammatory disease. However, the etiology and molecular mechanisms underlying BD development have not been thoroughly understood. Gene expression data for BD were obtained from the Gene Expression Omnibus database. We used robust rank aggregation (RRA) to identify differentially expressed genes (DEGs) between patients with BD and healthy controls. Gene ontology functional enrichment was used to investigate the potential functions of the DEGs. Protein–protein interaction (PPI) network analysis was performed to identify the hub genes. Receiver operating characteristic analyses were performed to investigate the value of hub genes in the diagnosis of BD. GSE17114 and GSE61399 datasets were included, comprising 32 patients with BD and 26 controls. The RRA integrated analysis identified 44 significant DEGs among the GSE17114 and GSE61399 CD4 + T lymphocytes. Functional enrichment analysis revealed that protein tyrosine/threonine phosphatase activity and immunoglobulin binding were enriched in BD. PPI analysis identified FCGR3B as a hub gene in the CD4 + T lymphocytes of BD patients. Our bioinformatic analysis identified new genetic features, which will enable further understanding of the pathogenesis of BD.

Major Histocompatibility Complex Class I-Related Chain A (MICA) Allelic Variants Associate With Susceptibility and Prognosis of Gastric Cancer

Gastric cancer (GC) is the fifth most prevalent type of cancer worldwide. Gastric tumor cells express MICA protein, a ligand to NKG2D receptor that triggers natural killer (NK) cells effector functions for early tumor elimination. MICA gene is highly polymorphic, thus originating alleles that encode protein variants with a controversial role in cancer. The main goal of this work was to study MICA gene polymorphisms and their relationship with the susceptibility and prognosis of GC. Fifty patients with GC and 50 healthy volunteers were included in this study. MICA alleles were identified using Sanger sequencing methods. The analysis of MICA gene sequence revealed 13 MICA sequences and 5 MICA-short tandem repeats (STR) alleles in the studied cohorts We identified MICA^*002 (^*A9) as the most frequent allele in both, patients and controls, followed by MICA^*008 allele (^*A5.1). MICA^*009/049 allele was significantly associated with increased risk of GC (OR: 5.11 [95% CI: 1.39–18.74], p = 0.014). The analysis of MICA-STR alleles revealed a higher frequency of MICA^*A5 in healthy individuals than GC patients (OR = 0.34 [95% CI: 0.12–0.98], p = 0.046). Survival analysis after gastrectomy showed that patients with MICA^*002/002 or MICA^*002/004 alleles had significantly higher survival rates than those patients bearing MICA^*002/008 (p = 0.014) or MICA^*002/009 (MICA^*002/049) alleles (p = 0.040). The presence of threonine in the position MICA-181 (MICA^*009/049 allele) was more frequent in GC patients than controls (p = 0.023). Molecular analysis of MICA-181 showed that the presence of threonine provides greater mobility to the protein than arginine in the same position (MICA^*004), which could explain, at least in part, some immune evasion mechanisms developed by the tumor. In conclusion, our findings suggest that the study of MICA alleles is crucial to search for new therapeutic approaches and may be useful for the evaluation of risk and prognosis of GC and personalized therapy.

Genetics of Behçet's Disease: Functional Genetic Analysis and Estimating Disease Heritability

Behçet's disease is a chronic multisystemic inflammatory disorder characterized by recurrent oral and genital ulcers. Although its etiology remains unclear, it is thought that both genetic and environmental factors contribute to the onset and progression of Behçet's disease. Here, we provide an updated view of the genetic landscape and architecture of Behçet's disease. Large-scale genetic studies performed to date revealed 21 genetic susceptibility loci associated with the disease at a GWAS level of significance (p-value = 5 × 10⁻⁸). We performed epigenetic pattern enrichment analysis in Behçet's disease associated loci, providing new insights into the molecular mechanisms underlying its pathophysiology. Our data suggest the crucial involvement of several immune cell types, including natural killer cells, monocytes, and B cells in the pathogenesis of the disease. Pathway enrichment analysis identified important biological processes involved. Using large-scale genetic data available from ~200 immune-related loci (Immunochip), we estimate Behçet's disease heritability to be at least 16%. We further used the same approach to estimate the heritability explained by the known Behçet's disease-associated loci, suggesting that they explain ~ 60% of the genetic component underlying Behçet's disease. These results indicate a significant role of non-genetic factors in causing Behçet's disease and that additional genetic variation influencing the risk of Behçet's disease remains to be identified. Finally, we calculated a cumulative genetic risk score across populations reinforcing the link between geographic variations in disease prevalence with its genetic component.

Association between Interleukin-10 Gene Polymorphisms and Behcet's Disease Susceptibility: Evidence from a Meta-Analysis

Epidemiological studies have demonstrated that interleukin-10 (IL-10) polymorphisms may be associated with the development of Behcet's disease (BD). However, the published results were inconsistent. Therefore, this meta-analysis was conducted to derive a more precise relationship between IL-10 polymorphisms and BD susceptibility. Online databases (PubMed, Embase, Science Citation Index (SCI), CNKI, and WanFang) were searched to identify eligible studies. Odds ratio (OR) and a 95% confidence interval (CI) were applied to assess the relationship strength between IL-10 -1082A>G (rs1800896), -819T>C (rs1800871), and -592A>C (rs1800872) polymorphisms and BD susceptibility. Publication bias, sensitivity, and cumulative analyses were conducted to measure the robustness of our findings. Finally, fifteen articles (36 independent case-control studies) involving 5,971 patients and 8,913 controls were included. Overall, significant associations between -819T>C polymorphisms and BD risk were observed in the total population (C vs. T: OR = 0.72, 95%CI = 0.67‐0.77, P < 0.01, I² = 36.6%; TC vs. TT: OR = 0.73, 95%CI = 0.66‐0.80, P < 0.01, I² = 23.0%; CC vs. TT: OR = 0.52, 95%CI = 0.39‐0.70, P < 0.01, I² = 53.7%; TC+CC vs. TT: OR = 0.67, 95%CI = 0.61‐0.71, P < 0.01, I² = 22.1%; and CC vs. TT+TC: OR = 0.66, 95%CI = 0.53‐0.82, P < 0.01, I² = 57.8%). Moreover, the IL-10 -592 A>C polymorphism and the ACC haplotype exhibited a significant, protective effect against BD susceptibility. In summary, our meta-analysis suggested that IL-10 gene polymorphisms may play a salient role for BD development.

Past and Present Behçet's Disease Animal Models

Behçet's disease (BD) is presumably an autoinflammatory disease of unknown etiology for which several animal models have been described over the years. Agents and methods used for the development of these models have ranged from the herpes simplex type one virus (hsv-1) pathogen to the use of transgenic mice. Other models have also been used to investigate a possible autoimmune component. Each model possesses its own unique set of benefits and shortcomings, with no one model fully being able to recapitulate the disease phenotype. Here, we review the proposed models and provide commentary on their effectiveness and usefulness in studying the disease.

Exploring the association of IL-10 polymorphisms in Behcet's disease: a systematic review and meta-analysis

Background

Polymorphisms in the interleukin-10 (IL-10) gene have been studied in various ethnic groups for possible association with Behçet’s disease (BD). This study aimed to perform a meta-analysis of eligible studies to calculate the association of IL-10 polymorphisms with BD.

A systematic literature search was carried out in PubMed, Embase, Web of Science, and Scopus databases to identify relevant publications, and extracted the respective results. Pooled odds ratio (OR) with 95% confidence interval (CI) was used to evaluate the power of association with a random-effects model.

Results

A total of 19 articles, consisting of 10,626 patients and 13,592 controls were included in the meta-analysis. The meta-analysis revealed significant associations in allelic and genotypic test models of − 819 (C vs. T: OR = 0.691, P < 0.001; CC vs. TT: OR = 0.466, P < 0.001; CC + CT vs. TT: OR = 0.692, P < 0.001; and CC vs. CT + TT: OR = 0.557, P < 0.001), − 592 (C vs. A: OR = 0.779, P = 0.002; CC + AA vs. AA: OR = 0.713, P = 0.021; and CA vs. AA: OR = 0.716, P = 0.016), rs1518111 (G vs. A: OR = 0.738, P < 0.001; GG vs. AA: OR = 0.570, P < 0.001; GG + AG vs. AA: OR = 0.697, P < 0.001; GG vs. GA + AA: OR = 0.701, P < 0.001; and AG vs. GG: OR = 0.786, P = 0.004) and rs1554286 (C vs. T: OR = 0.582, P < 0.001; CC vs. TT: OR = 0.508, P < 0.001; CC + CT vs. TT: OR = 0.605, P < 0.001; CC vs. CT + TT: OR = 0.665, P = 0.012; and CT vs. TT: OR = 0.646, P = 0.001). However, we failed to find any association between − 1082 polymorphism and susceptibility of BD.

Conclusion

This meta-analysis demonstrated that the interleukin-10 -819, − 596, rs1518111 and rs1554286 polymorphisms could be responsible against BD susceptibility, and should probably be regarded as a protective factor for Behçet’s disease.

Single-Nucleotide Polymorphisms in IL23R-IL12RB2 (rs1495965) Are Highly Prevalent in Patients with Behcet's Uveitis and Vary Between Populations

Purpose: To test the frequency of single-nucleotide polymorphisms in the IL-10, IL23R-IL12RB2 genes in patients with Behcet's uveitis. Methods: Blood samples were collected from 89 Israeli and Turkish patients, and from healthy control subjects of different origins. Genomic DNA was extracted from peripheral blood leukocytes and genotyped. Results: The risk allele, A, in rs1800871, of IL-10 gene was highly prevalent in Behcet's uveitis and healthy control samples alike; highest among the Turkish groups. Prevalence of G allele, in rs1495965, in the IL23R-IL12RB2 gene was high in Behcet's uveitis patients, and among healthy Turkish and Israelis of Middle Eastern origin, while lower among the other Israeli control group (77.9%, 78.9%, 27.8%, respectively, P < 0.001). Conclusion: Our findings highlight the differences between populations and may account for the increased prevalence of the disease among Turkish and Israelis of Middle Eastern origin. Further studies are required to map other healthy and affected populations.

A new MHC-linked susceptibility locus for primary Sjögren's syndrome: MICA

The association of primary Sjögren's syndrome (pSS) with Major Histocompatibility Complex (MHC) alleles is quintessential of MHC-disease associations. Indeed, although disease associations with classical HLA class I and II alleles/haplotypes are amply documented, further dissection is often prevented by the strong linkage disequilibrium across the entire MHC complex. Here we study the association of pSS, not with HLA genes, but with the non-conventional MHC encoded class I gene, MICA (MHC class I chain-related gene A). MICA is selectively expressed within epithelia, and is the major ligand for the activatory receptor, NKG2D, both attributes relevant to pSS' etiology. MICA-pSS association was studied in two independent (French and UK) cohorts representing a total of 959 cases and 1,043 controls. MICA*008 allele was shown to be significantly associated with pSS (pcor=2.61 × 10-35). A multivariate logistic regression showed that this association was independent of all major known MHC-linked risk loci/alleles, as well as other relevant candidate loci that are in linkage disequilibrium with MICA*008 i.e. HLA-B*08:01, rs3131619 (T), MICB*008, TNF308A, HLA-DRB1*03:01 and HLA-DRB1*15:01 (P = 1.84 × 10-04). Furthermore, independently of the MICA*008 allele, higher levels of soluble MICA proteins were detected in sera of pSS patients compared to healthy controls. This study hence defines MICA as a new, MHC-linked, yet HLA-independent, pSS risk locus and opens a new front in our understanding of the still enigmatic pathophysiology of this disease. The fact that the soluble MICA protein is further amplified in MICA*008 carrying individuals, might also be relevant in other auto-immune diseases and cancer.

Behçet's disease risk association fine-mapped on the IL23R-IL12RB2 intergenic region in Koreans

Background

Behçet’s disease (BD) susceptibility had been associated with single-nucleotide polymorphisms (SNPs) in IL23R–IL12RB2, IL10, STAT4, or ERAP1 locus in Japanese, Turkish, Chinese, and other populations, but not in a Korean genome-wide association study (GWAS). We aimed to fine-map BD risk association of these four loci using extensive imputation and additional genotyping for replication.

Methods

In the discovery phase, 369 patients with BD enrolled in the previous Korean GWAS and 2000 controls retrieved from a population-based cohort of healthy Koreans were imputed for their genotypes of all SNPs in the four loci using the Asian data of the 1000 Genomes Project as reference. For genotype imputation of ERAP1 SNPs, the adjacent ERAP2 SNPs were also covered. For the 10 most significantly associated SNPs (8 imputed and 2 GWAS-genotyped), an additional 84 patients with BD and 283 healthy controls were genotyped for replication. The results from the discovery and replication phases were pooled for meta-analysis using the Mantel-Haenszel test to estimate the odds ratio (OR) and 95% confidence interval (CI).

Results

An IL23R–IL12RB2 intergenic SNP rs1495965 was significantly associated with BD risk (OR (95% CI) = 1.5 (1.3, 1.7), P = 2.5 × 10⁻⁷) in the pooled meta-analysis of the discovery (1.4 (1.2, 1.7), P = 4.9 × 10⁻⁷) and replication (1.9 (1.3, 2.6), P = 6.0 × 10⁻⁴) phases. BD risk association was fine-mapped on the intergenic region rather than the two flanking genes, as rs1495966 and rs4655535, almost perfectly correlated with rs1495965 (r² = 0.99), were also located in the same intergenic region. Consistent with previous reports, the P values tended to be lower within IL23R than IL12RB2. On the other hand, several IL10 SNPs were suggested for association in the discovery phase but all failed in the replication phase. No SNP in ERAP1–ERAP2 and STAT4 was suggested even in the discovery phase.

Conclusions

BD susceptibility association was fine-mapped on the intergenic region between IL23R and IL12RB2 as marked by three correlated SNPs, rs1495965, rs1495966, and rs4655535.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-017-1435-5) contains supplementary material, which is available to authorized users.

The potential of genetically guided treatment in Behçet's disease

Continuous identification of specific targets and candidate genes together with novel approaches offers new promises for the future of gene therapy design in Behçet's disease (BD). Personalized medicine based on pharmacogenomics is being developed at the clinical stage to improve treatment response. Screening the whole gene and regulatory regions is important when searching for novel variants associated with such complex diseases. Different host genetic factors play significant roles in susceptibility to BD. Thus, identifying these genes responsible for susceptibility and resistance to BD may offer a notable contribution toward understanding its pathogenesis, and may lead to the development of novel prophylactic and treatment strategies. Evidenced-based treatment strategy is recommended for the management in BD patients. This review sheds light on the immunopathogenesis and pharmacogenetics of BD with special attention to the treatment targeting gene polymorphisms. In conclusion, the potential of genetically guided treatment in BD takes us back to the future for an accurate management strategy of this serious rheumatic disease. The ongoing discovery of pivotal genes related to the susceptibility, manifestations, disease activity and treatment options provide substantial hope to the reduced frequency of BD, effective control and improvement in the prognosis. Targeted gene therapy could be a leading option in the treatment armamentarium of BD.

Association between Functional MICA-TM and Behcet's Disease: A Systematic Review and Meta-analysis

The relationships between polymorphisms of the trans-membrane(TM) region located in the major histocompatibility complex (MHC) class I chain-related gene A (MICA) and Behcet's disease (BD) have been discussed previously, however, the results were contradictory. In this study, we thoroughly assess whether MICA-TM gene variants are associated with BD by means of a systematic review and meta-analysis. Our study focused on the effects of polymorphisms of MICA-A4, A5, A5.1, A6, and A9 from the included articles. Sixteen previous original publications representing 1,555 BD patients and 2,086 unrelated healthy controls analyzed the association of BD with MICA-TM gene polymorphisms. For the five alleles, MICA-A6 showed a strongly positive correlation with BD patients and could be viewed as an increased risk factor of BD (OR = 2.34, 95%CI: 2.02-2.70). Furthermore, MICA-A4, A5, A5.1, and A9 exhibited negative associations with BD (OR = 0.71, 95%CI: 0.58-0.86; OR = 0.75, 95%CI: 0.63-0.90; OR = 0.63, 95%CI: 0.44-0.91; OR = 0.70, 95%CI: 0.58-0.84, respectively). Our meta-analysis confirmed MICA-A6 could be responsible for BD in three ethnic regions and should probably be treated as a risk factor for BD. MICA-A4, A5, A5.1, and A9 could be regarded as protective factors, especially in the Middle East and East Asia.

Epigenomic functional characterization of genetic susceptibility variants in systemic vasculitis

Systemic vasculitides are poorly understood inflammatory diseases of the blood vessels that are frequently associated with significant organ damage. Genetic risk variants contribute to the susceptibility of vasculitis, but functional consequences of these genetic variants are largely unknown. Most genetic risk variants in immune-mediated diseases, including systemic vasculitis, are localized to non-coding genetic regions suggesting they might increase disease risk by influencing regulatory elements within the genome. Long range regulatory interactions pose an additional obstacle in localizing functional consequences associated with risk variants to specific genes or cell types. We used cell-type specific enrichment patterns of histone changes that mark poised, primed, and active enhancers, and DNase hypersensitivity to identify specific immune cells mediating genetic risk in vasculitis. Our data suggest that genetic risk variants in ANCA-associated vasculitis are significantly enriched in enhancer elements in Th17 cells, supporting a role for Th17 cells in this disease. Primed and active enhancer elements in B cells can be potentially affected by genetic risk variants associated with Kawasaki disease. Genetic risk in Behçet's disease and Takayasu arteritis might affect enhancer elements in multiple cell types, possibly explained by influencing enhancers in hematopoietic stem cells. Interestingly, our analyses indicate a role for B cells in Kawasaki disease, Behçet's disease, and Takayasu arteritis, and suggest that further work to characterize the involvement of B cells in these diseases is warranted.