key: cord-0978331-6axhqsg8 authors: Lin, Sheng-Nan; Rui, Jia; Chen, Qiu-Ping; Zhao, Bin; Yu, Shan-Shan; Li, Zhuo-Yang; Zhao, Ze-Yu; Wang, Yao; Zhu, Yuan-Zhao; Xu, Jing-Wen; Yang, Meng; Liu, Xing-Chun; Yang, Tian-Long; Luo, Li; Deng, Bin; Huang, Jie-Feng; Liu, Chan; Li, Pei-Hua; Liu, Wei-Kang; Xie, Fang; Chen, Yong; Su, Yan-Hua; Zhao, Ben-Hua; Chiang, Yi-Chen; Chen, Tian-Mu title: Effectiveness of potential antiviral treatments in COVID-19 transmission control: a modelling study date: 2021-04-19 journal: Infect Dis Poverty DOI: 10.1186/s40249-021-00835-2 sha: 29ddb661c7afd7feff966331fcc4fdeb9541a6ca doc_id: 978331 cord_uid: 6axhqsg8 BACKGROUND: Novel coronavirus disease 2019 (COVID-19) causes an immense disease burden. Although public health countermeasures effectively controlled the epidemic in China, non-pharmaceutical interventions can neither be maintained indefinitely nor conveniently implemented globally. Vaccination is mainly used to prevent COVID-19, and most current antiviral treatment evaluations focus on clinical efficacy. Therefore, we conducted population-based simulations to assess antiviral treatment effectiveness among different age groups based on its clinical efficacy. METHODS: We collected COVID-19 data of Wuhan City from published literature and established a database (from 2 December 2019 to 16 March 2020). We developed an age-specific model to evaluate the effectiveness of antiviral treatment in patients with COVID-19. Efficacy was divided into three types: (1) viral activity reduction, reflected as transmission rate decrease [reduction was set as v (0–0.8) to simulate hypothetical antiviral treatments]; (2) reduction in the duration time from symptom onset to patient recovery/removal, reflected as a 1/γ decrease (reduction was set as 1–3 days to simulate hypothetical or real-life antiviral treatments, and the time of asymptomatic was reduced by the same proportion); (3) fatality rate reduction in severely ill patients (f(c)) [reduction (z) was set as 0.3 to simulate real-life antiviral treatments]. The population was divided into four age groups (groups 1, 2, 3 and 4), which included those aged ≤ 14; 15–44; 45–64; and ≥ 65 years, respectively. Evaluation indices were based on outbreak duration, cumulative number of cases, total attack rate (TAR), peak date, number of peak cases, and case fatality rate (f). RESULTS: Comparing the simulation results of combination and single medication therapy s, all four age groups showed better results with combination medication. When 1/γ = 2 and v = 0.4, age group 2 had the highest TAR reduction rate (98.48%, 56.01–0.85%). When 1/γ = 2, z = 0.3, and v = 0.1, age group 1 had the highest reduction rate of f (83.08%, 0.71–0.12%). CONCLUSIONS: Antiviral treatments are more effective in COVID-19 transmission control than in mortality reduction. Overall, antiviral treatments were more effective in younger age groups, while older age groups showed higher COVID-19 prevalence and mortality. Therefore, physicians should pay more attention to prevention of viral spread and patients deaths when providing antiviral treatments to patients of older age groups. [Image: see text] SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s40249-021-00835-2. On January 30, 2020, the World Health Organization (WHO) declared novel coronavirus disease 2019 (COVID-19) a public health emergency of international concern [1] and a challenging public health crisis [2] . Many studies have reported that COVID-19 has caused an enormous disease burden globally [3, 4] , and it continues to spread vigorously in the United States, Brazil, and other countries [5] . Presently, international epidemic prevention and control strategies for COVID-19 include pharmaceutical interventions (PIs) (antiviral treatment and vaccination) and non-pharmaceutical interventions (NPIs) [isolation, wearing of masks, social distancing (closing of schools, cancellation of collective activities, and crowd gathering restrictions), and tourism restrictions]. The effectiveness of NPIs have been proven in many countries [6] [7] [8] [9] [10] [11] [12] [13] . China, especially, successfully controlled the first wave of the COVID-19 outbreak by strictly implementing various public health policies, including the timely isolation of cases, contact tracing, social distancing control, and entertainment activity restriction [12, [14] [15] [16] . Furthermore, Republic of Korea's timely public health measures also achieved remarkable results in epidemic control [17] . NPIs implemented were dependent on a variety of factors, such as social and economic systems and required good cooperation from the public. Further, 'lockdowns' produced substantial economic hardship [18] . Meanwhile, model analysis showed that characteristics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) allow it to stably coexist with people [19] . When SARS-CoV-2 reaches the endemic phase, its overall pathogenicity in the population may be equivalent to that of common cold viruses [19] . NPIs can neither be maintained indefinitely nor conveniently implemented in all countries. Furthermore, evidence has already indicated that these types of public health measures alone may be insufficient for controlling COVID-19 [20] . As work gradually resumes, and productivity increases, various provinces and cities in China are facing the potential risk of an additional wave of COVID-19 outbreak [21] . Since the virus continues to spread, PIs will be essential for preventing and controlling the spread of COVID-19 [14, 22] , thereby reducing social control dependence [23] . Recently, numerous studies have been published on the developmental process, safety, and efficacy of COVID-19 vaccines [24] [25] [26] [27] [28] [29] [30] . Moreover, our team has conducted relevant research on this subject [31] . However, the primary role of vaccination is to prevent COVID- 19 . Approximately 2000 ongoing trials on the efficacy of pharmacological therapy against SARS-CoV-2 infection have been registered in the WHO International Clinical Trials Registry Platform. Nevertheless, no specific drug has been confirmed to be effective [32] . A systematic review and network meta-analysis of the efficacy and safety of 31 existing drugs revealed that anti-inflammatory agents (tocilizumab, anakinra, and gamma globulin) and remdesivir might improve the prognosis of patients with COVID-19. Furthermore, hydroxychloroquine may be associated with cardiac and non-cardiac safety risks after clearing the virus [33] , and dexamethasone reduces the 28-day mortality risk in severely ill patients, especially in those receiving invasive ventilation [34] . The efficacy of antiviral treatments requires further verification, and most recent studies evaluating antiviral treatments focus on clinical efficacy specifically. Further, many studies have shown that the prevalence and case fatality rate of COVID-19 varies among different age groups [15, [35] [36] [37] [38] . Therefore, we aimed to evaluate the effectiveness of COVID-19 antiviral treatments among different age groups based on the hypothetical or real-life efficacy of antiviral treatments from the public health perspective. On the one hand, we hope to predict the public health effects of existing antiviral treatments; on the other hand, we hope to provide a model that can be directly applied to evaluate the effectiveness of novel antiviral treatments based on their clinical efficacy, when such antiviral treatments becomes available in the future. The following previously published data were collected from patients with COVID-19 in Wuhan from 2 December 2019 to 16 March 2020: age, clinical severity classification (mild, moderate, severe, and critical), date of onset, and date of death [31] . The study population was divided into four age groups based on previous study findings [31, 36, 39] [31] . An age-specific Susceptible-exposed-symptomaticasymptomatic-recovered/removed (SEIAR) model was developed previously [31, 39] . To implement the model, individuals were divided into the following five categories: susceptible (S), exposed (E), symptomatic (I), asymptomatic (A), and recovered/removed (R). The rate of infection transmissibility for each age group was estimated using the model, and the process of data fitting detailed previously [31] depended on the following assumptions: 1) S individuals in age group i (the range of i and j was 1-4, indicating different age groups) were infected by exposure to I and A individuals in age group i and other age groups. The transmission rates of S were β and κβ (0 < κ < 1) after effective contacts with I and A, respectively. 2) The transmission rate (β) was divided into two categories: β ii (within the age group i) and β ij (age groups i to j). 3) Parameter p was set as the proportion of A individuals, whereas incubation and latent periods of I and A were 1/ω and 1/ω′, respectively. 4) The times from categories I and A to category R were set as 1/γ and 1/γ′, respectively. 5) The case fatality rate was set as f for members of category I who died after infection. The flowchart of this model is presented in our previously published paper, as well as some parameter estimations, such as k, p, ω, ω′, γ, and γ′ [31] . The equations used in the current model are as follows: With subscripts i and j (i ≠ j) representing age groups 1-4. Based on the existing data stage grouping, we performed data fitting and calculated the dissemination capacity of the four age groups (stages 1, 2, 3, and 4 occurred from 2 December 2019 to 23 January 2020; 24 January to 2 February 2020; 3 February to 18 February 2020; and 19 February 2020 to 16 March 2020, respectively; the details of which are shown in our previous article [31] ). The results of data fitting are shown in Additional file 1: Fig. S1 and Additional file 3: Table S1 . Firstly, death only occurred in severely ill patients [15] , therefore, we classified all patients as either severely (severe and critical) or non-severely ill (mild and moderate), according to the COVID-19 clinical severity classification of each patient. Based on the existing age-specific model [39, 40] , we distinguished severely ill patients from category I patients. The framework of the SEIAR model for antiviral treatment is shown in Fig. 1 . In this model, the population was divided into six categories, as follows: susceptible (S), exposed (E), symptomatic (I), severely ill patients (C), asymptomatic (A), and recovered/removed (R). When implementing the model, the following assumptions were made: 1) S individuals in age group i (the range of i and j was 1-4, indicating different age groups) were infected by exposure to I and A individuals in age group i and other age groups. The transmission rates of S were β and κβ (0 < κ < 1) after effective contacts with I and A, respectively. 2) The transmission rate (β) was divided into two categories, as follows: β ii (within the age group i) and β ij (age group i to j), and the reduction ratio of β was set as v (the initial value of v was 0). 3) Parameter p was set as the proportion of A individuals, whereas the incubation and latent periods of I and A were 1/ω and 1/ω′, respectively. 4) The ratio of severely ill patients in age group i was set as q i , and the fatality rate of severely ill patients in age group i was set as f ci . The reduction ratio of f ci was set as z (initial value of z was 0), and the number of people who changed from to category R were set as 1/γ and 1/γ′, respectively. Therefore, the numbers of people who transitioned from I to R and A to R, respectively, were γI and γ′A at time t. The following equations were used in the model: Parameter values used in the model and the methods used for their determination are detailed in Table 1 . The model had 30 parameters (β ii , β ij , β jj , β ji , f ci , q i , κ, p, ω, ω′, γ and γ′). Fig. 1 The framework of Susceptible-Exposed-Symptomatic-Asymptomatic-Recovered/Removed (SEIAR) model 1) We used an age-specific model to fit data, and obtained results of four stages. The parameters of β that we used for the simulation were the results determined during stage 1 (Additional file 2: Fig. S2 ). 2) Based on the analysis of data, ratios of severely ill patients (q i ) were 3.25%, 10.99%, 19.14%, and 37.79% for age groups 1, 2, 3, and 4, respectively. Fatality rates for severely ill patients (f c ) of the groups were 5.00%, 5.16%, 18.20%, and 39.79%, respectively. 3) A previous study indicated that 4.11% and 6.30% of individuals would become infected after close contact with A and I patients, respectively [41] . Moreover, in the study, k was set as 0.65 [31] ; therefore, in this study, we assumed that the transmissibility of A infections was 0.65 times that of I infections. Further, it has been reported that the transmissibility of I is 3.9 times that of A, and an asymptomatic individual may infect 11 people [42] . Another study set κ as 1.0 [39] . Based on previous results, we set the range of k to 0-1. For instance, the proportion of A patients in the Diamond Princess cruise ship was 17.9% [95% confidence interval (CI): 15.5-20.2%] [43] and 20.75% in Ningbo City [41] . In addition, a study estimated an A ratio of 30.8% (95% CI: 7.7-53.8%) using a binomial distribution [44] , whereas another study reported that A patients constituted 5% to 28% of all patients with COVID-19 [35] . Furthermore, it was also reported that the A ratio could reach 78% [45] . Therefore, since we set the A proportion (p) to 0.36 previously (1.6-78%) [31] , we used the same value in the present study ( Fig. 2-a) . A previous study revealed that the incubation period early in the epidemics in Wuhan City was 4 days (interquartile range: 2-7 days) [46] , whereas other studies indicated that the incubation period in Wuhan [47] and Ningbo [41] raged from 2-14 days and 2-18 days, respectively. In addition, another study reported an incubation period of 5.1 days (95% CI: 4.5-5.8) [48] . Therefore, we set the incubation period to 5 days previously (ω = ω' = 0.2) based on the median incubation periods of multiple studies [31] . The same value was used in the present study. We show incubation periods reported in many different studies in Fig. 2c , and, based on current findings, determined that setting the incubation period as 5 days is representative of published data. 6) Various studies have reported the following durations from symptom onset to hospitalisation: 6.39 (range: 1.00-8.83), 7, 4-6, and 4.1-7.5 days [36, [49] [50] [51] . Moreover, right-truncated data indicated that the time from illness onset to hospitalisation ranged from 2.7 to 8.0 days [47] . Therefore, in our previous study, the infectious periods for I patients were set to 5 days (γ = 0.2) [31] . The same value, which was based on the median infectious period determined in multiple studies, was used in the present study. We show infectious periods of I determined previously in Fig. 2 -b. A prior study indicated that the median infectious period of 24 A patients was 9.5 days (range: 1-21 days) [52] . Therefore, we set γ' as 0.1 in the model [31] . Existing antiviral treatments and their clinical efficacies are shown in Table 2 . Based on current research, we divided the effectiveness of antiviral treatment into three types: (1) reducing transmission (β); (2) decreasing the infectiousness of I and A (1/γ and 1/γ'); and (3) reducing the fatality rate of severely ill patients (f c ). 1) A prior study showed that remdesivir efficiently inhibited viral infection in a human cell line (human liver cancer Huh-7 cells) sensitive to SARS-CoV-2; thus, remdesivir can reduce the infectivity of COVID-19 to an extent [46] . Similarly, researchers determined that chloroquine has antiviral activity and can synergistically enhance antiviral effects of remdesivir in vivo [46] . In the model, antiviral treatments capable of reducing viral activity were revealed by decreases in β. However, existing studies have not specified the extent to which antiviral treatments reduce β values, and therefore, the reduction of β was determined by adjusting the value of v (0-0.8) to simulate that of hypothetical antiviral treatments. 2) Previous studies have shown that remdesivir might not improve recovery or reduce mortality in patients with COVID-19 [53, 54] . However, a preliminary study in the United States, involving many hospitals globally, demonstrated that using remdesivir in hospitalised patients with COVID-19 receiving oxygen therapy shortened their recovery time from 15 to 11 days and reduced risk of respiratory infections [55] . Some clinical trials have shown that the recovery time of patients with COVID-19 who use remdesivir is 31% faster than that of patients provided a placebo [56] . For severely ill patients, combined use of ribavirin with interferon beta-1b and lopinavir/ ritonavir safely and efficiently shortens the duration and reduce symptoms of viral infection in patients with mild-to-moderate disease, with mild side effects [57] . In addition, other studies have shown that chloroquine phosphate can shorten the clinical course of COVID-19 [36] . In this model, antiviral treatments capable of reducing the duration needed to transition from categories I and A to category R were reflected in 1/γ and 1/γ′ decrease. The reduced duration needed to transition from category I was set to 1-3 days to simulate hypothetical or real-life antiviral treatments, and the duration from category A was reduced proportionally to 1/γ. 3) Dexamethasone use reduced the number of deaths among patients with COVID-19 using ventilators by 30% [34, 58] . In addition, the improvement effect of tocilizumab on severely or critical ill patients with COVID-19 has also been confirmed. Futher, interferon-alpha 2b was reported to affect survival in patients with COVID-19 [59] [60] [61] . Another research team successfully identified multiple highly active neutralising antibodies in the plasma of recovered patients [62] . Experimental treatments using plasma of recovered patients have also been shown to significantly effective in severely and critically ill patients [63] , but their widespread use is not feasible at present. In the model, antiviral treatments capable of reducing fatality rates of severely ill patients were reflected as decreases in f c (reduction was set as z). We set z as 0.3, to simulate real-life antiviral treatment. The evaluation indices were as follows: The values of these indices were the simulated number of cases provided different antiviral treatments rather than the actual number. Oster) was employed for curve fitting and simulation. We used the same simulation methods (Runge-Kutta fourth-order method with the tolerance set to 0.001) as described previously [64] [65] [66] [67] [68] [69] . Berkeley Madonna was used to adopt the curve fitting of the least root-meansquare deviation. Goodness of fit was judged by the coefficient of determination (R 2 ) value, which was calculated using IBM SPSS Statistics for Windows, version 21.0 (IBM Corp. Armonk, NY, USA) and Microsoft Office Excel 2016 (Microsoft, Redmond, WA, USA). In this study, two parameters were used to analyse the sensitivity of the model: k (0-1) and p (0.016-0.78), each was split into 1000 values according to its range. The mean and standard deviation (SD) were calculated for sensitivity analysis. This study evaluated the intervention effect of COVID-19 therapies by simulating several antiviral treatments with different efficacies. Overall, we observed age-dependent differences. Age group-specific results obtained without intervention indicated that an increase in age was associated with increased TAR and f, as well as advanced PD. However, with an increase in age, CNC, NPC, and OD initially increased and then decreased. Changes in β, γ andγ' did not vary. When β, γ, andγ' of each age group remained unchanged, changes in f c did not affect the trend and severity of the epidemic, but reduced f. However, by reducing β without changing f c , γ, and γ', we observed that CNC, NPC, and TAR continued to decrease, while OD was prolonged, and PD continued to be delayed. Although f continued to fall, β needed to be reduced by 70% in age group 1 and 80% in age groups 2, 3, and 4 to completely control the epidemic. When β and f c remained unchanged, γ and γ' increased, and CNC, NPC, and TAR continued to decrease; nevertheless, OD was prolonged, and PD continued to be delayed (Additional file 4: Table S2 , Additional file 5: Table S3, Additional file 6: Table S4 , Additional file 7: Table S5 ). When different antiviral treatments were used, the TAR reduction rate was higher than f for all age groups, and the reduction of different age groups varied. When comparing the reduction of TAR with f when only one type of antiviral treatment was used, we found that changes in β and γ had the greatest effects on TAR and f, respectively. When the value of v was 0.7, age group 2 had the highest TAR reduction rate (95.78%). However, when the value of 1/γ was 2, age group 1 had the highest f reduction rate (64.56%). Comparing the simulation results of combination and single medications, all four age groups showed better results with the combination medication. When 1/γ = 2, and v = 0.6 or 1/γ = 2, v = 0.6, and z = 0.3, age group 1 had the highest TAR reduction rate. However, When 1/γ = 2 and v = 0.4 or 1/γ = 2, v = 0.4, and z = 0.3, age group 2-4 had the highest TAR reduction rates. The TAR reduction rate of age group 2 was the highest among the 4 age groups considered (98.48%, from 56.01 to 0.85%) (Fig. 3b) , whereas the absolute TAR reduction value of age group 3 was highest among the 4 age groups (0.6070, from 0.6292 to 0.0222) (Fig. 3c) . When 1/γ = 2, z = 0.3, and v = 0.1, age group 1 had the highest f reduction rate. However, when v = 0.4, z = 0.3, and 1/γ = 2, age group 2-4 had the highest f reduction rate. The f reduction rate of age group 1 was the highest among all 4 age groups (83.08%, from 0.71 to 0.12%) (Fig. 4a) . whereas the absolute f reduction value of age group 4 was the highest (0.1451, from 0.2938 to 0.1487) (Fig. 4d) . Details are shown in Additional file 8: Table S6 , Additional file 9: Table S7 , Additional file 10: Table S8 , Additional file 11: Table S9 . In this study, we found that parameters k and p used in the model were included in the range of simulated mean ± SD values (Fig. 5) . These two parameters were not sensitive to the model. This study was based on our existing age-specific transmissibility model of SARS-CoV-2, which is the first model to quantify the transmissibility of COVID-19 within and between different age groups [39] . We pursued an innovative approach by incorporating multiple factors, such as the ratio of severely ill patients (q) and the fatality rate of severely ill patients (f c ) of each age group, to build an antiviral treatments evaluation model for a set of age groups. The model was used to assess intervention efficacy of different antiviral treatments by changing parameters as needed. However, most current COVID-19 pharmacological therapies remain at the clinical efficacy evaluation stage, and the reliability of antiviral treatments require further exploration. The focus of this study was to evaluate the effectiveness of antiviral treatments in different age groups from a public health perspective, using hypothetical or real-life antiviral treatments. Parameters of this study were as follows: β and disease natural history parameters. β parameters were obtained from the results of curve fitting. The date separating stages 1 and 2 was 23 January 2020, as Wuhan declared a lockdown on 22 January 2020. Thus, we can approximately regard the first stage as the non-intervention stage. It is universally accepted that β parameters of stage 1 represent the transmissibility of the disease in its natural state [31] . This is also the possible reason why the number of simulated cases exceeded the actual number of cases in Wuhan. Under this premise, we excluded the efficacy of other possible interventions when evaluating the efficacy of antiviral treatments. Moreover, results of R 2 calculations showed that the model fitted well. There are two sources of disease natural history parameters: the ratio of severely ill patients (q) and the fatality rate of severely ill patients (f c ), which were calculated using real data. Remaining parameters (κ, p, ω, ω′, γ, and γ′) were quoted from prior studies. Sensitivity analyses were performed for parameters of k and p, and other values were based on published findings. Median values of published findings were calculated. Some studies have reported variations in COVID-19 prevalence between different age groups [35, 36, 39] , akin to our study findings. Without intervention, groups 1 and 4 had the lowest and highest TAR values, respectively; these were similar to findings of a study conducted in Hungary [37] . Another three-age-group study using a generalised linear mixed model revealed that people aged ≥ 65 years had a higher risk of becoming infected with SARS-CoV-2 than those aged 15-64 years [70] . It was speculated that reasons for these age-based differences may be due to increased rates of underlying diseases observed in adults aged 65 years [38] that predisposed them to COVID-19 infection due to their low level of immunity. Differing fatality rates observed in severely ill patients did not have an impact on the epidemic. This may be due to the fact that most severely ill patients are in the intensive care units and reducing the fatality rate will not affect other populations. When antiviral treatments reduce β or γ and γ', each age group had a different degree of control. This may have been due to differences in disease transmissibility observed for each age group. For example, the age group 2 population was the largest group, which may suggest that the decrease in β had the greatest impact on the TAR in the age group. Moreover, we hypothesised that different age groups had different sensitivities to different antiviral treatments. However, the specific antiviral treatment sensitivity of each age group requires further pharmacological exploration. Furthermore, differences in autoimmunity may have an effect on the effectiveness of potential antiviral treatments. Some parameters were not included in the study, such as data associated with non-severe cases and that of immunity acquisition over time. Future scenario predictions may be inaccurate and unpredictable, and World Health Organization declares global emergency Current understanding of COVID-19 clinical course and investigational treatments Impact of the burden of COVID-19 in Italy: results of disabilityadjusted life years (DALYs) and productivity loss The burden of disease due to COVID-19 in korea using disability-adjusted life years WHO. Coronavirus disease (COVID-19 Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe Mathematical modeling of COVID-19: impact of non-pharmaceutical interventions in India The COVID-19 epidemic, its mortality, and the role of non-pharmaceutical interventions Effect of COVID-19 non-pharmaceutical interventions and the implications for human rights COVID-19 epidemic monitoring after non-pharmaceutical interventions: The use of time-varying reproduction number in a country with a large migrant population Centre for the Mathematical Modelling of Infectious Diseases C-wg. Effects of nonpharmaceutical interventions on COVID-19 cases, deaths, and demand for hospital services in the UK: a modelling study The end of "cordon sanitaire" in Wuhan: the role of non-pharmaceutical interventions against COVID-19 Impact of non-pharmaceutical interventions on the COVID-19 epidemic: a modelling study Active case finding with case management: the key to tackling the COVID-19 pandemic 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 Effect of non-pharmaceutical interventions to contain COVID-19 in China Transmission dynamics of the COVID-19 outbreak and effectiveness of government interventions: a data-driven analysis The South African response to the pandemic Immunological characteristics govern the transition of COVID-19 to endemicity COVID-19: lessons from SARS and MERS Effectiveness of interventions to control transmission of reemergent cases of COVID-19 COVID-19 vaccine: China's public good versus ' America First Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period Age and frailty in COVID-19 vaccine development Broad cross-national public support for accelerated COVID-19 vaccine trial designs The COVID-19 vaccine development landscape Evolution of the COVID-19 vaccine development landscape Postapproval vaccine safety surveillance for COVID-19 vaccines in the US Safety and efficacy of the Russian COVID-19 vaccine: more information needed Should we mandate a COVID-19 vaccine for children? The optimal vaccination strategy to control COVID-19: a modeling study based on the transmission scenario in Wuhan City Pharmacological interventions for COVID-19: protocol for a rapid living systematic review with network meta-analysis. medRxiv Comparative efficacy and safety of pharmacological interventions for the treatment of COVID-19: a systematic review and network meta-analysis of confounder-adjusted 20212 hospitalized patients Dexamethasone in hospitalized patients with COVID-19 Clinical and epidemiological features of 36 children with coronavirus disease 2019 (COVID-19) in Zhejiang, China: an observational cohort study Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia Early phase of the COVID-19 outbreak in Hungary and post-lockdown scenarios The first 100 cases of COVID-19 in a hospital in Madrid with a 2-month follow-up A five-compartment model of age-specific transmissibility of SARS-CoV-2 A mathematical model for simulating the phase-based transmissibility of a novel coronavirus Epidemiological characteristics of infection in COVID-19 close contacts in Ningbo city The relative transmissibility of asymptomatic COVID-19 infections among close contacts Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship Estimation of the asymptomatic ratio of novel coronavirus infections (COVID-19) Covid-19: four fifths of cases are asymptomatic, China figures indicate Clinical characteristics of coronavirus disease 2019 in China Incubation period and other epidemiological characteristics of 2019 novel coronavirus infections with right truncation: a statistical analysis of publicly available case data The incubation period of coronavirus disease 2019 (COVID-19) bion Epidemiological characteristics of the first 53 laboratory-confirmed cases of COVID-19 epidemic in Hong Kong Clinical features of patients infected with 2019 novel coronavirus in Wuhan Novel coronavirus outbreak in Wuhan, China, 2020: intense surveillance is vital for preventing sustained transmission in new locations Clinical characteristics of 24 asymptomatic infections with COVID-19 screened among close contacts in Nanjing Remdesivir in adults with severe COVID-19: a randomised, double-blind, placebo-controlled, multicentre trial Evaluation of the efficacy and safety of intravenous remdesivir in adult patients with severe COVID-19: study protocol for a phase 3 randomized, double-blind, placebocontrolled, multicentre trial Remdesivir for the treatment of COVID-19-preliminary report NIH clinical trial shows Remdesivir accelerates recovery from advanced COVID-19 Triple combination of interferon beta-1b, lopinavir-ritonavir, and ribavirin in the treatment of patients admitted to hospital with COVID-19: an open-label, randomised, phase 2 trial Association between administration of systemic corticosteroids and mortality among critically ill patients with COVID-19: a meta-analysis Effective treatment of severe COVID-19 patients with tocilizumab Impact of tocilizumab on mortality and/or invasive mechanical ventilation requirement in a cohort of 206 COVID-19 patients Tocilizumab therapy reduced intensive care unit admissions and/or mortality in COVID-19 patients Potent neutralizing antibodies against SARS-CoV-2 identified by high-throughput single-cell sequencing of convalescent patients' B cells Effect of convalescent plasma therapy on time to clinical improvement in patients with severe and life-threatening COVID-19: a randomized clinical trial Transmissibility of acute haemorrhagic conjunctivitis in small-scale outbreaks in Hunan Province Incidence dynamics and investigation of key interventions in a dengue outbreak in Ningbo City, China Estimating the transmissibility of hand, foot, and mouth disease by a dynamic model Seasonality of the transmissibility of hand, foot and mouth disease: a modelling study in Xiamen City Relative transmissibility of hand, foot and mouth disease from male to female individuals Evaluating the effects of common control measures for influenza A (H1N1) outbreak at school in China: a modeling study Changes in contact patterns shape the dynamics of the COVID-19 outbreak in China We wish to thank the Undergraduate Innovation Practice Platform of School of Public Health, Xiamen University for support. Fig. 3 Evaluation of the effectiveness of potential antiviral treatments (with highest TAR reduction in each age group). a ≤ 14 years; b15-44 years; c 45-64 years; d ≥ 65 years current calculation methods may overestimate or underestimate the effectiveness of potential antiviral treatments throughout the population considered. Antiviral treatments are more effective for COVID-19 transmission control than for mortality reduction. Overall, antiviral treatments were most effective when used to treat younger age groups, while older age groups showed higher disease prevalence and mortality. Therefore, older age groups require more attention with respect to the use of antiviral treatments in clinical practice. A: Asymptomatic; CNC: Cumulative number of cases; COVID-19: Novel coronavirus disease 2019; E: Exposed; f: Case fatality rate; f c : The fatality rate of severely ill patients; I: Symptomatic; NPC: Number of peak cases; OD: Outbreak duration; PD: Peak date; R: Recovered/removed; S: Susceptible; SARS-CoV-2: Severe acute respiratory syndrome coronavirus 2; SEIAR: Susceptible-exposed-symptomatic-asymptomatic-recovered/removed; TAR : Total attack rate; WHO: World Health Organisation. The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s40249-021-00835-2. Table S1 . The results of goodness of fit in four stages of four age groups. Table S2 . The effectiveness of potential antiviral treatments in group 1 (ages ≤ 14 years). Table S3 . The effectiveness of potential antiviral treatments in group 2 (ages 15-44 years). The datasets used and analysed during the current study are available from Miss Shengnan Lin (shengnanlin0228@163.com) on reasonable request. The study was approved by the Ethics Committee of School of Medicine, Xiamen University. Only broad information (such as the date of illness onset and age) of the cases were collected with no identifying patient information. Neither medical intervention nor biological samples were involved. Study procedures and results did not affect the clinical management of patients. The authors declare that they have no competing interests.