The genetic association of RUNX3 with ankylosing spondylitis can be explained by allele-specific effects on IRF4 recruitment that alter gene expression.

OBJECTIVES: To identify the functional basis for the genetic association of single nucleotide polymorphisms (SNP), upstream of the RUNX3 promoter, with ankylosing spondylitis (AS). METHODS: We performed conditional analysis of genetic association data and used ENCODE data on chromatin remodelling and transcription factor (TF) binding sites to identify the primary AS-associated regulatory SNP in the RUNX3 region. The functional effects of this SNP were tested in luciferase reporter assays. Its effects on TF binding were investigated by electrophoretic mobility gel shift assays and chromatin immunoprecipitation. RUNX3 mRNA levels were compared in primary CD8+ T cells of AS risk and protective genotypes by real-time PCR. RESULTS: The association of the RUNX3 SNP rs4648889 with AS (p<7.6×10(-14)) was robust to conditioning on all other SNPs in this region. We identified a 2 kb putative regulatory element, upstream of RUNX3, containing rs4648889. In reporter gene constructs, the protective rs4648889 'G' allele increased luciferase activity ninefold but significantly less activity (4.3-fold) was seen with the AS risk 'A' allele (p≤0.01). The binding of Jurkat or CD8+ T-cell nuclear extracts to the risk allele was decreased and IRF4 recruitment was reduced. The AS-risk allele also affected H3K4Me1 histone methylation and associated with an allele-specific reduction in RUNX3 mRNA (p<0.05). CONCLUSION: We identified a regulatory region upstream of RUNX3 that is modulated by rs4648889. The risk allele decreases TF binding (including IRF4) and reduces reporter activity and RUNX3 expression. These findings may have important implications for understanding the role of T cells and other immune cells in AS.

An ankylosing spondylitis-associated genetic variant in the IL23R-IL12RB2 intergenic region modulates enhancer activity and is associated with increased Th1-cell differentiation.

OBJECTIVES: To explore the functional basis for the association between ankylosing spondylitis (AS) and single-nucleotide polymorphisms (SNPs) in the IL23R-IL12RB2 intergenic region. METHODS: We performed conditional analysis on genetic association data and used epigenetic data on chromatin remodelling and transcription factor (TF) binding to identify the primary AS-associated IL23R-IL12RB2 intergenic SNP. Functional effects were tested in luciferase reporter assays in HEK293T cells and allele-specific TF binding was investigated by electrophoretic mobility gel shift assays. IL23R and IL12RB2 mRNA levels in CD4+ T cells were compared between cases homozygous for the AS-risk 'A' allele and the protective 'G' allele. The proportions of interleukin (IL)-17A+ and interferon (IFN)-γ+ CD4+ T-cells were measured by fluorescence-activated cell sorting and compared between these AS-risk and protective genotypes. RESULTS: Conditional analysis identified rs11209032 as the probable causal SNP within a 1.14 kb putative enhancer between IL23R and IL12RB2. Reduced luciferase activity was seen for the risk allele (p<0.001) and reduced H3K4me1 methylation observed in CD4+ T-cells from 'A/A' homozygotes (p=0.02). The binding of nuclear extract to the risk allele was decreased ∼3.5-fold compared with the protective allele (p<0.001). The proportion of IFN-γ+ CD4+ T-cells was increased in 'A/A' homozygotes (p=0.004), but neither IL23R nor IL12RB2 mRNA was affected. CONCLUSIONS: The rs11209032 SNP downstream of IL23R forms part of an enhancer, allelic variation of which may influence Th1-cell numbers. Homozygosity for the risk 'A' allele is associated with more IFN-γ-secreting (Th1) cells. Further work is necessary to explain the mechanisms for these important observations.

Structural changes of the brain in rheumatoid arthritis.

