Epigenetic and gene expression analysis of ankylosing spondylitis-associated loci implicate immune cells and the gut in the disease pathogenesis

Functional genomics implicates natural killer cells in the pathogenesis of ankylosing spondylitis

Objective

Multiple lines of evidence indicate that ankylosing spondylitis (AS) is a lymphocyte-driven disease. However, which lymphocyte populations are critical in AS pathogenesis is not known. In this study, we aimed to identify the key cell types mediating the genetic risk in AS using an unbiased functional genomics approach.

Methods

We integrated genome-wide association study (GWAS) data with epigenomic and transcriptomic datasets of human immune cells. To quantify enrichment of cell type-specific open chromatin or gene expression in AS risk loci, we used three published methods that have successfully identified relevant cell types in other diseases. We performed co-localization analyses between GWAS risk loci and genetic variants associated with gene expression (eQTL) to find putative target genes.

Results

Natural killer (NK) cell-specific open chromatin regions are significantly enriched in heritability for AS, compared to other immune cell types such as T cells, B cells, and monocytes. This finding was consistent between two AS GWAS. Using RNA-seq data, we validated that genes in AS risk loci are enriched in NK cell-specific gene expression. Using the human Space-Time Gut Cell Atlas, we also found significant upregulation of AS-associated genes predominantly in NK cells. Co-localization analysis revealed four AS risk loci affecting regulation of candidate target genes in NK cells: two known loci, ERAP1 and TNFRSF1A, and two under-studied loci, ENTR1 (aka SDCCAG3) and B3GNT2.

Conclusion

Our findings suggest that NK cells may play a crucial role in AS development and highlight four putative target genes for functional follow-up in NK cells.

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.

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⁻¹⁵) 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.

Single-cell RNA-seq reveals altered NK cell subsets and reduced levels of cytotoxic molecules in patients with ankylosing spondylitis

Ankylosing spondylitis (AS) is an autoimmune disease with unknown aetiology. To unravel the mechanisms mediating AS pathogenesis, we profiled peripheral blood mononuclear cells (PBMCs) from AS patients and healthy subjects using 10X single-cell RNA sequencing. The frequencies of immune cell subsets were evaluated by flow cytometry. NK cells were purified from PBMCs using isolation kit and were examined for gene expression by RT-qPCR. Plasma levels of cytolytic molecules were examined by enzyme-linked immunosorbent assay. Compared to healthy controls, AS patients showed a significant decrease in total NK cells as well as CD56^dim NK subset, whereas CD56^bright NK cells were increased. Additionally, impaired expression of cytotoxic genes in NK cells of AS patients was observed by bioinformatics algorithm and verified by RT-qPCR and flow cytometry. Consistent with changes in transcriptomics, we found decreased plasma levels of granzymes, but not granulysin, in AS patients. Furthermore, Pearson correlation analysis revealed a negative correlation between plasma GZMB levels and disease activity (r = -0.5275, p = 0.0358). No correlation was observed between plasma cytolytic molecules and biochemical indexes (ESR and CRP). Our findings uncover altered NK cell subsets and cytotoxic profiles in peripheral circulation of AS patients at single-cell resolution.

Differential expression of circular RNAs in plasma exosomes from patients with ankylosing spondylitis

Exosomes are being extensively studied as a source of valuable new biomarkers. The underlying mechanism of ankylosing spondylitis (AS) may include changes in the circular RNAs (circRNAs) of exosomes. However, there is a lack of reports on the role of exosomal cirRNAs in the plasma of patients with AS. We isolated exosomes from the plasma of patients with AS and healthy individuals (controls). Subsequently, we investigated the circRNA profiles of the exosomes via high-throughput RNA sequencing and identified 56 differentially expressed circRNAs in the exosomes of patients with AS compared with those of the healthy controls. Gene Ontology (GO) demonstrated that the differentially expressed circRNAs were mainly involved in the negative regulation of the activity of the transcription factor NF-κB and bone remodelling that is potentially related to AS. Kyoto Encyclopedia Genes and Genome (KEGG) demonstrated that the most highly AS-correlated pathways that were identified were "notch signalling pathway" and those primarily involved with "cholinergic synapse". Finally, we validated five differentially expressed circRNAs using quantitative real time PCR (qRT-PCR) and predicted their potential functions through the circRNA-miRNA-mRNA network.Our study is the first to report changes in exosomal circRNAs from plasma samples of patients with AS, and provides novel targets for further investigation of molecular mechanisms and potential intervention therapy targets of AS. This article is protected by copyright. All rights reserved.

