February 2018
Volume 59, Issue 2
Open Access
Genetics  |   February 2018
Association of Long Noncoding RNAs Polymorphisms With Ankylosing Spondylitis, Vogt-Koyanagi-Harada Disease, and Behcet's Disease
Author Affiliations & Notes
  • Yingying Yue
    The First Affiliated Hospital of Chongqing Medical University, Chongqing Key Laboratory of Ophthalmology, Chongqing Eye Institute, Chongqing, China
  • Jun Zhang
    The First Affiliated Hospital of Chongqing Medical University, Chongqing Key Laboratory of Ophthalmology, Chongqing Eye Institute, Chongqing, China
  • Lu Yang
    The First Affiliated Hospital of Chongqing Medical University, Chongqing Key Laboratory of Ophthalmology, Chongqing Eye Institute, Chongqing, China
  • Shengyun Liu
    Department of Rheumatology, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
  • Jian Qi
    The First Affiliated Hospital of Chongqing Medical University, Chongqing Key Laboratory of Ophthalmology, Chongqing Eye Institute, Chongqing, China
  • Qingfeng Cao
    The First Affiliated Hospital of Chongqing Medical University, Chongqing Key Laboratory of Ophthalmology, Chongqing Eye Institute, Chongqing, China
  • Chunjiang Zhou
    The First Affiliated Hospital of Chongqing Medical University, Chongqing Key Laboratory of Ophthalmology, Chongqing Eye Institute, Chongqing, China
  • Yao Wang
    The First Affiliated Hospital of Chongqing Medical University, Chongqing Key Laboratory of Ophthalmology, Chongqing Eye Institute, Chongqing, China
  • Aize Kijlstra
    University Eye Clinic Maastricht, Maastricht, The Netherlands
  • Peizeng Yang
    The First Affiliated Hospital of Chongqing Medical University, Chongqing Key Laboratory of Ophthalmology, Chongqing Eye Institute, Chongqing, China
    The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
  • Shengping Hou
    The First Affiliated Hospital of Chongqing Medical University, Chongqing Key Laboratory of Ophthalmology, Chongqing Eye Institute, Chongqing, China
  • Correspondence: Peizeng Yang, The First Affiliated Hospital of Chongqing Medical University, Youyi Road 1, Chongqing 400016, China; peizengycmu@126.com
  • Shengping Hou, The First Affiliated Hospital of Chongqing Medical University, Youyi Road 1, Chongqing 400016, China; sphou828@163.com
Investigative Ophthalmology & Visual Science February 2018, Vol.59, 1158-1166. doi:https://doi.org/10.1167/iovs.17-23247
  • Views
  • PDF
  • Share
  • Tools
    • Alerts
      ×
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      Yingying Yue, Jun Zhang, Lu Yang, Shengyun Liu, Jian Qi, Qingfeng Cao, Chunjiang Zhou, Yao Wang, Aize Kijlstra, Peizeng Yang, Shengping Hou; Association of Long Noncoding RNAs Polymorphisms With Ankylosing Spondylitis, Vogt-Koyanagi-Harada Disease, and Behcet's Disease. Invest. Ophthalmol. Vis. Sci. 2018;59(2):1158-1166. https://doi.org/10.1167/iovs.17-23247.

      Download citation file:


      © ARVO (1962-2015); The Authors (2016-present)

      ×
  • Supplements
Abstract

Purpose: Long noncoding RNAs (lncRNAs) are emerging as important regulators of inflammatory immune responses, whereby genetic variants may affect this biologic function. This study aimed to investigate the association of 110 single nucleotide polymorphisms (SNPs) of lncRNAs, known to be associated with autoimmune disease, in patients with ocular Vogt-Koyanagi-Harada (VKH) disease, Behcet's disease (BD), and acute anterior uveitis (AAU) with or without ankylosing spondylitis (AS).

Methods: A two-stage case-control study was performed on 1626 VKH patients, 384 BD patients, 624 AAU with AS, 751 AAU without AS, 720 AS without AAU, and 3305 healthy subjects. lncRNAs 110 SNPs were genotyped using MassARRAY System or TaqMan SNP assays. The gene expression and cytokine production were measured using real-time PCR or ELISA.

Results: The frequency of the C allele of rs4937362 in RP11-264E20.1 was markedly decreased in the AS without AAU group compared with controls (Combined P = 9.37 × 10−7, odds ratio [OR] = 0.73). An increased frequency of the A allele of rs6871626 between UBLCP1, IL12B, and LOC285627 was found in VKH patients compared with controls (Combined P = 1.88 × 10−4, OR = 1.19). UBLCP1, IL12B, and LOC285627 were expressed in human uveal tissues. Functional studies showed a decreased LOC285627 mRNA expression in peripheral blood mononuclear cells (PBMCs) and an increased IL-10 production in PBMCs following LPS stimulation in rs6871626 CC genotype carriers.

Conclusions: Our study is the first to show that rs4937362/RP11-264E20.1 is associated with AS and that rs6871626 is associated with VKH disease in Chinese Han. The protective rs6871626 genotype was shown to regulate the expression of LOC285627 and to increase the production of the anti-inflammatory cytokine IL-10.

