Primary biliary cirrhosis (PBC) is an uncommon autoimmune disease with a homogeneous clinical phenotype that reflects incomplete disease concordance in monozygotic (MZ) twins. We have taken advantage of a unique collection consisting of genomic DNA and mRNA from peripheral blood cells of female MZ twins (n = 3 sets) and sisters of similar age (n = 8 pairs) discordant for disease. We performed a genome-wide study to investigate differences in (i) DNA methylation (using a custom tiled four-plex array containing tiled 50-mers 19,084 randomly chosen methylation sites), (ii) copy number variation (CNV) (with a chip including markers derived from the 1000 Genomes Project, all three HapMap phases, and recently published studies), and/or (iii) gene expression (by whole-genome expression arrays). Based on the results obtained from these three approaches we utilized quantitative PCR to compare the expression of candidate genes. Importantly, our data support consistent differences in discordant twins and siblings for the (i) methylation profiles of 60 gene regions, (ii) CNV of 10 genes, and (iii) the expression of 2 interferon-dependent genes. Quantitative PCR analysis showed that 17 of these genes are differentially expressed in discordant sibling pairs. In conclusion, we report that MZ twins and sisters discordant for PBC manifest particular epigenetic differences and highlight the value of the epigenetic study of twins.

Systemic lupus erythematosus (SLE) is propelled by pathogenic autoantibody (AutoAb) and immune pathway dysregulation. Identifying populations at risk of reaching classified SLE is essential to curtail inflammatory damage. Lupus blood relatives (Rel) have an increased risk of developing SLE. We tested factors to identify Rel at risk of developing incomplete lupus (ILE) or classified SLE vs. clinically unaffected Rel and healthy controls (HC), drawing from two unique, well characterized lupus cohorts, the lupus autoimmunity in relatives (LAUREL) follow-up cohort, consisting of Rel meeting <4 ACR criteria at baseline, and the Lupus Family Registry and Repository (LFRR), made up of SLE patients, lupus Rel, and HC. Medical record review determined ACR SLE classification criteria; study participants completed the SLE portion of the connective tissue disease questionnaire (SLE-CSQ), type 2 symptom questions, and provided samples for assessment of serum SLE-associated AutoAb specificities and 52 plasma immune mediators. Elevated SLE-CSQ scores were associated with type 2 symptoms, ACR scores, and serology in both cohorts. Fatigue at BL was associated with transition to classified SLE in the LAUREL cohort (p≤0.01). Increased levels of BLyS and decreased levels of IL-10 were associated with type 2 symptoms (p<0.05). SLE-CSQ scores, ACR scores, and accumulated AutoAb specificities correlated with levels of multiple inflammatory immune mediators (p<0.05), including BLyS, IL-2Rα, stem cell factor (SCF), soluble TNF receptors, and Th-1 type mediators and chemokines. Transition to SLE was associated with increased levels of SCF (p<0.05). ILE Rel also had increased levels of TNF-α and IFN-γ, offset by increased levels of regulatory IL-10 and TGF-β (p<0.05). Clinically unaffected Rel (vs. HC) had higher SLE-CSQ scores (p<0.001), increased serology (p<0.05), and increased inflammatory mediator levels, offset by increased IL-10 and TGF-β (p<0.01). These findings suggest that Rel at highest risk of transitioning to classified SLE have increased inflammation coupled with decreased regulatory mediators. In contrast, clinically unaffected Rel and Rel with ILE demonstrate increased inflammation offset with increased immune regulation, intimating a window of opportunity for early intervention and enrollment in prevention trials.

Systemic lupus erythematosus (SLE) patients exist an imbalance between regulatory T (Treg) and T helper 17 cells (Th17), which might be contributed by defective immune regulation of bone marrow derived mesenchymal stem cells (BM-MSCs) from SLE patients. Our microRNA array analysis showed markedly down-regulated expression levels of microRNA let-7f in BM-MSCs from SLE patients compared to those from normal controls (NOR). To explore the role of let-7f in the disease pathogenesis, we showed that expression levels of let-7f in SLE BM-MSCs were negatively associated with SLE disease activity, and the predicted let-7 family targeted gene expression of interlukin-6 (IL-6) was significantly higher in BM-MSCs from SLE patients compared to normal controls (NOR). Transient transfection of BM-MSCs with let-7f mimics or inhibitors showed reduced levels of let-7f impaired the proliferation rate of BM-MSCs, BM-MSC-mediated downregulation of Th17 cells and upregulation of Treg cells, increased the apoptosis rate of BM-MSCs through targeting IL-6 and activating signal transducers and activators of transcription-3 (STAT3) pathway, but had no significant effect on the differentiation of Th1 and Th2. Our findings showed a key role of let-7f in the imbalance of Treg/Th17 mediated by SLE BM-MSCs, suggesting the potential of manipulating let-7f expression in BM-MSCs for treating SLE patients.

There is an urgent need to identify biomarkers for diagnosis and disease activity monitoring in rheumatoid arthritis (RA). We leveraged publicly available microarray gene expression data in the NCBI GEO database for whole blood (N=1,885) and synovial (N=284) tissues from RA patients and healthy controls. We developed a robust machine learning feature selection pipeline with validation on five independent datasets culminating in 13 genes: TNFAIP6, S100A8, TNFSF10, DRAM1, LY96, QPCT, KYNU, ENTPD1, CLIC1, ATP6V0E1, HSP90AB1, NCL and CIRBP which define the RA score and demonstrate its clinical utility: the score tracks the disease activity DAS28 (p = 7e-9), distinguishes osteoarthritis (OA) from RA (OR 0.57, p = 8e-10) and polyJIA from healthy controls (OR 1.15, p = 2e-4) and monitors treatment effect in RA (p = 2e-4). Finally, the immunoblotting analysis of six proteins on an independent cohort confirmed two proteins, TNFAIP6/TSG6 and HSP90AB1/HSP90.

Fibroblast-like synoviocytes (FLS) exhibit a unique aggressive phenotype in rheumatoid arthritis (RA). Increased FLS migration and subsequent invasion of the extracellular matrix are essential to joint destruction in RA. Our previous research reported that transcription factor SOX5 was highly expressed in RA-FLS. Here, the effects of SOX5 in RA-FLS migration and invasion will be investigated.