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

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.

HLA-B*51 allelic and carrier frequency in Kazakhstan: insights into Behçet's disease prevalence in Central Asia

Behçet's disease (BD) is a systemic vasculitis whose prevalence varies according to the ethnicity. HLA-B*51 is a well-known genetic factor predisposing to BD: populations with higher HLA-B*51 allelic frequencies are also those with a greater disease prevalence, as occurs in the Silk Road countries. However, the BD prevalence is unknown in many countries, including Kazakhstan and, in general, Central Asia, despite the location on or near the Silk Road. Due to the lack of clinical studies from this area, we analyzed the bone marrow donors registry in Kazakhstan, in order to assess the HLA-B*51 allelic and carrier frequencies, which resulted to be 8.14 and 15.28%, respectively. Through a comparison with the same parameters in other countries (by performing a data search in the medical literature and in the Allele Frequency Net Database), we might speculate a BD prevalence in Kazakhstan of 10-15/100.000 at least. This result should prompt a greater clinical consideration of this disease in Central Asia, as well as the publication of clinical data from this geographic area.

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.

A comprehensive overview on the genetics of Behçet's disease

Behçet's disease (BD) is a systemic and inflammatory disease, characterized mainly by recurrent oral and genital ulcers, eye involvement, and skin lesions. Although the exact etiopathogenesis of BD remains unrevealed, a bulk of studies have implicated the genetic contributing factors as critical players in disease predisposition. In countries along the Silk Road, human leukocyte antigen (HLA)-B51 has been reported as the strongest genetically associated factor for BD. Genome-wide association studies, local genetic polymorphism studies, and meta-analysis of combined data from Turkish, Iranian, and Japanese populations have also identified new genetic associations with BD. Among these, other HLA alleles such as HLA-B*15, HLA-B*27, HLA-B*57, and HLA-A*26 have been found as independent risk factors for BD, whereas HLA-B*49 and HLA-A*03 are independent protective alleles for BD. Moreover, other genes have also reached the genome-wide significance level of association with BD susceptibility, including IL10, IL23R-IL12RB2, IL12A, CCR1-CCR3, STAT4, TNFAIP3, ERAP1, KLRC4, and FUT2. Also, several rare nonsynonymous variants in TLR4, IL23R, NOD2, and MEFV genes have been reported to be involved in BD pathogenesis. According to genetic determinants in the loci outside the MHC region that are contributed to the host defense, immunity, and inflammation pathways, it is suggested that immune responses to the pathogen as an important environmental factor and mucosal immunity contribute to BD susceptibility.

Impact of systemic factors in shaping the periodontal microbiome

Since 2010, next-generation sequencing platforms have laid the foundation to an exciting phase of discovery in oral microbiology as it relates to oral and systemic health and disease. Next-generation sequencing has allowed large-scale oral microbial surveys, based on informative marker genes, such as 16S ribosomal RNA, community gene inventories (metagenomics), and functional analyses (metatranscriptomics), to be undertaken. More specifically, the availability of next-generation sequencing has also paved the way for studying, in greater depth and breadth, the effect of systemic factors on the periodontal microbiome. It was natural to investigate systemic diseases, such as diabetes, in such studies, along with systemic conditions or states, , pregnancy, menopause, stress, rheumatoid arthritis, and systemic lupus erythematosus. In addition, in recent years, the relevance of systemic "variables" (ie, factors that are not necessarily diseases or conditions, but may modulate the periodontal microbiome) has been explored in detail. These include ethnicity and genetics. In the present manuscript, we describe and elaborate on the new and confirmatory findings unveiled by next-generation sequencing as it pertains to systemic factors that may shape the periodontal microbiome. We also explore the systemic and mechanistic basis for such modulation and highlight the importance of those relationships in the management and treatment of patients.

Genetic Polymorphisms in Transforming Growth Factor-β, Interferon-γ and Interleukin-6 Genes and Susceptibility to Behcet's Disease in Saudi Population

Background

