Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 17 March 2023
Sec. Autoimmune and Autoinflammatory Disorders : Autoimmune Disorders
This article is part of the Research Topic Advances in understanding the pathogenesis of and designing therapies for connective tissue disease-associated interstitial lung disease View all 6 articles

Type I interferon score is associated with the severity and poor prognosis in anti-MDA5 antibody-positive dermatomyositis patients

Jinjing Qian&#x;Jinjing Qian1†Rui Li&#x;Rui Li1†Zhiwei ChenZhiwei Chen1Zehui Cao*&#x;Zehui Cao2*‡Liangjing Lu*&#x;Liangjing Lu1*‡Qiong Fu*&#x;Qiong Fu1*‡
  • 1Department of Rheumatology, Renji Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
  • 2Institute of Molecular Medicine, Renji Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China

Objectives: To investigate the clinical significance of the interferon (IFN) score, especially the IFN-I score, in patients with anti-melanoma differentiation-associated gene 5 (MDA5) antibody-positive dermatomyositis (anti-MDA5+ DM).

Methods: We enrolled 262 patients with different autoimmune diseases, including idiopathic inflammatory myopathy, systemic lupus erythematosus, rheumatoid arthritis, adult-onset Still’s disease, and Sjögren’s syndrome, as well as 58 healthy controls. Multiplex quantitative real-time polymerase chain reaction (RT-qPCR) using four TaqMan probes was used to evaluate type I IFN-stimulated genes (IFI44 and MX1), one type II IFN-stimulated gene (IRF1), and one internal control gene (HRPT1), which were used to determine the IFN-I score. The clinical features and disease activity index were compared between the high and low IFN-I score groups in 61 patients with anti-MDA5+ DM. The associations between laboratory findings and the predictive value of the baseline IFN-I score for mortality were analyzed.

Results: The IFN score was significantly higher in patients with anti-MDA5+ DM than in healthy controls. The IFN-I score was positively correlated with the serum IFN-α concentration, ferritin concentration, and Myositis Disease Activity Assessment Visual Analogue Scale (MYOACT) score. Compared with patients with a low IFN-I score, patients with a high IFN-I score showed a higher MYOACT score, C-reactive protein concentration, aspartate transaminase concentration, ferritin concentration, plasma cell percentage, and CD3+ T-cell percentage, as well as lower lymphocyte, natural killer cell, and monocyte counts. The 3-month survival rate was significantly lower in patients with an IFN-I score of >4.9 than in those with an IFN-I score of ≤4.9 (72.9% vs. 100%, respectively; P = 0.044).

Conclusion: The IFN score, especially the IFN-I score, measured by multiplex RT-qPCR is a valuable tool to monitor disease activity and predict mortality in patients with anti-MDA5+ DM.

1 Introduction

Inflammatory myopathies are a group of systemic autoimmune diseases that affect the muscles, skin, and other organs (1, 2). Anti-melanoma differentiation-associated gene 5 (MDA5) antibody is highly associated with a specific type of inflammatory myopathy known referred to as anti-MDA5+ dermatomyositis (anti-MDA5+ DM). MDA5, which is also named interferon induced with helicase C domain 1 (IFIH1), is a cytoplasmic sensor of viral RNA that activates a cascade of antiviral responses, including type I interferon (IFN) signaling and other proinflammatory cytokines (3, 4). Patients with anti-MDA5 antibody positivity are prone to mild or no muscle involvement and typical cutaneous manifestations, and they frequently develop rapidly progressive interstitial lung disease (RP-ILD) (5, 6). Of note, patients with RP-ILD usually do not respond to conventional immunosuppressive treatment. However, Janus kinase (JAK) inhibitors combined with glucocorticoids significantly improve the survival of patients with early-stage anti-MDA5+ DM with interstitial lung disease (7). However, the 6-month mortality rate remains high, especially in East Asian regions (5, 7).

The involvement of type 1 IFN in the pathogenesis of anti-MDA5+ DM has been proposed based on significant elevation in the expression of downstream stimulated genes in muscle, skin, lung, and peripheral blood. Several studies have demonstrated the presence of an elevated type I IFN signature in muscle biopsies from patients with DM, including anti-MDA5+ DM (8, 9). Furthermore, sarcoplasmic MxA expression detected by immunohistochemistry is a sensitive marker for diagnosing DM, and it reflects the severity of muscle involvement in patients with juvenile DM (1012). Similarities to muscle disease also exist in the skin. For example, the skin lesions of patients with DM show type I IFN-induced recruitment of CXCR3+ lymphocytes, and robust expression of the type I IFN marker MxA has been demonstrated in both blood vessels and interstitial fibroblasts in the skin of patients with anti-MDA5+ DM (1315). In the lungs of patients with anti-MDA5+ DM, increased proportions of interferon-stimulated gene ISG+ CD4+ T cells and ISG+ CD8+ T cells were detected, further highlighting overactivation of the type I IFN pathway in the affected lungs of patients with anti-MDA5+ DM (16). Compared with other cytokine-induced genes, such as the genes encoding tumor necrosis factor-α (TNF-α), interleukin (IL)-10, and IL-1β, the type I IFN signature in the blood is positively correlated with disease activity in individual patients with dermatomyositis (DM) or polymyositis (PM) during longitudinal follow-up according to a previous study (17). These findings suggest that the type I IFN signature might serve as a useful biomarker and could be used as a tool to monitor clinical efficacy outcomes, which is useful because the clinical management of patients with anti-MDA5+ DM is challenging.

Poor sensitivity and specificity have made enzyme-linked immunosorbent assay measurements of circulating IFN challenging (18). Despite the fact that RNA sequencing can yield a vast amount of data, it is time-consuming, which is not optimal in clinical practice where timely feedback is required to monitor disease progression. To set up a clinical routine test to measure the IFN score quickly and easily, we established a multiplex RT-qPCR assay to simultaneously measure four genes in a single reaction, which allows the two IFN families to be differentiated, including two type I IFN-stimulated genes, one type II IFN-stimulated gene, and one internal control gene.

