Skip to main content

ORIGINAL RESEARCH article

Front. Med., 18 August 2021
Sec. Infectious Diseases – Surveillance, Prevention and Treatment
This article is part of the Research Topic COVID-19: Integrating Artificial Intelligence, Data Science, Mathematics, Medicine and Public Health, Epidemiology, Neuroscience, and Biomedical Science in Pandemic Management View all 95 articles

Containing the Transmission of COVID-19: A Modeling Study in 160 Countries

\nYan Niu&#x;Yan Niu1Jia Rui&#x;Jia Rui2Qiupeng Wang,&#x;Qiupeng Wang2,3Wei ZhangWei Zhang2Zhiwei ChenZhiwei Chen4Fang XieFang Xie2Zeyu ZhaoZeyu Zhao2Shengnan LinShengnan Lin2Yuanzhao ZhuYuanzhao Zhu2Yao WangYao Wang2Jingwen XuJingwen Xu2Xingchun LiuXingchun Liu2Meng YangMeng Yang2Wei ZhengWei Zheng2Kaixin ChenKaixin Chen2Yilan XiaYilan Xia2Lijuan XuLijuan Xu2Shi ZhangShi Zhang2Rongrong JiRongrong Ji4Taisong JinTaisong Jin4Yong ChenYong Chen5Benhua ZhaoBenhua Zhao2Yanhua SuYanhua Su2Tie SongTie Song6Tianmu Chen
Tianmu Chen2*Guoqing Hu
Guoqing Hu7*
  • 1Chinese Center for Disease Control and Prevention, Beijing, China
  • 2State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Xiamen, China
  • 3School of Journalism and Communication, Peking University, Beijing, China
  • 4Media Analytics and Computing Lab, Department of Artificial Intelligence, School of Informatics, Xiamen University, Xiamen, China
  • 5Department of Stomatology, School of Medicine, Xiamen University, Xiamen, China
  • 6Guangdong Provincial Center for Disease Control and Prevention, Guangzhou, China
  • 7Department of Epidemiology and Health Statistics, Xiangya School of Public Health, Central South University, Changsha, China

Background: It is much valuable to evaluate the comparative effectiveness of the coronavirus disease 2019 (COVID-19) prevention and control in the non-pharmacological intervention phase of the pandemic across countries and identify useful experiences that could be generalized worldwide.

Methods: In this study, we developed a susceptible–exposure–infectious–asymptomatic–removed (SEIAR) model to fit the daily reported COVID-19 cases in 160 countries. The time-varying reproduction number (Rt) that was estimated through fitting the mathematical model was adopted to quantify the transmissibility. We defined a synthetic index (IAC) based on the value of Rt to reflect the national capability to control COVID-19.

Results: The goodness-of-fit tests showed that the SEIAR model fitted the data of the 160 countries well. At the beginning of the epidemic, the values of Rt of countries in the European region were generally higher than those in other regions. Among the 160 countries included in the study, all European countries had the ability to control the COVID-19 epidemic. The Western Pacific Region did best in continuous control of the epidemic, with a total of 73.76% of countries that can continuously control the COVID-19 epidemic, while only 43.63% of the countries in the European Region continuously controlled the epidemic, followed by the Region of Americas with 52.53% of countries, the Southeast Asian Region with 48% of countries, the African Region with 46.81% of countries, and the Eastern Mediterranean Region with 40.48% of countries.

Conclusion: Large variations in controlling the COVID-19 epidemic existed across countries. The world could benefit from the experience of some countries that demonstrated the highest containment capabilities.

Introduction

Coronavirus disease 2019 (COVID-19) is a novel infectious disease caused by severe acute respiratory syndrome coronavirus 2 virus (SARS-CoV-2) (1). On March 11, 2020, the World Health Organization (WHO) declared a global pandemic of COVID-19 (2). The COVID-19 pandemic has brought tremendous economic and medical burden, and has severely affected people's normal work and daily lives. In most countries, the health system faces collapse, and there are medical shortages of varying degrees (3). In countries affected by COVID-19, many factories and shops have been closed, and thus many people have lost their jobs and income. The International Labor Organization estimated that 25 million jobs were lost, which could lead to a loss of at least US$220 billion in income in developing countries (4). In 2020, the World Travel and Tourism Council warned that up to 1.67% of the world's gross domestic product could be lost due to the reduction in flight travel (5). Therefore, it is very important for the world to evaluate the comparative effectiveness of COVID-19 prevention and control across countries.