OBJECTIVE: To investigate whether structural changes are present in the cortical and subcortical gray matter of the brains of patients with rheumatoid arthritis (RA). METHODS: We used two surface-based style morphometry analysis programs and a voxel-based style analysis program to compare high-resolution structural magnetic resonance imaging data obtained for 31 RA patients and 25 age- and sex-matched healthy control subjects. RESULTS: We observed an increase in gray matter content in the basal ganglia of RA patients, mainly in the nucleus accumbens and caudate nucleus. There were no differences in the cortical gray matter. Moreover, patients had a smaller intracranial volume. CONCLUSION: Our results suggest that RA is associated with changes in the subcortical gray matter rather than with cortical gray matter atrophy. Since the basal ganglia play an important role in motor control as well as in pain processing and in modulating behavior in response to aversive stimuli, we suggest that these changes may result from altered motor control or prolonged pain processing. The differences in brain volume may reflect either generalized atrophy or differences in brain development.

Association of variants at 1q32 and STAT3 with ankylosing spondylitis suggests genetic overlap with Crohn's disease.

Ankylosing spondylitis (AS) is a common inflammatory arthritic condition. Overt inflammatory bowel disease (IBD) occurs in about 10% of AS patients, and in addition 70% of AS cases may have subclinical terminal ileitis. Spondyloarthritis is also common in IBD patients. We therefore tested Crohn's disease susceptibility genes for association with AS, aiming to identify pleiotropic genetic associations with both diseases. Genotyping was carried out using Sequenom and Applied Biosystems TaqMan and OpenArray technologies on 53 markers selected from 30 Crohn's disease associated genomic regions. We tested genotypes in a population of unrelated individual cases (n = 2,773) and controls (n = 2,215) of white European ancestry for association with AS. Statistical analysis was carried out using a Cochran-Armitage test for trend in PLINK. Strong association was detected at chr1q32 near KIF21B (rs11584383, P = 1.6 × 10(-10), odds ratio (OR) = 0.74, 95% CI:0.68-0.82). Association with disease was also detected for 2 variants within STAT3 (rs6503695, P = 4.6 × 10(-4). OR = 0.86 (95% CI:0.79-0.93); rs744166, P = 2.6 × 10(-5), OR = 0.84 (95% CI:0.77-0.91)). Association was confirmed for IL23R (rs11465804, P = 1.2 × 10(-5), OR = 0.65 (95% CI:0.54-0.79)), and further associations were detected for IL12B (rs10045431, P = 5.2 × 10(-5), OR = 0.83 (95% CI:0.76-0.91)), CDKAL1 (rs6908425, P = 1.1 × 10(-4), OR = 0.82 (95% CI:0.74-0.91)), LRRK2/MUC19 (rs11175593, P = 9.9 × 10(-5), OR = 1.92 (95% CI: 1.38-2.67)), and chr13q14 (rs3764147, P = 5.9 × 10(-4), OR = 1.19 (95% CI: 1.08-1.31)). Excluding cases with clinical IBD did not significantly affect these findings. This study identifies chr1q32 and STAT3 as ankylosing spondylitis susceptibility loci. It also further confirms association for IL23R and detects suggestive association with another 4 loci. STAT3 is a key signaling molecule within the Th17 lymphocyte differentiation pathway and further enhances the case for a major role of this T-lymphocyte subset in ankylosing spondylitis. Finally these findings suggest common aetiopathogenic pathways for AS and Crohn's disease and further highlight the involvement of common risk variants across multiple diseases.

Chromosome conformation capture approaches to investigate 3D genome architecture in Ankylosing Spondylitis.

Ankylosing Spondylitis (AS) is a chronic inflammatory arthritis of the spine exhibiting a strong genetic background. The mechanistic and functional understanding of the AS-associated genomic loci, identified with Genome Wide Association Studies (GWAS), remains challenging. Chromosome conformation capture (3C) and derivatives are recent techniques which are of great help in elucidating the spatial genome organization and of enormous support in uncover a mechanistic explanation for disease-associated genetic variants. The perturbation of three-dimensional (3D) genome hierarchy may lead to a plethora of human diseases, including rheumatological disorders. Here we illustrate the latest approaches and related findings on the field of genome organization, highlighting how the instability of 3D genome conformation may be among the causes of rheumatological disease phenotypes. We suggest a new perspective on the inclusive potential of a 3C approach to inform GWAS results in rheumatic diseases. 3D genome organization may ultimately lead to a more precise and comprehensive functional interpretation of AS association, which is the starting point for emerging and more specific therapies.

Comprehensive epigenomic profiling reveals the extent of disease-specific chromatin states and informs target discovery in ankylosing spondylitis.