By applying multiplex RT-qPCR, we assessed type I IFN activity by examining gene expression in peripheral blood mononuclear cells (PBMCs) obtained from patients with anti-MDA5+ DM to determine whether the IFN score could be used to differentiate anti-MDA5+ DM from other autoimmune diseases and different myopathy subtypes. Additionally, we examined whether the IFN score, especially the type I IFN score, correlates with disease activity in patients with anti-MDA5+ DM, along with its association with laboratory findings and ability to predict 3-month mortality.

2 Materials and methods

2.1 Patients, clinical features, and laboratory data

A total of 262 patients, including 61 patients with anti-MDA5+ DM and 201 patients with other autoimmune diseases, were recruited from February 2021 to September 2022 at the Department of Rheumatology, Renji Hospital, Shanghai, China. The study protocol was approved by the Ethics Committee of Renji Hospital (ID: 2013-126), Shanghai, China. Informed consent was obtained from the study participants. This study was conducted according to the Declaration of Helsinki. Patients with anti-MDA5+ DM fulfilled the modified Sontheimer’s criteria (19, 20) and were screened for myositis-specific autoantibodies (MSAs) and myositis-associated antibodies. Blood samples were collected at baseline, and the data were retrospectively obtained, including demographic information, laboratory data, clinical features, and measurements of idiopathic inflammatory myopathy (IIM) disease activity based on the Myositis Disease Activity Assessment Visual Analogue Scale (MYOACT) and The Cutaneous Dermatomyositis Disease Area and Severity Index (CDASI). Fifty-eight healthy volunteers with no known diseases were included as controls.

2.2 Detection of myositis-specific autoantibody

MSAs, including anti-MDA5 antibody, were analyzed using a commercially available immunoblot assay (EUROLINE Autoimmune Inflammatory Myopathies 16 Ag (IgG); Euroimmun, Lubeck, Germany).

2.3 PBMC isolation and RNA extraction

Whole blood was collected in ethylenediaminetetraacetic acid tubes (BD, 367863) and stored at 4°C until PBMC isolation. Serum samples were collected simultaneously. Human PBMCs were isolated using Lymphoprep (Stemcell, 07801), and erythrocytes were lysed with red blood cell lysis buffer (Tiangen, RT122-02). Following the manufacturer’s instructions, RNA was extracted from PBMCs using Direct-zol RNA Miniprep Kits (Zymo, R2052) and was used as a template for quantitative real-time polymerase chain reaction (RT-qPCR).

2.4 Real-time PCR analysis and scoring data

Multiplex Taqman-based RT-qPCR was performed to determine the expression of a panel of genes. HPRT1 was used as a housekeeping gene/internal control. The target genes were IFI44, Mx1/MxA, and IRF1. The relative expression of each target gene (T/R) was calculated using the 2−ΔCT method. The T/R of each target gene was then normalized as follows: (T/Rsubject − meanHC) ÷ SDHC, where HC is healthy control and SD is standard deviation. The IFN score was calculated as the mean of the normalized T/R of the three target genes. Scores higher than the mean of the HC plus two SDs were designated as high IFN scores; otherwise, the IFN score was considered low.

2.5 Statistical analysis

The Mann–Whitney U test with Bonferroni correction was performed to identify differences in variables between the high IFN-I score group and the low IFN-I score group. The chi-square test or Fisher’s exact test was used to compare categorical variables, as appropriate. Correlations between laboratory findings and disease activity indices were assessed using Spearman’s Rho. All graphs depict the mean ± SD, unless otherwise stated. All analyses were performed using SPSS (version 26.0), and a two-tailed P value of <0.05 was considered statistically significant.

3 Results

3.1 Multiplex RT-qPCR based on three IFN-stimulated genes revealed IFN pathway activation in autoimmune diseases

We collected peripheral blood samples from patients with five distinct autoimmune diseases, as well as from healthy controls. Among the 118 patients diagnosed with IIM, 61 had anti-MDA5+ DM, 20 had anti-MDA5 DM, 22 had anti-synthetase syndrome (ASS), six had polymyositis (PM), and nine had IMNM (immune-mediated necrotizing myopathy). There were 116 patients with systemic lupus erythematosus (SLE), 17 with rheumatoid arthritis (RA), five with adult-onset Still’s disease (AOSD), six with Sjögren’s syndrome (SS), and 58 healthy controls. Using RNA isolated from PBMCs as the template, multiplex RT-qPCR was performed to simultaneously quantify the relative expression of four genes, including two type I IFN-stimulated genes (IFI44 and MX1), one type II IFN-stimulated gene (IRF1), and one housekeeping gene (HPRT1) (Figure 1). After dividing the relative expression of each target gene by the mean of the normalized expression in HCs, the IFN score was calculated as the median of the relative expression of all contributing genes. The mean plus two SDs of the HC values was set as the cutoff for identifying a high IFN score. Therefore, patients with anti-MDA5+ DM were divided into two groups (the high and low IFN score groups). As shown in Figure 1, a significant clinical association was found between IFN-I score and MYOACT score.

FIGURE 1
www.frontiersin.org

Figure 1 Flowchart illustrating the overall experimental design of the study. Biorender and Adobe Illustrator were used to create the image.

3.2 IFN score in PBMC distinguishes anti-MDA5+DM from disease controls and health controls