Current research on COVID-19 mainly focuses on issues such as transmission mechanisms, pathogenic mechanisms, prevention and control measures, and infection control (6). Mathematical models of pandemics are useful for understanding the early transmission dynamics of infections, achieving pandemic prediction, and evaluating the effectiveness of control measures. They have been widely used in the prediction of the basic reproduction number (R0) of COVID-19, and reports of R0 for COVID-19 have ranged from 2.6 to 4.71 (612). One recent study suggested that the median R0 of COVID-19 may be up to 5.7 (13), which was much higher than the previously estimated maximum value of R0. The reported R0 varied across countries. The extent each country to control the epidemic of COVID-19 in the early stage significantly affect the future development of the global pandemic. Several models, such as the stochastic transmission model (14), and ordinary differential equation model (11, 12, 15, 16), have been used in the research on COVID-19; however, research has yet to be conducted at the global level.

This study developed an innovative mathematical model to assess the transmissibility and control capability of COVID-19 epidemic of 160 countries in the initial phase of global pandemic (before August 22, 2020). Considering that COVID-19 can cause a lot asymptomatic infections (A), we added “A” to the SEIR model and formulated a susceptible–exposure–infectious–asymptomatic–removed (SEIAR) model, which better characterized the transmission of COVID-19 compared to the SEIR model. Reproduction number was used to reflect the transmissibility of COVID-19. We also evaluated the ability of 160 countries to control the COVID-19 epidemic.

Methods

Study Design

The analysis process used in this study is shown in Figure 1. First, we collected COVID-19 data from each country in the six regions defined by the World Health Organization (WHO) from January 22, 2020, to August 22, 2020, and excluded countries (11 countries) with poor data quality. Second, we fit the mathematical model for the remaining 172 countries and excluded countries (12 countries) due to the poor goodness of fit (the p values of goodness of fit was <0.05). Finally, we calculated the transmissibility of COVID-19 and evaluated the prevention and control capacities of the remaining 160 countries.

FIGURE 1
www.frontiersin.org

Figure 1. Study design for analyzing the differences of capability to control COVID-19 among 160 countries.

The level of data monitoring varies across the six WHO regions. The African countries included in the study accounted for 74.47% of the total number of countries in the African region due to insufficient monitoring systems and economic backwardness (17). The remaining five regions were the Southeast Asian Region (100%), the Eastern Mediterranean Region (95.24%), the European Region (92.72%), the Region of Americas (88.57%), and the Western Pacific Region (86.67%). The three countries with the highest cumulative COVID-19 cases were the United States (5,667,112), Brazil (3,582,362), and India (3,044,940). Although some countries in the Western Pacific region, such as China and South Korea, were reported earlier among the countries with the highest cumulative total, the increase of number of cases slowed in the later period. The cumulative number of COVID-19 cases in the European region was between 50,000 and 1,000,000.

Data Sources

To determine the newly confirmed cases and cumulative deaths globally, we used data from the repository published by Johns Hopkins University Center for Systems Science and Engineering and the Esri Living Atlas team and Johns Hopkins University Applied Physics Lab. Then, we applied our codes to clean and process off-the-shelf data for model building. Our code repository will be published after the paper is accepted.

The Model of COVID-19

Considering that COVID-19 can cause asymptomatic infections, we developed a SEIAR model based on our previous work (15, 18). The SEIAR model involves five group of populations: susceptible individuals (S), individuals in incubation period (E), symptomatic individuals (I), asymptomatic individuals (A), and recovered individuals (R). A SEIAR model was developed based on the reported data of COVID-19. The model framework is shown in Figure 2.

FIGURE 2
www.frontiersin.org

Figure 2. Flowchart of the SEIAR model.

The model was based on the following assumptions:

1) Assuming that the infection rate coefficient after effective contact between S and I is β, and assuming that A is infectious, and the transmissibility is κ (0 < κ < 1) times that of I, then at time t, the number of new infections is βS (I + κA).

2) Assuming that the proportion of asymptomatic infection is p, and the incubation period and latent period are 1/ω and 1/ω′, respectively. At time t, the number of people who change from E to A and I are pω′E and (1–p) ωE, respectively.

