Functional Genomic Analysis of a RUNX3 Polymorphism Associated With Ankylosing Spondylitis

METTL14-m6A-FOXO3a axis regulates autophagy and inflammation in ankylosing spondylitis

The role of m6A in ankylosing spondylitis (AS) remains largely obscure. In this study, we found that m6A modification was decreased in T cells of AS, and the abnormal m6A modification was attributed to the downregulation of methyltransferase-like 14 (METTL14). METTL14 exerted a critical role in regulating autophagy activity and inflammation via targeting Forkhead box O3a (FOXO3a). Mechanistically, the loss of METTL14 decreased the expression of FOXO3a, leading to the damage of autophagic flux and the aggravation of inflammation. Inversely, the forced expression of METTL14 upregulated the expression of FOXO3a, thereby activating autophagy and alleviating inflammation. Furthermore, our results revealed that METTL14 targeted FOXO3a mRNA and regulated its expression and stability in a m6A-dependent manner. These findings uncovered the functional importance of m6A methylation mechanisms in the regulation of autophagy and inflammation, which expanded our understanding of this interaction and was critical for the development of therapeutic strategies for AS.

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.

CHD4 orchestrates the symphony of T and B lymphocytes development and a good mediator in preventing from autoimmune disease

Chromodomain helicase DNA binding protein 4 (CHD4) is an ATPase subunit of the nucleosome remodeling and deacetylation complex. It has been implicated in gene transcription, DNA damage repair, maintenance of genome stability, and chromatin assembly. Meanwhile, it is highly related to cell cycle progression and the proceeding of malignancy. Most of the previous studies were focused on the function of CHD4 with tumor cells, cancer stem cells, and cancer cells multidrug resistance. Recently, some researchers have explored the CHD4 functions on the development and differentiation of adaptive immune cells, such as T and B lymphocytes. In this review, we will discuss details of CHD4 in lymphocyte differentiation and development, as well as the critical role of CHD4 in the pathogenesis of the autoimmune disease.

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.

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.

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.