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

T cells in the pathogenesis of axial spondyloarthritis

Axial spondyloarthritis (axSpA) is the prototype of the spondyloarthritis spectrum. The involvement of T cells in its pathogenesis has long been suspected on the basis of the association with the major histocompatibility complex I molecule HLA-B27 and the pivotal role of interleukin 17 in the inflammatory mechanisms associated with the disease. Moreover, the presence of unconventional or "innate-like" T cells within the axial enthesis suggests an important role for these cells in the pathophysiology of the disease. In this review, we describe the characteristics and the interleukin 17 secretion capacity of the T-cell subsets identified in axSpA. We discuss the genetic and epigenetic mechanisms that support the alteration of T-cell functions and promote their activation in axSpA. We also discuss recent data on T cells that could explain the extra-articular manifestations of the SpA spectrum.

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.

Genetically encoded Runx3 and CD4+ intestinal epithelial lymphocyte deficiencies link SKG mouse and human predisposition to spondyloarthropathy

Disturbances in immune regulation, intestinal dysbiosis and inflammation characterize ankylosing spondylitis (AS), which is associated with RUNX3 loss-of-function variants. ZAP70^W163C mutant (SKG) mice have reduced ZAP70 signaling, spondyloarthritis and ileitis. In small intestine, Foxp3⁺ regulatory T cells (Treg) and CD4⁺CD8αα⁺TCRαβ⁺ intraepithelial lymphocytes (CD4-IEL) control inflammation. TGF-β and retinoic acid (RA)-producing dendritic cells and MHC-class II⁺ intestinal epithelial cells (IEC) are required for Treg and CD4-IEL differentiation from CD4⁺ conventional or Treg precursors, with upregulation of Runx3 and suppression of ThPOK. We show in SKG mouse ileum, that ZAP70^W163C or ZAP70 inhibition prevented CD4-IEL but not Treg differentiation, dysregulating Runx3 and ThPOK. TGF-β/RA-mediated CD4-IEL development, T-cell IFN-γ production, MHC class-II⁺ IEC, tissue-resident memory T-cell and Runx3-regulated genes were reduced. In AS intestine, CD4-IEL were decreased, while in AS blood CD4⁺CD8⁺ T cells were reduced and Treg increased. Thus, genetically-encoded TCR signaling dysfunction links intestinal T-cell immunodeficiency in mouse and human spondyloarthropathy.

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.

Methylation of RUNX3 Promoter 2 in the Whole Blood of Children with Ulcerative Colitis

Ulcerative colitis (UC) results from a complex interplay between the environment, gut microbiota, host genetics, and immunity. Runt-related transcription factor 3 (RUNX3) regulates Th1/Th2 balance and, thus, the synthesis of cytokines and inflammation. We aimed to analyze the dependence of RUNX3 promoter 2 (P2) methylation level on: age, sex, body mass index (BMI), C-reactive protein (CRP), serum albumin, disease duration, Pediatric Ulcerative Colitis Activity Index (PUCAI), the Paris classification, and exposure to medications. This multicenter, cross-sectional study recruited hospitalized children with UC. Methylation of RUNX3 P2 was measured with methylation-sensitive restriction enzymes in the whole blood DNA. Sixty-four children were enrolled, with a mean age of 14.5 ± 2.8 years. Half of them were female (51.6%), and the average BMI Z-score was −0.44 ± 1.14. The mean methylation of RUNX3 P2 was 54.1 ± 13.3%. The methylation level of RUNX3 P2 did not correlate with age, sex, nutritional status, CRP, albumin, PUCAI, or the extent of colitis (Paris E1–E4). RUNX3 P2 methylation did not differ between patients recruited within two and a half months of diagnosis and children who had UC for at least a year. Current or past exposure to biologics, immunosuppressants, or steroids was not associated with RUNX3 P2 methylation. Methylation of RUNX3 promoter 2 in whole blood DNA does not seem to be associated with the characteristics of UC in children.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

Monocytes and Macrophages in Spondyloarthritis: Functional Roles and Effects of Current Therapies

