University of California San Diego San Diego, CA, United States
Disclosure information not submitted.
Eunice Choi1, Camilla R. L. Machado2, David Boyle1, Wei Wang1 and Gary Firestein2, 1University of California San Diego, La Jolla, CA, 2University of California San Diego, San Diego, CA
Background/Purpose: Mechanisms responsible for the distribution and severity of joint involvement in RA are not known. To explore whether site-specific FLS biology might contribute to location-specific synovitis and explain the predilection for hand (wrist/metacarpal phalangeal joints) in RA, we generated transcriptomic and chromatin accessibility data from FLS to define the transcription factors (TFs) and pathways unique to each location.
Methods: FLS were derived from hand, knee and hip surgical samples. Whole genome ATAC-seq and RNA-seq were obtained from 10 hand, 10 hip, and 10 knee FLS under unstimulated (unstim) or stimulated conditions (TNF [50 ng/ml for 6 h]). Differentially expressed genes (DEGs) were identified using the Wilcoxon test. TF motifs were queried within promoter/enhancer open chromatin. The assays were integrated using Taiji (Nat Commun 2022;13:6221) and Personalized PageRanks were quantified. Cell growth was assayed using MTT in medium and PDGF-stimulated FLS on day 7.
Results: Multiple DEGs between joints were identified in their respective FLS (Table 1A). Hand vs knee or hand vs hip FLS comparisons revealed increased expression of genes associated with proliferation/inflammation (eg, TLR2, VCAM1, MMP13) and decreased cell-cycle inhibitors (eg, CDKN1A, CDKN1C, CDKN2D) in hand. Pathway analysis revealed greatest enrichment for hand vs hip and included proinflammatory (eg, TNFs bind physiological receptors) and proliferative (eg, Integrin cell surface interactions) pathways (Table 1A). Taiji distinguished joint-specific TFs and identified greater PageRanks for key TFs (eg, ATF2, CREB1) in hand (Fig 1, Table 1A). Pathway analysis revealed greatest enrichment for hand vs hip and included proliferation/inflammation pathways (eg, MAPK activation, PIP3 activates AKT signaling) (Table 1A). The predicted differences in proliferation were confirmed in a cell growth assay showing that hand FLS proliferation was greater than hip or knee after PDGF-stim (Hand/knee or hip ratio=1.4; P< 0.001). The joint-specific transcriptome differences for hand vs hip and knee were enhanced by TNF-stim. DEGs within a particular joint between unstim vs TNF-stim revealed greatest pathway enrichment in hand (Table 1B). TNF-stim hand FLS had unique induction of proliferation/inflammation genes (eg, TANK, KSR1, CDC14A) and pathways (eg, Diseases of signal transduction by growth factor receptors/second messengers, Regulation of necroptotic cell death), and these pathways were not observed in TNF-stim hip or knee. TNF-stim in hand FLS had increased PageRanks of key inflammatory TFs (eg, RELA, STAT3) and pathways (eg, MAPK activation, Signaling by ALK in cancer) that were not observed for hip/knee.
Conclusion: Integrated analysis of FLS identified joint-specific patterns of gene expression and TFs, including pathway enrichment in hand FLS associated with inflammation and proliferation. These differences are more prominent after TNF stimulation. Hand joints are commonly earlier and are more severe in RA. Distinctive joint-specific FLS biology associated with hand-specific pathways might explain the distribution and severity of joints involved in RA.
Fig 1. PageRanks of top differential TFs (Wilcoxon, top TFs selected by p-value for each pairwise comparison) for hand, knee and hip joint-derived FLS cultured in medium. Columns are colored by joint annotation. TF PageRank patterns are distinct for the 3 joint locations and involve key regulators of growth, differentiation and inflammation. Labeled TFs are TFs enriched in Reactome or Gene Ontology (GO) pathways. Unlabeled rows are not pathway enriched.
Table 1 A. Pairwise (hand vs hip, hand vs knee, hip vs knee) DEGs (Wilcoxon, p<0.05, abs log2FC > 0.58) were identified. Reactome and GO identified significantly enriched pathways using differentially expressed genes (p<0.05). Pairwise differential TFs were identified using the TF PageRank scores (Wilcoxon, p<0.05). Pathway enrichment for differential TFs were conducted as for the transcriptome. B. Responses to TNF stimulation were evaluated for each joint location and compared TNF-stim to unstim. DEGs between unstim and TNF-stim FLS for a particular joint location were identified (Wilcoxon, p<0.05, abs log2FC > 0.58). A similar process identified TNF-stim vs unstim differential TFs using the TF PageRank scores (Wilcoxon, p<0.05). Pathway enrichment for DEGs and differential TFs were conducted as part A.
E. Choi: None; C. R. L. Machado: Eli Lilly, 5; D. Boyle: None; W. Wang: None; G. Firestein: Eli Lilly, 5.