B Cell Involvement in the Pathogenesis of Ankylosing Spondylitis

Extensive research into ankylosing spondylitis (AS) has suggested the major role of genetics, immune reactions, and the joint-gut axis in its etiology, although an ultimate consensus does not yet exist. The available evidence indicates that both autoinflammation and T-cell-mediated autoimmune processes are actively involved in the disease process of AS. So far, B cells have received relatively little attention in AS pathogenesis; this is largely due to a lack of conventional disease-defining autoantibodies. However, against prevailing dogma, there is a growing body of evidence suggestive of B cell involvement. This is illustrated by disturbances in circulating B cell populations and the formation of auto-reactive and non-autoreactive antibodies, along with B cell infiltrates within the axial skeleton of AS patients. Furthermore, the depletion of B cells, using rituximab, displayed beneficial results in a subgroup of patients with AS. This review provides an overview of our current knowledge of B cells in AS, and discusses their potential role in its pathogenesis. An overarching picture portrays increased B cell activation in AS, although it is unclear whether B cells directly affect pathogenesis, or are merely bystanders in the disease process.

Dynamics of Adaptive Immune Cell and NK Cell Subsets in Patients With Ankylosing Spondylitis After IL-17A Inhibition by Secukinumab

Background: Anti-IL-17A therapy is generally effectively applied in patients with Ankylosing Spondylitis (AS) to achieve and maintain remission. However, the influence of anti-IL-17A on the composition of the immune system is not apparent. Our prospective study was to explore the changes in immune imbalance regarding T cell, B cell and natural killer (NK) cell subsets after secukinumab treatment in AS patients. Methods: Immune cell distribution of 43 AS patients treated with secukinumab for 12 weeks and 47 healthy controls (HC) were evaluated. Flow cytometry using monoclonal antibodies against 25 surface markers was accomplished to explore the frequencies of lineage subsets. The differences between HC, AS pre-treatment, and post-treatment were compared using the paired Wilcoxon test, Mann-Whitney U test, and ANOVA. Results: AS patients had altered immune cell distribution regarding T cell and B cell subsets. Apart from activated differentiation of CD4⁺ T cell, CD8⁺ T cell and B cell, higher levels of cytotoxic T (Tc) two cells and Tc17 cells were noted in AS patients. We confirmed that helper T (Th) one cell became decreased; however, Th17 cells and T follicular helper (Tfh) 17 cells went increased in AS. After 12 weeks of secukinumab therapy, CRP and ASDAS became significantly decreased, and meanwhile, the proportions of Th1 cells, Tfh17 cells and classic switched B cells were changed towards those of HC. A decreased CRP was positively correlated with a decrease in the frequency of naïve CD8⁺ T cells (p = 0.039) and B cells (p = 0.007) after secukinumab treatment. An elevated level of T cells at baseline was detected in patients who had a good response to secukinumab (p = 0.005). Conclusion: Our study confirmed that AS patients had significant multiple immune cell dysregulation. Anti-IL-17A therapy (Secukinumab) could reverse partial immune cell imbalance.

Dysregulation of ribosome-related genes in ankylosing spondylitis: a systems biology approach and experimental method

Background

Ankylosing spondylitis (AS) is an autoimmune rheumatic disease. Few candidate gene associations have been reported for AS and the current understanding of its pathogenesis remains still poor. Thus, the exact mechanism of AS is needed to urgently be disclosed. The purpose of this study was to identify candidate genes involving in AS disease.

Methods and results

GSE25101 publicly available microarray and GSE117769 RNA-seq datasets of AS patients were obtained for bioinformatics analyses. Gene set enrichment analysis showed that in the microarray dataset, the ribosome pathway was significantly up-regulated in AS compared with controls. Furthermore, some ribosomal components demonstrated overexpression in patients in the RNA-seq dataset. To confirm the findings, 20 AS patients and 20 matching controls were selected from the Rheumatology Research Center clinic, Shariati Hospital. PBMCs were separated from whole blood and RNA contents were extracted. Following the results of datasets analysis, the expression level of rRNA5.8S pseudogene, rRNA18S pseudogene, RPL23, RPL7, and RPL17 genes were measured through real-time PCR. Our findings showed dysregulation of rRNA5.8S and rRNA18S pseudogenes, and also the RPL17 gene in patients.

