key: cord-0735542-ytk0r4vp authors: Yan, Danying; Zhang, Xiaobao; Chen, Can; Jiang, Daixi; Liu, Xiaoxiao; Zhou, Yuqing; Huang, Chenyang; Zhou, Yiyi; Guan, Zhou; Ding, Cheng; Chen, Lu; Lan, Lei; Fu, Xiaofang; Wu, Jie; Li, Lanjuan; Yang, Shigui title: Characteristics of Viral Shedding Time in SARS-CoV-2 Infections: A Systematic Review and Meta-Analysis date: 2021-03-19 journal: Front Public Health DOI: 10.3389/fpubh.2021.652842 sha: 7176e3eea891315744a6f8d6e5886941edce7d3d doc_id: 735542 cord_uid: ytk0r4vp Background: The viral shedding time (VST) of SARS-CoV-2 mainly determines its transmission and duration of infectiousness. However, it was heterogeneous in the existing studies. Here, we performed a meta-analysis to comprehensively summarize the VST of SARS-CoV-2. Methods: We searched PubMed, Web of Science, MedRxiv, BioRxiv, CNKI, CSTJ, and Wanfang up to October 25, 2020, for studies that reported VSTs of SARS-CoV-2. Pooled estimates and 95% CIs for the VSTs were calculated using log-transformed data. The VSTs in SARS-CoV-2 infections based on different demographic and clinical characteristics, treatments and specimens were stratified by subgroup analysis. Results: A total of 35 studies involving 3,385 participants met the inclusion criteria. The pooled mean VST was 16.8 days (95% CI: 14.8–19.4, I(2) = 99.56%) in SARS-CoV-2 infections. The VST was significantly longer in symptomatic infections (19.7 days, 95% CI: 17.2–22.7, I(2) = 99.34%) than in asymptomatic infections (10.9 days, 95% CI: 8.3–14.3, I(2) = 98.89%) (P < 0.05). The VST was 23.2 days (95% CI: 19.0–28.4, I(2) = 99.24%) in adults, which was significantly longer than that in children (9.9 days, 95% CI: 8.1–12.2, I(2) = 85.74%) (P < 0.05). The VST was significantly longer in persons with chronic diseases (24.2 days, 95% CI: 19.2–30.2, I(2) = 84.07%) than in those without chronic diseases (11.5 days, 95% CI: 5.3–25.0, I(2) = 82.11%) (P < 0.05). Persons receiving corticosteroid treatment (28.3 days, 95% CI: 25.6–31.2, I(2) = 0.00%) had a longer VST than those without corticosteroid treatment (16.2 days, 95% CI: 11.5–22.5, I(2) = 92.27%) (P = 0.06). The VST was significantly longer in stool specimens (30.3 days, 95% CI: 23.1–39.2, I(2) = 92.09%) than in respiratory tract specimens (17.5 days, 95% CI: 14.9–20.6, I(2) = 99.67%) (P < 0.05). Conclusions: A longer VST was found in symptomatic infections, infected adults, persons with chronic diseases, and stool specimens. Coronaviruses (CoVs), belonging to Nidovirales order, have caused three global outbreaks in the past 20 years. The first epidemic was Severe Acute Respiratory Syndrome (SARS) caused by SARS-CoV-1 in 2003, the second outbreak was Middle East Respiratory Syndrome (MERS) caused by MERS-CoV in 2012, and the third and most recent pandemic was Coronavirus Disease 2019 (COVID-19) caused by SARS-CoV-2 (1, 2). As of February 4, 2021, more than 104 million cases of COVID-19 have been reported with over 2.2 million deaths globally (3) . Pulmonary clinical manifestations are the most common clinical presentations of COVID-19, such as fever, cough, shortness of breath, sputum production, respiratory failure and even acute respiratory distress syndrome (ARDS). Diarrhea, loss of smell or taste, and other extra-pulmonary clinical manifestations can also be found in some patients (4) (5) (6) (7) . Persons infected with SARS-CoV-2 with long viral shedding times (VSTs) have drawn considerable concern, which put greater challenges and difficulties on epidemic prevention and control (8) (9) (10) (11) . The VST is an important parameter for judging hospital discharge, discontinuation of quarantine and the effect of antiviral treatment for infectious diseases, which mainly determines disease transmission and the duration of infectiousness (12) . However, the characteristics of the VST in SARS-CoV-2 infections have not been well-clarified. Although there have been many studies on the VSTs of SARS-CoV-2, the results across studies so far have been heterogeneous (13, 14) . A meta-analysis performed by Muge Cevik found that the mean VST of SARS-CoV-2 in the upper respiratory tract, lower respiratory tract, stool and serum was 17.0, 14.6, 17.2, and 16.6 days, respectively (15) . However, a comprehensive summary of VSTs in SARS-CoV-2 infections with different demographic and clinical features is still lacking. Therefore, we performed a metaanalysis to estimate the mean VST in SARS-CoV-2 infections and explore the characteristics of VSTs in SARS-CoV-2 infections based on different demographic features, clinical characteristics, treatments and specimens. Our meta-analysis was strictly conducted in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses Protocols (PRISMA-P) guidelines (16) . The definition of the VST varied among the studies, so a unified definition was made. We defined the VST as the time from illness onset to viral shedding cessation. Illness onset was defined as the first appearance of the symptoms for symptomatic infections and the first positive RT-PCR results for asymptomatic infections. Viral shedding cessation referred to the occurrence of the last positive RT-PCR results or negative RT-PCR results. We searched PubMed, Web of Science, MedRxiv, BioRxiv, the China National Knowledge Infrastructure Database (CNKI), the In this systematic review with no study design limit, studies meeting the following inclusion criteria were eligible: (i) SARS-CoV-2 infections were based on positive RT-PCR results; (ii) the VSTs of SARS-CoV-2 infections including sample size, mean and standard deviation (SD) could be obtained directly from the original studies or by calculation; and (iii) the definition of the VST in the original studies was consistent with our definition. We excluded (i) duplicated data; (ii) case reports and case series with <5 participants due to reporting bias; and (iii) studies without original data (e.g., modeling studies and reviews). Studies presenting VSTs with medians and interquartile ranges (IQRs) or ranges were excluded to reduce the errors caused by data conversion. After removing duplicates, two reviewers (DY and XZ) independently performed the initial screening of titles and abstracts to exclude studies that clearly contained no data for VST of SARS-CoV-2. All retained full-text articles were scrutinized against the eligibility criteria by two independent reviewers (DY and XZ). Nine investigators (DY, XZ, CC, DJ, XL, YZ, CH, YZ, and ZG) participated in the data extraction. And data extraction from each study was performed by three independent investigators. Disagreements and uncertainties were consulted by SY to reach a consensus. The following data were extracted: basic information of the studies (first author, publication time, journal name, sample size), VSTs in SARS-CoV-2 infections based on sex, age (adult and child), infection status (symptomatic infection and asymptomatic infection), disease severity (severe infection and non-severe infection), treatments (corticosteroid treatment and antiviral therapy) and specimens (respiratory tract specimens (RTS), upper respiratory tract specimens (URTS), lower respiratory tract specimens (LRTS), stool and serum). The cutoff point for classifying adults and children was 18 years old. Asymptomatic infections referred to the absence of any clinical symptoms throughout the disease course. Non-severe infections included mild and moderate infections, and severe infections included severe and critical infections. Antiviral drugs included interferon, lopinavir/ritonavir, abidor, ribavirin, chloroquine and hydroxychloroquine. URTS included nasopharyngeal, oropharyngeal and oronasopharyngeal swabs, and LRTS included sputum and bronchoalveolar lavage fluid. The overall VST of the total participants was extracted to estimate the overall pooled VST. If studies did not report the overall VST of the total participants, the stratified VSTs were extracted to estimate the overall pooled VST. If a study included more than one independent study population, each population was extracted as a separate dataset in the metaanalysis. When the same study reported the VSTs of multiple specimens, the VSTs of the URTS were extracted to estimate the overall pooled VST, and the VSTs of other specimens were displayed in the subgroup analysis. In subgroup analysis, only studies having clear population characteristics were included in the corresponding subgroup, and studies having no clear information or mixed population group were included in the unclassified group. The scale recommended by the Agency for Healthcare Research and Quality was used to assess the quality of the included studies (17) . The scale consists of 11 items, and 1 point is given to each item when the conditions are met. It mainly focuses on information source, inclusion and exclusion criteria, study period, selection of participants, evaluation of subjective outcomes/components, quality assurance, possible confounding variables, handling of missing data, participants' response rates and completeness of data collection. According to the total score, the studies were divided into low-(0-3), medium-(4-7) and highquality (8-11) groups. EndNote (version X9) was used to manage the articles and citations. We first extracted the individual VSTs from the published articles and found that the distribution type of the VST was approximately in accordance with the log-normal distribution by using P-P plots (Supplementary Figure 1) . Then, we used the method developed by McAloon C (18) to transform the original VST data to make the data obey a normal distribution. We used random-effects model to perform the meta-analysis due to the high heterogeneity. Finally, we used the method developed by McAloon C (18) to back-transform the point estimates and their 95% confidence intervals (CIs). The I-squared (I 2 ) statistic was used to evaluate the heterogeneity among the studies. Metaregression was used to quantify the sources of heterogeneity and to explore the level of significance between subgroup comparisons. We did not assess publication bias because usual appraisal methods are uninformative when studies in the metaanalysis do not include a test of significance. The data cleaning and analysis were performed using the Microsoft Excel 2016 and R version 3.2.3. A total of 17,284 records were retrieved through a database search. The titles and abstracts of 11,911 records were screened after deleting duplicates, and then 526 records were selected for full-text review. Finally, 35 full texts met the inclusion criteria (Figure 1 ). This study included 35 observational studies and involved 3,385 individuals infected with SARS-CoV-2, of which 2,955 were symptomatic infections and 338 were asymptomatic infections ( Table 1) . According to the scale, 32 studies were of high quality, 3 studies were of medium quality and none were of low quality ( Table 1 and Supplementary Table 2 ). The initial pooled estimate of the log-transformed VST in SARS-CoV-2 infections was 2.82 (95% CI: 2.69-2.96) (Figure 2) (Figure 3 ). The mean VST was 19.4 days (95% CI: 9.5-39.4) in females and 11.9 days (95% CI: 8.4-16.9) in males. The VST was significantly shorter in the infected children (9.9 days, 95% CI: 8.1-12.2) than in the infected adults (23.2 days, 95% CI: 19.0-28.4) (P < 0.05). The VST of persons with chronic diseases was 24.2 days (95% CI: 19.2-30.2), which was significantly longer than that of persons without chronic diseases (11.5 days, 95% CI: 5.3-25.0) (P < 0.05) (Figure 3 ). In persons receiving corticosteroid treatment, the VST was 28.3 days (95% CI: 25.6-31.2), which was longer than that in those without corticosteroid treatment (16.2 days, 95% CI: 11.5-22.5). However, there was no statistically significant difference between them (P = 0.06). The VST was 17. (Figure 3) . Only one study reported the VSTs of 5 patients without antiviral therapy, and the result was 11.2 ± 5.2 days (33). Most studies (63%) reported the VSTs in the URTS. Among the different specimens, the mean VST was 17.5 days (95% CI: 14.9-20.6) in the RTS and 17.5 days (95% CI: 14.6-21.0) in the URTS. Compared with the RTS, a longer VST was found in the stool specimens (30.3 days, 95% CI: 23.1-39.2) (P < 0.05) (Figure 4) . No included study reported VSTs in LRTS or serum specimens. The univariate meta-regression model indicated that the mean age (R 2 = 35.28%, P < 0.05) and the proportion of the asymptomatic cases (R 2 = 22.64%, P < 0.05) could partly explain the overall heterogeneity. By introducing these two variables into the multivariate meta-regression model, nearly half of We performed a meta-analysis to clarify the characteristics of VSTs in SARS-CoV-2 infections, which was important for determining hospital discharge, discontinuation of quarantine and the effect of antiviral treatment for COVID-19. Compared with the meta-analysis conducted by Muge Cevik (15), our study not only estimated the VSTs in different specimens but also summarized the VSTs in SARS-CoV-2 infections based on different demographic features, clinical characteristics and treatments. Previous studies have shown that the basic reproduction number (R0) of SARS-CoV-2 is between 2 and 6.7, which indicates that SARS-CoV-2 is more infectious than SARS-CoV-1 and MERS-CoV (55) (56) (57) . In our study, we found that the mean VST of SARS-CoV-2 was 16.8 days (95% CI: 14. 8-19.4) , which was between that of SARS-CoV-1 (21.0 days) and MERS-CoV (13.2 days) (58, 59) . In addition to the VST, the viral load released is also important to evaluate the transmissibility. Some studies have found that the viral load of SARS-CoV-2 is highest during the 1st week after symptom onset and subsequently declines with time (60) (61) (62) . Based on the above analysis, from the perspective of epidemic prevention and control, strict precautions should be taken throughout the disease course, especially within 1 week after the onset of the disease. The duration of viral shedding is mainly related to the host immune status (63) . Persons with chronic diseases always have relatively low immunity, which might lead to longer viral shedding. In our study, we found that the VST of symptomatic infections was longer than that of asymptomatic infections. One reason is that virus clearance in asymptomatic individuals is indeed faster than that in symptomatic cases (38, 52) . Another reason was that the VST for asymptomatic infections was calculated from the first positive PCR results and depended mainly on close contact tracking investigations. These individuals might have begun viral shedding before the first positive PCR results, and were ignored due to the absence of clinical features. A higher proportion of asymptomatic infections and milder clinical symptoms were found in infected children compared with infected adults (64, 65) , which might explain the shorter VST of children. The VST is an important parameter for evaluating the effect of antiviral treatment for infectious diseases. Until now, there have been no specific antiviral drugs for COVID-19, and inhibiting the cytokine storm has been an important treatment for patients with severe COVID-19. Corticosteroids are used because of their rapid, powerful anti-inflammatory effects. In our study, we found that the patients who received corticosteroid treatment had longer VSTs, although no statistically significant difference was found. This phenomenon was also found in severe SARS and MERS, where high-dose corticosteroids could cause prolonged viral clearance, secondary infection and long-term complications (66) . Although corticosteroids can inhibit lung inflammation and alleviate possible immune-mediated pulmonary damage, it can also inhibit the systemic immune response dominated by T cell response, resulting in the delayed virus clearance (67) . This finding alerted us that high-dose corticosteroids might prolong VSTs in SARS-CoV-2 infections and that appropriate doses of corticosteroids should be used after weighing the advantages and disadvantages according to the patients' condition. The VST is also an important parameter for determining hospital discharge and discontinuation of quarantine. Two consecutive negative PCR results of RTS are one of the current criteria for hospital discharge or discontinuation of quarantine in China (68). The overexpression of ACE-2 in the gastrointestinal (GI) epithelial cells suggested the replication and shedding of SARS-CoV-2 in GI tract (69) . Similar to SARS-CoV-1 (59), the VST of SARS-CoV-2 in stool specimens was longer than that in RTS. One study suggested that the VST in stool specimens could be prolonged by 5 weeks after SARS-CoV-2 had turned negative in RTS (35) . Given that, the negative PCR results in RTS might not guarantee that patients no longer shed virus. Recently, several incidents of cold chain food polluted by SARS-CoV-2 have caused widespread concern by indicating that the virus could indeed infect individuals by polluting the environment. Considering the potential risk of oral-fecal transmission (70) and the long VST in stool specimens, more comprehensive protective measures should be taken for high-risk groups of oral-fecal transmission, such as GI endoscopy staff (2, 71) , and stool or anal swabs collection and testing staff. Our results might provide scientific support for the formulation of antiviral treatment and criteria for hospital discharge and discontinuation of quarantine for COVID-19, and help identify which patients need more attention and more effective preventive measures. Based on the mean VST of SARS-CoV-2 infections, hospitals could estimate the number of individuals with COVID-19 who can be admitted in a period of time, and reasonably allocate medical resources, such as the number of beds and medical staff. This study has several limitations. The mean age, disease severity, treatment regimens, underlying diseases and infection status of individuals infected with SARS-CoV-2 varied in the included studies, which might cause high statistical heterogeneity. In the multivariate meta-regression model, nearly half of the heterogeneity could be explained by mean age and the proportion of the asymptomatic cases (R 2 = 44.18%, P < 0.05). In some subgroup analyses, the number of included studies was small and most were case series with limited sample sizes, which might make the effect size of some outcomes insufficient. For example, the pooled mean VST in the stool specimens was based on estimates obtained in only two studies. More studies on the VST of SARS-CoV-2 are needed to provide further evidence. It would be better to incorporate as many studies as possible to obtain sufficient subgroup data and to ensure the homogeneity of the studies. Furthermore, the day of symptom onset for symptomatic infections depended on subjective memories and the day of the first positive RT-PCR results for asymptomatic infections relied mainly on close contact tracking investigations. If the individuals' recall was incorrect or close contact tracking investigations were not timely, these would cause the obtained VSTs to deviate from the real values. This study provided a comprehensive overview of VSTs in SARS-CoV-2 infections, which was important for determining hospital discharge, discontinuation of quarantine and the effect of antiviral treatment for COVID-19. The pooled mean VST was 16.8 days (95% CI: 14. 8-19.4 ) in SARS-CoV-2 infections. Due to the high infectivity of SARS-CoV-2, strict precautions should be taken to reduce the risk of disease transmission, especially for adults, persons with chronic diseases, symptomatic infections and persons with positive RT-PCR results in stool specimens, in whom longer VSTs were found. Given that highdose corticosteroids could alleviate possible immune-mediated pulmonary damage but might prolong VSTs in SARS-CoV-2 infections, corticosteroids should be used with caution after analyzing the risk of prolonged VST with reducing the disease severity. The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s. SY, LJL, and JW designed the study and revised the manuscript. DY and XZ independently performed the literature identification and uncertainties were consulted by SY to reach a consensus. DY, XZ, CC, DJ, XL, YZ, CH, YZ, and ZG extracted data. CD, LC, LL, and XF rechecked the data. DY carried out the data analysis. DY and XZ interpreted data and wrote the manuscript. All authors have read and approved the final version of the manuscript for submission. Clinical characteristics, diagnosis, and treatment of major coronavirus outbreaks COVID-19 and gastrointestinal endoscopies: current insights and emergent strategies COVID-19 extrapulmonary illness -special gastrointestinal and hepatic considerations Pulmonary and extra-pulmonary clinical manifestations of COVID-19 The association of "loss of smell" to COVID-19: a systematic review and meta-analysis Hyperlipasemia and potential pancreatic injury patterns in COVID-19: a marker of severity or innocent bystander? A special case of COVID-19 with long duration of viral shedding for 49 days. medRxiv Persistent viral RNA positivity during the recovery period of a patient with SARS-CoV-2 infection Prolonged shedding of severe acute respiratory syndrome coronavirus 2 in patients with COVID-19 Prolonged viral shedding in feces of pediatric patients with coronavirus disease 2019 MERS coronavirus: data gaps for laboratory preparedness Clinical characteristics and factors associated with long-term viral excretion in patients with severe acute respiratory syndrome coronavirus 2 infection: a single-center 28-day study Epidemiology and clinical course of COVID-19 in Shanghai SARS-CoV-2, SARS-CoV, and MERS-CoV viral load dynamics, duration of viral shedding, and infectiousness: a systematic review and meta-analysis The PRISMA extension statement for reporting of systematic reviews incorporating network meta-analyses of health care interventions: checklist and explanations Introduction from the agency for healthcare research and quality Incubation period of COVID-19: a rapid systematic review and meta-analysis of observational research A case series of children with 2019 novel coronavirus infection: clinical and epidemiological features A comparison between Chinese children infected with coronavirus disease-2019 and with severe acute respiratory syndrome A threegeneration family cluster with COVID-19 infection: should quarantine be prolonged? Public Health Asymptomatic infection and atypical manifestations of COVID-19: Comparison of viral shedding duration Clinical characteristics and outcomes of COVID-19 patients with gastrointestinal symptoms admitted to Jianghan Fangcang Shelter Hospital in Wuhan Clinical course and molecular viral shedding among asymptomatic and symptomatic patients with SARS-CoV-2 infection in a community treatment center in the Republic of Korea Clinical features of pediatric patients with coronavirus disease (COVID-19) Clinical progression of patients with COVID-19 in Shanghai Coronavirus disease 2019 pneumonia in immunosuppressed renal transplant recipients: a summary of 10 confirmed cases in Wuhan, China Digestive symptoms in COVID-19 patients with mild disease severity: clinical presentation, stool viral RNA testing, and outcomes Duration of SARS-CoV-2 viral RNA in asymptomatic carriers Effects of methylprednisolone use on viral genomic nucleic acid negative conversion and CT imaging lesion absorption in COVID-19 patients under 50 years old Effects of oral care on prolonged viral shedding in coronavirus disease 2019 (COVID-19). Special Care Dentistr Epidemiological and clinical characteristics of 26 asymptomatic severe acute respiratory syndrome coronavirus 2 carriers Epidemiological features and viral shedding in children with SARS-CoV-2 infection Nasopharyngeal SARS-CoV-2 viral RNA shedding in patients with diabetes mellitus Prolonged presence of SARS-CoV-2 viral RNA in faecal samples Prolonged shedding of SARS-CoV-2 in COVID-19 infected hemodialysis patients Viral transmission and clinical features in asymptomatic carriers of SARS-CoV-2 in Wuhan Early viral clearance and antibody kinetics of COVID-19 among asymptomatic carriers Familial clustering infection caused by asymptomatic COVID-19 (in Chinese) Clinical features of asymptomatic or subclinical COVID-19 in children (in Chinese) Influencing factors of negative conversion time of viral nucleic acid of nasopharyngeal swabs in patients with novel coronavirus infected pneumonia (in Chinese) Traditional Chinese medicine characteristics of 46 patients with COVID-19 (in Chinese) Clinical ananlysis of corona virus disease 2019: a report of 24 cases (in Chinese) Related factors of the conversion time of virus nucleic acid turning negative in patients with coronavirus disease 2019 and its effect on prognosis (in Chinese) Differences in laboratory indicators and prognosis between COVID-19 patients with or without diabetes mellitus (in Chinese) Retrospective analysis of viral nucleic acid test results of 116 fecal samples from COVID-19 patients (in Chinese) Clinical characteristics of 59 cases of COVID-19 in Guangzhou (in Chinese) Dynamic changes of IgM and IgG in confirmed COVID-19 cases: detection results analysis Nucleic acid detection analysis of 71 cases of COVID-19 with biological samples of diferent disease course(in Chinese) Clinical characteristics of 2019 novel coronavirus pneumonia in a mobile cabin hospital in Wuhan: analysis of 94 cases (in Chinese) Analysis of Factors influencing viral shedding time in severe and critical COVID-19 patients (in Chinese) Symptomatic infection is associated with prolonged duration of viral shedding in mild coronavirus disease 2019: a retrospective study of 110 children in Wuhan Influencing factors for positive duration of viral nucleic acid in patients with severe coronavirus disease 2019 (in Chinese) Chloroquine phosphate accelerates the conversion of nucleic acid to negative in 88 common COVID?19 patients (in Chinese) Estimation of the transmission risk of the 2019-nCoV and its implication for public health interventions Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: a data-driven analysis in the early phase of the outbreak Epidemiology and predictors of survival of MERS-CoV infections in Riyadh region Long-term SARS coronavirus excretion from patient cohort Temporal dynamics in viral shedding and transmissibility of COVID-19 Temporal profiles of viral load in posterior oropharyngeal saliva samples and serum antibody responses during infection by SARS-CoV-2: an observational cohort study Virological assessment of hospitalized patients with COVID-2019 Decreased B cells on admission associated with prolonged Viral RNA shedding from the respiratory tract in coronavirus disease 2019: a case-control study Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72,314 cases from the chinese center for disease control and prevention SARS-CoV-2 infection in children Clinical evidence does not support corticosteroid treatment for 2019-nCoV lung injury Effects of early corticosteroid treatment on plasma SARS-associated Coronavirus RNA concentrations in adult patients National Health Commission of the People's Republic of China. Diagnosis and Treatment for COVID-19 (Trial Version 8) Receptor recognition by the novel coronavirus from Wuhan: an analysis based on decade-long structural studies of SARS Coronavirus Current knowledge of COVID-19 and infection prevention and control strategies in healthcare settings: a global analysis Gastrointestinal endoscopy in the Era of COVID-19 The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh. 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.Copyright © 2021 Yan, Zhang, Chen, Jiang, Liu, Zhou, Huang, Zhou, Guan, Ding, Chen, Lan, Fu, Wu, Li and Yang. 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.