We compared the IFN scores among patients with different autoimmune diseases. Figure 2A shows that patients with IIM (P = 0.000), SLE (P = 0.000), and RA (P = 0.006) had significantly higher IFN scores than healthy controls. No significant differences were found between patients with AOSD or SS and patients in the other groups. Remarkably, IIM and SLE were the two types of autoimmune disease with the most highly upregulated IFN scores. As stated previously, anti-MDA5+ DM is a distinct subtype of IIM with a high mortality rate. Therefore, we compared the IFN scores of various subtypes of IIM with the IFN score of HCs. Anti-MDA5+ DM had the highest IFN score of all IIM subtypes. Patients with anti-MDA5+ DM had significantly higher IFN scores than patients with ASS (P = 0.000) and HCs (P = 0.000) (Figure 2B). Patients with anti-MDA5− DM also had a significantly higher IFN score than HCs (P = 0.000). Based on the analysis of the 20 patients of MDA5-DM, the study found that 15% (3 out of 20) exhibited anti-SAE1 antibodies. In addition, 30% (6 out of 20) of patients had anti-NXP2 antibodies, 30% of patients had anti-TIF1-γ antibodies,15% (3 out of 20) had anti-Mi-2 antibodies, and 10% (2 out of 20) had anti-Ro52 antibodies. No significant differences were found between the IFN scores of the ASS, PM, IMNM, and HC groups.

FIGURE 2
www.frontiersin.org

Figure 2 Baseline IFN scores in different autoimmune diseases and in healthy controls. (A). IFN scores in the IIM, SLE, RA, AOSD, SS, and HC groups. (B). IFN scores in the anti-MDA5+ DM, anti-MDA5− DM, ASS, PM, INIM, and HC groups. The Kruskal–Wallis test with Bonferroni correction was performed using SPSS software (version 23.0). The IFN score of each individual is represented as a single dot. Horizontal lines represent the mean ± SD. IFN, interferon; IIM, idiopathic inflammatory myopathy; SLE, systemic lupus erythematosus; RA, rheumatoid arthritis; AOSD, adult-onset Still’s disease; SS, Sjögren’s syndrome; HC, healthy control.

3.3 Characteristics of patients with anti-MDA5+DM

IFN-I has been recognized as a key player in the pathogenesis of anti-MDA5 + DM, with several studies highlighting its involvement in the disease. Therefore we decided to focus primarily on the IFN-I score in this study. The baseline characteristics of patients with anti-MDA5+ DM are listed in Table 1. On the basis of the IFN-I score, 51 of 61 patients with anti-MDA5+ DM were assigned to the high IFN-I score group, while 10 patients were assigned to the low IFN-I score group. A statistically significant difference in disease duration was found at 9.68 ± 13.36 months in the high IFN-I score group and at 49.80 ± 13.70 months in the low IFN-I score group (P = 0.007). In terms of age, gender, clinical features, and treatments, there were no significant differences between the high and low IFN-I score groups. There were 6 patients with anti-MDA5 antibodies who were diagnosed with classic DM and 56 patients with anti-MDA5 antibodies who were diagnosed with clinically amyopathic dermatomyositis (CADM).

TABLE 1
www.frontiersin.org

Table 1 Baseline characteristics of patients with anti-MDA5+ DM and comparisons between the high and low IFN-I score groups.

3.4 The IFN-I score is related to the clinical characteristics of patients with anti-MDA5+ DM

The Mann–Whitney U test was performed to understand whether the laboratory test results differed between the high and low IFN-I score groups. The serum IFN-α concentration (15.34 ± 25.19 pg/ml vs. 9.32 ± 18.75 pg/ml, P = 0.015), neutrophil-to-lymphocyte ratio (8.96 ± 6.67 vs. 4.63 ± 3.78, P = 0.018), C-reactive protein (CRP) concentration (4.30 ± 6.65 mg/L vs. 0.67 ± 0.94 mg/L, P = 0.008), aspartate transaminase (AST) concentration (96.75 ± 189.70 U/L vs. 27.00 ± 7.86 U/L, P = 0.002), and ferritin concentration (1225.30 ±1192.26 μg/L vs. 514.28 ± 584.03 μg/L, P = 0.041), as well as the percentage of peripheral plasma cells (7.74% ± 9.30% vs. 2.03% ± 2.85%, P = 0.004) and CD3+ T cells (65.12% ± 12.15% vs. 53.16% ± 19.00%, P = 0.039), were significantly higher in the high IFN-I score group than in the low IFN-I score group (Figures 3A–F, Figure S1C). In contrast, patients with anti-MDA5+ DM in the high IFN-I score group had lower lymphocyte (0.83 ± 0.54 cells/μL vs. 1.29 ± 0.68 cells/μL , P = 0.036) and absolute monocyte (0.45 ± 0.29 x109/L vs. 0.55 ± 0.13 x109/L, P = 0.027) counts than patients in the low IFN-I score group (Figures 3G, H). The natural killer cell count (92.77 ± 100.22 cells/μL vs. 353.91 ± 437.46 cells/μL, P =0.004) and the percentage of CD3− CD16+ CD56+ natural killer cells (11.59% ± 8.21% vs. 26.29% ± 21.09%, P = 0.029) were also lower in the high IFN-I score group (Figures S1A–B).

FIGURE 3
www.frontiersin.org

Figure 3 Comparison of laboratory data between patients with anti-MDA5+ DM in the high (n = 51) and low (n = 10) IFN-I score groups. (A) Serum IFN-α concentration, (B) NLR, (C) CRP concentration, (D) AST concentration, (E) ferritin concentration, (F) PC%, (G) lymphocyte count, and (H) absolute monocyte count. Comparisons were performed using the Mann–Whitney U test. The bars represent the mean ± SD. NLR, neutrophil-to-lymphocyte ratio; CRP, C‐reactive protein; AST, aspartate transaminase; PC, plasma cell.

3.5 Correlation of the IFN-I score with laboratory findings and disease activity in patients with anti-MDA5+DM