Spondyloarthritis (SpA) is a family of chronic inflammatory diseases, being the most prevalent ankylosing spondylitis (AS) and psoriatic arthritis (PsA). These diseases share genetic, clinical and immunological features, such as the implication of human leukocyte antigen (HLA) class I molecule 27 (HLA-B27), the inflammation of peripheral, spine and sacroiliac joints and the presence of extra-articular manifestations (psoriasis, anterior uveitis, enthesitis and inflammatory bowel disease). Monocytes and macrophages are essential cells of the innate immune system and are the first line of defence against external agents. In rheumatic diseases including SpA, the frequency and phenotypic and functional characteristics of both cell types are deregulated and are involved in the pathogenesis of these diseases. In fact, monocytes and macrophages play key roles in the inflammatory processes characteristics of SpA. The aim of this review is analysing the characteristics and functional roles of monocytes and macrophages in these diseases, as well as the impact of different current therapies on these cell types.

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.

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.

DNA Methylation Signature in Monozygotic Twins Discordant for Psoriatic Disease

Background: Psoriatic disease is a multifactorial inflammatory condition spanning from skin and nail psoriasis (Pso) to spine and joint involvement characterizing psoriatic arthritis (PsA). Monozygotic twins provide a model to investigate genetic, early life environmental exposure and stochastic influences to complex diseases, mainly mediated by epigenetics.

Methods: We performed a genome-wide DNA methylation study on whole blood of monozygotic twins from 7 pairs discordant for Pso/PsA using the Infinium Methylation EPIC array (Illumina). MeDiP—qPCR was used to confirm specific signals. Data were replicated in an independent cohort of seven patients with Pso/PsA and 3 healthy controls. Transcriptomic profiling was performed by RNAsequence on the same 7 monozygotic twin pairs.

Results: We identified 2,564 differentially methylated positions between psoriatic disease and controls, corresponding to 1,703 genes, 59% within gene bodies. There were 19 regions with at least two DMPs within 1 kb of distance and significant within-pair Δβ-values (p < 0.005), among them SNX25, BRG1 and SMAD3 genes, all involved in TGF-β signaling pathway, were identified. Co-expression analyses on transcriptome data identified IL-6/JAK/STAT3 and TNF-α pathways as important signaling axes involved in the disease, and they also suggested an altered glucose metabolism in patients’ immune cells, characteristic of pro-inflammatory T lymphocytes.

Conclusion: The study suggests the presence of an epigenetic signature in affected individuals, pointing to genes involved in immunological and inflammatory responses. This result is also supported by transcriptome data, that altogether suggest a higher activation state of the immune system, that could promote the disease status.

The interaction of Epstein-Barr virus encoded transcription factor EBNA2 with multiple sclerosis risk loci is dependent on the risk genotype

Background

Epstein-Barr virus (EBV) infection may be necessary for the development of Multiple sclerosis (MS). Earlier we had identified six MS risk loci that are co-located with binding sites for the EBV transcription factor Epstein-Barr Nuclear Antigen 2 (EBNA2) in EBV-infected B cells (lymphoblastoid cell lines – LCLs).

Methods

We used an allele-specific chromatin immunoprecipitation PCR assay to assess EBNA2 allelic preference. We treated LCLs with a peptide inhibitor of EBNA2 (EBNA2-TAT), reasoning that inhibiting EBNA2 function would alter gene expression at these loci if it was mediated by EBNA2.

Findings

We found that EBNA2 binding was dependent on the risk allele for five of these six MS risk loci (p < 0·05). Treatment with EBNA2-TAT significantly altered the expression of TRAF3 (p < 0·05), CD40 (p < 0·001), CLECL1 (p <0·0001), TNFAIP8 (p < 0·001) and TNFRSF1A (p < 0·001).

Interpretation

These data suggest that EBNA2 can enhance or reduce expression of the gene depending on the risk allele, likely promoting EBV infection. This is consistent with the concept that these MS risk loci affect MS risk through altering the response to EBNA2. Together with the extensive data indicating a pathogenic role for EBV in MS, this study supports targeting EBV and EBNA2 to reduce their effect on MS pathogenesis.

Funding

Funding was provided by grants from MS Research Australia, National Health and Medical Research Council of Australia, Australian Government Research Training Program, Multiple Sclerosis International Federation, Trish Multiple Sclerosis Research Foundation.

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.

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.

