BACKGROUND: To report baseline characteristics, patient reported outcomes and treatment of children with Juvenile Dermatomyositis (JDM) in the Childhood Arthritis and Rheumatology Research Alliance (CARRA) Registry. METHODS: Children newly diagnosed with JDM were enrolled in the CARRA Registry from 41 pediatric rheumatology centers. Baseline patient demographics, disease characteristics, assessments, patient reported outcome and treatments were recorded. RESULTS: In the first year, 119 JDM participants were enrolled. Most were female (63.4%), and white (72.3%) with a median diagnosis age 8.0 years (IQR 4.0-11.5), and median age of disease onset 7.0 years (IQR 3.5-7.5). They had characteristic rashes (92.4%), elevated muscle enzymes (83.2%), physician global score 4.0 (IQR 2.5-5.0) and manual muscle testing score 63.5 (IQR 51.0-75.0). Calcinosis (3.4%) and interstitial lung disease (< 1%) were uncommon. Myositis specific antibodies were measured and reported in nearly half of participants enrolled where anti-MJ followed by Anti-p155/140 were most common (11/49 and 7/53 respectively). Childhood Health Assessment Questionnaire (CHAQ) results showed mild-moderate disability (median 0.750, IQR 0.030-1.875), as did patient/parent global assessments of disease activity (median 3, patient IQR: 1.75-5.25; parent IQR: 1-7). Patient Reported Outcomes Measurement Information System (PROMIS®) Pediatric Global Health 7 scores, Pain Interference, Physical Function scores for Mobility, and Upper Extremity Function were commonly worse than 95% of the general pediatric population. CONCLUSIONS: In its inaugural year, 119 JDM patients were successfully enrolled in participapte in the New CARRA Registy. This registry will provide the necessary foundation to advance clinical research to improve outcomes using traditional measures and patient reported outcomes. With the CARRA biorepository, this infrastructure will enable future translational research. Together, these efforts may aid in future clinical trials, including comparative effectiveness trials.

Objective

We conducted a comprehensive gene expression meta-analysis in dermatomyositis (DM) muscle and skin tissues to identify shared disease-relevant genes and pathways across tissues.

Methods

Six publicly available data sets from DM muscle and two from skin were identified. Meta-analysis was performed by first processing data sets individually then cross-study normalization and merging creating tissue-specific gene expression matrices for subsequent analysis. Complementary single-gene and network analyses using Significance Analysis of Microarrays (SAM) and Weighted Gene Co-expression Network Analysis (WGCNA) were conducted to identify genes significantly associated with DM. Cell-type enrichment was performed using xCell.

Results

There were 544 differentially expressed genes (FC ≥ 1.3, q < 0.05) in muscle and 300 in skin. There were 94 shared upregulated genes across tissues enriched in type I and II interferon (IFN) signaling and major histocompatibility complex (MHC) class I antigen-processing pathways. In a network analysis, we identified eight significant gene modules in muscle and seven in skin. The most highly correlated modules were enriched in pathways consistent with the single-gene analysis. Additional pathways uncovered by WGCNA included T-cell activation and T-cell receptor signaling. In the cell-type enrichment analysis, both tissues were highly enriched in activated dendritic cells and M1 macrophages.

Conclusion

There is striking similarity in gene expression across DM target tissues with enrichment of type I and II IFN pathways, MHC class I antigen-processing, T-cell activation, and antigen-presenting cells. These results suggest IFN-γ may contribute to the global IFN signature in DM, and altered auto-antigen presentation through the class I MHC pathway may be important in disease pathogenesis.

Aging and cellular senescence are increasingly recognized as key contributors to pulmonary fibrosis. However, our understanding in the context of scleroderma-associated interstitial lung disease (SSc-ILD) is limited. To investigate, we leveraged previously established lung aging- and cell-specific senescence signatures to determine their presence and potential relevance to SSc-ILD. We performed a gene expression meta-analysis of lung tissues from 38 SSc-ILD and 18 healthy controls and found that markers (GDF15, COMP, and CDKN2A) and pathways (p53) of senescence were significantly increased in SSc-ILD. When probing the established aging and cellular senescence signatures, we found that epithelial and fibroblast senescence signatures had a 3.6- and 3.7-fold enrichment, respectively, in the lung tissue of SSc-ILD and that lung aging genes (CDKN2A, FRZB, PDE1A, and NAPI12) were increased in SSc-ILD. These signatures were also enriched in SSc skin and associated with degree of skin involvement (limited vs. diffuse cutaneous). To further support these findings, we examined telomere length (TL), a surrogate for aging, in the lung tissue and found that, independent of age, SSc-ILD had significantly shorter telomeres than controls in type II alveolar cells in the lung. TL in SSc-ILD was comparable to idiopathic pulmonary fibrosis, a disease of known aberrant aging. Taken together, this study provides novel insight into the possible mechanistic effects of accelerated aging and aberrant cellular senescence in SSc-ILD pathogenesis.

