Abstract

Behcet’s Disease (BD) is a multisystem chronic inflammatory disease. The pathology is believed to involve both genetic susceptibility and environmental factors. Hypomethylation leading to activation of interspersed repetitive sequences (IRSs) such as LINE-1 and Alu contributes to the pathologies of autoimmune diseases and cancer. Herein, the epigenetic changes of IRSs in BD were evaluated using combined bisulfite restriction analysis-interspersed repetitive sequences (COBRA-IRS). DNA from neutrophils and peripheral blood mononuclear cells (PBMCs) of BD patients with ocular involvement that were in active or inactive states and healthy controls were used to analyze LINE-1 and Alu methylation levels. For Alu sequences, significant differences were observed in the frequency of alleles between PBMCs of patients and controls (), and between inactive patients and controls (). For neutrophils, the frequency of was significantly higher between patients and controls () and between inactive patients and controls (). The partial methylation () frequencies of Alu between inactive patients and control samples also differed (). No statistically significant differences for LINE-1 were detected. Thus, changes in the methylation level of IRS elements might contribute to the pathogenesis of BD. The role of Alu transcripts in BD should be investigated further.

1. Introduction

Behcet’s Disease (BD) is a complex systemic inflammatory disorder, generally characterized by recurrent oral aphthous ulcers, genital ulcers, and uveitis [1]. However the clinical spectrum is wide and the manifestations of the disease, such as the involvement of the nervous and gastrointestinal systems, and vasculitis in large veins and arteries vary considerably depending on gender, individual differences, and ethnicity and can lead to mortality and organ loss in severe cases.

BD shares many similarities with autoinflammatory diseases, which comprise a group of disorders caused by genetic mutations in the components of the innate immune system. Among these similarities are the nonspecific inflammatory response, which manifests itself as flares and remissions, with the main involvement of neutrophils, and clinical findings, such as fever, increased acute phase proteins, and overexpression of proinflammatory cytokines, such as IL-1 B and TNF-alpha during the attacks [25]. Gene mutations in Familial Mediterranean fever, a prototypical autoinflammatory disease, have been found frequently in BD patients and are suggested to contribute to the severity of the disease [6]. BD also shows critical differences from classical autoimmune diseases, such as male dominance in severe disease [79], lack of association with autoimmune HLA class II haplotypes, and, more importantly, absence of disease-specific high titer autoantibodies or antigen-specific T cells [10].

The etiology of the disease is unknown; however, both genetic and environmental factors have been implicated in its pathogenesis. Occasional familial occurrence [11], genetic distribution along the ancient Silk Road, and an association with HLA-B51 are some of the factors pointing toward genetic involvement [12, 13]. Environmental conditions, such as bacterial or viral infections, are thought to trigger the disease in genetically susceptible individuals [10]. To date, HLA B-51 shows the strongest association with BD but accounts for less than 20% of the risk [14], which suggests the involvement of other genetic factors. Genome wide association studies revealed other candidate genes, such as IL-10, IL23R, STAT4, CCR1, and KLRC4, that could contribute to BD pathogenesis [15, 16]. One genome wide association study also suggested the epistasis between HLA-B51 and ERAP1 gene [17]. In another study, copy number variation in the DEFA1 defensin gene was associated with susceptibility to intestinal involvement in BD [18]. Other recent studies reported that more candidate gene polymorphisms involved in BD included ATG5, FAS, pre-miR-196a2, miR-182, and miR-146a [1923]. In 2014, a genome wide methylation array study in monocytes and CD4+ T lymphocytes revealed the role of epigenetics in BD pathogenesis. The authors identified abundant aberrant methylation patterns of cytoskeletal element genes in monocytes and CD4+ T lymphocytes as a major contributor to disease pathogenesis [24]. Importantly, it was reported that, after treatment, when the patients were in remission, their methylation patterns reversed back to the patterns seen in healthy controls, suggesting that a better understanding of epigenetic alterations might help us to find new disease markers and treatment options for BD patients with different symptoms. Active transcription factors and specific proteins that affect the binding of methyltransferases in BD likely determine the specific genes that are hypomethylated. However, it is clear that the regulation of methylation is defective in BD.

