Co-expression network analysis reveals the pivotal role of mitochondrial dysfunction and interferon signature in juvenile dermatomyositis

Juvenile idiopathic inflammatory myositis: an update on pathophysiology and clinical care

The childhood-onset or juvenile idiopathic inflammatory myopathies (JIIMs) are a heterogenous group of rare and serious autoimmune diseases of children and young people that predominantly affect the muscles and skin but can also involve other organs, including the lungs, gut, joints, heart and central nervous system. Different myositis-specific autoantibodies have been identified that are associated with different muscle biopsy features, as well as with different clinical characteristics, prognoses and treatment responses. Thus, myositis-specific autoantibodies can be used to subset JIIMs into sub-phenotypes; some of these sub-phenotypes parallel disease seen in adults, whereas others are distinct from adult-onset idiopathic inflammatory myopathies. Although treatments and management have much improved over the past decade, evidence is still lacking for many of the current treatments and few validated prognostic biomarkers are available with which to predict response to treatment, comorbidities (such as calcinosis) or outcome. Emerging data on the pathogenesis of the JIIMs are leading to proposals for new trials and tools for monitoring disease.

Role of CD14+ monocyte-derived oxidised mitochondrial DNA in the inflammatory interferon type 1 signature in juvenile dermatomyositis

OBJECTIVES

To define the host mechanisms contributing to the pathological interferon (IFN) type 1 signature in Juvenile dermatomyositis (JDM).

METHODS

RNA-sequencing was performed on CD4+, CD8+, CD14+ and CD19+ cells sorted from pretreatment and on-treatment JDM (pretreatment n=10, on-treatment n=11) and age/sex-matched child healthy-control (CHC n=4) peripheral blood mononuclear cell (PBMC). Mitochondrial morphology and superoxide were assessed by fluorescence microscopy, cellular metabolism by 13C glucose uptake assays, and oxidised mitochondrial DNA (oxmtDNA) content by dot-blot. Healthy-control PBMC and JDM pretreatment PBMC were cultured with IFN-α, oxmtDNA, cGAS-inhibitor, TLR-9 antagonist and/or n-acetyl cysteine (NAC). IFN-stimulated gene (ISGs) expression was measured by qPCR. Total numbers of patient and controls for functional experiments, JDM n=82, total CHC n=35.

RESULTS

Dysregulated mitochondrial-associated gene expression correlated with increased ISG expression in JDM CD14+ monocytes. Altered mitochondrial-associated gene expression was paralleled by altered mitochondrial biology, including 'megamitochondria', cellular metabolism and a decrease in gene expression of superoxide dismutase (SOD)1. This was associated with enhanced production of oxidised mitochondrial (oxmt)DNA. OxmtDNA induced ISG expression in healthy PBMC, which was blocked by targeting oxidative stress and intracellular nucleic acid sensing pathways. Complementary experiments showed that, under in vitro experimental conditions, targeting these pathways via the antioxidant drug NAC, TLR9 antagonist and to a lesser extent cGAS-inhibitor, suppressed ISG expression in pretreatment JDM PBMC.

CONCLUSIONS

These results describe a novel pathway where altered mitochondrial biology in JDM CD14+ monocytes lead to oxmtDNA production and stimulates ISG expression. Targeting this pathway has therapeutical potential in JDM and other IFN type 1-driven autoimmune diseases.

Respiratory chain dysfunction in perifascicular muscle fibres in patients with dermatomyositis is associated with mitochondrial DNA depletion

AIMS

Patients with dermatomyositis (DM) suffer from reduced aerobic metabolism contributing to impaired muscle function, which has been linked to cytochrome c oxidase (COX) deficiency in muscle tissue. This mitochondrial respiratory chain dysfunction is typically seen in perifascicular regions, which also show the most intense inflammatory reaction along with capillary loss and muscle fibre atrophy. The objective of this study was to investigate the pathobiology of the oxidative phosphorylation deficiency in DM.

METHODS

Muscle biopsy specimens with perifascicular COX deficiency from five juveniles and seven adults with DM were investigated. We combined immunohistochemical analyses of subunits in the respiratory chain including complex I (subunit NDUFB8), complex II (succinate dehydrogenase, subunit SDHB) and complex IV (COX, subunit MTCO1) with in situ hybridisation, next generation deep sequencing and quantitative polymerase chain reaction (PCR).

RESULTS

There was a profound deficiency of complexes I and IV in the perifascicular regions with enzyme histochemical COX deficiency, whereas succinate dehydrogenase activity and complex II were preserved. In situ hybridisation of mitochondrial RNA showed depletion of mitochondrial DNA (mtDNA) transcripts in the perifascicular regions. Analysis of mtDNA by next generation deep sequencing and quantitative PCR in affected muscle regions showed an overall reduction of mtDNA copy number particularly in the perifascicular regions.

CONCLUSION

The respiratory chain dysfunction in DM muscle is associated with mtDNA depletion causing deficiency of complexes I and IV, which are partially encoded by mtDNA, whereas complex II, which is entirely encoded by nuclear DNA, is preserved. The depletion of mtDNA indicates a perturbed replication of mtDNA explaining the muscle pathology and the disturbed aerobic metabolism.

Myocardial involvement in idiopathic inflammatory myopathies: a multi-center cross-sectional study in the CRDC-MYO Registry

OBJECTIVES

This study aimed to investigate the associated factors of myocardial involvements (MIs) in patients with idiopathic inflammatory myopathies (IIMs).