Waiting in the wings: RUNX3 reveals hidden depths of immune regulation with potential implications for inflammatory bowel disease

BACKGROUND

Complex interactions between the environment and the mucosal immune system underlie inflammatory bowel disease (IBD). The involved cytokine signalling pathways are modulated by a number of transcription factors, one of which is runt-related transcription factor 3 (RUNX3).

OBJECTIVE

To systematically review the immune roles of RUNX3 in immune regulation, with a focus on the context of IBD.

METHODS

Relevant articles and reviews were identified through a Scopus search in April 2020. Information was categorized by immune cell types, analysed and synthesized. IBD transcriptome data sets and FANTOM5 regulatory networks were processed in order to complement the literature review.

RESULTS

The available evidence on the immune roles of RUNX3 allowed for its description in twelve cell types: intraepithelial lymphocyte, Th1, Th2, Th17, Treg, double-positive T, cytotoxic T, B, dendritic, innate lymphoid, natural killer and macrophages. In the gut, the activity of RUNX3 is multifaceted and context-dependent: it may promote homeostasis or exacerbated reactions via cytokine signalling and regulation of receptor expression. RUNX3 is mostly engaged in pathways involving ThPOK, T-bet, IFN-γ, TGF-β/IL-2Rβ, GATA/CBF-β, SMAD/p300 and a number of miRNAs. RUNX3 targets relevant to IBD may include RAG1, OSM and IL-17B. Moreover, in IBD RUNX3 expression correlates positively with GZMM, and negatively with IFNAR1, whereas in controls, it strongly associates with TGFBR3.

CONCLUSIONS

Dysregulation of RUNX3, mostly in the form of deficiency, likely contributes to IBD pathogenesis. More clinical research is needed to examine RUNX3 in IBD.

Replication study and meta-analysis indicate a suggestive association of RUNX3 locus with primary biliary cholangitis

Susceptibility to primary biliary cholangitis (PBC) is in part genetically determined. In our previous PBC genome-wide association study (GWAS) in 1118 Han Chinese PBC and 4036 controls, we noted that multiple SNPs in the runt-related transcription factor 3 (RUNX3) regions showed a nominally significant association. The tag SNP rs7529070 was genotyped using a TaqMan assay in a separately collected 1435 PBC and 3205 controls. A meta-analysis with a combined 2553 PBC and 7241 controls showed that rs7529070 is still nominally associated with PBC (p = 1.7 × 10^-4, odds ratio (OR) = 1.18, 95% confidence interval (CI) = 1.08-1.28). Further analysis indicated that the risk allele of rs7529070 (G allele) is in complete linkage disequilibrium (LD) (r² = 1) with the G allele of rs4648889, which is known to be associated with increased RUNX3 expression. Bioinformatic analysis with existing expression data showed that the expression of RUNX3 is significantly increased in PBC patients (p = 0.001) and the expression level is correlated with disease severity. Consistently, we also found significantly increased RUNX3 expression (p < 0.01) in the livers of dnTGFβRII mice (a PBC mouse model). This study suggests that the RUNX3 locus may associate with PBC in Han Chinese.

The Genetics of Spondyloarthritis

The term spondyloarthritis (SpA) encompasses a group of chronic inflammatory diseases with common features in terms of clinical presentation and genetic predisposition. SpA is characterized by inflammation of the spine and peripheral joints, and is also be associated with extra-articular inflammatory manifestations such as psoriasis, uveitis, or inflammatory bowel disease (IBD). The etiology of SpA is not completely understood, but it is known to have a strong genetic component dominated by the human leukocyte antigen (HLA)-B27. In the last few years, our understanding of genetic susceptibility to SpA, particularly ankylosing spondylitis (AS), has greatly improved thanks to the findings derived from powered genome-wide association studies (GWAS) based on single nucleotide polymorphism (SNP) arrays. These studies have identified many candidate genes, therefore providing new potential directions in the exploration of disease mechanisms, especially with regard to the key role of the immune system in the pathogenesis of SpA. SpA is a complex disease where genetic variability, environmental factors, and random events interact to trigger pathological pathways. The aim of this review is to summarize current findings on the genetics of SpA, some of which might help to study new treatment approaches.