Almost half of the human genome comprises interspersed repetitive sequences (IRSs) [25], which can be divided into DNA transposons and retroelements. There are two types of retroelements, non-long terminal repeats (LTRs) and LTRs; non-LTR retroelements are further classified into long interspersed nuclear element (LINE) retroposons, which are represented by LINE-1 (20.1%), and short interspersed nuclear elements (SINEs), which are primarily represented by Alu sequences (13.1%) [26]. Most IRS sequences are densely methylated in normal somatic cells. In this state, they are mostly inactive and remain silent. However, hypomethylation and reactivation of these sequences are suggested to have several functional roles, such as controlling the activity of genes by regulating enhancers and repressors, or acting as an alternative promoter upon mobilization, which could lead to insertional mutations and chromosomal instability.

There is growing evidence for the contributions of epigenetic alterations to the pathogenesis of cancer and autoimmune diseases [27]. Recently, hypomethylation and reactivation of IRSs were shown to be important in autoimmune disorders such as systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), and psoriasis. In SLE, LINE-1, but not Alu, was shown to be significantly hypomethylated in CD4+ T lymphocytes and neutrophils of active patients compared with inactive patients and healthy controls [28, 29]. Similarly, significant hypomethylation of LINE-1, but not Alu, in the epidermis of psoriasis patients was reported [30]. In RA, abundant LINE-1 transcripts were detected in synovial fluids from patients [31], as well as hypomethylation in LINE-1 promoters in synovial fibroblasts [32].

In this report, we investigated for the first time the methylation in LINE-1 and Alu repetitive sequences from active and inactive Behcet’s patients compared with healthy controls.

2. Materials and Methods

2.1. Patients and Controls

Patients who fulfilled the International Study Group criteria for Behcet’s Disease were recruited from the ophthalmology clinic of the Ankara Education and Research Hospital. BD patients with an ongoing severe ocular involvement (uveitis, retinitis, papillitis, or vasculitis) were defined as the active group (; eight male, four female). BD patients that were free of any active ocular involvement for at least 3 months (confirmed by microscopic and imaging findings) and any extraocular involvement for at least 4 weeks were considered as the inactive group (; 12 men, five women). Fifteen ethnically matched healthy volunteers who were free of any acute or chronic immune-mediated illness (chronic infection, allergic and autoimmune diseases) at the time of sample collection were used as the control group (; eight women, seven men). None of the healthy controls has family history of BD and any immune-mediated diseases. The average and age range of the patients with active BD and inactive BD and healthy controls were 32.9 ± 9.0 (21–57), 36.8 ± 9.7 (18–49), and 33.8 ± 7.9 (20–48) years, respectively (mean ± SD). The details of the patient group are shown in Table 1. This study was approved by the Bogazici University Institutional Review Board for Research with Human Subjects. Informed consent was obtained from all patients and control subjects before entering the study.

2.2. Neutrophil and PBMC Isolation