We examined the correlations of the IFN-I score with the laboratory findings and disease activity in patients with anti-MDA5+ DM. Moderately positive correlations between the IFN-I score and serum IFN-α concentration (r = 0.335, P = 0.008), ferritin concentration (r = 0.302, P = 0.018), AST concentration (r = 0.343, P = 0.007), and percentage of peripheral plasma cells (r = 0.362, P = 0.004) were observed (Figures 4A, B, Figures S2A–B). In contrast, the IFN-I score was negatively correlated with the absolute natural killer cell count (r = −.0360, P = 0.004), the percentage of CD3-CD16+CD56+ natural killer cell (r = −.0364, P = 0.004), and the disease duration (r = −.0317, P = 0.013) (Figure S2C–E). The results show a positive association between the IFN-I score and laboratory findings. Overall, the IFN-I score was positively correlated with the MYOACT score (r = 0.426, P = 0.001) (Figure 4C). In terms of individual organ systems, the IFN-I score was moderately positively correlated with the MYOACT constitutional (r = 0.406, P = 0.001), MYOACT skeletal (r = 0.303, P = 0.018), MYOACT mucocutaneous (r = 0.378, P = 0.003), and severity according to the CDASI (r = 0.356, P = 0.005) (Figures 4D–G). These results imply that the IFN-I score is positively correlated with myositis disease activity in patients with anti-MDA5+ DM.

FIGURE 4
www.frontiersin.org

Figure 4 Correlation of the IFN-I score with disease activity parameters and clinical features. The IFN-I score in patients with anti-MDA5+ DM is associated with the concentrations of serum IFN-α (A) and ferritin (B). The IFN-I score also correlated with the MYOACT score (C), MYOACT constitutional (D), MYOACT skeletal (E), MYOACT mucocutaneous (F), and severity based on the CDASI (G). The data were analyzed using Spearman’s Rho.CDASI, cutaneous dermatomyositis disease area and severity index; MYOACT, Myositis Disease Activity Assessment Visual Analogue Scale.

3.6 A high IFN-I score identifies patients with active MDA5+DM and predicts mortality

We examined the difference in disease activity index between the high and low IFN-I score groups. Compared with the low IFN-I score group, the high IFN-I score group had a significantly higher MYOACT score (P=0.006), MYOACT constitutional (P = 0.001), MYOACT mucocutaneous (P = 0.001), and severity based on the CDASI (P = 0.007) (Figures 5A–D). In contrast, MYOACT skeletal, MYOACT pulmonary, and CDASI damage were unaffected (Figures S3A–C).

FIGURE 5
www.frontiersin.org

Figure 5  A high IFN-I score identifies patients with active anti-MDA5+ DM and predicts mortality. Grouped analysis of MYOACT scores and CDASI scores at baseline (A–F), (A) MYOACT score, (B) MYOACT constitutional, (C) MYOACT mucocutaneous, and (D) severity based on the CDASI. Patients were divided into two groups based on the IFN-I score (high vs. low). Comparisons were performed using the Mann–Whitney U test. (E, F) Correlations of the IFN-I score and the MYOACT score in 16 patients before and after treatment. (G) The Kaplan–Meier survival curve shows the survival rates stratified according to the IFN-I score. The high IFN-I score group had an IFN-I score of >4.9 (n = 48), whereas the low IFN-I score group had an IFN-I score of ≤4.9 (n = 13). Comparisons were performed using the log-rank test.

We analyzed the change in the IFN-I score in 16 patients during hospitalization. The results show that the IFN-I score paralleled the MYOACT disease activity index before (Figure 5E) and after (Figure 5F) treatment. During the observation period, 13 of 61 patients died of respiratory failure induced by interstitial lung disease. A significant difference was observed in the 3-month survival rate of patients with an IFN-I score of >4.9 and those with an IFN-I score of ≤4.9 (72.9% vs. 100%, respectively; P = 0.044) (Figure 5G). This suggests that a high IFN-I score of >4.9 may serve as a sensitive marker to predict mortality in patients with anti-MDA5+ DM.

4 Discussion

DM is a heterogeneous group of autoimmune diseases that can be subdivided based on the presence of various MSAs. Patients with anti-MDA5+ DM exhibit typical cutaneous manifestations and interstitial lung disease, but myositis is either less severe or non-existent in these patients. The poor prognosis results from the fact that traditional immunosuppressive treatment is frequently ineffective in patients with anti-MDA5+ DM with interstitial lung disease, especially rapidly progressive interstitial lung disease, which is a significant factor contributing to mortality. From the therapeutic perspective, additional simple and sensitive biomarkers should be identified to evaluate disease activity and forecast the disease course of patients with anti-MDA5+ DM. The selection of IFI44, Mx1, and IRF1 as targeted genes in this study was based on several considerations, including their previous identification as highly upregulated and specific interferon-stimulated genes (2123), the need to differentiate between IFN-I and IFN-II signaling pathways, and the importance of selecting genes with close expression and designing non-mutually exclusive primers for accurate and reliable detection using multiplex RT-qPCR assays. The application of multiplex RT-qPCR in this study of IFN-I scores in patients with anti-MDA5+ DM was innovative. The limitation of not being able to differentiate between the two IFN families was overcome by this approach, which covers both type I and type II IFN-stimulated genes, whereas previous gene sets only allowed differentiation into high and low IFN scores (24). As anticipated, the IFN-I score was increased in more than 80% of patients with anti-MDA5+ DM and was significantly higher in individuals with anti-MDA5+ DM than in patients with SLE, RA, AOSD, and SS. Meanwhile, 23% of patients with anti-MDA5+ DM exhibited elevated IFN-II scores and differences in CRP concentration, ALT concentration, and MYOACT skeletal were observed between the high IFN-II score and the low IFN-II score groups (Figure S4). Previous research has suggested that increased expression of IFN-II may contribute to the development and severity of MDA5+DM (25). However, the number of significant variations associated with IFN-II score was fewer than expected. This may be because IFN-II has a local effect on the tissue microenvironment and could not be easily detected in blood circulation. Thus, we primarily focused on analyzing the IFN-I score, as it appeared to be more strongly associated with practical parameters than the IFN-II score. The IFN-I score reflected the disease severity because it was significantly associated with ferritin concentration, AST concentration, and MYOACT scores. During the course of hospitalization, the IFN-I score was associated with disease activity. Notably, anti-MDA5 antibody-positive patients with an IFN-I score of >4.9 were likely to have a worse prognosis than those with an IFN-I score of ≤4.9, which may aid in the clinical care of these patients. Thus, patients with an IFN-I score of >4.9 may respond to anti-type I IFN therapy, such as JAK inhibitors. These results imply that the IFN-I score may contribute to the pathogenesis of anti-MDA5+ DM and could potentially serve as a diagnostic tool for this disease. However, it’s important to note that further research is needed to fully understand the role of IFN-II in this disease, as our findings highlight the importance of both IFN-I and IFN-II in anti-MDA5+ DM.