Behcet’s disease (BD) is a complex, inflammatory, immune-mediated multi-systemic disease of unknown etiology. Cytokines play major roles in the pathophysiology of BD, and its production may be affected by polymorphism in cytokine genes. Hence, the present study was planned to investigate any possible association between the polymorphism in TGF-β, IFN-γ, and IL-6 genes and BD in the Saudi population.

Materials and Methods

The present study includes 79 BD cases and 117 age-matched controls. Polymerase chain reaction with restriction fragment length polymorphism (PCR-RFLP) and amplification refractory mutation systems (ARMS) PCR methods were used for polymorphic analysis.

Results

The heterozygous (CT) and (CT+TT) combined genotypes of TGF-β (509-C/T), heterozygous (AT), variant (TT), and (AT+TT) combined genotypes of IFN-γ (874-A/T) were significantly (P<0.05) associated with BD in the Saudi population. No significant differences were observed for IL-6 (174-G/C) genotypes and alleles between BD cases and controls. Gender does not reflect any significant genotypic and allelic association with males and females.

Conclusion

CT genotype of TGF-β, and AT and TT genotypes of IFN-γ could be associative genetic risk factors for BD in Saudis. Regulatory region polymorphism in cytokines gene can increase inflammation and deregulated immune response, which could be risk factor for BD.

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.

Distribution of rs17482078 and rs27044 ERAP1 polymorphisms in a group of Italian Behçet's syndrome patients: a preliminary case-control study

The Endoplasmic reticulum aminopeptidase protein 1 (ERAP1) trims N-terminal amino acids from epitope precursors for Major Histocompatibility Complex class I presentation. Genome-wide association studies demonstrated that ERAP1 gene single nucleotide polymorphisms (SNPs) are associated with Behçet's syndrome (BS). This study was conducted on the two most consistently BS-associated ERAP1 polymorphisms, rs17482078 (NG_027839.1:g.35983G>A) and rs27044 (NG_027839.1:g.35997C>G) to analyse their distribution in 55 Italian BS patients and 65 ethnically matched controls (healthy controls, HC) and to test their association with BS risk. SNPs were detected by isolation, amplification of genomic DNA and direct sequencing. SNPs functional effects were predicted by bioinformatics software. The odds ratio (OR) with 95% confidence intervals was calculated to assess the strength of BS association for genotypes and alleles, also validated by logistic regression (LR). LR was used to test the association between both SNPs and patients HLA genetic data. Bonferroni correction was also applied. Comparing patients and controls, we found a significant higher frequency of rs17482078 A allele (32.73% BS vs 17.69% HC, p = 0.007) and AA genotype (18.18% BS vs 0% HC; p = 0.0003) and rs27044 G allele (63.64% BS vs 46.92% HC; p = 0.0096) in BS group after Bonferroni correction. No association was found between HLA-B*51 and both ERAP1 SNPs. Although preliminary, our data show a stronger association of rs17482078 with BS compared to rs27044 by means of case-control genetic analysis and bioinformatics prediction of protein structure change. A larger series of patients and controls is required to confirm our preliminary findings.

New Insights into Behçet's Syndrome Metabolic Reprogramming: Citrate Pathway Dysregulation

To date, a major research effort on Behçet's syndrome (BS) has been concentrated on immunological aspects. Little is known about the metabolic reprogramming in BS. Citrate is an intermediary metabolite synthesized in mitochondria, and when transported into the cytosol by the mitochondrial citrate carrier—SLC25A1-encoded protein—it is cleaved into acetyl-CoA and oxaloacetate by ATP citrate lyase (ACLY). In induced macrophages, mitochondrial citrate is necessary for the production of inflammatory mediators. The aim of our study was to evaluate SLC25A1 and ACLY expression levels in BS patients. Following a power analysis undertaken on few random samples, the number of enrolled patients was set. Thirty-nine consecutive BS patients fulfilling ISG criteria, and 21 healthy controls suitable for age and sex were recruited. BS patients were divided into two groups according to the presence (active) or absence (inactive) of clinical manifestations. Real-time PCR experiments were performed on PBMCs to quantify SLC25A1 and ACLY mRNA levels. Data processing through the Kruskal-Wallis test and Dunn's multiple comparison test as post hoc showed higher SLC25A1 and ACLY mRNA levels in BS patients compared to those in healthy controls. Therefore, SLC25A1 and ACLY upregulation suggests that metabolic reprogramming in BS involves the citrate pathway dysregulation.