Conclusion

Considering that genes involved in ribosome biogenesis contributed to some AS-associated biological processes as well as diseases that have comorbidities with AS, our results might advance our understanding of the pathological mechanisms of ankylosing spondylitis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12891-021-04662-2.

TNF-α inhibitor therapy can improve the immune imbalance of CD4+ T cells and negative regulatory cells but not CD8+ T cells in ankylosing spondylitis

BACKGROUND

Studies into ankylosing spondylitis (AS) and its relationship with immune imbalance are controversial, and the correlation between the efficacy of TNF-α inhibitor and changes in immune imbalance is unclear.

METHODS

A total of 40 immune cells were tested with flow cytometry, and the results of 105 healthy control (HC) subjects, 177 active-stage AS patients, and 23 AS cases before and after 12 weeks of TNF-α inhibitor therapy (Anbainuo) were analyzed.

RESULTS

Compared with the HC group, the proportion of immune cells, such as naïve and central memory CD4+T cells, in AS increased (P < 0.0001), but effector memory and terminally differentiated CD4+T cells were decreased (P < 0.01 and 0.0001, respectively). Naïve, central memory, and effector memory CD8+T cells were increased (P < 0.0001, 0.001, and 0.01, respectively), but terminally differentiated CD8+T cells were decreased (P < 0.0001). Th1 cells (helper T cells-1), Tfh1 cells (follicular helper T cells-1), Tc1 cells (cytotoxic T cells-1), and Tregs (regulatory T cells) were lower (P < 0.01, 0.05, 0.0001, and 0.001, respectively), but Th17 cells, Tfh17 cells, and Tc cells were higher (P < 0.001, 0.0001, and 0.001, respectively). The proportions of total B cells and class-switched B cells were increased (P < 0.05), but non-switched B cells, plasma cells, memory B cells, and immature Bregs (regulatory B cells) were lower (P < 0.01, 0.0001, 0.0001, and 0.0001, respectively). After Anbainuo therapy, the percentage of naïve CD4+ T cells had decreased (P < 0.05) but Tregs and B10 cells (IL-10-producing regulatory B cells) had increased (P < 0.01 and 0.05, respectively), and the increase in Tregs was positively correlated with the decrease in C-reactive protein (CRP) (r = 0.489, P = 0.018).

CONCLUSIONS

We found that active-stage AS patients have an immunity imbalance of frequency involving multiple types of immune cells, including CD4+T cells, CD8+T cells, Th cells, Tfh cells, Tc cells, Tregs, Bregs, and B cells. TNF-α inhibitor Anbainuo can not only help to inhibit disease activity but can also improve the immune imbalance of CD4+ T cells and negative regulatory cells in frequency. But CD8+ T cells have not been rescued.

HLA associations in inflammatory arthritis: emerging mechanisms and clinical implications

Our understanding of the mechanisms underlying HLA associations with inflammatory arthritis continues to evolve. Disease associations have been refined, and interactions of HLA genotype with other genes and environmental risk factors in determining disease risk have been identified. This Review provides basic information on the genetics and molecular function of HLA molecules, as well as general features of HLA associations with disease. Evidence is discussed regarding the various peptide-dependent and peptide-independent mechanisms by which HLA alleles might contribute to the pathogenesis of three types of inflammatory arthritis: rheumatoid arthritis, spondyloarthritis and systemic juvenile idiopathic arthritis. Also discussed are HLA allelic associations that shed light on the genetic heterogeneity of inflammatory arthritides and on the relationships between adult and paediatric forms of arthritis. Clinical implications range from improved diagnosis and outcome prediction to the possibility of using HLA associations in developing personalized strategies for the treatment and prevention of these diseases.

Shotgun metagenomics reveals an enrichment of potentially cross-reactive bacterial epitopes in ankylosing spondylitis patients, as well as the effects of TNFi therapy upon microbiome composition

Objectives

Diverse evidence including clinical, genetic and microbiome studies support a major role of the gut microbiome in the common immune-mediated arthropathy, ankylosing spondylitis (AS). We set out to (1) further define the key microbial characteristics driving disease, and (2) examine the effects of tumour necrosis factor-inhibitor (TNFi) therapy upon the microbiome.

Methods