Revisiting the gut-joint axis: links between gut inflammation and spondyloarthritis

Gut inflammation is strongly associated with spondyloarthritis (SpA), as exemplified by the high prevalence of inflammatory bowel disease (IBD) and the even higher occurrence of subclinical gut inflammation in patients with SpA. The gut-joint axis of inflammation in SpA is further reinforced by similarities in immunopathogenesis at both anatomical sites and by the clinical success of therapies blocking TNF and IL-23 in IBD and in some forms of SpA. Many genetic risk factors are shared between SpA and IBD, and changes in the composition of gut microbiota are seen in both diseases. Current dogma is that inflammation in SpA initiates in the gut and leads to joint inflammation; however, although conceptually attractive, some research does not support this causal relationship. For example, therapies targeting IL-17A are efficacious in the joint but not the gut, and interfering with gut trafficking by targeting molecules such as α4β7 in IBD can lead to onset or flares of SpA. Several important knowledge gaps remain that must be addressed in future studies. Determining the true nature of the gut-joint axis has real-world implications for the treatment of patients with co-incident IBD and SpA and for the repurposing of therapeutics from one disease to the other.

Potential T Cell-Intrinsic Regulatory Roles for IRF5 via Cytokine Modulation in T Helper Subset Differentiation and Function

Interferon Regulatory Factor 5 (IRF5) is one of nine members of the IRF family of transcription factors. Although initially discovered as a key regulator of the type I interferon and pro-inflammatory cytokine arm of the innate immune response, IRF5 has now been found to also mediate pathways involved in cell growth and differentiation, apoptosis, metabolic homeostasis and tumor suppression. Hyperactivation of IRF5 has been implicated in numerous autoimmune diseases, chief among them systemic lupus erythematosus (SLE). SLE is a heterogeneous autoimmune disease in which patients often share similar characteristics in terms of autoantibody production and strong genetic risk factors, yet also possess unique disease signatures. IRF5 pathogenic alleles contribute one of the strongest risk factors for SLE disease development. Multiple models of murine lupus have shown that loss of Irf5 is protective against disease development. In an attempt to elucidate the regulatory role(s) of IRF5 in driving SLE pathogenesis, labs have begun to examine the function of IRF5 in several immune cell types, including B cells, macrophages, and dendritic cells. A somewhat untouched area of research on IRF5 is in T cells, even though Irf5 knockout mice were reported to have skewing of T cell subsets from T helper 1 (Th1) and T helper 17 (Th17) toward T helper 2 (Th2), indicating a potential role for IRF5 in T cell regulation. However, most studies attributed this T cell phenotype in Irf5 knockout mice to dysregulation of antigen presenting cell function rather than an intrinsic role for IRF5 in T cells. In this review, we offer a different interpretation of the literature. The role of IRF5 in T cells, specifically its control of T cell effector polarization and the resultant T cell-mediated cytokine production, has yet to be elucidated. A strong understanding of the regulatory role(s) of this key transcription factor in T cells is necessary for us to grasp the full picture of the complex pathogenesis of autoimmune diseases like SLE.

Analysis of 47 Non-MHC Ankylosing Spondylitis Susceptibility Loci Regarding Associated Variants across Whites and Han Chinese

OBJECTIVE

To present a systematic evaluation of 47 non-MHC ankylosing spondylitis (AS) susceptibility loci that have been initially discovered through white genome-wide association studies in Han Chinese.

METHODS

Originally, 10,743 samples representing north and south Chinese in 4 datasets were obtained. After data quality control and imputation, metaanalysis results of 94,621 variants within 47 loci were extracted. Four ERAP1 single-nucleotide polymorphisms (SNP) and HLA-B27 tag SNP rs13202464 were used for interaction analysis. Population-attributable risk percentages of AS-associated variants were compared. Functional annotations of AS-associated variants were conducted using HaploReg, RegulomeDB, and rVarBase databases.

RESULTS