3) Assuming that the time interval from the onset of the case I to the first diagnosis is 1/γ, then at time t, the number of people who change from I to R is γI. Because there are deaths in the reported data, we set f as the case fatality rate, and then the number of deaths at time t is fI.

4) Assuming that the infectious period of the recessive infected person A is 1/γ′, then at time t, the number of people who change from A to R is γ′A.

Based on these assumptions, the equations of the SEIAR model are presented as follows:

dSdt=-βS(I+κA)dEdt=βS(I+κA)-pωE-(1-p)ωEdIdt=(1-p)ωE-γI-fIdAdt=pωE-γAdRdt=γI+γA

Parameter Estimation

The values of the model parameters and their methods are listed in Table 1. There were 8 parameters in the model: β, κ, p, 1/ω, 1/ω′, f, 1/γ, and 1/γ′.

TABLE 1
www.frontiersin.org

Table 1. Parameter estimation.

1) The collected data of reported COVID-19 cases were employed to fit the SEIAR model to calculate β in each region. In this study, the nodes were segmented according to the inflection point of the epidemic curve shown by the date of onset, and the models were fitted to estimate the β values at different time periods, namely, β1, β2,…,βn.

2) We assumed that the transmissibility of asymptomatic infection was one times that of symptomatic infection (κ = 1) from information in a previous study (19).

3) Because free-access data of asymptomatic infections were unavailable for each country, we determined the parameter p based on the published literatures. Studies showed that the values of p ranged from 16.1 to 50.6% (20, 21). In this study, we set a mean value of 0.2 in the model. We also performed a sensitivity analysis with the p value of 0.1, 0.2, …, and 0.9, respectively, to assess the uncertainty of our results related to with the parameter p.

4) Based on previous studies (22, 23), we set the incubation period of COVID-19 to 5 days, and the latent period to 3 days in this study; therefore, 1/ω = 0.2, 1/ω′ = 0.3333.

5) Parameter f was the fatality rate of COVID-19 that was calculated using actual data.

6) Based on a previous study (12), we set the infectious period of symptomatic infection as 5 days and the infectious period of asymptomatic infection as 10 days; therefore, 1/γ = 0.2 and 1/γ′ = 0.1.

Transmissibility Evaluation Index

In this study, the population was not completely susceptible and artificially adopted some prevention and control measures, so we chose the time-varying reproduction number (Rt) to calculate transmissibility. The calculation formula was as follows:

Rt=βS(1-pγ+κpγ)

The Rt value at the beginning of the epidemic is defined as the basic reproduction number (R0), which was the initial Rt value of each country. When countries have adopted control measures for COVID-19, Rt was defined as the effective reproduction number (Reff).

A synthetic index (IAC) was constructed to evaluate the capability to control COVID-19 in each country. The calculation formula was as follows:

IAC=ITn+IWNi+ICj

In the equation, ITn was defined as the grade of time it took to make COVID-19 under control; n represented the level of the time required to control the COVID-19 epidemic. We divided the time into five grades (≤ 30, 31–60, 61–90, 91–120, 121–150, and 151–180 days) and gave a value of n = 1, 2, 3, 4, 5, and 6, respectively. The specific values are shown in Table 2.

TABLE 2
www.frontiersin.org

Table 2. The specific values of ITn, IWNi, and ICj.

IWNi is the grade of wave number of COVID-19; i represents the number of COVID-19 epidemic fluctuation in each country. When the number of fluctuation is 1, i = 1. When the number is 2, 3, and 4, the value of i is 2, 3, and 4 respectively. The specific values are shown in Table 2.

ICj is the grade of countries where COVID-19 has been brought under control (Reff <1) as of August 22; j represents whether the country under study has the ability to control the COVID-19 epidemic. The condition of whether it can be controlled is whether the country can reduce Reff to below 1. When the country is able to control COVID-19, j = 1, and when it is unable to control it, j = 0. The specific values are shown in Table 2.

Simulation Methods and Statistical Analysis

Berkeley Madonna 8.3.18 software (developed by Robert Macey and George Oster of the University of California, Berkeley at Berkeley. Copyright©1993-2001 Robert I. Macey & George F. Oster) was used for the curve fitting. The fourth-order Runge–Kutta method, with tolerance set at 0.001, was used to perform the curve fitting. The coefficient of determination (R2) was used to assess the goodness of fit. SPSS 13.0 software (IBM Corp., Armonk, NY, USA) was used to calculate the R2.