Ankylosing spondylitis (AS) is a common, highly heritable inflammatory arthritis characterized by enthesitis of the spine and sacroiliac joints. Genome-wide association studies (GWASs) have revealed more than 100 genetic associations whose functional effects remain largely unresolved. Here, we present a comprehensive transcriptomic and epigenomic map of disease-relevant blood immune cell subsets from AS patients and healthy controls. We find that, while CD14+ monocytes and CD4+ and CD8+ T cells show disease-specific differences at the RNA level, epigenomic differences are only apparent upon multi-omics integration. The latter reveals enrichment at disease-associated loci in monocytes. We link putative functional SNPs to genes using high-resolution Capture-C at 10 loci, including PTGER4 and ETS1, and show how disease-specific functional genomic data can be integrated with GWASs to enhance therapeutic target discovery. This study combines epigenetic and transcriptional analysis with GWASs to identify disease-relevant cell types and gene regulation of likely pathogenic relevance and prioritize drug targets.

Investigating the genetic association between ERAP1 and ankylosing spondylitis.

A strong association between ERAP1 and ankylosing spondylitis (AS) was recently identified by the Wellcome Trust Case Control Consortium and the Australo-Anglo-American Spondylitis Consortium (WTCCC-TASC) study. ERAP1 is highly polymorphic with strong linkage disequilibrium evident across the gene. We therefore conducted a series of experiments to try to identify the primary genetic association(s) with ERAP1. We replicated the original associations in an independent set of 730 patients and 1021 controls, resequenced ERAP1 to define the full extent of coding polymorphisms and tested all variants in additional association studies. The genetic association with ERAP1 was independently confirmed; the strongest association was with rs30187 in the replication set (P = 3.4 x 10(-3)). When the data were combined with the original WTCCC-TASC study the strongest association was with rs27044 (P = 1.1 x 10(-9)). We identified 33 sequence polymorphisms in ERAP1, including three novel and eight known non-synonymous polymorphisms. We report several new associations between AS and polymorphisms distributed across ERAP1 from the extended case-control study, the most significant of which was with rs27434 (P = 4.7 x 10(-7)). Regression analysis failed to identify a primary association clearly; we therefore used data from HapMap to impute genotypes for an additional 205 non-coding SNPs located within and adjacent to ERAP1. A number of highly significant associations (P < 5 x 10(-9)) were identified in regulatory sequences which are good candidates for causing susceptibility to AS, possibly by regulating ERAP1 expression.

Perspectives on the Genetic Associations of Ankylosing Spondylitis.

Ankylosing spondylitis (AS) is a common form of inflammatory spinal arthritis with a complex polygenic aetiology. Genome-wide association studies have identified more than 100 loci, including some involved in antigen presentation (HLA-B27, ERAP1, and ERAP2), some in Th17 responses (IL6R, IL23R, TYK2, and STAT3), and others in macrophages and T-cells (IL7R, CSF2, RUNX3, and GPR65). Such observations have already helped identify potential new therapies targeting IL-17 and GM-CSF. Most AS genetic associations are not in protein-coding sequences but lie in intergenic regions where their direct relationship to particular genes is difficult to assess. They most likely reflect functional polymorphisms concerned with cell type-specific regulation of gene expression. Clarifying the nature of these associations should help to understand the pathogenic pathways involved in AS better and suggest potential cellular and molecular targets for drug therapy. However, even identifying the precise mechanisms behind the extremely strong HLA-B27 association with AS has so far proved elusive. Polygenic risk scores (using all the known genetic associations with AS) can be effective for the diagnosis of AS, particularly where there is a relatively high pre-test probability of AS. Genetic prediction of disease outcomes and response to biologics is not currently practicable.

Disruption of c-MYC Binding and Chromosomal Looping Involving Genetic Variants Associated With Ankylosing Spondylitis Upstream of the RUNX3 Promoter.