The results of this study were in accordance with those of a previous study by Higgs et al., who showed that the type I IFN signature can be found in the blood of patients with SLE, RA, PM, Systemic Sclerosis (SSc), and DM (26). Compared with ASS and other myopathies, DM is typically associated with a strong type I IFN signature (27, 28). These characteristics were corroborated in the present study by the existence of the IFN-I signature in the blood of patients with anti-MDA5+ DM.

Earlier investigators compared the IFN scores between disease types. In the current investigation, in addition to evaluating the IFN score between selected diseases, patients with anti-MDA5+ DM were separated into two groups (the high and low IFN-I score groups). Intriguingly, the percentage of plasma cells in the high IFN-I score group was considerably higher than in the low IFN-I score group. This result may be because RNA-containing immune complexes generated by MDA5 and anti-MDA5 are potent inducers of IFN-α through Toll-like receptor-7, and type I IFN promotes B-cell differentiation into plasma cells (2932). In a multivariate logistic regression model developed by Xu et al., anti-MDA5 positivity, a CRP concentration of >50 μg/L, and a lymphocyte count of <500/μL predicted rapidly progressive interstitial lung disease with an overall accuracy of 90% (33). In a recent study, the lymphocyte and monocyte counts of the non-survivor group in the first month were significantly lower than the survivor group (34). In this study, the high IFN-I score group also had a higher CRP concentration and decreased lymphocyte and monocyte counts. These findings suggest that in clinical practice, patients with a reasonably high IFN-I score, either alone or in combination with other serum indicators, should be prioritized for higher mortality.

The IFN-I score has been used in the past to quantify disease activity. The type I IFN signature score in the blood is strongly associated with disease activity in DM and PM during baseline evaluation and longitudinal follow-up according to previous studies (17, 26, 27, 35). The disease activity index, as well as the ferritin, serum IFN-α, and AST concentrations, was positively correlated with the IFN-I score in the present study. IFN-α promotes the synthesis of ferritin (36). Horai et al. discovered that anti-MDA5+ individuals with interstitial lung disease had higher IFN-α and ferritin concentrations than anti-MDA5− patients with interstitial lung disease (37). According to Lian et al. and Gono et al., the serum ferritin concentration is significantly higher in patients with anti-MDA5+ DM and predicts the disease severity and prognosis of rapidly progressive interstitial lung disease (3840). Ferritin was also elevated in the high IFN-I score group, indicating that it could be useful to monitor disease activity and assist in stratifying high-risk groups. A statistically significant difference in disease duration was found between high and low IFN-I score groups suggesting treatment over a longer period could lead to decreasing levels of biomarkers. The lack of significant differences in the use of GC, the dosage of GC, and several immunosuppressants between high and low IFN-I groups as shown in Table 1 suggest that medication use may not be a major confounding factor in the cross-sectional study analysis of the IFN-I score. It is possible that low IFN-I score patients have milder conditions and therefore longer survival time, which could contribute to a longer follow-up period while patients with a shorter disease duration may have a more acute and severe form of the disease. However, long duration of treatment and high doses of GCs and/or immunosuppressants might affect the levels of cytokines and gene expression, which can in turn influence the levels of IFN-I scores. Exploring and studying the use of a normal or decreasing IFN-I score as a target to be achieved with drug therapy could be a promising area of research to improve the treatment of anti-MDA5-positive dermatomyositis.

Multiple markers are associated with the prognosis of anti-MDA5+ DM, including anti-MDA5 titers (41, 42), the presence of anti-Ro52 antibody (43), lactate dehydrogenase (18), ferritin (5, 24, 26), KL-6 (27), the proportion of CD4+ CXCR4+ T cells (28), and a high proportion of ISG15+ CD8+ T cells (16). The presence of interstitial lung disease, particularly rapidly progressive interstitial lung disease, poses significant obstacles to the prognosis of patients with anti-MDA5+ DM. Overactivation of the type I IFN pathway is apparent in the affected lungs of patients with anti-MDA5+ DM according to a previous study (16). Type I IFN recruits CX3CR1+ M2 macrophages in the lungs by inducing the secretion of CX3CL (44, 45), and M2 macrophages are responsible for producing tumor growth factor-β to promote pulmonary fibrosis (4648). However, little previous evidence supports the idea that the IFN score in patients with anti-MDA5+ DM is directly associated with a poor prognosis. We propose that an IFN-I score of >4.9 could be useful to predict outcomes and assess treatment efficacy, especially in patients treated with targeted synthetic disease-modifying antirheumatic drugs, thereby addressing a previously unfulfilled need. The lack of a strong association between the IFN-I score and the MYOACT pulmonary score in patients with MDA5+ DM suggests that the IFN-I score may not be a reliable indicator of lung involvement in these patients. It is possible that the IFN-I score reflects a more systemic immune response rather than a localized response in the lungs. However, A high IFN-I score is still considered an indicator of a high-risk immune and inflammatory state that is linked to disease severity and poor prognosis. Therefore, patients with a high IFN-I score should still receive clinical attention and aggressive treatment, as this may be an early warning sign before the onset of significant pulmonary symptoms.

