Abstract
Recent studies have provided insights into the autoinflammation triggered by severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) infection, which is associated with high mortality of coronavirus disease 2019 (COVID-19). Striking similarities has been noted between COVID-19 and anti-melanoma differentiation-associated gene 5 (MDA5) antibody (Ab)-related dermatomyositis (DM), implying a shared autoinflammatory aberrance. However, it is unclear whether anti-MDA5 Ab is present in COVID-19 and correlates with the severity and adverse outcome of COVID-19 patients. Here, we found that the positive rate of anti-MDA5 Ab in patients with COVID-19 was 48.2% and the anti-MDA5 Ab positive patients tended to develop severe disease (88.6% vs 66.9%, P<0.0001). In particular, the titer of anti-MDA5 Ab was increased in the non-survivals (5.95±5.16 vs 8.22±6.64, P=0.030) and the positive rate was also higher than that in the survivals (23.5% vs 12.0%, P=0.012). Regarding to severe COVID-19 patients, we found that high titer of anti-MDA5 Ab (≥10.0 U/mL) was more prevalent in the non-survivals (31.2% vs 14.0%, P=0.006). Moreover, early profiling of anti-MDA5 Ab could distinguish severe patients from those with non-severe ones. Overall, our data reveal that anti-MDA5 Ab is prevalent in the COVID-19 patients and high titer of this antibody is correlated with severe disease and unfavorable outcomes.
Coronavirus Disease 2019 (COVID-19), caused by highly contagious severe acute respiratory syndrome coronavirus 2 (SARS-CoV2), has become a pandemic involving more than 12 million cases globally by July 2020.1 The average mortality is estimated to be 1%,2 but can raise up to 62% in critically ill patients, mostly due to acute respiratory distress syndrome (ARDS).3 Therefore, reducing mortality of severe COVID-19 patients has become an urgent task in this battle. Unfortunately, few antiviral agents have been proved to be effective enough to treat the disease, and whether to use corticosteroids and other immunomodulatory drugs remains controversial.4
Accumulating evidence has demonstrated that a bunch of cytokines (e.g. IL-2, IL-6, IL-7, IL-10, G-CSF, IP10, MCP1, TNFα, etc.) are elevated in the blood of patients with severe COVID-19, resembling cytokine storm syndrome. The COVID-19 patients with hypercytokinemia and hyperferritinemia were more likely to exhibit extensive lung damage and ARDS.5-7 In addition, high prevalence of antinuclear antibodies (35.6%) and lupus anticoagulant (46.6%), along with antineutrophil cytoplasmic antibodies (6.6%) and anti-Ro antibody (4.4%), were also identified in hospitalized patients with COVID-19.8 Thus, hypothesis that SARS-CoV2 might trigger autoimmune and/or autoinflammatory aberrance in genetically predisposed subjects has been raised.9 It is reasonable that all patients with severe COVID-19 should be screened for hyperinflammation to identify the subgroup of patients for whom adjunctive immunosuppression therapy would improve mortality.4
Striking similarities has been noted between multifaceted features of COVID-19 and a rare autoimmune disease, the anti-melanoma-differentiation-associated gene 5 (MDA5) antibody (Ab)-related dermatomyositis (DM).10,11 Both diseases can involve the lungs, skin12,13 and muscles.14 Anti-MDA5 Ab-related DM patients often develop rapid progressive interstitial lung disease (RP-ILD), whose prognosis is disappointingly poor with more than half of patients dying within 3 months of disease onset.15 The initial radiological features of lung in Anti-MDA5 Ab-related DM patients are mainly subpleural ground-glass opacities or mixed with consolidation and signs of ARDS, which resemble severe and critical COVID-19.5,16 Besides, similar skin eruptions have been reported in both diseases.13 Furthermore, serum cytokine profiles are also similar in these two conditions, as serum levels of ferritin, IL-6, IL-8, and IL-10 usually were elevated in patients with RP-ILD secondary to Anti-MDA5 Ab-related DM.17 The similarity of these two diseases implies a shared underlying autoinflammatory/autoimmune mechanisms. To date, there is no report on whether anti-MDA5 Ab also exists in patients with COVID-19.
MDA5 is a well-known cytoplasmic sensor for viral RNA and its expression is induced by RNA viruses. This activates the expression of antiviral type I and III interferons (IFNs) with inflammatory cytokines. Correspondingly, IFN signaling can induce the expression of MDA5. SARS-CoV2 infection has been reported to trigger the expression of MDA5.18,19 In addition, MDA5 is involved in pathogenesis of several autoimmune disorders as well,20 such as systemic lupus erythematosus21,22, multiple sclerosis23, and even type 1 diabetes24,25. Nevertheless, it remains unclear whether the anti-MDA5 Ab plays a role in the pathophysiology of COVID-19 or whether it correlates with the disease severity. Some researchers have called for screening the anti-MDA5 Ab in severe COVID-19 patients.10,11
In the present study, we report the presence of anti-MDA5 Ab in patients with SARS-CoV2 infection and address its correlation with the clinical severity and outcomes of COVID-19 by measuring the anti-MDA5 Ab.
Results
Patients description
A total of 274 patients were included in this study (Figure 1a; Table 1). Of these patients, 230 cases from Wuhan Jinyintan hospital, the time admitted to hospital is from Dec 1, 2019 to Jan 27, 2020. 24 cases from Harbin infectious disease hospital, the time admitted to hospital is from Apr 2 to Apr 19. 20 cases from Beijing Ditan hospital, the time admitted to hospital is from Jan 20, 2020 to Jan 26, 2020. The median age was 56 years (IQR, 45-65 years), and 159 (58.0%) patients were male. The average disease course from onset of symptoms to discharge was 22.8±9.6 days. According to the definition of disease severity, 212 (77.4%) patients were classified as severe disease. Nearly half of the patients (n=119, 43.4%) had underlying chronic diseases, including hypertension, coronary arterial disease, chronic lung disease, and diabetes mellitus. On admission, 43 (15.7%) patients were complicated with respiratory failure, shock or other organ dysfunctions. Thirty-one (11.3%) were transferred to intensive care unit during their hospital stay, and 48 (17.5%) patients died.
Anti-MDA5 Ab is identified in the plasma of patients with COVID-19
To determine the presence of anti-MDA5 Ab in COVID-19 patients, ELISA analysis was employed to test the plasma collected from COVID-19 patients. We found that the titer of anti-MDA5 Ab was increased in these examined samples as compared to normal controls (1.85±0.67 vs 6.60±5.50, P<0.0001) (Figure 1b). The plasma from five patients of anti-MDA5 Ab-related DM were used as positive controls (Supplementary Figure 1). Based on the cut-off value (5.0 U/mL) of the anti-MDA5 Ab ELISA Kit, we noticed that the positive rate of anti-MDA5 Ab was also higher in patients with COVID-19 than that in normal controls (48.2%, 132/274, P<0.0001) (Figure 1c). These data were further validated by immunoblots in selected COVID-19 plasma samples. To this aim, we firstly performed MDA5 overexpression in 293T cells as shown by Western blotting analysis (Figure 1d). The plasma of anti-MDA5 Ab-related DM patients included in the ELISA were also confirmed by Western blots (Figure 1e). Next, a total of 17 plasma samples of COVID-19 were conducted Western blotting analysis, which included five non-severe and 12 severe COVID-19 patients. These data showed that the anti-MDA5 Ab were detected in these examined samples as well (Figures 1f and 1g).
Altogether, our findings indicate that SARS-CoV2 infection leads to an increased anti-MDA5 Ab titer in patients with COVID-19.
COVID-19 patients with positive anti-MDA5 Ab tend to exhibit severe disease
Anti-MDA5 Ab is first identified in DM and correlated with the status of this disease, which promoted us to further investigate whether the titer of anti-MDA5 Ab was related to the clinical severity of patients with COVID-19. To this end, the confirmed COVID-19 patients were stratified into two groups based on the cutoff value, that is, the anti-MDA5 Ab negative (<5.0 U/mL) and positive (≥5.0 U/mL). We found that the percentage of severe COVID-19 patients was much higher in anti-MDA5 Ab positive group than that in negative group (88.6% vs 66.9%, P<0.0001) (Figure 2a). The survival rate of anti-MDA5 Ab positive patients with COVID-19 is much lower compared with the negative group (76.5% vs 88.0%, P=0.012) (Figure 2b). As expected, the COVID-19 patients with positive anti-MDA5 Ab tended to have much longer disease course at discharge and develop respiratory failure, shock and other organ dysfunction (Figures 2c and 2d). These data indicated that COVID-19 patients with increased anti-MDA5 Ab titer exhibited severe clinical performance. Next, a univariate analysis was employed to investigate the correlation between anti-MDA5 Ab and other COVID-19 prognostic factors (Table 2). We found that the titer of anti-MDA5 Ab was positively correlated with the age of COVID-19 patients (Figure 2f). We also noticed that COVID-19 patients with positive anti-MDA5 Ab depicted decreased lymphocytes, increased neutrophils (Figures 2f and 2g). The levels of albumin were found to reduce in anti-MDA5 Ab positive patients compared with the negative (Figure 2h). The ratio of neutrophils versus lymphocytes (NLR) or CRP versus albumin (CAR) was much higher in anti-MDA5 Ab positive samples than that in the negative, indicating much severer inflammatory damage (Figures 2j and 2k). No significant difference was observed in Creatine Kinase (CK), lactate dehydrogenase (LDH), ferritin, and CRP (Table 2).
Taken together, our findings suggest that anti-MDA5 Ab is positively correlated with the clinical severity of COVID-19 patients.
The correlation between anti-MDA5 Ab and COVID-19 outcomes
As mentioned above, COVID-19 patients with positive anti-MDA5 Ab tended to develop severe disease. The titer of anti-MDA5 Ab was higher in severe COVID-19 patients as compared to the non-severe diseases (Figure 3a). The positive rate of this antibody was also higher in COVID-19 patients with severe performance (Figure 3b). We also observed that the titer of anti-MDA5 Ab depicted a significant increase in COVID-19 patients suffering from the chronic comorbidities, for instance, hypertension, diabetes, and cardiovascular disease (Figure 3c). An augment of this antibody was noticed in the plasma samples of COVID-19 patients suffering from shock, respiratory or other organ failure (Figure 3d).
We next investigated whether the titer of anti-MDA5 Ab was correlated with the outcome of COVID-19 patients. To this end, a further comparison of anti-MDA5 Ab was employed in the survival COVID-19 patients and the non-survivals, which indicated that the titer of anti-MDA5 Ab was upregulated in dead COVID-19 patients as compared to the survival (Figure 3e). Its positive rate was lower in the survivals with COVID-19 than that in the non-survivals (Figure 3f). These data suggested that anti-MDA5 Ab had the potential to serve as a prognostic factor for COVID-19. Consistent with published predictive factors for COVID-19 outcomes, we found that the levels of LDH, ferritin, and CRP were significantly decreased in the non-survivals as compared to that in the survivals, and the number of lymphocytes was also markedly reduced in the non-survivals (Figure 3g; Table 1).
We further performed a comparison of the anti-MDA5 Ab in COVID-19 patients with non-severe, severe performance and those deceased. The titer of anti-MDA5 Ab and positive rate were increased in severe and deceased patients compared with the non-severe ones (Figures 3h and 3i). Although both of the titer and positive rate of anti-MDA5 Ab depicted a moderate increase in the deceased patients as compared to the severe ones, no significant difference was observed between these two clusters (Figures 3h and 3i). In addition, we addressed the difference between the survivals and non-survivals in severe COVID-19 patients using 2-fold cut-off value based on the ELISA kit and found that the percentage of COVID-19 patients with high titer of anti-MDA5 Ab (≥ 10.0 U/mL) was elevated in the non-survivals than that in the survivals (Figure 3j).
Taken together, our data indicate that anti-MDA5 Ab is a marker for prognosis of COVID-19 patients and severe COVID-19 patients with high titer of anti-MDA5 Ab tend to have increased mortality.
Early profile of anti-MDA5 Ab distinguishes the prognosis of non-severe and severe COVID-19
Since the alteration of anti-MDA5 Ab titer is correlated with the activity and outcome of DM, we asked whether the change of anti-MDA5 Ab was associated with the clinical features of COVID-19. To this end, a cross-sectional analysis was employed using the titer of anti-MDA5 Ab achieved from the whole disease course. A total of 273 cases was stratified into three clusters based on the weeks following symptoms onset (WFSO) and shown as follows, WFSO-1, WFSO-2, and WFSO-3 (Figure 4a). A longitudinal profiling of anti-MDA5 Ab from 23 cases of severe COVID-19 patients was also determined at four sequential time-points (Figure 4a). A significant increase of the positive rate of anti-MDA5 Ab was observed in the samples from WFSO-2 and WFSO-3, as compared to WFSO-1, although no difference of the anti-MDA5 Ab titer was noticed in these three clusters (Figures 4b and 4c). These data revealed that the dynamic alteration of anti-MDA5 Ab might be various in the disease course of COVID-19 patients with diverse clinical performance. To test this idea, we compared anti-MDA5 Ab at three intervals as stated above in non-severe and severe patients, respectively. Interestingly, the titer of anti-MDA5 Ab in non-severe patients with COVID-19 was significantly increased at WFSO-2 as compared to that in WFSO-1, and then decreased at WFSO-3 (Figure 4d,). Similar result was found in the positive rate of anti-MDA5 Ab (Figure 4e). On the contrary, COVID-19 patients with severe performance exhibited high titer of anti-MDA5 Ab at the disease onset (WFSO-1) and then decreased at WFSO-2 and -3 (Figure 4f). However, no significant alteration was observed in the positive rate of this antibody at all three intervals, indicating that high titer of anti-MDA5 was preserved in the disease course of severe COVID-19 (Figure 4g).
We further determined the titer of anti-MDA5 Ab in sequential samples from severe COVID-19 disease as shown in Figure 4a. The titer of anti-MDA5 Ab in patient #1, #2, and #3 depicted a similar alteration compared with that in the cross-sectional analysis (Figures 4h and 4i). Next, the titer of anti-MDA5 Ab was examined in the samples collected at the WFSO-2, that is, the days following symptoms onset (DFSO) 8-14. We found that the titer and positive rate of anti-MDA5 Ab remained substantial (Figures 4j, 4k, and 4l; Supplementary Figure 2).
Collectively, our data indicate that COVID-19 patients with high titer of anti-MDA5 Ab initially tend to develop severe disease.
Discussion
The present study, for the first time, identified and confirmed the prevalence of anti-MDA5 Ab in patients with COVID-19 by both ELISA and Western blots. We also demonstrated that the positive rate and titer of anti-MDA5 Ab was associated with the clinical severity and outcomes of COVID-19. In severe COVID-19 patients, we found that high titer of anti-MDA5 Ab (≥10.0 U/mL) was more prevalent in non-survival patients. Moreover, early profile of anti-MDA5 Ab could distinguish severe patients from non-severe ones. Our study provides the evidence that early screening of anti-MDA5 Ab might help identify high risk population and predict the outcome of patients with COVID-19.
MDA5 is a intracellular RNA-specific helicase that belongs to a family of retinoic acid-inducible gene I-like receptors.26 It is triggered to oligomerize upon binding RNA via its helicase domain, which transmits a signal through its homotypic-interacting caspase recruitment domain (CARD) to induce oligomerization of the adaptor mitochondrial antiviral signaling (MAVS) protein and then activates the production of type I and III IFNs, which plays a key role in anti-viral immune response. 27,28 Thus, MDA5 is a crucial antiviral factor and has been previously reported to involve in SARS-CoV, MERS-CoV and SARS-CoV2 infections.18, 29,30 Interestingly, MDA5 is also involved in several autoimmune disorders such as anti-MDA5 Ab-related DM. Therefore, it is not surprising that COVID-19 and anti-MDA5 Ab-related DM share similar features of hyperinflammation and multi-systemic manifestations, especially rapid progressive interstitial lung disease that results in ARDS and death. In this study, we determined anti-MDA5 Ab in as many as 48.2% patients with COVID-19.
Our study revealed a positive correlation between the anti-MDA5 Ab and the severity of COVID-19, and high titer of anti-MDA5 Ab was associated with higher mortality in severe COVID-19 patients. Similar observation was reported in anti-MDA5 Ab-related DM patients.31 However, the titer of this antibody is even higher in anti-MDA5 Ab-related DM than that in COVID-19. This may indicate that high titer of anti-MDA5 Ab probably is related to an uncontrolled autoinflammation and autoimmune response to SARS-CoV2 infection in genetically predisposed hosts. Furthermore, our study also demonstrated that elder age, chronic comorbidities, lymphocytopenia, hypoalbuminemia, hyperferritinemia, increased D-dimer and CRP levels were more prevalent in COVID-19 patients with organ dysfunction and the mortality was comparatively high, which has been reported in previous studies and implies a dysregulation of inflamation.5-7,32-34
It has been reported that the change of anti-MDA5 Ab titer correlates with disease activity and predicts treatment response and disease outcome in patients with DM and rapidly progressive interstitial lung disease. 31 Our data also indicated that the dynamic titer alteration of anti-MDA5 Ab clearly varied in COVID-19 patients with diverse clinical severities. In the non-severe patients, the titer of anti-MDA5 Ab is upregulated in week 2 after symptom onset and then decreased, suggesting that the IFNs-MDA5 circuit is under fine-tuning regulation and the immune homeostasis is preserved in the total process of SARS-CoV2 infection (Figures 4m and 4n). However, in the severe COVID-19 patients, the titer of anti-MDA5 Ab boosts up in the 1st week after symptom onset and subsequently remains at a high positivity although a decreased titer is observed at weeks 2, 3, and 4 (Figures 4m and 4n). These data further supported that the MDA5 signaling might be persistently over-activated in severe COVID-19 patients. These findings also suggest that early screening and serially monitoring of anti-MDA5 Ab titer has the potential to predict the disease progression of COVID-19.
Several studies have already shown effectiveness of tocilizumab (IL-6 receptor blockade)35, ruxolitinib (JAK inhibitor)36 and tacrolimus37 in inhibiting SARS-CoV2 replication, improving the chest tomography or facilitating clinical improvement. Recently, dexamethasone has also been reported to improve survival in severe COVID-19 patients as well.38 Our findings provide supportive evidence that anti-inflammation and immunosuppressive therapy might be compromising strategy for the treatment of COVID-19, especially in those with high titer of anti-MDA5 Ab.
There are several limitations in our study. Firstly, since MDA5 is validated as a general sensor for diverse RNA viruses, no evidence has addressed whether anti-MDA5 Ab is present in the infection of other RNA viruses, for instance, influenza virus, enterovirus, and other coronaviruses. Therefore, the specificity of anti-MDA5 Ab in COVID-19 need to be further investigated in other RNA viruses. Secondly, due to limited sample size and endpoint events, we were not able to further evaluate whether anti-MDA5 Ab is an independent predictive factor for the death in COVID-19 or could be included in a risk stratification model. Thirdly, all patients were from China and it is not clear whether patients with other genetic backgrounds would have same results. Our findings are to be validated in a larger population of different ethnicities in future.
Overall, we, for the first time, revealed that anti-MDA5 Ab is present in patients with COVID-19 and correlates with severe disease and poor outcomes. Early screening and serially monitoring of anti-MDA5 Ab titer has the potential to predict the disease progression of COVID-19.
Methods
Study design and population
This retrospective study included three cohorts of adult patients (≥18 years old) from Jinyintan Hospital (Wuhan, China), Beijing Ditan Hospital (Beijing, China), and Heilongjiang Infectious Disease Hospital (Harbin, China), who were hospitalized from Dec 1, 2019 to Apr 19, 2020. All patients who were diagnosed with COVID-19 according to the Protocol for Prevention and Control of COVID-19 (Edition 7) promulgated by National Health Commission of China.39 All patients with COVID-19 were tested positive for SARS-CoV2 by use of quantitative, real-time PCR technology on samples from the upper airway specimens (pharyngeal or nasal swabs, nasopharyngeal secretions).
The plasma of patients with COVID-19 were collected within 24 hours after admission, and stored in -80 degree refrigerator. The plasma of five patients with anti-MDA5 Ab-related DM were provided by the Biobank of Myositis Registry of Department of Rheumatology, Peking Union Medical College Hospital. All of the 5 DM patients were diagnosed based on the criteria of Bohan and Peter.40,41
The study was approved by the Research Ethics Committee of the participating hospitals and the ethical board of the Institute of Pathogen Biology, Chinese Academy of Medical Sciences. The requirement for informed consent was waived by the Ethics Commission of the designated hospitals for emerging infectious diseases as described previously5.
Data collection
We extracted demographic, clinical, laboratory, treatment, and outcome data from medical and nursing records using standardized data collection forms (a revised version of case record form for severe acute respiratory infection shared by WHO and the International Severe Acute Respiratory and Emerging Infection Consortium). All data were checked by two investigators (QW and CL) and a third researcher (BC) adjudicated any difference in interpretation between the two primary reviewers.
According to the clinical classification of COVID-19 by the Protocol for Prevention and Control of cases of COVID-19 (Edition 7),39 we divided the patients into two groups on hospitalization: (1) Severe group, the patients fulfilled the diagnostic criteria of severe and critical cases, who meets any one of follows: i) Respiratory rate > 30/min, ii) Pulse oxygen saturation < 93%, iii) Oxygenation index < 300 mmHg, or iv) respiratory failure or other organ dysfunction requiring transmission to intensive care unit. (2) Non-severe group, the patients’ severity was mild or moderate that didn’t meet the above criteria.
Western Blots
293T cell line was obtained from the American Type Culture Collection (ATCC) and grown in DMEM with 10% FBS (Hyclone) at 37°C in 5% CO2 cell culture incubator. This cell line was tested 1 month before the experiment by methods of morphology check by microscopy, growth curve analysis, and mycoplasma detection according to the ATCC cell line verification test recommendations. The transfection of the plasmid expressing human MDA5 cDNA was performed using Lipofectamine 2000 transfection reagent (Invitrogen) according to the manufacturer′s instruction. Western blotting of proteins was performed as described previously42. The antibodies used included those against Flag and MDA5 were purchased from Sigma-Aldrich Co. and the antibody against β-actin were obtained from Abcam Co.
ELISA
IgG against MDA5 were detected in plasma samples using Anti-MDA5 Ab ELISA Kits (Medical & Biological Laboratories Co.), according to the manufacturer’s instructions. Briefly, during plasma samples incubation, anti-MDA5 antibodies in the plasma react with the immobilized human MDA5 antigen coated in the Microwells before. The test result is determined photometrically by measuring the absorbance (wave length: 450nm) and plotting the results performed on microplate reader (SpectraMax M5, MD, USA). The formula used for the unit value calculation is: Unit value (U/ml) =(A450 < Sample > - A450 < Calibrator 1> / (A450 < Calibrator 2 > - A450 < Calibrator 1 >)×100, which refers to the manufacturer’s instructions. The unit value ≥ 5.0 U/mL is considered positive and the unit value ≥ 10.0 U/mL is defined as high titer of anti-MDA5 Ab. The plasma from 5 cases of anti-MDA5 Ab-related DM were used as positive control.
Statistical analysis
For the detection of anti-MDA5 Ab, each experiment was repeated 3 times. Unpaired, two-sided Mann-Whitney U-test was performed to compare two groups unless otherwise indicated (X2 test). For the clinical analysis of anti-MDA5 Ab, descriptive statistics (percentages, means, standard deviations [SDs], medians, interquartile [IQR]) were provided for describe baseline demographic and clinical characteristics. The comparison of demographic, clinical, laboratory characteristics and outcomes across anti-MDA5 Ab positive/negative and survival/non-survival subgroups was performed by the Chi-squared tests or analysis of variance as appropriate. All statistical analyses were performed using SPSS 16.0 software (SPSS Inc., Chicago, IL, USA). P-values <0.05 were considered statistically significant.
Data Availability
In this retrospective study of adult inpatients in three hospitals in China, we identified anti-MDA5 Ab in 48.2% of patients with COVID-19. The anti-MDA5 Ab positive patients tended to represent with severe disease (88.6% vs 66.9%, P<0.0001), and the titer of anti-MDA5 Ab was significantly elevated in the non-survivals (P=0.030) than that in the survivals. We also demonstrated that early profiling of anti-MDA5 Ab could distinguish severe patients from those with non-severe ones. Overall, our data suggest that early screening and serially monitoring of anti-MDA5 Ab might help identify severe COVID-19 patients with poor outcomes, or those who would probably benefit from anti-inflammation and immunosuppressive therapy. The findings are yet to be validated in a larger population of different ethnicities in future.
Author contribution
All data were checked by two investigators (QW and CL) and a third researcher (BC)
Competing Interests statement
All the authors have no disclosures.
Acknowledgements
We thank Dr. Jinmin Peng and Jingjing Bai for providing assistance with statistical analysis and preparing samples of anti-MDA5 Ab-related DM. We are grateful to clinicians of Wuhan Jinyintan Hospital for sample collection and Dr. colleagues of Hubei Provincial Center for Disease Control and Prevention for sample transportation. This study was supported by National Key R&D Program of China (2020YFA0707600) and CAMS Innovation Fund for Medical Sciences (2018-I2M-1-003, 2019-I2M-1-006, 2019-I2M-2-008), National Science Grant for Distinguished Young Scholars (81425001/H0104), The Beijing Science and Technology Project (D151100002115004)