We revealed 16 AS-associated variants with nominal evidence in Han Chinese, including rs10865331 (p = 6.30 × 10^-10), rs10050860 (p = 4.09 × 10^-5) and rs8070463 (p = 1.03 × 10^-4). Potential susceptible SNP within these 47 loci were also identified, such as rs13024541 (2p15), rs17401719 (5q15), and rs62074054 (17q21). Epistatic interactions between 3 ERAP1 SNP (rs17401719, rs30187, and rs10050860) and HLA-B27 were confirmed. Among the 16 AS-associated variants, rs30187 showed weaker risk effect, while rs10050860 and rs12504282 seemed to attribute more risk in Han Chinese than in whites. Further genomic annotation pinpointed 35 candidate functional SNP, especially in the 2p15, ERAP1, and NPEPPS-TBKBP1 regions.

CONCLUSION

Our results provided a detailed spectrum of all the reported non-MHC AS susceptibility loci in Han Chinese, which comprehensively exhibited the ethnic heterogeneity of AS susceptibility and highlighted that 2p15, ERAP1, and NPEPPS-TBKBP1 regions may play a critical role in AS pathogenesis across diverse populations.

RUNX3 Polymorphisms Affect the Risk of Ankylosing Spondylitis

Background

We aimed to assess the potential association of runt-related transcription factor 3 (RUNX3) gene variants with ankylosing spondylitis (AS) susceptibility among Chinese Han people.

Material/Methods

Genotyping for RUNX3 variants was accomplished through polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) in 115 AS patients and 102 healthy controls. Genotypes distributions of the polymorphisms in controls was assessed for their deviation from Hardy-Weinberg equilibrium (HWE). Moreover, odds ratio (OR) with 95% confidence interval (95%CI) was achieved using chi-square analysis to evaluate AS risk related to RUNX3 polymorphisms. Additionally, logistic regression analysis produced adjusted OR values.

Results

Genotypes distribution of rs760805 and rs11249206 polymorphisms conformed to HWE in the control group (P>0.05). TT genotype of rs760805 appeared more frequently among AS cases than in controls (P=0.033), indicating its significant association with increased risk of AS onset (OR=2.309, 95%CI=1.069–4.892). The carriage of T allele in rs760805 also heightened AS incidence, in comparison to A allele (OR=1.578, 95%CI=1.075–2.316, P=0.020). Moreover, the carriage of AT+TT genotype in rs760805 and TT genotype in rs11249206 obviously increased risk of AS onset (OR=2.585, 95%CI=1.062–6.288).

Conclusions

RUNX3 rs760805 polymorphism can contribute to AS incidence in Chinese Han people. The interaction of the 2 polymorphisms may be a risk factor for AS.

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.

A Functional Polymorphism-Mediated Disruption of EGR1/ADAM10 Pathway Confers the Risk of Sepsis Progression

Increasing evidence has indicated that single nucleotide polymorphisms (SNPs) are related to the susceptibility of sepsis and might provide potential evidence for the mechanisms of sepsis. Our recent preliminary study showed that the ADAM10 genetic polymorphism was clinically associated with the development of sepsis, and little is known about the underlying mechanism. The aim of this study was to confirm the association between the ADAM10 promoter rs653765 G→A polymorphism and the progression of sepsis and to discover the underlying mechanism. Clinical data showed that the rs653765 G→A polymorphism was positively correlated with the development of sepsis, as evidenced by a multiple-center case-control association study with a large sample size, and showed that EGR1 and ADAM10 levels were associated well with the different subtypes of sepsis patients. In vitro results demonstrated that the rs653765 G→A variants could functionally modulate ADAM10 promoter activity by altering the binding of the EGR1 transcription factor (TF) to the ADAM10 promoter, affecting the transcription and translation of the ADAM10 gene. Electrophoretic mobility shift assay (EMSA) followed by chromatin immunoprecipitation (ChIP) assay indicated the direct interaction. Functional studies further identified that the EGR1/ADAM10 pathway is important for the inflammatory response. EGR1 intervention in vivo decreased host proinflammatory cytokine secretion and rescued the survival and tissue injury of the mouse endotoxemia model.IMPORTANCE Sepsis is characterized as life-threatening organ dysfunction, with unacceptably high mortality. Evidence has indicated that functional SNPs within inflammatory genes are associated with susceptibility, progression, and prognosis of sepsis. These mechanisms on which these susceptible sites depended often suggest the key pathogenesis and potential targets in sepsis. In the present study, we confirmed that a functional variant acts as an important genetic factor that confers the progression of sepsis in a large sample size and in multiple centers and revealed that the variants modulate the EGR1/ADAM10 pathway and influence the severity of sepsis. We believe that we provide an important insight into this new pathway involving the regulation of inflammatory process of sepsis based on the clinical genetic evidence, which will enhance the understanding of nosogenesis of sepsis and provide the potential target for inflammation-related diseases.