This study had several limitations. First, anti-MDA5+ DM is relatively uncommon, and thus the research was limited by the small sample size. Second, the study was performed at a single medical center in Shanghai, China. Therefore, external validation and large-scale studies are required to validate the findings. Third, pulmonary function testing was not performed because the overall medical condition of some patients was pretty severe, and these patients could therefore not undergo pulmonary function testing. Finally, due to their severe medical condition, some patients were already being treated with high glucocorticoid doses or JAK inhibitors at the time of admission, which is likely to interfere with the type I IFN pathway and can cause alterations in gene expression. Hence, large-scale studies in patients with the new-onset disease are required.

5 Conclusion

The present study offers a novel approach to comprehending the dysregulated type I IFN production in patients with anti-MDA5+ DM and understanding how it relates to the disease severity. In the future, dissecting the cellular and molecular mechanisms underpinning this process will be essential. Precise risk assessment based on clinical phenotypes and IFN scores to help implement early targeted and personalized therapy in high-risk patients may significantly improve the long-term prognosis of patients with anti-MDA5+ DM.

Equations

Normalized T/Rsubject=(T/RsubjectmeanHC)÷SDHC

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving human participants were reviewed and approved by The Ethics Committee of Renji Hospital (ID: 2013-126), Shanghai, China. The patients/participants provided their written informed consent to participate in this study.

Author contributions

QF, LL and ZHC supervised the study and revised the manuscript. JQ performed the experiments and drafted the manuscript. RL and ZWC collected and analyzed clinical data. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China [81974251], and Shanghai Hospital Development Center, Joint Research of New Advanced Technology Project [SHDC12018106].

Acknowledgments

We thank all the patients for participating in this study. We thank Emily Woodhouse, PhD, from Edanz for editing the English text of a draft of this manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1151695/full#supplementary-material

References

1. Lundberg IE, Fujimoto M, Vencovsky J, Aggarwal R, Holmqvist M, Christopher-Stine L, et al. Idiopathic inflammatory myopathies. Nat Rev Dis Primers (2021) 7(1):86. doi: 10.1038/s41572-021-00321-x

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Dalakas MC. Inflammatory muscle diseases. N Engl J Med (2015) 372(18):1734–47. doi: 10.1056/NEJMra1402225

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Liu G, Lee JH, Parker ZM, Acharya D, Chiang JJ, van Gent M, et al. Isg15-dependent activation of the sensor Mda5 is antagonized by the sars-Cov-2 papain-like protease to evade host innate immunity. Nat Microbiol (2021) 6(4):467–78. doi: 10.1038/s41564-021-00884-1

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Wu XM, Zhang J, Li PW, Hu YW, Cao L, Ouyang S, et al. Nod1 promotes antiviral signaling by binding viral rna and regulating the interaction of Mda5 and mavs. J Immunol (2020) 204(8):2216–31. doi: 10.4049/jimmunol.1900667

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Kurasawa K, Arai S, Namiki Y, Tanaka A, Takamura Y, Owada T, et al. Tofacitinib for refractory interstitial lung diseases in anti-melanoma differentiation-associated 5 gene antibody-positive dermatomyositis. Rheumatol (Oxford) (2018) 57(12):2114–9. doi: 10.1093/rheumatology/key188

CrossRef Full Text | Google Scholar

6. Wu W, Guo L, Fu Y, Wang K, Zhang D, Xu W, et al. Interstitial lung disease in anti-Mda5 positive dermatomyositis. Clin Rev Allergy Immunol (2021) 60(2):293–304. doi: 10.1007/s12016-020-08822-5

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Chen Z, Wang X, Ye S. Tofacitinib in amyopathic dermatomyositis–associated interstitial lung disease. New Engl J Med (2019) 381(3):291–3. doi: 10.1056/NEJMc1900045

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Moneta GM, Pires Marafon D, Marasco E, Rosina S, Verardo M, Fiorillo C, et al. Muscle expression of type I and type ii interferons is increased in juvenile dermatomyositis and related to clinical and histologic features. Arthritis Rheumatol (2019) 71(6):1011–21. doi: 10.1002/art.40800

PubMed Abstract | CrossRef Full Text | Google Scholar

9. O'Connor KA, Abbott KA, Sabin B, Kuroda M, Pachman LM. Mxa gene expression in juvenile dermatomyositis peripheral blood monon uclear cells: Association with muscle involvement. Clin Immunol (2006) 120(3):319–25. doi: 10.1016/j.clim.2006.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Uruha A, Nishikawa A, Tsuburaya RS, Hamanaka K, Kuwana M, Watanabe Y, et al. Sarcoplasmic mxa expression: A valuable marker of dermatomyositis. Neurology (2017) 88(5):493–500. doi: 10.1212/wnl.0000000000003568

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Greenberg SA, Pinkus JL, Pinkus GS, Burleson T, Sanoudou D, Tawil R, et al. Interferon-α/β-Mediated innate immune mechanisms in dermatomyositis. Ann Neurol (2005) 57(5):664–78. doi: 10.1002/ana.20464

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Preuße C, Allenbach Y, Hoffmann O, Goebel HH, Pehl D, Radke J, et al. Differential roles of hypoxia and innate immunity in juvenile and adult dermatomyositis. Acta Neuropathol Commun (2016) 4(1):45. doi: 10.1186/s40478-016-0308-5

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Wenzel J, Scheler M, Bieber T, Tüting T. Evidence for a role of type I interferons in the pathogenesis of dermatomyositis. Br J Dermatol (2005) 153(2):462–3. doi: 10.1111/j.1365-2133.2005.06786.x

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Wenzel J, Schmidt R, Proelss J, Zahn S, Bieber T, Tüting T. Type I interferon-associated skin recruitment of Cxcr3+ lymphocytes in dermatomyositis. Clin Exp Dermatol (2006) 31(4):576–82. doi: 10.1111/j.1365-2230.2006.02150.x

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Ono N, Kai K, Maruyama A, Sakai M, Sadanaga Y, Koarada S, et al. The relationship between type 1 ifn and vasculopathy in anti-Mda5 antibody-positive dermatomyositis patients. Rheumatology (2018) 58(5):786–91. doi: 10.1093/rheumatology/key386