HLA class I variation in Iranian Lur and Kurd populations: high haplotype and allotype diversity with an abundance of KIR ligands

HLA-A, -B and -C alleles of 285 individuals, representing three Iranian Lur populations and one Iranian Kurd population were sequenced completely, yielding human leukocyte antigen (HLA) class I genotypes at high resolution and filling four fields of the official HLA nomenclature. Each population has 87-99 alleles, evenly distributed between the three HLA class I genes, 145 alleles being identified in total. These alleles were already known, named and deposited in the HLA database. The alleles form 316 different HLA A-B-C haplotypes, with each population having between 80 and 112 haplotypes. The four Iranian populations form a related group that is distinguished from other populations, including other Iranians. All four KIR ligands - the A3/11, Bw4, C1 and C2 epitopes - are well represented, particularly Bw4, which is carried by three high-frequency allotypes: HLA-A*24:02, HLA-A*32:01 and HLA-B*51:01. In the Lur and Kurd populations, between 82% and 94% of individuals have the Bw4 epitope, the ligand for KIR3DL1. HLA-B*51:01 is likely of Neandertal origin and associated with Behcet's disease, also known as the Silk Road disease. The Lordegan Lur have the highest frequency of HLA-B*51:01 in the world. This allele is present on 46 Lur and Kurd haplotypes. Present at lower frequency is HLA-B*51:08, which is also associated with Behcet's disease. In the four Iranian populations, 31 haplotypes encode both Bw4(+) HLA-A and Bw4(+) HLA-B, a dual combination of Bw4 epitopes that is relatively rare in other populations, worldwide. This study both demonstrates and emphasizes the value of studying HLA class I polymorphism at highest resolution in anthropologically well-defined populations.

Genetic Association of HLA-A*26, -A*31, and -B*51 with Behcet's Disease in Saudi Patients

BACKGROUND

HLA-B*51 has been universally associated with Behcet’s disease (BD) susceptibility, while different alleles of HLA-A have also been identified as independent BD susceptibility loci in various ethnic populations. The objective of this study was to investigate associations of HLA-A and -B alleles with BD in Saudi patients.

MATERIALS AND METHODS

Genotyping for HLA-A and HLA-B was performed using HLA genotyping kit (Lab type^(R) SSO) in 120 Saudi subjects, including 60 BD patients and 60 matched healthy controls.

RESULTS

Our results revealed that frequencies of HLA-A*26, -A*31, and -B*51 were significantly higher in BD patients than in controls, suggesting that HLA-A*26, -A*31, and -B*51 are associated with BD. The frequency of HLA-B*15 was significantly lower in BD patients than in controls. Stratification of genotyping results into active and nonactive forms of BD revealed that the frequency of HLA-A*31 was significantly higher in the nonactive form than in the active form of BD, while there was no significant difference in the distribution of other alleles between the two forms of BD.

CONCLUSION

This study suggests that HLA-A*26, -A*31, and -B*51 are associated with susceptibility risk to BD, while HLA-B*15 may be protective in Saudi patients. However, larger scale studies are needed to confirm these findings.

Sequencing of 16S rRNA reveals a distinct salivary microbiome signature in Behçet's disease