Ankylosing spondylitis: etiology, pathogenesis, and treatments

Ankylosing spondylitis (AS), a common type of spondyloarthropathy, is a chronic inflammatory autoimmune disease that mainly affects spine joints, causing severe, chronic pain; additionally, in more advanced cases, it can cause spine fusion. Significant progress in its pathophysiology and treatment has been achieved in the last decade. Immune cells and innate cytokines have been suggested to be crucial in the pathogenesis of AS, especially human leukocyte antigen (HLA)‑B27 and the interleukin‑23/17 axis. However, the pathogenesis of AS remains unclear. The current study reviewed the etiology and pathogenesis of AS, including genome-wide association studies and cytokine pathways. This study also summarized the current pharmaceutical and surgical treatment with a discussion of future potential therapies.

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.

Genetics and Functional Genomics of Spondyloarthritis

Spondyloarthritis (SpA) is a chronic inflammatory disorder with high heritability but with complex genetics. It encompasses several entities that share common clinical features. Most of the genetic studies in SpA have been restricted to ankylosing spondylitis (AS), the prototypical form of SpA. However, there is growing evidence of shared genetic background between all the SpA subtypes and also with some other immune-mediated diseases. The most important part of SpA heritability comes from the HLA-B27 allele in the major histocompatibility complex (MHC) that explains around 25% of the attributable heredity. Several other loci outside of the MHC have been shown to be involved in the disease. However, all these non-MHC loci explain only a small additional fraction of disease predisposition. Thus, a substantial fraction of SpA genetic basis remains poorly understood. Gene expression profiling is a complementary approach to elucidate the underlying mechanisms and pathways that drive the disease. Several expression profiling studies have been undertaken in SpA. However, results have been quite disappointing with little overlap between the studies largely due to the small sample sizes, resulting in limited power to discover small effects. In this review, we summarize current knowledge on genetic findings concerning SpA and we describe strategic approaches for identification of additional variants, with a focus on rare variants in familial forms. We also provide an overview of gene expression studies in SpA and discuss the possibilities offered by high-throughput RNA sequencing technologies, in particular in sorted cells. Finally, issues in establishing molecular mechanisms underlying genetic association hits and potential translational applications will be addressed.

The Early Phases of Ankylosing Spondylitis: Emerging Insights From Clinical and Basic Science

In our paper, we discuss how the early phases of ankylosing spondylitis (AS) are linked to peri-firbocartilagenous osteitis in the sacroiliac joint and entheseal bone related anchoring sites. This skeletal proclivity is linked to an abnormal immunological response to skeletal biomechanical stress and associated microdamage. A key event in the early stages of AS appears to be the association with subclinical Crohn's-like colitis with this gut inflammation being pivotal to the osteitis reaction. Whether this osteitis is consequent to non-specific intestinal innate immune activation or adaptive immune responses against specific microbiotal or self-antigens is unknown. Recurrent iritis is an HLA-B27 associated feature that may predate AS and pursues a course independent of joint involvement, and points toward the pivotal role of organ specific immunology over generalized systemic immune responses in disease expression. Human genetics and animal model studies strongly incriminate the IL23/17 axis and TNF-α in disease pathogenesis. Preliminary work shows a strong convergence of innate immune cells including type 3 innate lymphoid-cells (ILC3) and γδ T-cells in skin, gut, entheseal, and eye inflammation. Despite the HLA-B27 association, the role of adaptive immunity, especially CD8+ T-cells mediated responses remains unproven and alternative theories have been proposed. The emerging non-dependence of axial inflammation on IL-23 but dependence on IL-17A is an unexpected new twist that awaits full explanation. In this mini-review, we discuss the key events in the early stages of human AS from clinical and basic science aspects, which could be crucial for attempted disease prevention studies in at risk subjects.