CrossRef Full Text | Google Scholar

16. Ye Y, Chen Z, Jiang S, Jia F, Li T, Lu X, et al. Single-cell profiling reveals distinct adaptive immune hallmarks in Mda5+ dermatomyositis with therapeutic implications. Nat Commun (2022) 13(1):6458. doi: 10.1038/s41467-022-34145-4

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Greenberg SA, Higgs BW, Morehouse C, Walsh RJ, Kong SW, Brohawn P, et al. Relationship between disease activity and type 1 interferon- and other cytokine-inducible gene expression in blood in dermatomyositis and polymyositis. Genes Immun (2012) 13(3):207–13. doi: 10.1038/gene.2011.61

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Chasset F, Mathian A, Dorgham K, Ribi C, Trendelenburg M, Huynh-Do U, et al. Serum interferon-α levels and ifn type I-stimulated genes score perform equally to assess systemic lupus erythematosus disease activity. Ann Rheum Dis (2022) 81(6):901–3. doi: 10.1136/annrheumdis-2021-221835

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Sontheimer RD. Would a new name hasten the acceptance of amyopathic dermatomyositis (Dermatomyositis siné myositis) as a distinctive subset within the idiopathic inflammatory dermatomyopathies spectrum of clinical illness? J Am Acad Dermatol (2002) 46(4):626–36. doi: 10.1067/mjd.2002.120621

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Ye S, Chen XX, Lu XY, Wu MF, Deng Y, Huang WQ, et al. Adult clinically amyopathic dermatomyositis with rapid progressive interstitial lung disease: A retrospective cohort study. Clin Rheumatol (2007) 26(10):1647–54. doi: 10.1007/s10067-007-0562-9

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Baechler EC, Bauer JW, Slattery CA, Ortmann WA, Espe KJ, Novitzke J, et al. An interferon signature in the peripheral blood of dermatomyositis patients is associated with disease activity. Mol Med (2007) 13(1-2):59–68. doi: 10.2119/2006-00085

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Huard C, Gullà SV, Bennett DV, Coyle AJ, Vleugels RA, Greenberg SA. Correlation of cutaneous disease activity with type 1 interferon gene signature and interferon β in dermatomyositis. Br J Dermatol (2017) 176(5):1224–30. doi: 10.1111/bjd.15006

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Rosain J, Neehus AL, Manry J, Yang R, Le Pen J, Daher W, et al. Human Irf1 governs macrophagic ifn-Γ immunity to mycobacteria. Cell (2023) 186(3):621–45.e33. doi: 10.1016/j.cell.2022.12.038

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Nocturne G, Mariette X. Interferon signature in systemic autoimmune diseases: What does it mean? RMD Open (2022) 8(2):e002687. doi: 10.1136/rmdopen-2022-002687

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Gono T, Kaneko H, Kawaguchi Y, Hanaoka M, Kataoka S, Kuwana M, et al. Cytokine profiles in polymyositis and dermatomyositis complicated by rapidly progressive or chronic interstitial lung disease. Rheumatology (2014) 53(12):2196–203. doi: 10.1093/rheumatology/keu258

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Higgs BW, Liu Z, White B, Zhu W, White WI, Morehouse C, et al. Patients with systemic lupus erythematosus, myositis, rheumatoid arthritis and scleroderma share activation of a common type I interferon pathway. Ann Rheum Dis (2011) 70(11):2029–36. doi: 10.1136/ard.2011.150326

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Pinal-Fernandez I, Casal-Dominguez M, Derfoul A, Pak K, Plotz P, Miller FW, et al. Identification of distinctive interferon gene signatures in different types of myositis. Neurology (2019) 93(12):e1193–e204. doi: 10.1212/wnl.0000000000008128

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Rigolet M, Hou C, Baba Amer Y, Aouizerate J, Periou B, Gherardi RK, et al. Distinct interferon signatures stratify inflammatory and dysimmune myopathies. RMD Open (2019) 5(1):e000811. doi: 10.1136/rmdopen-2018-000811

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Bekeredjian-Ding IB, Wagner M, Hornung V, Giese T, Schnurr M, Endres S, et al. Plasmacytoid dendritic cells control Tlr7 sensitivity of naive b cells Via type I ifn. J Immunol (2005) 174(7):4043–50. doi: 10.4049/jimmunol.174.7.4043

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Wang K, Zhao J, Wu W, Xu W, Sun S, Chen Z, et al. Rna-containing immune complexes formed by anti-melanoma differentiation associated gene 5 autoantibody are potent inducers of ifn-α. Front Immunol (2021) 12:743704. doi: 10.3389/fimmu.2021.743704

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Barnas JL, Albrecht J, Meednu N, Alzamareh DF, Baker C, McDavid A, et al. B cell activation and plasma cell differentiation are promoted by ifn-Λ in systemic lupus erythematosus. J Immunol (2021) 207(11):2660–72. doi: 10.4049/jimmunol.2100339

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Jego G, Palucka AK, Blanck JP, Chalouni C, Pascual V, Banchereau J. Plasmacytoid dendritic cells induce plasma cell differentiation through type I interferon and interleukin 6. Immunity (2003) 19(2):225–34. doi: 10.1016/s1074-7613(03)00208-5

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Xu Y, Yang CS, Li YJ, Liu XD, Wang JN, Zhao Q, et al. Predictive factors of rapidly progressive-interstitial lung disease in patients with clinically amyopathic dermatomyositis. Clin Rheumatol (2016) 35(1):113–6. doi: 10.1007/s10067-015-3139-z

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Lv X, Jin Y, Zhang D, Li Y, Fu Y, Wang S, et al. Low circulating monocytes is in parallel with lymphopenia which predicts poor outcome in anti-melanoma differentiation-associated gene 5 antibody-positive dermatomyositis-associated interstitial lung disease. Front Med (Lausanne) (2021) 8:808875. doi: 10.3389/fmed.2021.808875

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Walsh RJ, Kong SW, Yao Y, Jallal B, Kiener PA, Pinkus JL, et al. Type I interferon-inducible gene expression in blood is present and reflects disease activity in dermatomyositis and polymyositis. Arthritis Rheum (2007) 56(11):3784–92. doi: 10.1002/art.22928

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Stam TC, Swaak AJ, Kruit WH, Eggermont AM. Regulation of ferritin: A specific role for interferon-alpha (Ifn-alpha)? the acute phase response in patients treated with ifn-Alpha-2b. Eur J Clin Invest (2002) 32(Suppl 1):79–83. doi: 10.1046/j.1365-2362.2002.0320s1079.x

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Horai Y, Koga T, Fujikawa K, Takatani A, Nishino A, Nakashima Y, et al. Serum interferon-α is a useful biomarker in patients with anti-melanoma differentiation-associated gene 5 (Mda5) antibody-positive dermatomyositis. Mod Rheumatol (2015) 25(1):85–9. doi: 10.3109/14397595.2014.900843

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Lian X, Zou J, Guo Q, Chen S, Lu L, Wang R, et al. Mortality risk prediction in amyopathic dermatomyositis associated with interstitial lung disease: The flair model. Chest (2020) 158(4):1535–45. doi: 10.1016/j.chest.2020.04.057

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Gono T, Kawaguchi Y, Hara M, Masuda I, Katsumata Y, Shinozaki M, et al. Increased ferritin predicts development and severity of acute interstitial lung disease as a complication of dermatomyositis. Rheumatol (Oxford) (2010) 49(7):1354–60. doi: 10.1093/rheumatology/keq073

