Abstract:Introduction: GWAS have identified multiple regions that confer risk for juvenile idiopathic arthritis (JIA). However, identifying the single nucleotide polymorphisms (SNPs) that drive disease risk is impeded by the SNPs' that identify risk loci being in linkage disequilibrium (LD) with hundreds of other SNPs. Since the causal SNPs remain unknown, it is difficult to identify target genes and use genetic information to inform patient care. We used genotyping and functional data in primary human monocytes/macrophages to nominate disease-driving SNPs on JIA risk haplotypes and identify their likely target genes. Methods: We identified JIA risk haplotypes using Immunochip data from Hinks et al (Nature Gen 2013) and the meta-analysis from McIntosh et al (Arthritis Rheum 2017). We used genotyping data from 3,939 children with JIA and 14,412 healthy controls to identify SNPs that: (1) were situated within open chromatin in multiple immune cell types and (2) were more common in children with JIA than the controls (p< 0.05). We intersected the chosen SNPs (n=846) with regions of bi-directional transcription initiation characteristic of non-coding regulatory regions detected using dREG to analyze GRO-seq data. Finally, we used MicroC data to identify gene promoters interacting with the regulatory regions harboring the candidate causal SNPs. Results: We identified 190 SNPs that overlap with dREG peaks in monocytes and126 SNPs that overlap with dREG peaks in macrophages. Of these SNPs, 101 were situated within dREG peaks in both monocytes and macrophages, suggesting that these SNPs exert their effects independent of the cellular activation state. MicroC data in monocytes identified 20 genes/transcripts whose promoters interact with the enhancers harboring the SNPs of interest. Conclusion: SNPs in JIA risk regions that are candidate causal variants can be further screened using functional data such as GRO-seq. This process identifies a finite number of candidate causal SNPs, the majority of which are likely to exert their biological effects independent of cellular activation state in monocytes. Three-dimensional chromatin data generated with MicroC identifies genes likely to be influenced by these SNPs. These studies demonstrate the importance of investigations into the role of innate immunity in JIA.

Complex regulatory interactions at GDF5 shape joint morphology and osteoarthritis disease risk

Age-dependent genetic regulation of osteoarthritis: independent effects of immune system genes

A comprehensive meta-analysis of association between genetic variants of GDF5 and osteoarthritis of the knee, hip and hand

Large-scale Analysis of Association Between GDF5 and FRZB Variants and Osteoarthritis of the Hip, Knee, and Hand.

Genetics of osteoarthritis

A meta-analysis of European and Asian cohorts reveals a global role of a functional SNP in the 5' UTR of GDF5 with osteoarthritis susceptibility.

A common cis- regulatory variant impacts normal-range and disease-associated human facial shape through regulation of PKDCC during chondrogenesis

A functional polymorphism in the 5′ UTR of GDF5 is associated with susceptibility to osteoarthritis

Genetic and Epigenetic Fine‐Tuning of TGFB1 Expression Within the Human Osteoarthritic Joint

Comprehensive profiling of transcriptional regulation in cartilage reveals pathogenesis of osteoarthritis

Systematic Localization of Common Disease-Associated Variation in Regulatory DNA

Osteoarthritis as an Enhanceropathy: Gene Regulation in Complex Musculoskeletal Disease

Juvenile idiopathic arthritis-associated genetic loci exhibit spatially constrained gene regulatory effects across multiple tissues and immune cell types

Response eQTLs, chromatin accessibility, and 3D chromatin structure in chondrocytes provide mechanistic insight into osteoarthritis risk

Variant to gene mapping for carpal tunnel syndrome risk loci implicates skeletal muscle regulatory elements

Interplay between genetics and epigenetics in osteoarthritis

Regulatory elements in SEM1-DLX5-DLX6 (7q21.3) locus contribute to genetic control of coronal nonsyndromic craniosynostosis and bone density-related traits

Genome-wide enhancer maps link risk variants to disease genes

Regulatory Element-Based Prediction Identifies New Susceptibility Regulatory Variants for Osteoporosis.

Using Functional Genomic Data in Monocytes/Macrophages and Genotyping to Nominate Disease-Driving Single Nucleotide Polymorphisms and Target Genes in Juvenile Idiopathic Arthritis