Identification of key genes in Ankylosing spondylitis

Ankylosing spondylitis (AS) is a common form of seronegative spondyloarthritis. We had identified differentially expressed genes (DEGs) in AS based on our previous RNA-sequencing results. Our study aimed to identify key genes in AS by integrated microarray analysis. In this present study, we identified 1328 DEGs between AS and normal control by using integrated analysis of two datasets derived from the Gene Expression Omnibus (GEO) database. Functional annotation of DEGs were performed. Pathways in cancer, Pancreatic cancer and Natural killer cell mediated cytotoxicity were significantly enriched pathways for DEGs. Based on the shared DEGs of AS in both integrated analysis and our previous RNA-sequencing results, we constructed the protein and protein interaction (PPI) network. BIRC2, MAPILC3A and MAGED1 were hub proteins. Validation of gene expression by qRT-PCR were performed and the results were consistent with our integrated analysis generally. Our finding provided new clues for understanding the mechanism of AS.

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.

ERAP1/ERAP2 and RUNX3 polymorphisms are not associated with ankylosing spondylitis susceptibility in Chinese Han

Previous studies show that endoplasmic reticulum‐associated aminopeptidase (ERAP1/ERAP2) and runt‐related transcription factor 3 (RUNX3) gene polymorphisms are associated with AS (ankylosing spondylitis) in European Caucasians. However, contradictory results were reported in different Asian populations. The purpose of this study was to determine whether eleven candidate single nucleotide polymorphisms (SNPs) in ERAP1/ERAP2 and six in RUNX3 genes confer susceptibility to AS with or without acute anterior uveitis (AAU) [AS⁺AAU⁺ or AS⁺AAU^–] in Chinese Han. Therefore, a case–control association study was performed in 882 AS⁺AAU^–, 884 AS⁺AAU⁺ and 1727 healthy controls. Genotyping was performed using the iPLEXGold genotyping assay. A meta‐analysis was performed to assess the association of polymorphisms of ERAP1 with AS susceptibility in Asian populations. No association was found between SNPs of ERAP1/ERAP2/RUNX3 and susceptibility of AS with or without AAU. A case–control study between patients with human leucocyte antigen HLA‐B27‐positive and healthy controls also failed to demonstrate an association of the tested SNP with AS with or without AAU. Moreover, a meta‐analysis showed that there was no association of rs30187, rs27037, rs27980, rs27434 and rs27582 in ERAP1 with AS in Chinese Han. Taken together, 17 SNPs in ERAP1/ERAP2 and RUNX3 genes did not confer disease susceptibility to AS in Chinese Han.

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.

Association Between a CCL17 Genetic Variant and Risk of Coronary Artery Disease in a Chinese Han Population

BACKGROUND

In the present study we investigated the effects of genetic variations in the C-C motif chemokine ligand 17 (CCL17) gene on serum CCL17 levels and risk of coronary artery disease (CAD).

METHODS AND RESULTS

A case-control study was conducted to determine causal inferences amongCCL17single-nucleotide polymorphisms (SNPs), serum CCL17 levels, and risk of CAD. Luciferase assays, electrophoretic mobility shift assays (EMSA), and allele-specific quantitative chromatin immunoprecipitation (ChIP) assays were used to assess the function of the SNPs. In all, 947 participants (794 with CAD, 153 without CAD) were included in the study. The T allele in rs223828, located in intron of theCCL17gene, was associated with increased serum CCL17 levels as well as increased CAD risk. A causal inference test using mediation analysis suggested that rs223828 had a significant indirect casual effect on the increased risk of CAD mediated via serum CCL17 levels. Luciferase assays confirmed that the rs223828T allele enhancesCCL17promoter activity. Protein-DNA binding studies using EMSA and allele-specific quantitative ChIP assays indicated preferential activator protein-1 (AP-1) complex formation and recruitment with the rs223828 T allele compared with the C allele.

CONCLUSIONS

We propose that theCCL17SNP rs223828 is associated with increased risk of CAD, and that this site may be a potential AP-1 binding site.