Blood samples were collected into heparinized-Vacutainer tubes (Becton, Dickinson and Company, Plymouth, UK) and cell isolation was performed using density gradient centrifugation. In brief, the collected blood was carefully layered on top of lymphoprep (http://AXIS-SHIELD-PoC.com/, Oslo, Norway) at a blood : lymphoprep ratio of 2 : 1. The peripheral blood mononuclear cell (PBMC) layer was collected and washed once with Roswell Park Memorial Institute (RPMI) medium and stored at −80°C for DNA isolation. Remaining polymorphonuclear granulocytes (PMNs) or neutrophils were diluted 1 : 1 with RPMI and layered carefully on top of polymorphprep (http://AXIS-SHIELD-PoC.com/, Oslo, Norway) at a 1 : 1 ratio. The PMN fraction obtained from the gradient centrifugation was collected and washed once with RPMI. The red blood cells were removed from the PMN fraction using hypotonic lysis buffer. PMNs were collected and centrifuged and the pellet was stored at −80°C for further experiments. The purity of the PMNs was confirmed by Giemsa-Wright staining.

2.3. Combined Bisulfite Restriction Analyses (COBRA)

DNA was extracted from collected cells using a Qiagen EZ DNA isolation kit (Qiagen, Valencia, CA, USA), according to the manufacturer’s protocol. An EZ DNA methylation kit (Zymo Research, Orange, CA, USA) was used for the bisulfite conversion of 500 ng of the DNA, following the protocol supplied by the manufacturer. Global assessment of methylation levels in LINE-1 and Alu sequences was performed using the combined bisulfite restriction analysis-interspersed repetitive sequences (COBRA-IRS) technique, as previously described [2830] This technique enabled us to detect methylation levels in two CpG sites in LINE-1 and Alu repetitive elements among thousands of copies present in the human genome. The following primers were used for amplification of the bisulfite treated DNA: Alu forward 5′ GGY GYG GTG GTT TAY GTT TGT AA 3′ and Alu reverse, 5′ CTAACTTTTTATATTTTTAATAAAA ACRAAATTTCACCA 3′; LINE-1 forward 5′ GTTAAAGAAAGGGGTGAYGGT 3′ and LINE-1 reverse 5′ AATACRCCRTTTC TTAAA CCRATCTA-3′. The 92 bp PCR product for LINE-1 was obtained using the following cycle conditions: initial denaturation at 95°C for 5 min; 35 cycles of 95°C for 30 s, 55°C for 30 s, 72°C for 45 s; and a final extension at 72°C for 7 min. The 133 bp Alu PCR product was obtained by using the following cycle conditions: initial denaturation at 95°C for 5 min; 40 cycles of 95°C for 30 s, 57°C for 30 s, and 72°C for 45 s; and a final extension at 72°C for 7 min. After amplification, the LINE-1 PCR products were digested with two units of TaqI (Fermentas International, Canada) and eight units of TasI (Fermentas International) restriction enzymes. The Alu PCR products were digested with four units of TaqI. Each digestion reaction was incubated overnight at 65°C and separation was performed the following day on 8% Tris-borate-EDTA (TBE) polyacrylamide gels. The gel was stained with SYBR green nucleic acid stain (Invitrogen, Carlsbad, CA, USA). The intensity of the staining of the DNA fragments was measured using a phosphorimager equipped with ImageQuant software (Molecular Dynamics, GE Healthcare, Slough, UK). The samples were normalized against the data obtained from DNA samples of HeLa, Jurkat, and Daudi cell lines to eliminate interassay variation between each experiment.

2.4. LINE-1 Methylation Analysis

The intensity of the fragments generated by COBRA from LINE-1 depended on the methylation state of the two CpG dinucleotides and were represented by the following unique bands: fully methylated () 50 bp; partial methylation () 18 bp and () 92 bp; and hypomethylation () 60 bp (Figure 1(a)). To calculate the percentage methylation for each of the four possible states, initially, the intensity of each band was divided by the length (bp) of the double-stranded DNA: = intensity of the 92 bp fragment/92, = intensity of the 60 bp fragment/56, = intensity of the 50 fragment/48, = intensity of the 42 bp fragment/40, = intensity of the 32 bp fragment/28, and . Note that the length of the double-stranded DNA was calculated due to the overhanging end from restriction enzymes. In addition, please also note that since we could not see 18 bp band clearly, we calculate intensity of 18 bp using instead. The values obtained were then applied to the following formula; , = , , and , % Methylation () = (/()) = .

2.5. Alu Element Methylation Analysis

The COBRA method used in this study is based on the detection of methylation status of two CpG dinucleotides in the 133 bp Alu amplicon by the restriction enzyme TaqI, as shown in Figure 1(b). The different methylation states of these two CpG sites were analyzed in four groups. Fully methylated loci (), represented by a 32 bp fragment; unmethylated loci () represented by a 133 bp fragment; and partially methylated loci, (), represented by a 90 bp fragment, and (), represented by a 75 bp fragment (Figure 1(b)).

Enzymatic digestion of 133 bp COBRA-Alu products with the TaqI enzyme generated bands of 133, 90, 75, 58, 43, and 32 bp in length, with different intensities based on the methylation status of the two CpG dinucleotides. The frequency of each pattern was calculated according to the following formulas. First, the intensity of each band was divided by the length (bp) of the double-stranded DNA: = intensity of 133 bp fragment/131, = intensity of 58 bp fragment/56, = intensity of 75 bp fragment/73, = intensity of 90 bp fragment/88, = intensity of 43 bp fragment/41, and = intensity of 32 bp fragment/30. Next, the percentage of hypermethylated loci () = , the percentage of partially methylated loci () = , () = , and the percentage of hypomethylated loci () = , the percentage of methylated loci, and % Methylation () = (/( + )) = were calculated.

2.6. Statistical Analyses

Independent sample -tests (two-tailed) were applied to compare various methylation patterns of LINE-1 and Alu among active, inactive, and healthy controls. In some groups where the data were not in normal distribution, we used Mann-Whitney test instead. A value of <0.05 was considered statistically significant. Calculations were performed using GraphPad v5.0 (San Diego, CA) and SPSS software version 15.0 (SPSS Inc., Chicago, IL).

3. Results

3.1. LINE-1 Methylation Analysis

The differences in the frequency of methylation in genomic DNA obtained from PBMCs and neutrophils of patients with active and inactive BD as well as the controls were compared. The averages of the frequency of methylation are shown in Table 2. There were no statistically significant differences in the overall methylation between BD patients and healthy controls or between the active and inactive groups in both PBMC and neutrophil subsets. When we compared the differences between each methylation pattern (i.e., , , , and ), we also found no significant differences among the groups in the cell types.

3.2. Alu Methylation Analysis

The frequency of Alu methylation and paired comparisons in the PBMC and neutrophil cell types are shown in Table 3. There were no statistically significant differences in overall methylation levels of PBMCs of active patients when compared with controls and inactive patients, as well as between controls and inactive patients. However, the frequency of in PBMCs was significantly higher in BD patients compared with controls (), particularly for the inactive patients group compared with the controls (). In the neutrophil population, no statistically significant differences were observed between the groups in terms of the overall methylation. However, the level of was higher, indicating a hypomethylated status of Alu sequences in neutrophils of BD patients compared with those of healthy controls (). This finding was more striking in the inactive group compared with controls (). When we analyzed the pattern, a significant difference was also found between inactive patients and controls ().

4. Discussion

Hypomethylation of IRSs has been investigated extensively and shown to have a role in pathogenesis of complex diseases such as cancer and autoimmune related conditions [27, 33, 34]. The mechanism of how IRS hypomethylation contributes to the pathogenesis of autoimmune diseases is unclear; however, it has been suggested that it could cause aberrant regulation of neighboring genes that have a role in the regulation of the immune response [35, 36], as well as causing nonspecific autoimmune activation via the presence of viral dsRNA or ssRNA transcripts in the circulation [3739]. In this study, we investigated the role of global hypomethylation in BD and found that levels of overall methylation were not statistically different for Alu and LINE-1 among the active, inactive, and healthy control groups. However, we observed statistically significant differences in the frequency of and partial methylation of Alu in BD compared with healthy controls in both PBMCs and neutrophils. This hypomethylation status was more significant in inactive patients compared with healthy controls. Interestingly, this finding is quite different from our previous finding in SLE [28, 29].

The mechanisms driving demethylation of IRSs are complex and not fully understood; however, age-dependent and age-independent mechanisms accounting for the hypomethylation of IRSs have been reported. Both age-related (possibly caused by malnutrition and/or oxidative stress) and unrelated loss of methylation events have been proposed to occur in Alu elements, whereas hypomethylation of LINE-1 sequences was reported to be age-independent [40, 41]. LINE-1, but not Alu, hypomethylation was reported in SLE and several other autoimmune diseases. Contrastingly, for many different types of cancer, both Alu and LINE-1 hypomethylation play a role in carcinogenesis [33, 34]. In particular, LINE-1 methylation has been suggested as a poor prognostic factor for several cancer types [42]. In this study, we observed a significant increase in the frequency of the allele in both PBMCs and neutrophils in BD patients compared with controls. In addition, the frequency of partial methylation was higher in the controls compared with inactive patients in DNA obtained from neutrophils.

Although earlier reports proposed that Alu hypomethylation is age-related [40], the decrease was only significant for the 34–68-year age interval. The average ages of our control and patient cohort were comparable; therefore, the changes detected are unlikely to be related to the aging process. These findings indicated that there might be an increased tendency for hypomethylation of Alu sequences in BD patients. At present, it is unclear why there should be more change between the inactive group and the controls; however, previously, hypomethylation was found in inactive cells that represented a primed stage, suggesting that the cells were more ready to respond [29].

Interestingly, a recent report observed that the accumulation of Alu transcripts in retinal-pigmented epithelial cells induced NLRP3 activation through mitochondrial oxygen species generation, resulting in production of the proinflammatory cytokine IL-18, which contributed to the pathogenesis of age-related macular degeneration (AMD) [43]. Recently, another study showed that iron accumulation is responsible for Alu transcript accumulation and contributed to the inflammatory phenotype in AMD [44]. Thus, Alu transcript overexpression might induce inflammatory phenotypes in other diseases with similar etiologies. The BD patients who participated in this study were admitted to the eye clinic with severe eye complications, such as uveitis, in addition to oral and urogenital ulcers. An earlier report suggested that at low concentrations Alu transcripts were activators of dsRNA-dependent protein kinase R, which has an important role in the antiviral response, such as the activation of NF-KB, one of the key molecules in the initiation of the inflammatory response. It is not clear whether the slight hypomethylation of Alu seen in BD patients has any effect on the pathogenesis of BD. It should be kept in mind that BD is complex and manifests with different complications. Studies on BD have shown that different inflammatory cytokine profiles were present in different BD subtypes, such as dominant eye involvement, and vascular or neurological involvement [45]. It is possible that Alu-dependent inflammatory pathways play a role in a certain subset of BD profiles, but not in others. It should be noted that the sample size in this study was rather limited and the difference was marginal. Further analysis in larger patient groups, as well as studies on different patient profiles, will be useful to delineate whether Alu hypomethylation contributes to BD pathogenesis.

5. Conclusions

In conclusion, we demonstrated that there is an increase in the frequency of unmethylated () Alu alleles in PBMCs and neutrophils of inactive BD patients, while there was no significant difference in terms of hypomethylation between active BD samples and controls. This is contrary to the findings for autoimmune diseases such as SLE, RA, and psoriasis, in which LINE-1, but not Alu, is significantly hypomethylated in patients compared with healthy controls. This result might be explained in part by the different pathogenesis between BD, which is more similar to the autoinflammatory disease, and other classic autoimmune diseases. Further study of the epigenetic alterations, and the role and regulation of Alu transcripts in BD and its clinical significance, should be pursued to gain a better understanding of the disease. Such research could lead to the discovery of uncharacterized mechanisms that could be useful as diagnostic biomarkers and therapeutic targets.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to acknowledge the Grant for International Research Integration, Chula Research Scholar, Ratchadaphiseksomphot Endowment Fund, and Center of Excellence in Immunology and Immune-mediated Diseases. Şahru Yüksel was supported by the Rachadapisaek Sompote Postdoctoral Fund, Chulalongkorn University. They would like to thank Associate Professor Mayda Gursel and her team from Middle East Technical University for providing the infrastructure and technical support during the initial processing of the cells.