Juvenile dermatomyositis (JDM) is a rare autoimmune condition with insufficient biomarkers and treatments, in part, due to incomplete knowledge of the cell types mediating disease. We investigated immunophenotypes and cell-specific genes associated with disease activity using multiplexed RNA and protein single-cell sequencing applied to PBMCs from 4 treatment-naïve JDM (TN-JDM) subjects at baseline, 2, 4, and 6 months post-treatment and 4 subjects with inactive disease on treatment. Analysis of 55,564 cells revealed separate clustering of TN-JDM cells within monocyte, NK, CD8⁺ effector T and naïve B populations. The proportion of CD16⁺ monocytes was reduced in TN-JDM, and naïve B cells and CD4⁺ Tregs were expanded. Cell-type differential gene expression analysis and hierarchical clustering identified a pan-cell-type IFN gene signature over-expressed in TN-JDM in all cell types and correlated with disease activity most strongly in cytotoxic cell types. TN-JDM CD16⁺ monocytes expressed the highest IFN gene score and were highly skewed toward an inflammatory and antigen-presenting phenotype at both the transcriptomic and proteomic levels. A transitional B cell population with a distinct transcriptomic signature was expanded in TN-JDM and characterized by higher CD24 and CD5 proteins and less CD39, an immunoregulatory protein. This data provides new insights into JDM immune dysregulation at cellular resolution and serves as a novel resource for myositis investigators.

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.

Juvenile dermatomyositis (JDM) is one of several childhood-onset autoimmune disorders characterized by a type I IFN response and autoantibodies. Treatment options are limited due to an incomplete understanding of how the disease emerges from dysregulated cell states across the immune system. We therefore investigated the blood of patients with JDM at different stages of disease activity using single-cell transcriptomics paired with surface protein expression. By immunophenotyping peripheral blood mononuclear cells, we observed skewing of the B cell compartment toward an immature naive state as a hallmark of JDM at diagnosis. Furthermore, we find that these changes in B cells are paralleled by T cell signatures suggestive of Th2-mediated inflammation that persist despite disease quiescence. We applied network analysis to reveal that hyperactivation of the type I IFN response in all immune populations is coordinated with previously masked cell states including dysfunctional protein processing in CD4+ T cells and regulation of cell death programming in NK cells, CD8+ T cells, and γδ T cells. Together, these findings unveil the coordinated immune dysregulation underpinning JDM and provide insight into strategies for restoring balance in immune function.

The Muller F element (4.2 Mb, ~80 protein-coding genes) is an unusual autosome of Drosophila melanogaster; it is mostly heterochromatic with a low recombination rate. To investigate how these properties impact the evolution of repeats and genes, we manually improved the sequence and annotated the genes on the D. erecta, D. mojavensis, and D. grimshawi F elements and euchromatic domains from the Muller D element. We find that F elements have greater transposon density (25-50%) than euchromatic reference regions (3-11%). Among the F elements, D. grimshawi has the lowest transposon density (particularly DINE-1: 2% vs. 11-27%). F element genes have larger coding spans, more coding exons, larger introns, and lower codon bias. Comparison of the Effective Number of Codons with the Codon Adaptation Index shows that, in contrast to the other species, codon bias in D. grimshawi F element genes can be attributed primarily to selection instead of mutational biases, suggesting that density and types of transposons affect the degree of local heterochromatin formation. F element genes have lower estimated DNA melting temperatures than D element genes, potentially facilitating transcription through heterochromatin. Most F element genes (~90%) have remained on that element, but the F element has smaller syntenic blocks than genome averages (3.4-3.6 vs. 8.4-8.8 genes per block), indicating greater rates of inversion despite lower rates of recombination. Overall, the F element has maintained characteristics that are distinct from other autosomes in the Drosophila lineage, illuminating the constraints imposed by a heterochromatic milieu.