Results

Epidemiological Characteristics of the COVID-19 Pandemic of 184 Countries

Figure 3 shows the spatial and temporal distribution of the cumulative incidence of COVID-19 in 184 countries since January 31, 2020. By 24:00 on January 31, 2020, the cumulative number of COVID-19 cases in China exceeded 10,000, and there were only a few confirmed cases in other countries. On February 29, the cumulative number of confirmed cases in China exceeded 70,000, and more than 10 cases occurred in some countries in the Western Pacific, Europe, and Canada in the Americas, of which the number of cases in Italy and South Korea exceeded 500 and 3,000, respectively. By March 11, the pandemic had spread globally. The cumulative number of confirmed cases in Italy, France, Spain, Germany, Iran, and other places exceeded 1,000, while the growing rate of China slowed down at this time. By March 31, the COVID-19 pandemic became very serious. The cumulative number of confirmed cases in the United States and Italy exceeded 100,000 cases. By April 21, the cumulative number of confirmed cases in the United States exceeded 800,000, and the cumulative number of newly confirmed cases in France, Spain, Germany, and the United Kingdom exceeded 100,000.

FIGURE 3
www.frontiersin.org

Figure 3. The reported cumulative number of cases of COVID-19 in 160 countries as of August 22, 2020. The map depicted in these figures were taken from Wikimedia Commons (http://commons.wikimedia.org/wiki/Main_Page).

Model Fitting

The goodness-of-fit tests showed that the SEIAR model fitted well for the COVID-19 data of most countries (Figure 4). The R2, which is determined by the number of confirmed cases and the number of fitted new cases, was >0.5 for 112 countries, p < 0.05. The R2 of the Southeast Asian Region, Eastern Mediterranean Region, and European Region were much higher than those of the other regions. However, the SEIAR model did not fit the data of 12 countries. Thus, we excluded the 12 countries from further data analyses.

FIGURE 4
www.frontiersin.org

Figure 4. The time-varying reproduction number (Rt), the coefficient of determination (R2), and p of model fitting results.

Transmissibility of the COVID-19 Pandemic of 160 Countries

As shown in Figure 5, the R0 in the early stage of the pandemic was clearly higher than the Reff value after the control was carried out in the later period, and the Reff value in a few countries increased in the later period. As shown in Figure 4A, in the early stage of the pandemic, the European Region had the strongest transmissibility of COVID-19, where the R0 of many countries exceeds 10.

FIGURE 5
www.frontiersin.org

Figure 5. The spatial distribution of basic reproduction number (R0) and the effective reproduction number (Reff) of COVID-19 in different countries. The map depicted in these figures were taken from Wikimedia Commons (http://commons.wikimedia.org/wiki/Main_Page).

Figure 4 shows a clear relationship between the reproduction number (Rt), R2, and P in 160 countries in the six regions. All regions, except for the African region, began to spread COVID-19 on January 22. Among them, the largest Rt range were in the Western Pacific Region and the Region of Americas, and for each, the maximum value was more than 20. Comparted with other regions, the Rt value in the Region of Americas was larger. Secondly, the Rt ranges of the European Region and the African Region were largest, followed by the Eastern Mediterranean Region. The region with the smallest Rt range was the Southeast Asian Region.

Overall, the Rt value of COVID-19 in most countries showed a trend from strong to weak (e.g., China in the Western Pacific Region), but there are also some countries for which the Rt value increased in the later period (e.g., France in the European Region). However, there are still some countries showing further upward trends, such as Sweden, which implemented the policy of “herd immunity”. Of all the regions, changes in the Rt of various countries in the African region are more complicated, and the changes in the Rt values of the other regions are obviously greater. In the Region of Americas, the majority of countries' Rt at the beginning of the pandemic were concentrated and were above 2.0. In terms of overall changes in Rt trends, prevention and control measures were not as effective in the two regions, compared with the other regions. In these two regions, COVID-19 was not effectively controlled in the early stage, making them epicenters of the COVID-19 pandemic today.

Ability to Control COVID-19 Epidemic of 160 Countries

Figure 6A shows that from January 22, 2020, to August 11, 2020, there was great variation in the number of days in the period in which Rt of COVID-19 was <1 among the 160 countries. Compared with Figure 6B, the longer the period in which Rt was 1, the shorter the time needed to control COVID-19. For example, in China in the Western Pacific region, where Rt was <1 for more than 180 days, it took <20 days to reduce R0 below 1. At the time of the study, many countries in the African Region were not under control, and it took a relatively long time to control COVID-19 in the Regions of Americas.

FIGURE 6
www.frontiersin.org

Figure 6. The spatial distribution of ITn and ICj in 160 countries. ITn is defined as the grade of time it will take to bring COVID-19 under control, and n represents the level of the time required to control the COVID-19 epidemic. ICj is the grade of countries where COVID-19 has been brought under control (Reff <1) as of August, 22, 2020; j represents whether the country under study has the ability to control the COVID-19 epidemic. The map depicted in these figures were taken from Wikimedia Commons (http://commons.wikimedia.org/wiki/Main_Page).

Figure 7A shows that among the 160 countries included in the study, those in the European Region had the strongest national control capability. All 51 European countries included in the study were able to control COVID-19, and the European Region had the highest rate of all the regions (100%), followed by the Region of Americas (85.71%), the Eastern Mediterranean Region (80.95%), the Western Pacific Region (80.87%), the Southeast Asian Region (80%), and the African Region (74.47%). Among all the countries in the European Region, 76.4% were able to control COVID-19 within 60 days. However, in terms of continuous control, the Western Pacific Region did relatively better than the others: 73.76% of its countries on August 22, 2020, at the time of this study. This region was followed by the Region of Americas (52.53%), the Southeast Asian Region (48%), the African Region (46.81%), the European Region (43.63%), and the Eastern Mediterranean Region (40.48%). In countries where the COVID-19 epidemic was continuously controlled, the wave number (IWN) varied among the six regions. In the Western Pacific Region and the Region of Americas, the IWN was up to 2. Control of COVID-19 was relatively stable. In the European Region, Eastern Mediterranean Region, Southeast Asian Region, and African region, the IWN was frequently and relatively unstable.

FIGURE 7
www.frontiersin.org

Figure 7. The transmissibility and control effect classification flow chart to control COVID-19 in 160 countries. (A) shows the classification process of 160 countries; (B) shows the time each region needed to bring COVID-19 under control.

We further evaluated the ability to control COVID-19 (IAC) in 160 countries. At the most powerful level (IAC = 15), six regions had countries with control of COVID-19. Although the European Region and the Western Pacific Region each had three countries, in terms of proportion, 5.88% of the European Region and 23.08% of the Western Pacific Region were able to control COVID-19. The Region of Americas had two countries, and the remaining three continents each had one country. Seven countries, including Morocco, Syrian Arab Republic, Libya, Philippines, India, Indonesia, and Mexico, had the worst epidemic control effect (IAC = 0), and the Rt was <1 for all of these countries throughout the study. From Figure 8, we can see that the countries in the Western Pacific Region were mainly distributed in the front position. Although many countries in the European Region did better, some countries did worse. The Eastern Mediterranean Region and the Southeast Asian Region were relatively evenly distributed in each interval, the countries of the Region of Americas are concentrated at the middle level, and the countries of the African Region are concentrated at the lower level.

FIGURE 8
www.frontiersin.org

Figure 8. The ranking of IAC among160 countries.

Sensitivity Analysis

The six countries randomly selected (South Africa, USA, China, Iran, the UK, and India) located in six WHO regions, respectively. The results of the sensitivity analysis showed that p was significantly affected (see Figure 9). The results of sensitivity analysis of countries are consistent in six WHO regions, the larger the value of p, the smaller the value of Rt. We set the p value of each country to 0.2 for each country, which may result in a larger or smaller Rt value deviating from the actual situation. The range of Rt values corresponding to different p values in six countries was 1.12–3.45, 1.24–3.97, 4.27–8.78, 5.44–13.77, 4.52–12.19, and 0.99–3.29, respectively. We set the same p value in the SEIAR model for all countries to ensure the comparability of results across 160 countries.

FIGURE 9
www.frontiersin.org

Figure 9. Estimated Rt value based on different p in six selected countries in six WHO region.

Discussion

We constructed a theoretical epidemiological model, used Rt to reflect the transmissibility of COVID-19 in each country, and explained the differences in the cumulative incidence of COVID-19 across the six regions. Changes in Rt between January 22 and August 22 indicated a slight or moderate effect of local epidemic prevention and control measures. Based on the indicators above, we further constructed the IAC to evaluate the comprehensive effectiveness of prevention and control measures in 160 countries.

Model Fitting

From Figure 4, we can see that the model fitted the data of most countries, suggesting that the SEIAR model was useful. Therefore, the SEIAR model that we choose was applicable to the data of these countries. The R2 of the Southeast Asian Region, the Eastern Mediterranean Region, and the European Region were much higher than those of other regions. This finding may be due to the early outbreak of COVID-19 in these areas, the large number of infected people, or the high detection rate. The reported data from these areas are more in line with the actual data, thus improving the model's fitness. Some countries with low incidences had lower R2 (e.g., Equatorial Guinea, Guinea, Cameroon, Mongolia, and Liechtenstein). This finding may be related to factors such as the economy of these countries, the level of testing, or the lack of public awareness. There were no obvious epidemiological characteristics in the reported data, which makes poor model fitting results.

Effect of Prevention and Control Measures Among Six Regions

Coronavirus disease spread regionally in a natural transmission state. However, the estimates of Rt at the beginning of the pandemic may depend on numerous biological, socio-behavioral, and environmental factors (24). The duration and severity of each pandemic phase can vary depending on the characteristics of COVID-19 and the public health response. Italy and Iran were the original hot spots outside China, from which new infections were disseminated to nearby countries. There have been higher magnitudes of COVID-19 infection in Europe and America than in Asia, whereas Europe and America show similar strengths of infection (25). The prevention and control measures adopted and implemented by local governments of various countries in the early stages helped to a certain extent (26, 27). Globally, the number of new cases per day tends to be flat. However, there were still some countries without any reduction in transmissibility after interventions, which suggests that prevention and control measures may not have played important roles in the early stages of the pandemic. Such containment is unlikely to be achievable in most countries around the world, although China succeeded in containing the spread of the pandemic for 2 months by carrying out timely and effective measures (28). Some of the literature has shown that without non-pharmaceutical interventions (including intensive case and contact tracking, isolation of moderately ill patients in containment centers, social distancing, and shutting down public life of a whole province and many major cities outside Hubei), the COVID-19 cases in China would likely have shown a 67-fold increase (interquartile range from 44 to 94) by February 29, 2020 (29). However, it is worth noting that the intervention effect of the African Region was not as good as that of other regions because of its weak healthcare system, which urgently needs to be improved.

The lack of effective interventions is a major driver of the pandemic. Moreover, the fact that the virus is an unknown, new mutant type may be a major factor for the current overwhelming increase in cases despite the preemptively implemented precautions and countermeasures taken by at-risk countries (6). In general, although countries in the Western Pacific region may not be able to control COVID-19 as well as those in the European Region, countries in the European Region have relatively weak continuous response and prevention and control, which is also one of the reasons for the repeated severe epidemics.

Of note, there are some limitations in our study, such as setting of p value as 0.2 in 160 countries, but we set the p value in the SEIAR model to the same value to ensure comparability among 160 countries. In practice, some asymptomatic patients may die because of the COVID-19 infection. This issue was not considered in our mathematical model and may probably affect the results slightly. However, it would not have a significant impact on our results considering that the mortality of asymptomatic COVID-19 infection was much lower than that of symptomatic infection.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author Contributions

GH, TC, JR, and YN made substantial contributions to conception and design. QW, ZC, FX, WZha, ZZ, SL, and YW collected the data. JR, YN, YZ, MY, JX, XL, and WZhe conceived the experiments. JR, FX, QW, KC, YX, LX, SZ, TJ, BZ, YS, and TS conducted the experiments and analyzed the results. JR and QW wrote the manuscript. TC, GH, BZ, and YS revised it critically for important intellectual content. All authors approved the final manuscript and agreed to be accountable for all aspects of the work.

Funding

This study was partly supported by the Bill & Melinda Gates Foundation (INV-005834). The funders had no role in the study design, data collection and analysis, and decision to publish.

Conflict of Interest

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

Publisher's Note

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

Supplementary Material

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

References

1. Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, et al. Early transmission dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. N Engl J Med. (2020) 1199–1207. doi: 10.1056/NEJMoa2001316

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Zhao S, Lin Q, Ran J, Musa SS, Yang G, Wang W, et al. 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. Int J Infect Dis. (2020) 92:214–7. doi: 10.1101/2020.01.23.916395

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Chattu VK, Yaya S Emerging infectious diseases and outbreaks: implications for women's reproductive health and rights in resource-poor settings. Reprod Health. (2020) 17:43. doi: 10.1186/s12978-020-0899-y

PubMed Abstract | CrossRef Full Text | Google Scholar

4. I.L. Organization. Almost 25 million jobs could be lost worldwide as a result of COVID-19, Says ILO. (2020). Available online at: https://www.ilo.org/hanoi/Informationresources/Publicinformation/Pressreleases/WCMS_738992/lang–en/index.htm

5. Iacus SM, Natale F, Santamaria C, Spyratos S, Vespe M. Estimating and projecting air passenger traffic during the COVID-19 coronavirus outbreak and its socio-economic impact. Safety Science. (2020) 129:104791. doi: 10.1016/j.ssci.2020.104791

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Abduljalil JM, Abduljalil BM. Epidemiology, genome, and clinical features of the pandemic SARS-CoV-2: a recent view. New Microbes New Infect. (2020) 35:100672. doi: 10.1016/j.nmni.2020.100672

PubMed Abstract | CrossRef Full Text | Google Scholar

7. C.A, Imai N, Dorigatti I, Baguelin M, Donnelly CA, Riley S, et al. Report 3: Transmissibility of 2019-nCoV. (2020). Available online at: https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/covid-19/report-3-transmissibility-of-covid-19/

Google Scholar

8. Read JM, Bridgen JR, Cummings DA, Ho A, Jewell CP. Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions. Philos. Trans. R. Soc. (2020) 365. doi: 10.1098/rstb.2020.0265

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Riou J, Althaus CL. Pattern of early human-to-human transmission of Wuhan 2019 novel coronavirus (2019-nCoV), December 2019 to January 2020. Eurosurveillance. (2020) 25:7–11. doi: 10.2807/1560-7917.ES.2020.25.4.2000058

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Shen M, Peng Z, Xiao Y, Zhang L. Modelling the epidemic trend of the 2019 novel coronavirus outbreak in China. The Innovation. (2020) 1:100048. doi: 10.1016/j.xinn.2020.100048

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Wu JT, Leung K, Leung GM Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan China: a modelling study. Lancet. (2020) 395:689–97. doi: 10.1016/S0140-6736(20)30260-9

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Luo L, Guo Z, Lei Z, Hu Q, Chen M, Chen F, et al. Epidemiology of tsutsugamushi disease and its relationship with meteorological factors in Xiamen city, China. PLoS Negl Trop Dis. (2020) 14:e0008772. doi: 10.1371/journal.pntd.0008772

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Sanche S, Lin YT, Xu C, Romero-Severson E, Hengartner N, Ke R High contagiousness and rapid spread of severe acute respiratory syndrome coronavirus 2. Emerg Infect Dis. (2020) 26:1470–7. doi: 10.3201/eid2607.200282

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Kucharski AJ, Russell TW, Diamond C, Liu Y, Edmunds J, Funk S, et al. Early dynamics of transmission and control of COVID-19: a mathematical modelling study. Lancet Infect Dis. (2020) 20:553–8. doi: 10.1016/S1473-3099(20)30144-4

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Zhao ZY, Zhu YZ, Xu JW, Hu SX, Hu QQ, Lei Z, et al. A five-compartment model of age-specific transmissibility of SARS-CoV-2. Infect Dis Poverty. (2020) 9:117. doi: 10.1186/s40249-020-00735-x

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Lin SN, Rui J, Chen QP, Zhao B, Yu SS, Li ZY, et al. Effectiveness of potential antiviral treatments in COVID-19 transmission control: a modelling study. Infect Dis Poverty. (2021) 10:53. doi: 10.1186/s40249-021-00835-2

PubMed Abstract | CrossRef Full Text | Google Scholar

17. U.N.D. PROGRAMME, Global Human Development Indicators. Available online at: http://hdr.undp.org/en/2020-report

18. Chen TM, Rui J, Wang QP, Zhao ZY, Cui JA, Yin L. A mathematical model for simulating the phase-based transmissibility of a novel coronavirus. Infect Dis Poverty. (2020) 9:24. doi: 10.1186/s40249-020-00640-3

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Xuying Lao LL, Lei Z, Fang T, Chen Y, Liu Y, Ding K, et al. Epidemiological characteristics and the effectiveness of countermeasures to control coronavirus disease 2019 in Ningbo City, China. Sci Rep. (2021) 11:9545. doi: 10.1038/s41598-021-88473-4

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Mizumoto K, Kagaya K, Zarebski A, Chowell G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Eurosurveillance. (2020) 25:2000180. doi: 10.2807/1560-7917.ES.2020.25.10.2000180

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Miura F, Matsuyama R, Nishiura H. Estimating the asymptomatic ratio of norovirus infection during foodborne outbreaks with laboratory testing in Japan. J Epidemiol. (2018) 28:382–7. doi: 10.2188/jea.JE20170040

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Lauer SA, Grantz KH, Bi Q, Jones FK, Zheng Q, Meredith HR, et al. The incubation period of coronavirus disease 2019 (COVID-19) from publicly reported confirmed cases: estimation and application. Ann Intern Med. (2020) 172:577–82. doi: 10.7326/M20-0504

PubMed Abstract | CrossRef Full Text | Google Scholar

23. C.C.f.D.C.a. Prevention., Investigation and management guide for close contacts of novel coronavirus pneumonia case (Trial Version), 2020.

24. Delamater PL, Street EL, Leslie TF, Yang YT, Jacobsen KH. Complexity of the Basic Reproduction Number (R0). Emerg Infect Dis. (2019) 25:1–4. doi: 10.3201/eid2501.171901

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Acter T, Uddin N, Das J, Akhter A, Choudhury TR, Kim S. Evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as coronavirus disease 2019 (COVID-19) pandemic: A global health emergency. Sci Total Environ. (2020) 730:138996. doi: 10.1016/j.scitotenv.2020.138996

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Lai S, Ruktanonchai NW, Zhou L, Prosper O, Luo W, Floyd JR, et al. Tatem, Effect of non-pharmaceutical interventions to contain COVID-19 in China. Nature. (2020) 585:410–3. doi: 10.1038/s41586-020-2293-x

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Salzberger B, Gluck T, Ehrenstein B. Successful containment of COVID-19: the WHO-Report on the COVID-19 outbreak in China. Infection. (2020) 48:151–3. doi: 10.1007/s15010-020-01409-4

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Verity R, Okell LC, Dorigatti I, Winskill P, Whittaker C, Imai N, et al. Estimates of the severity of coronavirus disease 2019: a model-based analysis. Lancet Infect Dis. (2020). 20:669–77. doi: 10.1016/S1473-3099(20)30243-7

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Costs vs benefits of tackling COVID-19. PharmacoEcon Outcomes News. (2020) p. 852:1. doi: 10.1007/s40274-020-6760-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: COVID-19, transmissibility, mathematical model, the effective reproduction number, epidemic

Citation: Niu Y, Rui J, Wang Q, Zhang W, Chen Z, Xie F, Zhao Z, Lin S, Zhu Y, Wang Y, Xu J, Liu X, Yang M, Zheng W, Chen K, Xia Y, Xu L, Zhang S, Ji R, Jin T, Chen Y, Zhao B, Su Y, Song T, Chen T and Hu G (2021) Containing the Transmission of COVID-19: A Modeling Study in 160 Countries. Front. Med. 8:701836. doi: 10.3389/fmed.2021.701836

Received: 28 April 2021; Accepted: 20 July 2021;
Published: 18 August 2021.

Edited by:

Reza Lashgari, Shahid Beheshti University, Iran

Reviewed by:

Pankaj Bhardwaj, All India Institute of Medical Sciences, Jodhpur, India
Abed Khorasani, Kerman Medical University, Iran
Hadi Choubdar, Shahid Beheshti University of Medical Sciences, Iran

Copyright © 2021 Niu, Rui, Wang, Zhang, Chen, Xie, Zhao, Lin, Zhu, Wang, Xu, Liu, Yang, Zheng, Chen, Xia, Xu, Zhang, Ji, Jin, Chen, Zhao, Su, Song, Chen and Hu. 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: Guoqing Hu, huguoqing009@gmail.com; Tianmu Chen, chentianmu@xmu.edu.cn; 13698665@qq.com

These authors have contributed equally to this work

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