Background: Ankylosing Spondylitis (AS) is a common form of inflammatory spinal arthritis with a complex aetiology and high heritability, involving more than 100 genetic associations. These include several AS-associated single nucleotide polymorphisms (SNPs) upstream of RUNX3, which encodes the multifunctional RUNT-related transcription factor (TF) 3. The lead associated SNP rs6600247 (p = 2.6 × 10-15) lies ∼13kb upstream of the RUNX3 promoter adjacent to a c-MYC TF binding-site. The effect of rs6600247 genotype on DNA binding and chromosome looping were investigated by electrophoretic mobility gel shift assays (EMSA), Western blotting-EMSA (WEMSA) and Chromosome Conformation Capture (3C). Results: Interrogation of ENCODE published data showed open chromatin in the region overlapping rs6600247 in primary human CD14+ monocytes, in contrast to the Jurkat T cell line or primary human T-cells. The rs6600247 AS-risk allele is predicted to specifically disrupt a c-MYC binding-site. Using a 50bp DNA probe spanning rs6600247 we consistently observed reduced binding to the AS-risk "C" allele of both purified c-MYC protein and nuclear extracts (NE) from monocyte-like U937 cells. WEMSA on U937 NE and purified c-MYC protein confirmed these differences (n = 3; p < 0.05). 3C experiments demonstrated negligible interaction between the region encompassing rs6600247 and the RUNX3 promoter. A stronger interaction frequency was demonstrated between the RUNX3 promoter and the previously characterised AS-associated SNP rs4648889. Conclusion: The lead SNP rs6600247, located in an enhancer-like region upstream of the RUNX3 promoter, modulates c-MYC binding. However, the region encompassing rs6600247 has rather limited physical interaction with the promoter of RUNX3. In contrast a clear chromatin looping event between the region encompassing rs4648889 and the RUNX3 promoter was observed. These data provide further evidence for complexity in the regulatory elements upstream of the RUNX3 promoter and the involvement of RUNX3 transcriptional regulation in AS.

Functional Genomic Analysis of a RUNX3 Polymorphism Associated With Ankylosing Spondylitis.

OBJECTIVE: To investigate the functional consequences of the single-nucleotide polymorphism rs4648889 in a putative enhancer upstream of the RUNX3 promoter associated with susceptibility to ankylosing spondylitis (AS). METHODS: Using nuclear extracts from Jurkat cells and primary human CD8+ T cells, the effects of rs4648889 on allele-specific transcription factor (TF) binding were investigated by DNA pull-down assay and quantitative mass spectrometry (qMS), with validation by electrophoretic mobility shift assay (EMSA), Western blotting of the pulled-down eluates, and chromatin immunoprecipitation (ChIP)-quantitative polymerase chain reaction (qPCR) analysis. Further functional effects were tested by small interfering RNA knockdown of the gene for interferon regulatory factor 5 (IRF5), followed by reverse transcription-qPCR (RT-qPCR) and enzyme-linked immunosorbent assay (ELISA) to measure the levels of IFNγ messenger RNA (mRNA) and protein, respectively. RESULTS: In nuclear extracts from CD8+ T cells, results of qMS showed that relative TF binding to the AS-risk A allele of rs4648889 was increased 3.7-fold (P < 0.03) for Ikaros family zinc-finger protein 3 (IKZF3; Aiolos) and components of the NuRD complex, including chromodomain helicase DNA binding protein 4 (CHD4) (3.6-fold increase; P < 0.05) and retinoblastoma binding protein 4 (RBBP4) (4.1-fold increase; P < 0.03). In contrast, IRF5 bound significantly more to the AS-protective G allele compared to the AS-risk A allele (fold change 8.2; P = 0.003). Validation with Western blotting, EMSA, and ChIP-qPCR confirmed the differential allelic binding of IKZF3, CHD4, RBBP4, and IRF5. Silencing of IRF5 in CD8+ T cells increased the levels of IFNγ mRNA as measured by RT-qPCR (P = 0.03) and IFNγ protein as measured by ELISA (P = 0.02). CONCLUSION: These findings suggest that the association of rs4648889 with AS reflects allele-specific binding of this enhancer-like region to certain TFs, including IRF5, IKZF3, and members of the NuRD complex. IRF5 may have crucial influences on the functions of CD8+ lymphocytes, a finding that could reveal new therapeutic targets for the management of AS.

The chromosome 16q region associated with ankylosing spondylitis includes the candidate gene tumour necrosis factor receptor type 1-associated death domain (TRADD).