The stools from a case–control cohort of 250 Han-Chinese subjects underwent shotgun metagenomic sequencing. All subjects were genotyped using the Illumina CoreExome SNP microarray.

Results

Previous reports of gut dysbiosis in AS were reconfirmed and several notable bacterial species and functional categories were differentially abundant. TNFi therapy was correlated with a restoration the perturbed microbiome observed in untreated AS cases to that of healthy controls, including several important bacterial species that have been previously associated with AS and other related diseases. Enrichment of bacterial peptides homologous to HLA-B27-presented epitopes was observed in the stools of patients with AS, suggesting that either HLA-B27 fails to clear these or that they are involved in driving HLA-B27-associated immune reactions. TNFi therapy largely restored the perturbed microbiome observed in untreated AS cases to that of healthy controls, including several important bacterial species that have been previously associated with AS and other related diseases. TNFi therapy of patients with AS was also associated with a reduction of potentially arthritogenic bacterial peptides, relative to untreated patients.

Conclusion

These findings emphasise the key role that the gut microbiome plays in driving the pathogenesis of AS and highlight potential therapeutic and/or preventative targets.

IL-17 in the immunopathogenesis of spondyloarthritis

Spondyloarthritis (SpA) is a term that refers to a group of inflammatory diseases that includes psoriatic arthritis, axial SpA and nonradiographic axial SpA, reactive arthritis, enteropathic arthritis and undifferentiated SpA. The disease subtypes share clinical and immunological features, including joint inflammation (peripheral and axial skeleton); skin, gut and eye manifestations; and the absence of diagnostic autoantibodies (seronegative). The diseases also share genetic factors. The aetiology of SpA is still the subject of research by many groups worldwide. Evidence from genetic, experimental and clinical studies has accumulated to indicate a clear role for the IL-17 pathway in the pathogenesis of SpA. The IL-17 family consists of IL-17A, IL-17B, IL-17C, IL-17D, IL-17E and IL-17F, of which IL-17A is the best studied. IL-17A is a pro-inflammatory cytokine that also has the capacity to promote angiogenesis and osteoclastogenesis. Of the six family members, IL-17A has the strongest homology with IL-17F. In this Review, we discuss how IL-17A and IL-17F and their cellular sources might contribute to the immunopathology of SpA.

Genome-wide DNA methylation analysis in ankylosing spondylitis identifies HLA-B*27 dependent and independent DNA methylation changes in whole blood

BACKGROUND AND OBJECTIVE

Ankylosing spondylitis is a chronic inflammatory disease characterized by inflammation of the sacroiliac joints and the spine that can lead to significant pain, immobility, and disability. The etiology and pathogenesis of ankylosing spondylitis are incompletely understood, though most patients carry the HLA-B*27 allele. The objective of this study was to evaluate DNA methylation changes in ankylosing spondylitis with the goal of revealing novel mechanistic insights into this disease.

METHODS

Genome-wide DNA methylation analysis was performed in whole blood DNA samples using the Infinium MethylationEPIC array in patients with ankylosing spondylitis compared to age, sex, and race matched patients with osteoarthritis as a non-inflammatory disease control. We studied 24 patients with ankylosing spondylitis, including 12 patients who carry HLA-B*27 and 12 patients who are HLA-B*27 negative. DNA methylation analysis was performed with adjustment for blood cell composition in each sample.

RESULTS

We identified a total of 67 differentially methylated sites between ankylosing spondylitis patients and osteoarthritis controls. Hypermethylated genes found included GTPase-related genes, while hypomethylated genes included HCP5, which encodes a lncRNA within the MHC region, previously associated with genetic risk for psoriasis and toxic epidermal necrolysis. Carrying HLA-B*27 was associated with robust hypomethylation of HCP5, tubulin folding cofactor A (TBCA) and phospholipase D Family Member 6 (PLD6) in ankylosing spondylitis patients. Hypomethylation within HCP5 involves a CpG site that contains a single nucleotide polymorphism in linkage disequilibrium with HLA-B*27 and that controls DNA methylation at this locus in an allele-specific manner.

CONCLUSIONS

A genome-wide DNA methylation analysis in ankylosing spondylitis identified DNA methylation patterns that could provide potential novel insights into this disease. Our findings suggest that HLA-B*27 might play a role in ankylosing spondylitis in part through inducing epigenetic dysregulation.

RUNX3 and T-Bet in Immunopathogenesis of Ankylosing Spondylitis-Novel Targets for Therapy?