The role of histone modification and a regulatory single-nucleotide polymorphism (rs2071166) in the Cx43 promoter in patients with TOF

Abnormal level of Cx43 expression could result in CHD. Epigenetic modification and disease-associated, non-coding SNPs might influence gene transcription and expression. Our study aimed to determine the role of histone modification and an rSNP (rs2071166) in the Cx43 promoter in patients with TOF. Our results indicate that H3K18ac bind to Cx43 promoter and that their levels are reduced in TOF patients relative to controls. The relationship between the non-coding SNP in the Cx43 gene and TOF patients was evaluated in 158 patients and 300 controls. The C allele of rs2071166 was confirmed to result in an increased risk of TOF (OR = 1.586, 95%CI 1.149–2.189). Individuals with the CC genotype at rs2071166 also showed a significant susceptibility to TOF (OR = 2.961, 95%CI 1.452–6.038). The mRNA level in TOF who were CC genotype was lower than that in patients with the AA/AC genotype. Functional analysis in cells and transgenic zebrafish models showed that rs2071166 decreased the activity of the promoter and could block the interaction between RXRα and RARE. This is the first study to illustrate that epigenetic modification and an rSNP in the Cx43 promoter region play a critical role in TOF by impacting the transcriptional activity and expression level of Cx43.

Genetics and the Causes of Ankylosing Spondylitis

Ankylosing spondylitis (AS) is a common inflammatory arthritis in which genetic factors are the primary determinants of disease risk and severity. Substantial progress has been made in identifying genetic pathways involved in the disease, and in translating those discoveries to drug discovery programs. Recently discovered novel disease pathways include those involved in control of DNA methylation, bacterial sensing, and mucosal immunity. Additional pathways are likely to be identified as a higher proportion of the genetic risk of AS is determined.

A Variant in RUNX3 Is Associated with the Risk of Ankylosing Spondylitis in Koreans

Ankylosing spondylitis (AS) is a chronic autoinflammatory disease that affects the spine and sacroiliac joints. Regarding its etiology, although HLA-B27 is known to be the strongest genetic factor of AS, much evidence suggests the potential contribution of non-MHC genes to the susceptibility to AS. Most of these non-MHC genes have been discovered in non-Asian populations; however, just some of them have been validated in Koreans. In this study, we aimed to identify additional AS-associated single-nucleotide polymorphism (SNP) candidates by replicating the candidate SNPs in Korean AS patients and healthy controls. For this, we selected three SNPs (rs11249215 in RUNX3, rs6556416 in IL12B, and rs8070463 in TBKBP1), which were previously reported as risk factors of AS but have not been studied in Koreans, and performed genotyping assays using a total of 1138 Korean samples (572 AS patients and 566 healthy controls). Of the three SNP candidates, one SNP in RUNX3 (rs11249215) was significantly associated with the risk of AS (odds ratio, 1.31; 95% confidence interval, 1.02 to 1.68, p = 0.03). These results will be helpful in elucidating the pathogenesis of AS and may be useful for developing AS risk prediction models in Koreans.

Runx3 in Immunity, Inflammation and Cancer

In this chapter we summarize the pros and cons of the notion that Runx3 is a major tumor suppressor gene (TSG). Inactivation of TSGs in normal cells provides a viability/growth advantage that contributes cell-autonomously to cancer. More than a decade ago it was suggested that RUNX3 is involved in gastric cancer development, a postulate extended later to other epithelial cancers portraying RUNX3 as a major TSG. However, evidence that Runx3 is not expressed in normal gastric and other epithelia has challenged the RUNX3-TSG paradigm. In contrast, RUNX3 is overexpressed in a significant fraction of tumor cells in various human epithelial cancers and its overexpression in pancreatic cancer cells promotes their migration, anchorage-independent growth and metastatic potential. Moreover, recent high-throughput quantitative genome-wide studies on thousands of human samples of various tumors and new investigations of the role of Runx3 in mouse cancer models have unequivocally demonstrated that RUNX3 is not a bona fide cell-autonomous TSG. Importantly, accumulating data demonstrated that RUNX3 functions in control of immunity and inflammation, thereby indirectly influencing epithelial tumor development.