OBJECTIVE: To replicate and refine the reported association of ankylosing spondylitis (AS) with two non-synonymous single nucleotide polymorphisms (nsSNPs) on chromosome 16q22.1. METHODS: Firstly, 730 independent UK patients with AS were genotyped for rs9939768 and rs6979 and allele frequencies were compared with 2879 previously typed historic disease controls. Secondly, the two data sets were combined in meta-analyses. Finally, 5 tagging SNPs, located between rs9939768 and rs6979, were analysed in 1604 cases and 1020 controls. RESULTS: The association of rs6979 with AS was replicated, p=0.03, OR=1.14 (95% CI 1.01 to 1.28), and a trend for association with rs9939768 detected, p=0.06, OR=1.25 (95% CI 0.99 to 1.57). Meta-analyses revealed association of both SNPs with AS, p=0.0008, OR=1.31 (95% CI 1.12 to 1.54) and p=0.0009, OR=1.15 (95% CI 1.06 to 1.23) for rs9939768 and rs6979, respectively. New associations with rs9033 and rs868213 (p=0.00002, OR=1.23 (95% CI 1.12 to 1.36) and p=0.00002 OR=1.45 (95% CI 1.22 to 1.72), respectively, were identified. CONCLUSIONS: The region on chromosome 16 that has been replicated in the present work is interesting as the highly plausible candidate gene, tumour necrosis factor receptor type 1 (TNFR1)-associated death domain (TRADD), is located between rs9033 and rs868213. It will require additional work to identify the primary genetic association(s) with AS.

Differences in intracellular localisation of ANKH mutants that relate to mechanisms of calcium pyrophosphate deposition disease and craniometaphyseal dysplasia.

ANKH mutations are associated with calcium pyrophosphate deposition disease and craniometaphyseal dysplasia. This study investigated the effects of these ANKH mutants on cellular localisation and associated biochemistry. We generated four ANKH overexpression-plasmids containing either calcium pyrophosphate deposition disease or craniometaphyseal dysplasia linked mutations: P5L, E490del and S375del, G389R. They were transfected into CH-8 articular chondrocytes and HEK293 cells. The ANKH mutants dynamic differential localisations were imaged and we investigated the interactions with the autophagy marker LC3. Extracellular inorganic pyrophosphate, mineralization, ENPP1 activity expression of ENPP1, TNAP and PIT-1 were measured. P5L delayed cell membrane localisation but once recruited into the membrane it increased extracellular inorganic pyrophosphate, mineralization, and ENPP1 activity. E490del remained mostly cytoplasmic, forming punctate co-localisations with LC3, increased mineralization, ENPP1 and ENPP1 activity with an initial but unsustained increase in TNAP and PIT-1. S375del trended to decrease extracellular inorganic pyrophosphate, increase mineralization. G389R delayed cell membrane localisation, trended to decrease extracellular inorganic pyrophosphate, increased mineralization and co-localised with LC3. Our results demonstrate a link between pathological localisation of ANKH mutants with different degrees in mineralization. Furthermore, mutant ANKH functions are related to synthesis of defective proteins, inorganic pyrophosphate transport, ENPP1 activity and expression of ENPP1, TNAP and PIT-1.

The IL-17/IL-23 Axis and Its Genetic Contribution to Psoriatic Arthritis.

Psoriatic arthritis (PsA) is a chronic inflammatory disease belonging to the family of spondyloarthropathies (SpA). PsA commonly aggravates psoriasis of the skin and frequently manifests as an oligoarthritis with axial skeletal involvement and extraarticular manifestations including dactylitis, enthesitis, and uveitis. The weight of genetic predisposition to psoriasis and PsA is illustrated by the concordance rates in monozygotic twins which clearly demonstrate that genomics is insufficient to induce the clinical phenotype. The association of PsA with several single nucleotide polymorphisms (SNPs) at the IL23R locus and the involvement of Th17 cells in the immunopathogenesis of PsA clearly put the IL-23/IL-17 axis in the spotlight. The IL-23 and IL-17 cytokines have a pivotal role in the chronic inflammation of the synovium in PsA and are also prominent in the skin lesions of those with PsA. In this review, we focus on the genetic association of the IL-23/IL-17 axis with PsA and the contribution of these master cytokines in the pathophysiology of the disease, highlighting the main cell types incriminated in PsA and their specific role in the peripheral blood, lesional skin and joints of patients. We then provide an overview of the approved biologic drugs targeting the IL-23/IL-17 axis and discuss the advantages of genetic stratification to enhance personalized therapies in PsA.

Genomewide Association Study of Acute Anterior Uveitis Identifies New Susceptibility Loci.

Purpose: Acute anterior uveitis (AAU) is a common intraocular inflammatory disease. AAU occurs in 30% to 50% of patients with ankylosing spondylitis (AS), and both conditions are strongly associated with human leukocyte antigen (HLA)-B27, implying a shared etiology. This study aims to apply genomewide association study (GWAS) to characterize the genetic associations of AAU and their relationship to the genetics of AS. Methods: We undertook the GWAS analyses in 2752 patients with AS with AAU (cases) and 3836 patients with AS without AAU (controls). There were 7,436,415 single-nucleotide polymorphisms (SNPs) available after SNP microarray genotyping, imputation, and quality-control filtering. Results: We identified one locus associated with AAU at genomewide significance: rs9378248 (P = 2.69 × 10-8, odds ratio [OR] = 0.78), lying close to HLA-B. Suggestive association was observed at 11 additional loci, including previously reported AS loci ERAP1 (rs27529, P = 2.19 × 10-7, OR = 1.22) and NOS2 (rs2274894, P = 8.22 × 10-7, OR = 0.83). Multiple novel suggestive associations were also identified, including MERTK (rs10171979, P = 2.56 × 10-6, OR = 1.20), KIFAP3 (rs508063, P = 5.64 × 10-7, OR = 1.20), CLCN7 (rs67412457, P = 1.33 × 10-6, OR = 1.25), ACAA2 (rs9947182, P = 9.70 × 10-7, OR = 1.37), and 5 intergenic loci. The SNP-based heritability is approximately 0.5 for AS alone, and is much higher (approximately 0.7) for AS with AAU. Consistent with the high heritability, a genomewide polygenic risk score shows strong power in identifying individuals at high risk of either AS with AAU or AS alone. Conclusions: We report here the first GWAS for AAU and identify new susceptibility loci. Our findings confirm the strong overlap in etiopathogenesis of AAU with AS, and also provide new insights into the genetic basis of AAU.

Genetics in ankylosing spondylitis - Current state of the art and translation into clinical outcomes.

Ankylosing spondylitis (AS) is the prototypic form of axial spondyloarthritis (axSpA). It is highly heritable, with studies conducted in twins and in unrelated cases and controls showing that the heritability for AS is much higher than those for inflammatory bowel disease or rheumatoid arthritis. To date, 116 loci have been identified, contributing to 28% of the genetic variation in the disease. These loci provide important clues into pathogenic pathways in the disease that have led to therapeutic advances such as the repositioning of IL-17 inhibitors in the disease. Much more research is currently required to determine the functional mechanisms by which the genetic associations operate, from which it is likely that novel therapeutic approaches will be developed.

Genetics of ankylosing spondylitis--insights into pathogenesis.

Ankylosing spondylitis (AS), an immune-mediated arthritis, is the prototypic member of a group of conditions known as spondyloarthropathies that also includes reactive arthritis, psoriatic arthritis and enteropathic arthritis. Patients with these conditions share a clinical predisposition for spinal and pelvic joint dysfunction, as well as genetic associations, notably with HLA-B(*)27. Spondyloarthropathies are characterized by histopathological inflammation in entheses (regions of high mechanical stress where tendons and ligaments insert into bone) and in the subchondral bone marrow, and by abnormal osteoproliferation at involved sites. The association of AS with HLA-B(*)27, first described >40 years ago, led to hope that the cause of the disease would be rapidly established. However, even though many theories have been advanced to explain how HLA-B(*)27 is involved in AS, no consensus about the answers to this question has been reached, and no successful treatments have yet been developed that target HLA-B27 or its functional pathways. Over the past decade, rapid progress has been made in discovering further genetic associations with AS that have shed new light on the aetiopathogenesis of the disease. Some of these discoveries have driven translational ideas, such as the repurposing of therapeutics targeting the cytokines IL-12 and IL-23 and other factors downstream of this pathway. AS provides an excellent example of how hypothesis-free research can lead to major advances in understanding pathogenesis and to the development of innovative therapeutic strategies.