Susceptibility to ankylosing spondylitis (AS) is polygenic with more than 100 genes identified to date. These include HLA-B27 and the aminopeptidases (ERAP1, ERAP2, and LNPEPS), which are involved in antigen processing and presentation to T-cells, and several genes (IL23R, IL6R, STAT3, JAK2, IL1R1/2, IL12B, and IL7R) involved in IL23 driven pathways of inflammation. AS is also strongly associated with polymorphisms in two transcription factors, RUNX3 and T-bet (encoded by TBX21), which are important in T-cell development and function. The influence of these genes on the pathogenesis of AS and their potential for identifying drug targets is discussed here.

Association of High-Sensitivity Troponin T With Left Ventricular Dysfunction in Ankylosing Spondylitis

BACKGROUND

Ankylosing spondylitis (AS) is a systemic inflammatory disease, and cardiac dysfunction has not been clearly described clinically. High-sensitivity cardiac troponin T (hs-cTnT) is a noninvasive marker for subclinical myocardial injury.

OBJECTIVE

In this study, we aimed to investigate any relationship between hs-cTnT and left ventricular (LV) function evaluated via tissue Doppler imaging in AS patients with no known cardiac risk factor.

METHODS

Our study used a cross-sectional case protocol design and was conducted between January 2016 and June 2016. In total, 40 AS patients (17 females and 23 males) were age and sex matched with healthy volunteers (20 females and 20 males) and enlisted for this study. Detailed transthoracic echocardiography was performed, and tissue Doppler imaging was used to assess systolic and diastolic functions. High-sensitivity cardiac troponin T levels were measured and compared between 2 groups.

RESULTS

Compared with control subjects, AS patients had lower early (Em)/late (Am) diastolic myocardial velocities, mitral annular plane systolic excursion, and end-diastolic distance from the mitral annulus to the LV apex. Conversely, they had greater systolic myocardial velocity (Sm), isovolumetric relaxation time, and displacement index (p < 0.001, for all). Higher hs-cTnT levels were measured in AS patients (0.45 ± 0.22 vs. 1.11 ± 0.27, p < 0.001), and multivariate logistic regression analyses revealed that hs-cTnT was an independent predictor of LV diastolic dysfunction in AS patients.

CONCLUSIONS

These data show that AS patients had impaired LV functions and increased hs-cTnT levels. Tissue Doppler imaging may be a useful tool for detection of early functional LV abnormalities, and hs-cTnT may be valuable biomarker of diastolic LV dysfunction in AS patients.

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.

Evidence for a second ankylosing spondylitis-associated RUNX3 regulatory polymorphism

Objectives

To explore the functions of RUNX3 single nucleotide polymorphisms (SNPs) associated with ankylosing spondylitis (AS).

Methods

Individual SNP associations were evaluated in 4230 UK cases. Their effects on transcription factor (TF) binding, transcription regulation, chromatin modifications, gene expression and gene interactions were tested by database interrogation, luciferase reporter assays, electrophoretic mobility gel shifts, chromatin immunoprecipitation and real-time PCR.

Results

We confirmed the independent association of AS with rs4265380, which was robust (P=4.7×10⁻⁶) to conditioning on another nearby AS-associated RUNX3 SNP (rs4648889). A RUNX3 haplotype incorporating both SNPs was strongly associated with AS (OR 6.2; 95% CI 3.1 to 13.2, P=1.4×10⁻⁸). In a large UK cohort, rs4265380 is associated with leucocyte counts (including monocytes). RUNX3 expression is lower in AS peripheral blood mononuclear cells than healthy controls (P<0.002), independent of rs4265380 genotype. Enhancer function for this RUNX3 region was suggested by increased luciferase activity (approximately tenfold; P=0.005) for reporter constructs containing rs4265380. In monocytes, there was differential allelic binding of nuclear protein extracts to a 50 bp DNA probe containing rs4265380 that was strongly augmented by lipopolysaccharide activation. TF binding also included the histone modifier p300. There was enrichment for histone modifications associated with active enhancer elements (H3K27Ac and H3K79Me2) that may be allele dependent. Hi-C database interrogation showed chromosome interactions of RUNX3 bait with the nearby RP4-799D16.1 lincRNA.

Conclusions

The association of AS with this RUNX3 regulatory region involves at least two SNPs apparently operating in different cell types. Monocytes may be potential therapeutic targets in AS.