METHODS

In this multi-center cross-sectional study, 1946 patients with IIMs were enrolled from Chinese Rheumatism Data Center-Myositis Registry (CRDC-MYO). A total of 108 (5.5%) patients were identified with MIs, including congestive heart failure (n = 67, 62.0%), and severe arrhythmias (n = 61, 56.5%). The other 1838 IIM patients without IMs were set as the control group. Clinical features were collected including age, gender, comorbidities, clinical symptoms, clinical signs of both IIMs and MIs, lab findings including myositis-specific antibodies (MSAs) and myositis-associated antibodies (MAAs), echocardiogram, and radiological exams. Multivariate logistic analysis was used to explore independent associated factors of MIs in patients with IIMs.

RESULTS

Several independent associated factors were identified in multi-variate logistic regression, including positivity for anti-mitochondrial antibody-subtype 2 (AMA-M2) (OR 5.194, 95% CI 2.509-10.753, P < 0.001), elevation of creatine kinase (CK) (OR 2.611, 95% CI 1.312-5.198, P = 0.006), elevation of C-reactive protein (CRP) (OR 2.150, 95% CI 1.211-3.818, P = 0.001), and pulmonary hypertension (OR 4.165, 95% CI 1.765-9.882, P = 0.009). AMA-M2 and pulmonary hypertension were the most consistent associated factors in the polymyositis subgroup and the dermatomyositis/clinically amyopathic dermatomyositis subgroup.

CONCLUSIONS

MIs are rare but serious complication of IIMs could lead to congestive heart failure and severe arrhythmias. IIM patients with AMA-M2 positivity, elevation of CK and CRP, and pulmonary hypertension are more likely to develop MI complications. Key Points • This study investigated the independent associated factors for clinically significant myocardial involvements among idiopathic inflammatory myopathies in a large-scale, nation-wide multi-center cross-sectional study.

Identification of diagnostic genes and vital microRNAs involved in rheumatoid arthritis: based on data mining and experimental verification

Background

The pathogenesis of rheumatoid arthritis (RA) is complex. This study aimed to identify diagnostic biomarkers and transcriptional regulators that underlie RA based on bioinformatics analysis and experimental verification.

Material and Methods

We applied weighted gene co-expression network analysis (WGCNA) to analyze dataset GSE55457 and obtained the key module most relevant to the RA phenotype. We then conducted gene function annotation, gene set enrichment analysis (GSEA) and immunocytes quantitative analysis (CIBERSORT). Moreover, the intersection of differentially expressed genes (DEGs) and genes within the key module were entered into the STRING database to construct an interaction network and to mine hub genes. We predicted microRNA (miRNA) using a web-based tool (miRDB). Finally, hub genes and vital miRNAs were validated with independent GEO datasets, RT-qPCR and Western blot.

Results

A total of 367 DEGs were characterized by differential expression analysis. The WGCNA method divided genes into 14 modules, and we focused on the turquoise module containing 845 genes. Gene function annotation and GSEA suggested that immune response and inflammatory signaling pathways are the molecular mechanisms behind RA. Nine hub genes were screened from the network and seven vital regulators were obtained using miRNA prediction. CIBERSORT analysis identified five cell types enriched in RA samples, which were closely related to the expression of hub genes. Through ROC curve and RT-qPCR validation, we confirmed five genes that were specific for RA, including CCL25, CXCL9, CXCL10, CXCL11, and CXCL13. Moreover, we selected a representative gene (CXCL10) for Western blot validation. Vital miRNAs verification showed that only the differences in has-miR-573 and has-miR-34a were statistically significant.

Conclusion

Our study reveals diagnostic genes and vital microRNAs highly related to RA, which could help improve our understanding of the molecular mechanisms underlying the disorder and provide theoretical support for the future exploration of innovative therapeutic approaches.

Weighted gene co-expression network analysis and drug-gene interaction bioinformatics uncover key genes associated with various presentations of malaria infection in African children and major drug candidates

Malaria is a fatal parasitic disease with unelucidated pathogenetic mechanism. Herein, we aimed to uncover genes associated with different clinical aspects of malaria based on the GSE1124 dataset that is publicly accessible by using WGCNA. We obtained 16 co-expression modules and their correlations with clinical features. Using the MCODE tool, we identified THEM4, STYX, VPS36, LCOR, KIAA1143, EEA1, RAPGEF6, LOC439994, ZBTB33, PTPN22, ESCO1, and KLF3 as hub genes positively associated with Plasmodium falciparum infection (ASPF). These hub genes were involved in the biological processes of endosomal transport, regulation of natural killer cell proliferation, and KEGG pathways of endocytosis and fatty acid elongation. For the purple module negatively correlated with ASPF, we identified 19 hub genes that were involved in the biological processes of positive regulation of cellular protein catabolic process and KEGG pathways of other glycan degradation. For the salmon module positively correlated with severe malaria anemia (SMA), we identified 17 hub genes that were among those driving the biological processes of positive regulation of erythrocyte differentiation. For the brown module positively correlated with cerebral malaria (CM), we identified eight hub genes and these genes participated in phagolysosome assembly and positive regulation of exosomal secretion, and animal mitophagy pathway. For the tan module negatively correlated with CM, we identified four hub genes that were involved in CD8-positive, alpha-beta T cell differentiation and notching signaling pathway. These findings may provide new insights into the pathogenesis of malaria and help define new diagnostic and therapeutic approaches for malaria patients.