- 1Department of Ultrasound Medicine, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 2Hubei Province Key Laboratory of Molecular Imaging, Wuhan, China
- 3Department of Gynecology and Obstetrics, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 4British Heart Foundation Centre of Research Excellence, King's College London, London, United Kingdom
Background: Lung injury is a common condition among hospitalized patients with coronavirus disease 2019 (COVID-19). However, whether lung ultrasound (LUS) score predicts all-cause mortality in patients with COVID-19 is unknown. The aim of the present study was to explore the predictive value of lung ultrasound score for mortality in patients with COVID-19.
Methods: Patients with COVID-19 who underwent lung ultrasound were prospectively enrolled from three hospitals in Wuhan, China between February 2020 and March 2020. Demographic, clinical, and laboratory data were collected from digital patient records. Lung ultrasound scores were analyzed offline by two observers. Primary outcome was in-hospital mortality.
Results: Of the 402 patients, 318 (79.1%) had abnormal lung ultrasound. Compared with survivors (n = 360), non-survivors (n = 42) presented with more B2 lines, pleural line abnormalities, pulmonary consolidation, and pleural effusion (all p < 0.05). Moreover, non-survivors had higher global and anterolateral lung ultrasound score than survivors. In the receiver operating characteristic analysis, areas under the curve were 0.936 and 0.913 for global and anterolateral lung ultrasound score, respectively. A cutoff value of 15 for global lung ultrasound score had a sensitivity of 92.9% and specificity of 85.3%, and 9 for anterolateral score had a sensitivity of 88.1% and specificity of 83.3% for prediction of death. Kaplan–Meier analysis showed that both global and anterolateral scores were strong predictors of death (both p < 0.001). Multivariate Cox regression analysis showed that global lung ultrasound score was an independent predictor (hazard ratio, 1.08; 95% confidence interval, 1.01–1.16; p = 0.03) of death together with age, male sex, C-reactive protein, and creatine kinase-myocardial band.
Conclusion: Lung ultrasound score as a semiquantitative tool can be easily measured by bedside lung ultrasound. It is a powerful predictor of in-hospital mortality and may play a crucial role in risk stratification of patients with COVID-19.
Background
The coronavirus disease 2019 (COVID-19) is a newly recognized infectious disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Although chest computed tomography (CT) has been regarded as an important diagnostic tool for COVID-19 diagnosis (1), it is limited by high cost, radiation exposure, infection control challenges, and lack of continuous monitoring, particularly for critically ill patients (2). Lung ultrasound (LUS), with the advantage of being non-invasive, low cost, and radiation free, has been increasingly used as a bed-side tool for evaluation and monitoring of lung diseases, particularly in the intensive care unit (ICU) (2, 3). It was found to have high accuracy in diagnosing viral community-acquired pneumonia with 94% sensitivity and 89% specificity for the detection of viral pneumonia in symptomatic patients (4). Global LUS score, a semiquantitative numerical score of lung aeration across 12 lung regions, has been shown as a useful tool to diagnose acute respiratory distress syndrome (ARDS) (5).
We therefore hypothesized that LUS score may play an important role in detecting lung lesions and optimizing risk stratification in patients with COVID-19. To test this hypothesis, LUS images in patients prospectively recruited from three hospitals in Wuhan, China were analyzed to evaluate the prognostic value of LUS score for in-hospital mortality in patients with COVID-19.
Materials and Methods
Patient Population
Patients with confirmed COVID-19 who underwent lung ultrasound were consecutively recruited from the West Branch of Wuhan Union Hospital, Cancer Centre of Union Hospital, and Jianghan Mobile Cabin Hospital Wuhan, China between February 6, 2020 and March 15, 2020. The study was approved by the ethics committee, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology (No. 20200021). Written informed consent was waived because of the unprecedented nature of COVID-19 pandemic.
Inclusion criteria were age ≥18 years and confirmed COVID-19. Exclusion criteria were incomplete image acquisition, missing clinical data, and cardiac failure causing cardiogenic pulmonary oedema.
Demographic, clinical history, comorbidities, laboratory data, and outcomes of all patients were obtained from electronic medical records (Dthealth Medical Systems CO, Tianjin, China). Primary outcome was all-cause mortality. All patients were followed up until April 7, 2020 when the last patient in the study was discharged.
Lung Ultrasound
LUS examinations were performed by nine qualified ultrasound doctors using Mindray M9 potable ultrasound machines (Mindray Bio-medical electronics Co, Shenzhen, China) with 1- to 5- MHz convex probes. LUS consisted of 12 different regions (two anterior, two lateral, and two posterior thoracic regions) (Supplementary Figure 1) as previously described (6). All video files were recorded in a hospital local archive and were interpreted and scored offline by two experienced observers within 24 h of LUS examinations who were blinded to the clinical data and outcomes. In case of disagreement between observers, the two observers agreed by consensus on the LUS score.
Examples of ultrasound findings including the patterns of B lines, consolidations, pleural line abnormalities, pleural effusion, and the lesion distribution are shown in Figure 1.
Figure 1. Ultrasonographic features and lung ultrasound (LUS) score in patients with coronavirus disease 2019 (COVID-19). (A) Normal: the presence of A lines beyond the pleural line characterizes mornal pulmonay aeration, LUS score: 0. (B) B1 line: the presence of multiple vertical B lines (comet tails) with well-defined spacing regularly spaced B lines 7 mm apart, LUS score: 1. (C) B2 line: the presence of coalescent B lines <3 mm apart, LUS score: 2. (D) Lung consolidation: the presence of a tissue pattern (yellow arrowhead), LUS score: 3. (E) Pleural effusion at costophrenic angle (red arrowhead).
Lung Ultrasound Score
LUS score was determined based on four lung patterns (Supplementary Table 1): N = 0, B1 = 1, B2 = 2, and C = 3 as described previously (7):
a. N pattern—normal aeration: A lines or <3 isolated B lines;
b. B1 pattern—moderate loss of lung aeration: a clear number of multiple visible B lines with horizontal spacing between adjacent B lines ≤7 mm (B1 lines);
c. B2 pattern—severe loss of lung aeration: multiple B lines fused together with horizontal spacing between adjacent B lines ≤3 mm, including “white lung” (B2 lines); and
d. C pattern—complete loss of aeration: pulmonary consolidation, presence of tissue pattern accompanied by static or dynamic air bronchograms.
Global LUS score was calculated by summing the scores of all 12 lung regions (ranging from 0 to 36). An adjusted composite score, antero-lateral score, was also derived by summing the anterior and lateral regional scores (range from 0 to 24) (5, 7).
Repeatability and Reproducibility of Lung Ultrasound Score
Intra- and interobserver variability of global LUS score was assessed in 30 randomly selected subjects by repeat measurements on the same images 1 month apart by two observers. Bland–Altman plots were produced.
Statistical Analysis
Demographic, clinical, and outcome variables were presented as percentages for categorical variables and as medians with interquartile ranges (IQRs) for continuous variables. The Mann–Whitney U-test was used to compare LUS scores between survivors and non-survivors.
Receiver operating characteristic (ROC) curves for death were drafted for global and anterolateral score. The area under the receiver operating characteristic curves (AUCs) was calculated to determine the diagnostic accuracy for death. The optimal cutoffs were determined as the highest Youden's index (sensitivity + specificity – 1).
Kaplan–Meier curves were used to examine cumulative death rate, and differences between groups were tested using a log rank test. Univariate and multivariate Cox regression analysis was performed to identify potential predictors of death. Multivariate models were constructed to assess the prognostic utility of global and anterolateral scores, incorporating covariables that were significant (p < 0.05) in the univariate analysis. All statistical analyses were performed using SPSS version 25 (SPSS Inc. Chicago, Illinois).
Results
Patient Characteristics
A total of 407 patients with COVID-19 meeting the inclusion criteria were recruited, of whom 5 were excluded due to suboptimal LUS images (n = 3) and congestive heart failure (n = 2) (Figure 2). Four hundred two patients were included in the final analysis, of whom 42 died with median time to death 21 (IQR, 14–29) days. Cause of death was recorded as multiorgan failure (42.9%), respiratory failure (26.1%), cardiac (9.5%), septic shock (9.5%), unknown (7.1%), and stroke (4.8%). Baseline characteristics are summarized in Table 1. Non-survivors were older and more male gender compared to survivors. There was a higher prevalence of preexisting conditions including hypertension, coronary heart disease (CHD), and malignancy in non-survivors compared to survivors.
Laboratory Findings
Laboratory data on hospital admission are summarized in Table 2. Overall, non-survivors had significant worse laboratory results, including increased white blood cell count, prothrombin time, activated partial thromboplastin time, D-dimer, aspartate aminotransferase, total bilirubin, blood urea nitrogen, creatinine, hypersensitive troponin I (hs-Tnl), lactate dehydrogenase (LDH), creatine kinase–myocardial band (CK-MB), hypersensitive C-reactive protein (CRP), and procalcitonin and decreased lymphocyte count, platelet count, hemoglobin, total protein, and albumin (all p < 0.05) compared to survivors. Patients with a higher global LUS score (>15) had significant worse laboratory results, in particular, significantly increased D-dimer and CRP compared to those with a global LUS score <15.
Lung Ultrasound Findings and Lung Ultrasound Score
Lung ultrasound was performed within a median of 3 (IQR, 2–5) days from hospital admission. Lung ultrasound findings are shown in Table 3. Eighty-four patients (20.9%) had normal LUS. The presence of B lines was the most common finding (318/402, 79.1%), followed by pleural line abnormalities (137/402, 31.8%) and consolidation (117/402, 25.6%). Pleural effusions were detected in 36 (8.2%) patients. Compared to survivors, non-survivors were more likely to have B2 lines, pleural line abnormalities, pulmonary consolidation, and pleural effusion, but there was no difference in the presence of B1 lines. All non-survivors had bilateral involvement. Survivors had significantly lower global and anterior–lateral LUS scores compared to non-survivors (Figure 3). Findings of each of 12 lung regions are shown in Supplementary Figure 2. Regional LUS scores including anterior, lateral, and posterior scores are presented in Supplementary Figure 3. Bland–Altman plots for intra- and interobserver variability of global LUS score are shown in Supplementary Figure 4. All repeated measures were within 1.96 × standard deviation of the mean, which suggested a good reproducibility of global LUS score.
Figure 3. Box plots representation of (A) global and (B) anterolateral lung ultrasound (LUS) scores in survivors and non-survivors.
Prediction of Mortality by LUS Global and Anterolateral Score
After a median of 27 (IQR, 20–39) days of follow-up, 42 patients died. ROC curve analyses of global and anterolateral LUS score for predicting mortality are shown in Figure 4. The area under the curve were 0.936 and 0.913 for global and anterolateral LUS score, respectively. A cutoff value of 15 for global LUS score had a sensitivity of 92.9% and specificity of 85.3% for prediction of death, and a cutoff value of 9 for anterolateral LUS score had a sensitivity of 88.1% and specificity of 83.3%. Clinical characteristics and laboratory findings dichotomized according to global LUS score optimal value of 15 are shown in Tables 1, 2.
Figure 4. Receiver operating characteristic (ROC) curves of global and anterolateral lung ultrasound (LUS) score for prediction of death.
Kaplan–Meier analysis showed that both global and anterolateral LUS scores were strong predictors of death (Figure 5). When global LUS score was >15, 37/92 (40.2%) patients died compared to only 5/310 (1.6%) death in those with a global LUS score <15. When patients were dichotomized by anterolateral LUS score of 9, there were 36/97 (37.1%) deaths in patients with a high score compared to 6/305 (2.0%) deaths in those with a low anterolateral score.
Figure 5. Kaplan–Meier curves of (A) global lung ultrasound (LUS) score with optimal cutoff value of 15 and (B) anterolateral LUS score with optimal cutoff value of 9 for prediction of in-hospital mortality.
On univariate Cox regression analysis, age, male gender, malignancy, CHD, CRP, hs-cTnl, CK-MB, D-dimer, global LUS score, and anterolateral LUS score were significantly associated with mortality (Table 4). In multivariate model 1, considering global LUS score together with other significant predictors in the univariate model, age, male sex, CRP, CK-MB, and global LUS score [hazard ratio (HR), 1.08; 95%CI, 1.01–1.16, p = 0.03) remained as a significant predictor. In multivariate model 2, when anterolateral LUS score was tested with other variables, the predictive power of anterolateral LUS score did not remain significant.
Discussion
Our data suggested that global LUS score was a predictor of in-hospital mortality independent of age, gender, comorbidities, and biochemical markers and was superior to LUS anterolateral score. The optimal threshold of 15 for global LUS score and 9 for anterolateral LUS score were in line with those derived from previous investigations (5, 8). These findings supported the clinical utility of LUS in patients with COVID-19 (7, 9) given its ease of use at point of care, low cost, lack of radiation exposure, and ready combination with other components of critical care ultrasonography (10, 11).
LUS features in patients with COVID-19 in our study manifested as multiple lesions, various types of B lines, irregularly pleural lines, and subpleural consolidations. B lines presented in 79.1% patients. B2 lines and consolidations were more common in non-survivors than in survivors. Pleural effusion, pleural thickening, and pneumothorax were less common in COVID-19 patients, which were consistent with the latest autopsy report (12) that COVID-19 patients presented with acute interstitial lung disease.
Bass et al. showed that LUS had high sensitivity for detection of interstitial and alveolar–interstitial lung disease with peripheral distribution (13). Consistent with these features, our findings suggested that global LUS score was highly predictive of death in COVID-19 and independent of other previously identified predictors. Non-survivors in our study were older and more male with higher prevalence of preexisting conditions including hypertension, CHD, and malignancy and higher levels of cardiac injury and systematic inflammation markers than survivors, which were in consistence with previous studies (14).
Another interesting finding of our study was that when the posterior regions were excluded, the predictive power of anterolateral LUS score disappeared in the multivariate cox regression model. This finding was consistent with chest CT findings that the most commonly involved lung segments in patients with COVID-19 were the dorsal segment of the right lower lobe, the posterior basal segment of the right lower lobe, the lateral basal segment of the right lower lobe, and the dorsal segment and the posterior basal segment of the left lower lobe (15). Despite some studies showing that the posterior regions had the lowest diagnostic accuracy (5), scores from these regions could play an important role in risk stratification. In the present study, lung lesions were mainly located in the right lateral inferior area (39.8%), left lateral inferior area (38.1%), left posterior inferior area (38.8%), and right posterior inferior area (37.3%) (the lower posterior and lateral segments of the lungs). This finding also supported that the potential benefit of prone position in patients affected with COVID-19 acute respiratory distress syndrome (ARDS) due to a more even distribution of the gas–tissue ratios along the dependent–non-dependent axis and a more homogeneous distribution of lung stress and strain (16).
Although anterolateral LUS score had less predive power compared to global LUS score, it may still play an important role particularly in patients on ICU.
Clinical Implications
COVID-19 as a global pandemic imposes a huge burden on medical systems. Early quantification of patients with severe lung involvement may be critical for optimization of treatment and management. LUS as a non-invasive and cost-effective diagnostic tool can be performed rapidly, particularly in ICU. Severe studies have also demonstrated that echocardiography is a crucial tool in detecting cardiovascular complications (in particular on assessment of left and right ventricular function) and predicts poor prognosis in patients COVID-19 (17, 18). Combining LUS with echocardiography may add additional value to identify patients at higher risk of poor outcomes.
Limitations
Our study has several limitations. First, mortality rate was relatively low, which limits the strength of our conclusion. Low mortality rate may be due to the fact that majority of patients in the present study were not in ICU, while this rate was similar to previously published data (19). Second, the follow-up period was relatively short, as majority of patients were discharged within 28 days from admission.
Although our findings suggested that LUS may add additional value in risk stratification, the strength of our conclusion may be limited by the nature of an observational study. There are several other limitations of LUS that cannot be ignored such as the requirement of special training to perform high-quality LUS, lack of evidence-based guidelines, the high risk of infection when performing LUS examination in patients with COVID-19 (20, 21).
Patients included in this study were recruited from three hardest-hit hospitals in Wuhan, and these patients may not represent the population in other areas. Finally, the relationship between LUS and lung CT was not explored, as the majority of patients did not have lung CT due to limited availabilities and the nature of infectious disease.
Conclusion
Global LUS score as a semiquantitative measure of lung conditions is a powerful predictor of in-hospital mortality in patients with COVID-19 and may add additional value in patient monitoring and risk stratification.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by Union Hospital, Tongji Medical College, Huazhong University of Science and Technology. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
ZS, HW, HG, and MX conceived and designed the study. ZZ and JL contributed to the literature search. ZS, CC, SQ, YS, YD, WZ, MY, and LJ contributed to data collection. ZZ, CY, and HG contributed to data analysis. JW, YY, QL, and HG contributed to data interpretation. YX and RW contributed to the figures. ZS, HG, ZZ, and JL drafted the article. All authors contributed to the article and approved the submitted version.
Funding
We acknowledge the financial support from the National Natural Science Foundation of China (grant 81922033 to Li Zhang; grant 81727805 to MX; and grant 81701716 to ZS) and National Institute for Health Research, UK (ICA-CL-2018-04-ST2-012) to HG and by British Heart Foundation, UK (PG/19/23/34259) to HG.
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.
Acknowledgments
We are grateful for all our colleagues, especially Dr. Li Zhang and Hongliang Yuan, for their support to the present study. We are also grateful to the many frontline medical staffs for their dedication in the face of this outbreak despite the potential threat to their own lives and the lives of their families.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2021.633539/full#supplementary-material
References
1. Shi H, Han X, Jiang N, Cao Y, Alwalid O, Gu J, et al. Radiological findings from 81 patients with COVID-19 pneumonia in Wuhan, China: a descriptive study. Lancet Infect Dis. (2020) 20:425–34. doi: 10.1016/S1473-3099(20)30086-4
2. Lepri G, Orlandi M, Lazzeri C, Bruni C, Hughes M, Bonizzoli M, et al. The emerging role of lung ultrasound in COVID-19 pneumonia. Eur J Rheumatol. (2020) 7:129–33. doi: 10.5152/eurjrheum.2020.2063
3. Volpicelli G, Lamorte A, Villen T. What's new in lung ultrasound during the COVID-19 pandemic. Intensive Care Med. (2020) 46:1445–48. doi: 10.1007/s00134-020-06048-9
4. Zhang YK, Li J, Yang JP, Zhan Y, Chen J. Lung ultrasonography for the diagnosis of 11 patients with acute respiratory distress syndrome due to bird flu H7N9 infection. Virol J. (2015) 12:176. doi: 10.1186/s12985-015-0406-1
5. Pisani L, Vercesi V, van Tongeren PSI, Lagrand WK, Leopold SJ, Huson MAM, et al. The diagnostic accuracy for ARDS of global versus regional lung ultrasound scores - a post hoc analysis of an observational study in invasively ventilated ICU patients. Intensive Care Med Exp. (2019) 7:44. doi: 10.1186/s40635-019-0241-6
6. Bouhemad B, Mongodi S, Via G, Rouquette I. Ultrasound for “lung monitoring” of ventilated patients. Anesthesiology. (2015) 122:437–47. doi: 10.1097/ALN.0000000000000558
7. Volpicelli G, Elbarbary M, Blaivas M, Lichtenstein DA, Mathis G, Kirkpatrick AW, et al. International evidence-based recommendations for point-of-care lung ultrasound. Intensive Care Med. (2012) 38:577–91. doi: 10.1007/s00134-012-2513-4
8. Tierney DM, Boland LL, Overgaard JD, Huelster JS, Jorgenson A, Normington JP, et al. Pulmonary ultrasound scoring system for intubated critically ill patients and its association with clinical metrics and mortality: a prospective cohort study. J Clin Ultrasound. (2018) 46:14–22. doi: 10.1002/jcu.22526
9. Nazerian P, Volpicelli G, Vanni S, Gigli C, Betti L, Bartolucci M, et al. Accuracy of lung ultrasound for the diagnosis of consolidations when compared to chest computed tomography. Am J Emerg Med. (2015) 33:620–5. doi: 10.1016/j.ajem.2015.01.035
10. Peng QY, Wang XT, Zhang LN. Findings of lung ultrasonography of novel corona virus pneumonia during the 2019-2020 epidemic. Intensive Care Med. (2020) 46:849–50. doi: 10.1007/s00134-020-05996-6
11. Vetrugno L, Bove T, Orso D, Barbariol F, Bassi F, Boero E, et al. Our Italian experience using lung ultrasound for identification, grading and serial follow-up of severity of lung involvement for management of patients with COVID-19. Echocardiography. (2020) 37:625–27. doi: 10.1111/echo.14664
12. Xu Z, Shi L, Wang Y, Zhang J, Huang L, Zhang C, et al. Pathological findings of COVID-19 associated with acute respiratory distress syndrome. Lancet Respir Med. (2020) 8:420–22. doi: 10.1016/S2213-2600(20)30076-X
13. Bass CM, Sajed DR, Adedipe AA, West TE. Pulmonary ultrasound and pulse oximetry versus chest radiography and arterial blood gas analysis for the diagnosis of acute respiratory distress syndrome: a pilot study. Crit Care. (2015) 19:282. doi: 10.1186/s13054-015-0995-5
14. Luo X, Zhou W, Yan X, Guo T, Wang B, Xia H, et al. Prognostic value of C-reactive protein in patients with COVID-19. Clin Infect Dis. (2020) 71:2174–9. doi: 10.1093/cid/ciaa641
15. Wu J, Wu X, Zeng W, Guo D, Fang Z, Chen L, et al. Chest CT findings in patients with coronavirus disease 2019 and its relationship with clinical features. Invest Radiol. (2020) 55:257–61. doi: 10.1097/RLI.0000000000000670
16. Guérin C, Albert RK, Beitler J, Gattinoni L, Jaber S, Marini JJ, et al. Prone position in ARDS patients: why, when, how and for whom. Intensive Care Med. (2020) 46:2385–96. doi: 10.1007/s00134-020-06306-w
17. Moody WE, Mahmoud-Elsayed HM, Senior J, Gul U, Khan-Kheil AM, Horne S, et al. Impact of right ventricular dysfunction on mortality in patients hospitalized with COVID-19, according to race. CJC Open. (2021) 3:91–100. doi: 10.1016/j.cjco.2020.09.016
18. Cameli M, Pastore MC, Soliman Aboumarie H, Mandoli GE, D'Ascenzi F, et al. Usefulness of echocardiography to detect cardiac involvement in COVID-19 patients. Echocardiography. (2020) 37:1278–86. doi: 10.1111/echo.14779
19. Chen N, Zhou M, Dong X, Qu J, Gong F, Han Y, et al. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study. Lancet. (2020) 395:507–13. doi: 10.1016/S0140-6736(20)30211-7
20. Gargani L, Soliman-Aboumarie H, Volpicelli G, Corradi F, Pastore MC, Cameli M. Why, when, and how to use lung ultrasound during the COVID-19 pandemic: enthusiasm and caution. Eur Heart J Cardiovasc Imaging. (2020) 21:941–8. doi: 10.1093/ehjci/jeaa163
Keywords: COVID-19, SARS-CoV-2, lung ultrasound score, mortality, prognosis
Citation: Sun Z, Zhang Z, Liu J, Song Y, Qiao S, Duan Y, Cao H, Xie Y, Wang R, Zhang W, You M, Yu C, Ji L, Cao C, Wang J, Yang Y, Lv Q, Wang H, Gu H and Xie M (2021) Lung Ultrasound Score as a Predictor of Mortality in Patients With COVID-19. Front. Cardiovasc. Med. 8:633539. doi: 10.3389/fcvm.2021.633539
Received: 25 November 2020; Accepted: 22 March 2021;
Published: 25 May 2021.
Edited by:
Matteo Cameli, University of Siena, ItalyReviewed by:
Maria Concetta Pastore, Università del Piemonte Orientale, ItalyCiro Santoro, Federico II University Hospital, Italy
Copyright © 2021 Sun, Zhang, Liu, Song, Qiao, Duan, Cao, Xie, Wang, Zhang, You, Yu, Ji, Cao, Wang, Yang, Lv, Wang, Gu and Xie. 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: Haotian Gu, aGFvdGlhbi5ndSYjeDAwMDQwO2tjbC5hYy51aw==; Hongbo Wang, ZHJ3YW5naGI2OSYjeDAwMDQwO3NpbmEuY29t; Mingxing Xie, eGllbXgmI3gwMDA0MDtodXN0LmVkdS5jbg==
†These authors share first authorship
‡These authors share senior authorship