Behçet's disease (BD) is characterized by recurrent oro-genital ulcers, mucocutaneous lesions, and serious organ involvement. We investigated the salivary microbiome in BD using high-throughput sequencing of the 16S rRNA V4 region. Stimulated saliva samples were collected from 31 BD patients and 15 healthy controls, and in 9 BD patients, a second saliva sample was collected following dental and periodontal treatment. Sequence analysis identified a total of 908 operational taxonomic units (OTUs) present across all samples. Patients had a microbial community structure that is significantly less diverse than healthy controls. The most overabundant species in BD was Haemophilus parainfluenzae, while the most depleted included Alloprevotella rava and species in the genus Leptotrichia. Periodontal treatment improved oral health indices in BD but had no short-term effect on bacterial community structure. Neither the BD-associated genetic risk locus within the HLA-B/MICA region nor being on immunosuppressive medications explained the differences between patients and controls.

Brief report: association of CCR1, KLRC4, IL12A-AS1, STAT4, and ERAP1 With Behçet's disease in Iranians

OBJECTIVE

To independently replicate the top findings from 4 published genome-wide association studies (GWAS) of susceptibility genes in Behçet's disease (BD).

METHODS

We tested 14 single-nucleotide polymorphisms (SNPs) in 13 genomic loci (excluding the major histocompatibility complex [MHC], IL10, and IL23R-IL12RB2, which have already been associated with BD in Iranians) for allelic and genotypic associations with BD in 973 patients and 828 controls from Iran and performed meta-analyses of the significantly associated markers.

RESULTS

Six SNPs (in decreasing order of significance, rs7616215 located 38 kb downstream of CCR1, rs2617170 [p.Asn104Ser] in KLRC4, rs17810546 in IL12A-AS1, rs7574070 in STAT4, and rs10050860 [p.Asp575Asn] and rs13154629 in ERAP1) were nominally associated with BD in both allelic association tests (5.05 × 10(-9) ≤ Pallele ≤ 7.55 × 10(-3) ) and sex-adjusted genotypic association tests (6.01 × 10(-9) ≤ adjusted P value ≤ 1.30 × 10(-2) ). For all 6 SNPs tested by meta-analysis (Pmeta ), the association with BD was strengthened, because the direction and magnitude of association were similar across populations (e.g., for rs7574070, odds ratio [OR] for A allele 1.29 [95% confidence interval (95% CI) 1.21-1.37], Pmeta  = 2.34 × 10(-16) ; for rs7616215, OR for C allele 0.70 [95% CI 0.65-0.76], Pmeta  = 1.54 × 10(-19) ; for rs17810546, OR for A allele 0.60 [95% CI 0.52-0.70], Pmeta  = 6.34 × 10(-11) ; for rs2617170, OR for T allele 0.76 [95% CI 0.70-0.81], Pmeta  = 2.75 × 10(-14) ; for rs13154629, OR for TT genotype 2.76 [95% CI 2.01-3.80], Pmeta  = 3.57 × 10(-10) ).

CONCLUSION

This study reinforces the notion that CCR1, KLRC4, IL12A-AS1, STAT4, and ERAP1 are bona fide susceptibility genes for BD, in addition to the MHC, IL10, and IL23R-IL12RB2 loci. Future genetic and functional studies are now warranted to uncover the roles of these genes in the pathogenesis of BD.

Associations of MICA Polymorphisms with Inflammatory Rheumatic Diseases

Inflammatory rheumatic diseases are characterized by inflammation resulting from the immune dysregulation that usually attacks joints, skin and internal organs. Many of them are considered as complex disease that may be predisposed by multiple genes and/or genetic loci, and triggered by environmental factors such as microbiome and cellular stress. The major histocompatibility complex class I chain-related gene A (MICA) is a highly polymorphic gene that encodes protein variants expressed under cellular stress conditions, and these MICA variants play important roles in immune activation and surveillance. Recently, accumulating evidences from both genetic and functional studies have suggested that MICA polymorphisms may be associated with various rheumatic diseases, and the expression of MICA variants may attribute to the altered immune responses in the diseases. The objective of this review is to discuss potential genetic associations and pathological relevance of MICA in inflammatory rheumatic diseases that may help us to understand pathogenesis contributing to the development of these diseases.

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.