CrossRef Full Text | Google Scholar

40. Gono T, Kawaguchi Y, Satoh T, Kuwana M, Katsumata Y, Takagi K, et al. Clinical manifestation and prognostic factor in anti-melanoma differentiation-associated gene 5 antibody-associated interstitial lung disease as a complication of dermatomyositis. Rheumatol (Oxford) (2010) 49(9):1713–9. doi: 10.1093/rheumatology/keq149

CrossRef Full Text | Google Scholar

41. Bodewes ILA, Al-Ali S, van Helden-Meeuwsen CG, Maria NI, Tarn J, Lendrem DW, et al. Systemic interferon type I and type ii signatures in primary sjogren's syndrome reveal differences in biological disease activity. Rheumatol (Oxford) (2018) 57(5):921–30. doi: 10.1093/rheumatology/kex490

CrossRef Full Text | Google Scholar

42. Crow MK. Type I interferon in the pathogenesis of lupus. J Immunol (2014) 192(12):5459–68. doi: 10.4049/jimmunol.1002795

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Lambers WM, de Leeuw K, Doornbos-van der Meer B, Diercks GFH, Bootsma H, Westra J. Interferon score is increased in incomplete systemic lupus erythematosus and correlates with myxovirus-resistance protein a in blood and skin. Arthritis Res Ther (2019) 21(1):260. doi: 10.1186/s13075-019-2034-4

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Nakano M, Fujii T, Hashimoto M, Yukawa N, Yoshifuji H, Ohmura K, et al. Type I interferon induces Cx3cl1 (Fractalkine) and Ccl5 (Rantes) production in human pulmonary vascular endothelial cells. Clin Exp Immunol (2012) 170(1):94–100. doi: 10.1111/j.1365-2249.2012.04638.x

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Ishida Y, Kimura A, Nosaka M, Kuninaka Y, Hemmi H, Sasaki I, et al. Essential involvement of the Cx3cl1-Cx3cr1 axis in bleomycin-induced pulmonary fibrosis Via regulation of fibrocyte and M2 macrophage migration. Sci Rep (2017) 7(1):16833. doi: 10.1038/s41598-017-17007-8

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Zhang L, Wang Y, Wu G, Xiong W, Gu W, Wang CY. Macrophages: Friend or foe in idiopathic pulmonary fibrosis? Respir Res (2018) 19(1):170. doi: 10.1186/s12931-018-0864-2

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Gono T, Miyake K, Kawaguchi Y, Kaneko H, Shinozaki M, Yamanaka H. Hyperferritinaemia and macrophage activation in a patient with interstitial lung disease with clinically amyopathic dm. Rheumatol (Oxford) (2012) 51(7):1336–8. doi: 10.1093/rheumatology/kes012

CrossRef Full Text | Google Scholar

48. Horiike Y, Suzuki Y, Fujisawa T, Yasui H, Karayama M, Hozumi H, et al. Successful classification of macrophage-mannose receptor Cd206 in severity of anti-Mda5 antibody positive dermatomyositis associated ild. Rheumatol (Oxford) (2019) 58(12):2143–52. doi: 10.1093/rheumatology/kez185

CrossRef Full Text | Google Scholar

Keywords: dermatomyositis, anti-MDA5 antibody, interstitial lung disease (ILD), IFN signature, multiplex RT-qPCR assay

Citation: Qian J, Li R, Chen Z, Cao Z, Lu L and Fu Q (2023) Type I interferon score is associated with the severity and poor prognosis in anti-MDA5 antibody-positive dermatomyositis patients. Front. Immunol. 14:1151695. doi: 10.3389/fimmu.2023.1151695

Received: 26 January 2023; Accepted: 07 March 2023;
Published: 17 March 2023.

Edited by:

Xiaoxia Zhu, Fudan University, China

Reviewed by:

Wenfeng Tan, Nanjing Medical University, China
Takahisa Gono, Nippon Medical School, Japan

Copyright © 2023 Qian, Li, Chen, Cao, Lu and Fu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Qiong Fu, ZnVxaW9uZzVAMTYzLmNvbQ==; Liangjing Lu, bHVfbGlhbmdqaW5nQDE2My5jb20=; Zehui Cao, emluY2ltbUBvdXRsb29rLmNvbQ==

†These authors share first authorship

These authors contributed equally to this work and share senior authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.