Uveitis is an inflammation of the uveal tract and can lead to vision loss. According to the anatomic location, uveitis can be divided into panuveitis, anterior uveitis, intermediate uveitis, and posterior uveitis.1 Acute anterior uveitis (AAU) is a sudden-onset uveitis and mostly presents in a unilateral fashion. Almost 50% of anterior uveitis patients carry the human leukocyte antigen B27 (HLA-B27).2 Ankylosing spondylitis (AS) is an inflammatory disease that also shows a strong association with HLA-B27.3 Approximately one-third of HLA-B27-positive AS patients also suffer from AAU.4 In addition to AAU, Vogt-Koyanagi-Harada (VKH) disease and Behcet's disease (BD) are the two most general uveitis entities seen in China.5 BD is an inflammatory disease accompanied by the following manifestations: oral aphthae, skin lesions, and genital ulcers. VKH disease is an autoimmune disease characterized by a bilateral granulomatous uveitis that is often associated with systemic symptoms, including dysacusis, alopecia, poliosis, and vitiligo.6 
The exact etiology and pathogenesis of these diseases remain unclear. Many studies have shown that genetic factors and the dysfunction of immune responses play an important role in these diseases. HLA-DR4 and HLA-DRw53 are associated with VKH disease, and HLA-B51 is associated with BD.79 Additionally, many non-HLA genes, including IL17F, IL-23A, TNFAIP3, C3, IL-12B, MIF, TRAF5, and TRAF3IP2, are associated with both BD and VKH disease.10 An abnormal immune response is also involved in these three diseases. The imbalance in T-cell homeostasis, including cytotoxic T cells, Th1cells, Treg cells, and Th17 cells, have been implicated in the pathogenesis of both BD and VKH disease.10,11 Genome-wide association studies (GWAS) showed that loci at 1p31.2 and 10q21.3 were associated with the susceptibility to VKH disease. GWAS also identified the variants in the IL10, MHC class I, and IL23R/IL12RB2 regions were associated with the development of BD.1214 Both AAU and AS are strongly associated with HLA-B27 as well as with a variety on non-HLA genes, such as FOXO1, GATA-3, and FOXP3.15,16 However, some of the associated regions were located on intergenic regions, suggesting that noncoding RNAs (ncRNAs) may be involved in the development of these diseases. 
The mammalian genome has a large number of RNAs that do not encode proteins, which are called ncRNAs.17,18 According to the transcript size, ncRNAs can be divided into many types, including microRNAs (miRNAs) and long ncRNAs (lncRNAs).19 miRNAs range from 19 to 25 nucleotides, and have shown to play important roles in various human immune diseases, including VKH disease, BD, AAU, systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), multiple sclerosis (MS), systemic sclerosis, and primary Sjögren's syndrome.2032 miR-146a had an effect on cytokine production of IL-8, TNF-α, IL-17, and IL-1β in VKH disease.33 miR-155 regulated the mRNA expression of Ets-1 and the Th17 immune response in active BD.34 Copy number variants of miR-143, miR-146a, miR-9-3, miR-205, miR-301a, and miR-23a were shown to be associated with the susceptibility to AAU.30 lncRNAs are a group of noncoding potential transcripts that are longer than 200 nucleotides.19 Over the past few decades, lncRNAs have been considered as transcriptional “noise.” Recent studies have, however, shown that lncRNAs have an important role in various biological and physiological processes, including transcription, chromatin modification, and posttranscriptional processing.17,35,36 A number of studies have suggested that lncRNAs are also involved in the development of autoimmunity.19,3743 lncRNAs can regulate the development of various immune cells, such as B and T lymphocytes, macrophages, and dendritic cells (DCs).44,45 Recent studies showed that several lncRNAs are especially expressed in cells responsible for the innate and adaptive immune response.4446 Examples include LncRNA-CD244, lincRNA-Cox2, THRIL, TH2-LCR, and lnc-DC, which have been shown to be expressed in T cells, macrophages, and DCs.45 These lncRNAs regulate the transcription of immune-response genes and the production of related cytokines, such as IFN-γ, TNF-α, IL-13, IL-5, IL-4, and IL-12, which have all been implicated in the pathogenesis of a variety of uveitis entities.37,4650 
A number of studies have shown that the pathogenesis of autoimmune diseases, such as SLE, RA, and MS, were associated with lncRNAs.19,3743 Lnc-USP50–2, lnc-ZNF354A-1, lnc-FRG2C-3, and lnc-LIN54–1 have been shown to be involved in the development of AS.51 
Up to now, several lncRNAs have been found to be implicated with ocular disorders, such as ocular tumors, glaucoma, diabetic retinopathy, proliferative vitreoretinopathy, and corneal vascularization.52 As mentioned above, lncRNAs play a role in inflammation, and this regulatory control may be affected by genetic variants, which in turn might explain why a certain variant predisposes or protects against immune disease. Up to now, the relationship between lncRNAs and uveitis has yet to be reported, and the current study was therefore performed to address this issue. The results showed that rs6871626 was associated with VKH disease and SNP rs4937362 in RP11-264E20.1 with the group of AS patients without AAU. Preliminary experiments were performed to investigate the biological function of LncRNA genotypes. 
Materials and Methods
Study Population
A total of 1626 VKH patients, 384 patients with BD, 624 AAU+AS+ patients, 751 AAU+AS patients, 720 AS+AAU patients, and 3305 unrelated healthy controls were recruited from The First Affiliated Hospital of Chongqing Medical University (Chongqing, China). All patients and controls were Chinese Han. The study was performed in two stages, whereby patients and controls for the first stage came from Southwest China, including Chongqing and Sichuan. These included patients with BD (n = 384), VKH (n = 384), AAU with AS (n = 384), AAU without AS (n = 384), and AS without AAU (n = 384), and healthy controls (n = 768). A confirmatory second-stage study was performed with a cohort of patients coming from North China (Hebei, Beijing), Central China (Henan, Hubei), and East China (Anhui, Shandong, and Zhejiang). This group included 1242 VKH patients, 240 AAU with AS patients, 367 AAU without AS patients, 336 AS without AAU patients, and 2537 controls. The nonparametric Mann-Whitney U test or χ2 test was used to compare the differences of the clinical characteristics between the first and second disease groups. However, there was no difference between the first- and the second-stage groups (P > 0.05). All uveitis patients were seen at the Department of Ophthalmology, and AS patients were seen by our rheumatology department. The diagnosis of VKH disease and BD was made strictly according to the International Study Group's criteria for VKH disease and BD, respectively.53,54 The diagnosis of AAU was made according to the clinical features, including eye pain, the loss of vision, ciliary congestion, keratic precipitates, anterior chamber flare, iris abnormalities, and the appearance of inflammatory cells in the anterior chamber. Additionally, AAU patients show a sudden onset, as well as recurrent episodes, with a duration of no more than 3 months.55 Optical coherence tomography and ultrasound biomicroscopy were helpful in the diagnosis of AAU.56 The diagnosis of AS was made according to the 1984 modified criteria of New York.57 HLA-B27 was tested in 83% of the AS patients. All participants enrolled in this study signed a written informed consent voluntarily. All procedures of this study followed the principles of the Declaration of Helsinki and were approved by the ethical committee of Chongqing Medical University (Permit Number: 2009-201008). 
The Selection of Single Nucleotide Polymorphisms (SNPs)
The selection of lncRNAs SNPs was based on three databases, including SNP@lincTFBS (http://210.46.85.180:8080/SNP_linc_tfbs/, in the public domain), lncRNASNP-human (http://bioinfo.life.hust.edu.cn/lncRNASNP/, in the public domain), and lincSNP (http://210.46.80.146/lincsnp/, in the public domain). Finally, 110 SNPs, shown previously to be associated with autoimmune disease, were selected for this study (Supplementary Table S1). Linkage disequilibrium was tested between SNPs in case a positive association was found. 
The Extraction of DNA and Genotyping
Genomic DNA from peripheral blood samples was extracted by the QIAamp DNA Mini blood kit (QIAGEN, Valencia, CA, USA) based on the manufacturer's protocols. The extracted DNA was diluted and then stored at −80°C until used. 
The genotypes of lncRNAs SNPs were examined by MassARRAY System (Sequenom, San Diego, CA, USA) or TaqMan SNP assays. The primers used to genotype were made according to MassARRAY Assay design software. All SNPs were genotyped using Sequenom MassARRAY system (Sequenom) based on the manufacturer's manuals in the first stage of our study except rs1991866. Rs1991866 was genotyped by TaqMan SNP Genotyping Assay (Applied Biosystems, Foster City, CA, USA) on the Applied Biosystems 7500 Real-Time PCR system according to the manufacturer's manuals. In the second-stage study, selected SNPs, including rs6871626, rs793108, rs4937362, rs73013527, and rs2156698, were genotyped by TaqMan SNP Genotyping Assay according to the manufacturer's manuals. The results of genotyping assay data were calculated by TYPER software version 4.0 (Sequenom, San Diego, CA, USA) or TaqMan Genotyper Software. 
Cell Preparation, RNA Extraction, and Real-Time Quantitative PCR
Fresh venous blood of healthy subjects was used to isolate peripheral blood mononuclear cells (PBMCs) using Ficoll-Hypaque density gradient centrifugation. Total RNA was extracted from PBMCs using TRIzol reagent (Invitrogen, Carlsbad, CA, USA). The Takara transcriptase kit (Takara, Dalian, China) was used for reverse transcription to synthesize cDNA in accordance with the manufacturer's manuals. Power SYBR Green PCR Master MIX (Applied Biosystems, Warrington, UK) was used to detect the mRNA expression of IL12B, UBLCP1, LOC285627, ETS1, FLI1, and RP11–264E20.1 by real-time quantitative PCR with the ABI 7500 System (Applied Biosystems). The primers of the forward and reverse were as follows: IL12B: 5′-GCTGCTTCACAAAAAGGAAGATG-3′ and 5′-ACCAGCAGGTGAAACGTCCA-3′; UBLCP1: 5′-TCTCGCAGAGTGAAAGAGTACA-3′ and 5′-CTCTGCACAAGACCTGTGGT-3′; ETS1: 5′-GCAGCCAGTCATCTTTCAACAGCC-3′ and 5′- TCAGCACGGTCCCGCACATA-3′; FLI1: 5′-CCAAAGTGCACGGCAAAAGA-3′ and 5′-GGCATGGTAGGAAGGCATGT-3′; RP11–264E20.1: 5′-TGTGACTGTTGAGGGCTGTC-3′ and 5′-CAAGCGCAAAGGGATGTCAC-3′; LOC285627: 5′-CTACATTGCTCCCACCACTC-3′ and 5′-TAGTCGTAACATGGGGCTCT-3′; β-actin: 5′- GGATGCAGAAGGAGATCACTG-3′ and 5′-CGATCCACACGGAGTACTTG-3′. Relative expression levels of genes were calculated using the 2−ΔΔCt method. 
Cell Culture and ELISA
PBMCs isolated from peripheral blood were seeded at 1  ×  106 cells per well in 24-well plates in RPMI 1640 medium containing 10% fetal bovine serum, 100 U/mL penicillin, and 100 μg/mL streptomycin. The supernatants obtained after 24 hours of lipopolysaccharide (LPS) stimulation (100 ng/mL; Sigma-Aldrich Corp., St. Louis, MO, USA) were tested for the concentration of a number of cytokines, including IL-10, IL-1β, IL-17, IL-8, IFN-γ, TNF-α, IL-6, and MCP-1 using human Duoset ELISA development kits (R&D Systems, Minneapolis, MN, USA). 
The Expression of IL12B, LOC285627, and UBLCP1 in the Human Uveal Tissues
Human ocular tissue, including iris, ciliary body, and choroid, was dissected from a Han Chinese donor eye (male, 19 years old, died of leucocythemia) and TRIzol reagent was used to extract RNA. The mRNA expression of UBLCP1, IL12B, and LOC285627 was examined according to the manufacturer's manuals. Amplified production was stained by GoldView (SBS Genetech Beijing, China) and visualized on 3% agarose gels using a camera system (ChemiDoc XRS, Bio-Rad, Hercules, CA, USA). 
Statistical Analysis
The Hardy-Weinberg equilibrium in healthy subjects of all detected SNPs was tested by SHEsis software (Shanghai Jiao Tong University, Shangahi, China). The χ2 test was used to compare the differences of genotype or allele frequency between patients and healthy controls using SPSS 17.0 (SPSS, Inc., Chicago, IL, USA). The obtained P value was adjusted using Bonferroni correction according to the number of analyses. P < 0.05/110 = 4.55 × 10−4 was considered as statistically significant. The nonparametric Mann-Whitney U test or independent samples T test was used to compare the expression of IL12B, UBLCP1, ETS1, FLI1, LOC285627, RP11–264E20.1, and cytokine production (IL-10, IL-1β, IL-17, IL-8, IFN-γ, TNF-α, IL-6, MCP-1) in the various genotype groups. 
Results
Clinical Characteristics of Participating Uveitis and AS Patients
The details of the clinical characteristics, age, and gender of the participants included for the two stages of the study are shown in Table 1. No statistical difference was observed concerning disease manifestations between the two groups. 
Table 1
 
Clinical Characteristics of Participants Enrolled in This Study
Table 1
 
Clinical Characteristics of Participants Enrolled in This Study
The Frequencies of Genotypes and Alleles of Involved lncRNAs SNPs in the First-Phase Study
The frequencies of genotypes and alleles of 110 SNPs of lncRNAs were analyzed in the first-stage discovery cohort that included patients with BD (n = 384), VKH (n = 384), AAU with AS (n = 384), AAU without AS (n = 384), and AS without AAU (n = 384) and were compared with data obtained from healthy controls (n = 768) in the first-stage study. All participants for the first-stage study came from Southwest China, including Chongqing and Sichuan. The results showed a significantly increased frequency of the rs6871626 CA genotype, A allele, and a decreased frequency of the CC genotype (P = 6.04 ×  10−5, odds ratio [OR] = 1.66; P = 5.77 ×  10−5, OR = 1.45; P = 2.13 ×  10−6, OR = 0.54, respectively) in VKH patients as compared with controls (Table 2). We also found increased frequencies for the T allele of rs793108 (Supplementary Table S2), an increased frequency for the TT genotype, and a decreased frequency of the C allele of RP11-264E20.1/rs4937362 in the AS without AAU group as compared with controls (P = 1.28 ×  10−4, OR = 1.42; P = 3.66 × 10−5, OR = 1.69, P = 2.11 × 10−5, OR = 0.67, respectively) (Table 3). No association was observed for the other patient groups tested (Supplementary Tables S3S5). 
Table 2
 
The Association of rs6871626 With VKH Disease
Table 2
 
The Association of rs6871626 With VKH Disease
Table 3
 
The Association of rs4937362, rs2156698, and rs73013527 With AS+AAU
Table 3
 
The Association of rs4937362, rs2156698, and rs73013527 With AS+AAU
The Frequencies of Genotype and Allele of Involved lncRNAs SNPs in the Second and Combined Study
To further verify the first-stage results, we performed a confirmatory study of the disease-associated SNPs rs6871626, rs793108, and rs4937362 using an independent cohort including 1242 VKH patients, 240 AAU with AS patients, 367 AAU without AS patients, 336 AS without AAU patients, and 2537 controls from North China (Hebei, Beijing), Central China (Henan, Hubei), and East China (Anhui, Shandong, and Zhejiang). The results confirmed that the frequency of the rs6871626 CC genotype was decreased and A allele was increased in VKH patients as compared with controls (P = 2.10 × 10−2, OR = 0.85, P = 4.20 × 10−2, OR = 1.11, respectively) (Table 2). Additionally, the frequencies of RP11-264E20.1/rs4937362 TT genotype was increased and C allele was decreased in the AS without AAU patients as compared with controls (P = 2.00 × 10−3, OR = 1.45, P = 1.00 × 10−3, OR = 0.74, respectively) (Table 3). The combined results showed an increased frequency in the RP11-264E20.1/rs4937362 TT genotype and a decreased frequency in the C allele in AS+AAU patients as compared with healthy controls (P = 2.29 × 10−6, OR = 1.49; P = 9.37 × 10−7, OR = 0.73, respectively) (Table 3). In addition, the results showed a decreased frequency of the rs6871626 CC genotype and an increased frequency of the A allele in VKH patients as compared with controls (P = 1.59 × 10−5, OR = 0.76; P = 1.88 × 10−4, OR = 1.19, respectively) (Table 2). To verify the association results of rs4937362, a linkage disequilibrium (LD) block was constructed using Haploview software (Whitehead Institute for Biomedical Research Cambridge, MA, USA) (Supplementary Fig. S1) and the other two SNPs rs2156698 and rs73013527 linked with rs4937362 (r2 = 0.89) were chosen to confirm the association. The results showed a consistent association of rs2156698 and rs73013527 with AS+AAU (P = 4.38 × 10−5, OR = 0.67; P = 1.37 × 10−4, OR = 0.69, respectively) (Table 3). However, no association was found for SNP rs793108 (Supplementary Table S2). 
The Effect of rs6871626, rs4937362 Genotypes on IL12B, UBLCP1, LOC285627, ETS1, FLI1, and RP11-264E20.1
The aforementioned results showed the association of rs6871626 with VKH disease and rs4937362 with AS. Additional experiments were performed to elucidate the functional role of rs6871626 and rs4937362 on gene expression. Because immunosuppressive drug treatment may affect gene expression results of patients, normal genotyped controls were collected for the expression study. The results showed that individuals with the rs6871626/CC genotype had a decreased LOC285627 mRNA expression as compared with AC/AA individuals (Fig. 1a, P = 0.032). However, we did not find any differences between IL12B, UBLCP1, RP11–264E20.1, ETS1, FLI1 mRNA expression and the various genotypes of rs6871626 or rs4937362 (Figs. 1b, 1c; Supplementary Figs. S2a–c). 
Figure 1
 
The influence of diverse genotypes of rs6871626 on the expression of LOC285627 (CC: n = 30, AC+AA: n = 24), IL12B (CC: n = 15, AC+AA: n = 16), UBLCP1 (CC: n = 18, AC+AA: n = 20), in PBMCs of healthy individuals. The mRNA expression levels of LOC285627 (a), IL12B (b), and UBLCP1 (c) were examined in PBMCs of healthy individuals with a variety of genotypes of rs6871626.
Figure 1
 
The influence of diverse genotypes of rs6871626 on the expression of LOC285627 (CC: n = 30, AC+AA: n = 24), IL12B (CC: n = 15, AC+AA: n = 16), UBLCP1 (CC: n = 18, AC+AA: n = 20), in PBMCs of healthy individuals. The mRNA expression levels of LOC285627 (a), IL12B (b), and UBLCP1 (c) were examined in PBMCs of healthy individuals with a variety of genotypes of rs6871626.
The Influence of rs6871626 and rs4937362 Genotypes on Cytokine Production
The secretion of cytokines was examined in supernatants of PBMC cultures from healthy individuals following a 24-hour in vitro stimulation with LPS. The concentrations of IL-10, IL-1β, IL-17, IL-8, IFN-γ, TNF-α, IL-6, and MCP-1 were detected using ELISA (Fig. 2, Supplementary Fig. S3). An increased IL-10 production was found in CC carriers as compared with AA carriers (P = 0.016) (Fig. 2a). No association of the other cytokines with rs6871626 and rs4937362 genotypes was found (Figs. 2b–h, Supplementary Fig. S3). 
Figure 2
 
The effect of rs6871626 on cytokine production. PBMCs, collected from healthy individuals with different rs6871626 genotypes, were stimulated with LPS to examine the level of IL-10 (a), IL-1β (b), IL-17 (c), IL-8 (d), IFN-γ (e), TNF-α (f), IL-6 (g), and MCP-1 (h) (CC: n = 20–37, AC: n = 28–35, AA: n = 6–7).
Figure 2
 
The effect of rs6871626 on cytokine production. PBMCs, collected from healthy individuals with different rs6871626 genotypes, were stimulated with LPS to examine the level of IL-10 (a), IL-1β (b), IL-17 (c), IL-8 (d), IFN-γ (e), TNF-α (f), IL-6 (g), and MCP-1 (h) (CC: n = 20–37, AC: n = 28–35, AA: n = 6–7).
The Expression of IL12B, LOC285627, and UBLCP1 in Human Uveal Tissues
The expression of IL12B, LOC285627, and UBLCP1 was examined in human uveal tissues. The results showed that all three genes were expressed in human iris and ciliary body. Moreover, LOC285627 and UBLCP1 were found also to be expressed in the human choroid (Fig. 3). 
Figure 3
 
The mRNA expression of IL12B, LOC285627, UBLCP1, and β-actin in human uveal tissues (iris, ciliary body, and choroid) from a 19-year-old male who died from leucocythemia.
Figure 3
 
The mRNA expression of IL12B, LOC285627, UBLCP1, and β-actin in human uveal tissues (iris, ciliary body, and choroid) from a 19-year-old male who died from leucocythemia.
Discussion
The current study is the first to show an association of the lncRNAs SNP rs4937362 in RP11-264E20.1 in patients with AS who are not affected with uveitis. In uveitis patients we found that the SNP rs6871626 near the LOC285627, IL12B, and UBLCP1 genes is associated with VKH disease. LOC285627, IL12B, and UBLCP1 were expressed in human uveal tissues. Our study adds to the existing knowledge that has revealed that lncRNAs play important roles in the development of autoimmunity.3743,51 
Functional studies showed a decreased mRNA expression of LOC285627 in PBMCs in carriers of the CC genotype (VKH protective genotype) of rs6871626 as compared with the carriers of other genotypes and a significantly increased IL-10 production by LPS-stimulated PBMCs in rs6871626 carriers of the CC genotype as compared with carriers of the AA genotype. The increased production of an anti-inflammatory cytokine such as IL-10 in CC genotype carriers might provide an explanation for its protective role in disease susceptibility. 
SNP rs6871626 is located between UBLCP1, IL12B, and LOC285627 at 5q33.3.58 Several studies showed that rs6871626 was associated with leprosy, Crohn's disease, Takayasu arteritis, and ulcerative colitis.5860 The current study showed that rs6871626 was associated with VKH disease but not with AS. Our study is in agreement with a previous study that showed that rs6871626 was not associated with the susceptibility to AS in mainland Han Chinese.61 The reason why rs6871626 is only associated with VKH disease but not with the other uveitis groups tested (AAU or BD) is not clear, but may be due to the different etiology between these uveitis entities. BD is an inflammatory disorder with a complex genetic background characterized by uveitis, oral aphthae, skin lesions, and genital ulcers, whereas VKH disease is a systemic autoimmune disease that targets melanocytes not only in the uvea, but also elsewhere in the body, resulting in multiple systemic abnormalities leading to dysacusis, alopecia, poliosis, and vitiligo.62,63 
As mentioned above, rs4937362 was associated with AS in our Chinese Han population but not with the uveitis entities tested. The SNP rs4937362 is located at 11q24.3 in lncRNA RP11-264E20.1 and has been shown to be strongly associated with follicular lymphoma (P = 6.76 × 10−11).64 To further validate our findings, the other SNPs rs2156698 and rs73013527, which are in strong LD with rs4937362, were also genotyped and the results were consistent with the rs4937362 genotyping data. Intriguingly, rs4937362 locates at approximately 35 kb upstream of ETS1, which plays an important role in autoimmune diseases, such as SLE, psoriasis, RA, MS, pediatric uveitis, and hematopoietic tumors.65,66 These findings suggest that RP11–264E20.1 rs4937362 possibility affects disease susceptibility by regulating ETS1. It is not clear why rs4937362 associated with AS+AAU but not in AS patients with AAU. Further studies are needed to investigate whether the manifestations of AS in these two patient groups may differ from one another. 
Our study also has several limitations. We only examined the association of lncRNAs SNPs with uveitis or AS in a Chinese Han population, and further studies in other ethnic populations are needed to generalize our findings. This study only included 110 SNPs of lncRNAs identified to be associated with immune-mediated diseases by earlier studies. and it is possible that we may have missed potential not yet discovered lncRNA risk or protective SNPs. In view of the fact that these patients are often on a treatment regimen with immunosuppressive drugs and that it was difficult to collect sufficient patient samples due to the low frequency of the rs6871626 AA genotype (12.4%) and rs4937362 CC genotype (8.4%) in patients, we only examined the expression of cytokines in genotyped healthy controls. Further studies are needed to validate the role of rs6871626 and rs4937362 in uveitis as well as in AS patients. 
In conclusion, this study shows novel associations between several lncRNAs SNPs with VKH disease and AS, supporting the role of these neglected gene regions in autoimmune disease. 
Acknowledgments
Supported by National Key R&D Program of China (2016YFC0904000), National Natural Science Foundation Project (81522013, 81400389), Natural Science Foundation Major International (Regional) Joint Research Project (81320108009, 81720108009), Chongqing Outstanding Youth Grant (cstc2014jcyjjq10005), National Key Clinical Specialties Construction Program of China, and Chongqing Key Laboratory of Ophthalmology (CSTC, 2008CA5003). 
Disclosure: Y. Yue, None; J. Zhang, None; L. Yang, None; S. Liu, None; J. Qi, None; Q. Cao, None; C. Zhou, None; Y. Wang, None; A. Kijlstra, None; P. Yang, None; S. Hou, None 
References
Suttorp-Schulten MS, Rothova A. The possible impact of uveitis in blindness: a literature survey. Br J Ophthalmol. 1996; 80: 844–848.
Chang JH, McCluskey PJ, Wakefield D. Acute anterior uveitis and HLA-B27. Surv Ophthalmol. 2005; 50: 364–388.
Gouveia EB, Elmann D, Morales MS. Ankylosing spondylitis and uveitis: overview. Rev Bras Reumatol. 2012; 52: 742–756.
Feltkamp TE, Ringrose JH. Acute anterior uveitis and spondyloarthropathies. Curr Opin Rheumatol. 1998; 10: 314–318.
Yang P, Zhang Z, Zhou H, et al. Clinical patterns and characteristics of uveitis in a tertiary center for uveitis in China. Curr Eye Res. 2005; 30: 943–948.
Yang P, Ren Y, Li B, Fang W, Meng Q, Kijlstra A. Clinical characteristics of Vogt-Koyanagi-Harada syndrome in Chinese patients. Ophthalmology. 2007; 114: 606–614.
Ohno S. Immunological aspects of Behcet's and Vogt-Koyanagi-Harada's diseases. Trans Ophthalmol Soc U K. 1981; 101: 335–341.
Zhao M, Jiang Y, Abrahams IW. Association of HLA antigens with Vogt-Koyanagi-Harada syndrome in a Han Chinese population. Arch Ophthalmol. 1991; 109: 368–370.
Gul A, Ohno S. HLA-B*51 and Behcet Disease. Ocul Immunol Inflamm. 2012; 20: 37–43.
Du L, Kijlstra A, Yang P. Vogt-Koyanagi-Harada disease: novel insights into pathophysiology, diagnosis and treatment. Prog Retin Eye Res. 2016; 52: 84–111.
Kim J, Park JA, Lee EY, Lee YJ, Song YW, Lee EB. Imbalance of Th17 to Th1 cells in Behcet's disease. Clin Exp Rheumatol. 2010; 28: S16–S19.
Hou S, Du L, Lei B, et al. Genome-wide association analysis of Vogt-Koyanagi-Harada syndrome identifies two new susceptibility loci at 1p31.2 and 10q21.3. Nat Genet. 2014; 46: 1007–1011.
Remmers EF, Cosan F, Kirino Y, et al. Genome-wide association study identifies variants in the MHC class I, IL10, and IL23R-IL12RB2 regions associated with Behcet's disease. Nat Genet. 2010; 42: 698–702.
Mizuki N, Meguro A, Ota M, et al. Genome-wide association studies identify IL23R-IL12RB2 and IL10 as Behcet's disease susceptibility loci. Nat Genet. 2010; 42: 703–706.
Yu H, Liu Y, Zhang L, et al. FoxO1 gene confers genetic predisposition to acute anterior uveitis with ankylosing spondylitis. Invest Ophthalmol Vis Sci. 2014; 55: 7970–7974.
Bai L, Liu Y, Hou S, Liao D, Kijlstra A, Yang P. Association of T-Bet, GATA-3, RORC, and FOXP3 copy number variations with acute anterior uveitis with or without ankylosing spondylitis in Chinese Han. Invest Ophthalmol Vis Sci. 2016; 57: 1847–1852.
Mercer TR, Dinger ME, Mattick JS. Long non-coding RNAs: insights into functions. Nat Rev Genet. 2009; 10: 155–159.
Saha P, Verma S, Pathak RU, Mishra RK. Long noncoding RNAs in mammalian development and diseases. Adv Exp Med Biol. 2017; 1008: 155–198.
Wu GC, Pan HF, Leng RX, et al. Emerging role of long noncoding RNAs in autoimmune diseases. Autoimmun Rev. 2015; 14: 798–805.
Iorio MV, Visone R, Di Leva G, et al. MicroRNA signatures in human ovarian cancer. Cancer Res. 2007; 67: 8699–8707.
Singh RP, Massachi I, Manickavel S, et al. The role of miRNA in inflammation and autoimmunity. Autoimmun Rev. 2013; 12: 1160–1165.
Jimenez SA, Piera-Velazquez S. Potential role of human-specific genes, human-specific microRNAs and human-specific non-coding regulatory RNAs in the pathogenesis of systemic sclerosis and Sjogren's syndrome. Autoimmun Rev. 2013; 12: 1046–1051.
Chang R, Yi S, Tan X, et al. MicroRNA-20a-5p suppresses IL-17 production by targeting OSM and CCL1 in patients with Vogt-Koyanagi-Harada disease. Br J Ophthalmol. 2017; 102: 282–290.
Hou S, Ye Z, Liao D, et al. miR-23a, miR-146a and miR-301a confer predisposition to Vogt-Koyanagi-Harada syndrome but not to Behcet's disease. Sci Rep. 2016; 6: 20057.
Yu H, Liu Y, Bai L, Kijlstra A, Yang P. Predisposition to Behcet's disease and VKH syndrome by genetic variants of miR-182. J Mol Med. 2014; 92: 961–967.
Fang J, Hou S, Xiang Q, et al. Polymorphisms in genetics of vitamin D metabolism confer susceptibility to ocular Behcet disease in a Chinese Han population. Am J Ophthalmol. 2014; 157: 488–494.e6.
Qi J, Hou S, Zhang Q, et al. A functional variant of pre-miRNA-196a2 confers risk for Behcet's disease but not for Vogt-Koyanagi-Harada syndrome or AAU in ankylosing spondylitis. Hum Genet. 2013; 132: 1395–1404.
Woo MY, Yun SJ, Cho O, Kim K, Lee ES, Park S. MicroRNAs differentially expressed in Behcet disease are involved in interleukin-6 production. J Inflamm (Lond). 2016; 13: 22.
Erre GL, Piga M, Carru C, et al. Global microRNA profiling of peripheral blood mononuclear cells in patients with Behcet's disease. Clin Exp Rheumatol. 2015; 33: S72–S79.
Yang L, Du L, Yue Y, et al. miRNA copy number variants confer susceptibility to acute anterior uveitis with or without ankylosing spondylitis. Invest Ophthalmol Vis Sci. 2017; 58: 1991–2001.
Alipour S, Nouri M, Sakhinia E, et al. Epigenetic alterations in chronic disease focusing on Behcet's disease: Review. Biomed Pharmacother. 2017; 91: 526–533.
Ugurel E, Sehitoglu E, Tuzun E, Kurtuncu M, Coban A, Vural B. Increased complexin-1 and decreased miR-185 expression levels in Behcet's disease with and without neurological involvement. Neurol Sci. 2016; 37: 411–416.
Zhou Q, Hou S, Liang L, et al. MicroRNA-146a and Ets-1 gene polymorphisms in ocular Behcet's disease and Vogt-Koyanagi-Harada syndrome. Ann Rheum Dis. 2014; 73: 170–176.
Na SY, Park MJ, Park S, Lee ES. MicroRNA-155 regulates the Th17 immune response by targeting Ets-1 in Behcet's disease. Clin Exp Rheumatol. 2016; 34: S56–S63.
Aryal B, Rotllan N, Fernandez-Hernando C. Noncoding RNAs and atherosclerosis. Curr Atheroscler Rep. 2014; 16: 407.
Esteller M. Non-coding RNAs in human disease. Nat Rev Genet. 2011; 12: 861–874.
Wang P, Xue Y, Han Y, et al. The STAT3-binding long noncoding RNA lnc-DC controls human dendritic cell differentiation. Science. 2014; 344: 310–313.
Zhang X, Lian Z, Padden C, et al. A myelopoiesis-associated regulatory intergenic noncoding RNA transcript within the human HOXA cluster. Blood. 2009; 113: 2526–2534.
Williams GT, Mourtada-Maarabouni M, Farzaneh F. A critical role for non-coding RNA GAS5 in growth arrest and rapamycin inhibition in human T-lymphocytes. Biochem Soc Trans. 2011; 39: 482–486.
Verma-Gaur J, Torkamani A, Schaffer L, Head SR, Schork NJ, Feeney AJ. Noncoding transcription within the Igh distal V(H) region at PAIR elements affects the 3D structure of the Igh locus in pro-B cells. Proc Natl Acad U S A. 2012; 109: 17004–17009.
Long H, Yin H, Wang L, Gershwin ME, Lu Q. The critical role of epigenetics in systemic lupus erythematosus and autoimmunity. J Autoimmun. 2016; 74: 118–138.
Song J, Kim D, Han J, Kim Y, Lee M, Jin EJ. PBMC and exosome-derived Hotair is a critical regulator and potent marker for rheumatoid arthritis. Clin Exp Med. 2015; 15: 121–126.
Vigneau S, Rohrlich PS, Brahic M, Bureau JF. Tmevpg1, a candidate gene for the control of Theiler's virus persistence, could be implicated in the regulation of gamma interferon. J Virol. 2003; 77: 5632–5638.
Atianand MK, Caffrey DR, Fitzgerald KA. Immunobiology of long noncoding RNAs. Ann Rev Immunol. 2017; 35: 177–198.
Chen YG, Satpathy AT, Chang HY. Gene regulation in the immune system by long noncoding RNAs. Nat Immunol. 2017; 18: 962–972.
Spurlock CFIII, Tossberg JT, Guo Y, Collier SP, Crooke PSIII, Aune TM. Expression and functions of long noncoding RNAs during human T helper cell differentiation. Nat Commun. 2015; 6: 6932.
Kucuksezer UC, Aktas-Cetin E, Bilgic-Gazioglu S, Tugal-Tutkun I, Gul A, Deniz G. Natural killer cells dominate a Th-1 polarized response in Behcet's disease patients with uveitis. Clin Exp Rheumatol. 2015; 33: S24–S29.
Ooi KG, Galatowicz G, Towler HM, Lightman SL, Calder VL. Multiplex cytokine detection versus ELISA for aqueous humor: IL-5, IL-10, and IFNgamma profiles in uveitis. Invest Ophthalmol Vis Sci. 2006; 47: 272–277.
Ilott NE, Heward JA, Roux B, et al. Long non-coding RNAs and enhancer RNAs regulate the lipopolysaccharide-induced inflammatory response in human monocytes. Nat Commun. 2014; 5: 3979.
Wang Y, Zhong H, Xie X, et al. Long noncoding RNA derived from CD244 signaling epigenetically controls CD8+ T-cell immune responses in tuberculosis infection. Proc Natl Acad U S A. 2015; 112: E3883–E3892.
Xie Z, Li J, Wang P, et al. Differential expression profiles of long noncoding RNA and mRNA of osteogenically differentiated mesenchymal stem cells in ankylosing spondylitis. J Rheumatol. 2016; 43: 1523–1531.
Li F, Wen X, Zhang H, Fan X. Novel insights into the role of long noncoding RNA in ocular diseases. Int J Mol Sci. 2016; 17: 478.
Read RW, Holland GN, Rao NA, et al. Revised diagnostic criteria for Vogt-Koyanagi-Harada disease: report of an international committee on nomenclature. Am J Ophthalmol. 2001; 131: 647–652.
International Study Group for Behcet's Disease. Criteria for diagnosis of Behcet's disease. Lancet. 1990; 335: 1078–1080.
Jabs DA, Nussenblatt RB, Rosenbaum JT; for the Standardization of Uveitis Nomenclature Working Group. Standardization of uveitis nomenclature for reporting clinical data. Results of the First International Workshop. Am J Ophthalmol. 2005; 140: 509–516.
Yang P, Wan W, Du L, et al. Clinical features of HLA-B27-positive acute anterior uveitis with or without ankylosing spondylitis in a Chinese cohort. Br J Ophthalmol. 2017; 102: 215–219.
van der Linden S, Valkenburg HA, Cats A. Evaluation of diagnostic criteria for ankylosing spondylitis. A proposal for modification of the New York criteria. Arthritis Rheum. 1984; 27: 361–368.
Liu H, Irwanto A, Tian H, et al. Identification of IL18RAP/IL18R1 and IL12B as leprosy risk genes demonstrates shared pathogenesis between inflammation and infectious diseases. Am J Hum Genet. 2012; 91: 935–941.
Terao C, Yoshifuji H, Kimura A, et al. Two susceptibility loci to Takayasu arteritis reveal a synergistic role of the IL12B and HLA-B regions in a Japanese population. Am J Hum Genet. 2013; 93: 289–297.
Anderson CA, Boucher G, Lees CW, et al. Meta-analysis identifies 29 additional ulcerative colitis risk loci, increasing the number of confirmed associations to 47. Nat Genet. 2011; 43: 246–252.
Zhang L, Fan D, Liu L, et al. Association study of IL-12B polymorphisms susceptibility with ankylosing spondylitis in mainland Han Population. PLoS One. 2015; 10: e0130982.
Gul A. Behcet's disease as an autoinflammatory disorder. Curr Drug Targets Inflamm Allergy. 2005; 4: 81–83.
Lavezzo MM, Sakata VM, Morita C, et al. Vogt-Koyanagi-Harada disease: review of a rare autoimmune disease targeting antigens of melanocytes. Orphanet J Rare Dis. 2016; 11: 29.
Skibola CF, Berndt SI, Vijai J, et al. Genome-wide association study identifies five susceptibility loci for follicular lymphoma outside the HLA region. Am J Hum Genet. 2014; 95: 462–471.
Wei L, Zhou Q, Hou S, et al. MicroRNA-146a and Ets-1 gene polymorphisms are associated with pediatric uveitis. PLoS One. 2014; 9: e91199.
Garrett-Sinha LA. Review of Ets1 structure, function, and roles in immunity. Cell Mol Life Sci. 2013; 70: 3375–3390.
Figure 1
 
The influence of diverse genotypes of rs6871626 on the expression of LOC285627 (CC: n = 30, AC+AA: n = 24), IL12B (CC: n = 15, AC+AA: n = 16), UBLCP1 (CC: n = 18, AC+AA: n = 20), in PBMCs of healthy individuals. The mRNA expression levels of LOC285627 (a), IL12B (b), and UBLCP1 (c) were examined in PBMCs of healthy individuals with a variety of genotypes of rs6871626.
Figure 1
 
The influence of diverse genotypes of rs6871626 on the expression of LOC285627 (CC: n = 30, AC+AA: n = 24), IL12B (CC: n = 15, AC+AA: n = 16), UBLCP1 (CC: n = 18, AC+AA: n = 20), in PBMCs of healthy individuals. The mRNA expression levels of LOC285627 (a), IL12B (b), and UBLCP1 (c) were examined in PBMCs of healthy individuals with a variety of genotypes of rs6871626.
Figure 2
 
The effect of rs6871626 on cytokine production. PBMCs, collected from healthy individuals with different rs6871626 genotypes, were stimulated with LPS to examine the level of IL-10 (a), IL-1β (b), IL-17 (c), IL-8 (d), IFN-γ (e), TNF-α (f), IL-6 (g), and MCP-1 (h) (CC: n = 20–37, AC: n = 28–35, AA: n = 6–7).
Figure 2
 
The effect of rs6871626 on cytokine production. PBMCs, collected from healthy individuals with different rs6871626 genotypes, were stimulated with LPS to examine the level of IL-10 (a), IL-1β (b), IL-17 (c), IL-8 (d), IFN-γ (e), TNF-α (f), IL-6 (g), and MCP-1 (h) (CC: n = 20–37, AC: n = 28–35, AA: n = 6–7).
Figure 3
 
The mRNA expression of IL12B, LOC285627, UBLCP1, and β-actin in human uveal tissues (iris, ciliary body, and choroid) from a 19-year-old male who died from leucocythemia.
Figure 3
 
The mRNA expression of IL12B, LOC285627, UBLCP1, and β-actin in human uveal tissues (iris, ciliary body, and choroid) from a 19-year-old male who died from leucocythemia.
Table 1
 
Clinical Characteristics of Participants Enrolled in This Study
Table 1
 
Clinical Characteristics of Participants Enrolled in This Study
Table 2
 
The Association of rs6871626 With VKH Disease
Table 2
 
The Association of rs6871626 With VKH Disease
Table 3
 
The Association of rs4937362, rs2156698, and rs73013527 With AS+AAU
Table 3
 
The Association of rs4937362, rs2156698, and rs73013527 With AS+AAU
Supplement 1
Supplement 2
×
×

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.

×