- 1Department of Epidemiology and Health Statistics, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu, China
- 2Sichuan Center for Disease Control and Prevention, Chengdu, China
- 3Thomas E. Starzl Transplantation Institute, University of Pittsburgh Medical Center, Pittsburgh, PA, United States
Introduction: As of June 7, 2021, the outbreak of Coronavirus Disease 2019 (COVID-19) has spread to more than 200 countries. The global number of reported cases is more than 172.9 million, with more than 3.7 million deaths, and the number of infected individuals is still growing rapidly. Consequently, events and activities around the world were canceled or postponed, and the preparation for sporting events were greatly challenged. Under such circumstances, about 11,000 athletes from ~206 countries are arriving in Tokyo for the 32nd Summer Olympic Games. Therefore, it is urgently necessary to assess the occurrence and spread risk of COVID-19 for the Games.
Objectives: To explore effective prevention and control measures for COVID-19 in large international events through simulations of different interventions according to risk assessment.
Methods: We used a random model to calculate the number of initial infected patients and used Poisson distribution to determine the number of initial infected patients based on the number of countries involved. Furthermore, to simulate the COVID-19 transmission, the susceptible-exposed-symptomatic-asymptomatic-recovered-hospitalized (SEIARH) model was established based on the susceptible-exposed-infectious-recovered (SEIR) mathematical model of epidemic diseases. According to risk assessment indicators produced by different scenarios of the simulated interventions, the risk of COVID-19 transmission in Tokyo Olympic Games was assessed.
Results: The current COVID-19 prevention measures proposed by the Japan Olympic Committee need to be enhanced. And large-scale vaccination will effectively control the spread of COVID-19. When the protective efficacy of vaccines is 78.1% or 89.8%, and if the vaccination rate of athletes reaches 80%, an epidemic prevention barrier can be established.
Introduction
The 32nd Summer Olympic Games (Games of the XXII Olympiad) will be held from July 23 to August 8, 2021 in Tokyo, Japan (1). However, athletes and spectators of the Olympics Games are easily exposed to the Coronavirus Disease 2019 (COVID-19). As of June 7, 2021, the outbreak of COVID-19 has spread to more than 200 countries (2). The global epidemic situation remains grim, with more than 172.9 million global cumulative reported cases and more than 3.7 million deaths (3), while the number of infected individuals is still growing rapidly. Consequently, numerous major events and activities around the world have been canceled or postponed, and preparations for all parts of upcoming sport events were greatly challenged. Under such circumstances, about 11,000 athletes from ~206 countries and regions are arriving in Tokyo for the 32nd Summer Olympic during the summer of 2021 (4), along with coaches, referees, and associated International Sport Organization officials, which will definitely increase the possibility of infectious disease outbreaks and transmissions. Although the International Olympic Committee (IOC) has stated that overseas audiences will not be allowed to enter Japan to view the Tokyo Olympic Games (5), it is far from enough to ensure the safety of the Olympic Games in terms of COVID-19. Therefore, it is urgently necessary to commence an evaluation for the risk of COVID-19 under different prevention measures.
So far the COVID-19 prevention measures proposed by the Japan Olympic Committee (JOC) include the following aspects (6, 7):
1. Keep a minimum of two meters from athletes at all times. Keep a minimum of one meter from others. All Games participants must minimize contact within one meter of Games participants who have already been in Japan for more than 14 days, and Japanese residents.
2. During your stay in Japan, you will be expected to limit your activities to what is required to carry out your role.
3. All participants are required to take two COVID-19 tests before their flight to Japan. All other Games participants will be tested daily for three days after their arrival. After the first 3 days and throughout their stay, they will be tested regularly, based on the operational nature of their role and level of contact with athletes.
Thomas Bach, the president of IOC, said that if a vaccine becomes available in time for the July 23–August 8 Games in 2021, the IOC would foot the bill (8). However, since vaccination is voluntary for athletes, different vaccine coverage may produce different effects on an outbreak in the Olympic Village. In light of how this has not been studied so far, we used different vaccine coverage statistics to simulate the transmission of COVID-19 in this study.
In this study, to assess the risk of COVID-19 during the Tokyo Olympic Games, a simulation study based on the transmission dynamic model was carried out. Firstly, we collected the number of athletes from different countries participating in the Tokyo Olympic Games, the current COVID-19 infection probability of each country, and the transmission parameters of the COVID-19 model. Secondly, utilizing the initial number of asymptomatic infections, the number of contacts and other aspects, the susceptible-exposed-symptomatic-asymptomatic-recovered-hospitalized (SEIARH) model was established. Thirdly, in order to carry out risk assessments, the secondary infectors on the 17th day and peak hour of onset were calculated. Through realizing a comparison of the expected risks of COVID-19 under different prevention strategies, this study provided quantitative reference evidence regarding the formulation of COVID-19 prevention and control programs for the Tokyo Olympic Games. The specific process of analysis is shown in Figure 1.
Materials and Methods
The methods of this study consisted of three parts. Firstly, to determine the number of initial infected patients based on the number of countries coming to Japan, Poisson distribution was used. Secondly, based on the COVID-19 transmission mechanism, the SEIARH model was established to simulate COVID-19 transmission. And then we analyzed the robustness of the SEIARH model. Thirdly, according to the results of the SEIARH model, the risk assessment indicators (secondary infectors on the 17th day, peak hour of onset/d) were calculated.
Data Collection and Preparation
National infection probability πi(i = 1, 2, ⋯ , n) was approximated according to public data (3):
where Ii is the number of infections in country i; NC_i is the total population number of country i.
As the Olympic Games are held every 4 years, little change existed in the total number of the adjacent two sessions (9). Besides, the number of athletes was not available before the announcement of IOC, so we used the number of athletes announced in the last Olympic Games as a reference (10). Furthermore, we set the initial asymptomatic infected population A (0) as 10 because the previous imported asymptomatic infected cases did not exceed 10 at one time.
For the parameters of the SEIARH model, we referred to the published COVID-19 classic retrospective study (11–18).
Determining the Number of Initial Infected Patients
We assumed the presence of COVID-19 infectors at entry who were not identified by entry quarantine or health check-up. Here we leveraged the probability of infection πi to approximate the number of overseas import infectors. Furthermore, since the number of participants and the spread of COVID-19 vary across different countries, the initial infected patients referred to unidentified infectors among athletes, coaches, referees, officials, and others who entered Japan for the Tokyo Olympics Games.
The number of infected patients assumed from different regions i(i = 1, 2, ⋯ , n) is Xi. It can be approximately seen that the Poisson distribution of compliance parameter λi = Niπi, which is in the form of:
where Ni is the number of athletes, coaches, referees, and officials from different regions; πi is the infection probability of immigrants from different regions. The infected input probability Pi_imported of a person from region i was:
Furthermore, the total number of initial infected patients could be represented as:
Simulating the Transmission of COVID-19 During the Tokyo Olympic Games
The mathematical models of infectious diseases can be classified into two types by the level of data unit, i.e., the micro-dynamic and macro-dynamic models. The former type of model includes scale-free networks, small-world networks, and so on (19, 20), aiming to show the transmission process of the disease based on the individual level. Thus, it requires a large amount of high-quality personal data. However, since the purpose of this study was to simulate the transmission of COVID-19 on the population level, it was unnecessary and inefficient to collect personal data for model building. Hence, this made the micro-dynamic model unsuitable for our simulation study. On the contrary, the macro-dynamic model provides us with tools for monitoring population flow among different health statuses. For example, the traditional susceptible-exposed-infected-recovered (SEIR) model classifies the population into four categories by its name (21–24). However, for the specific purpose of simulating COVID-19 transmission in this study, neither the difference between asymptomatic and symptomatic patients nor the requirement of isolation treatment was reflected by the SEIR model. To this end, we added both the A (asymptomatic) and H (hospitalized) components to the SEIR model to establish the SEIARH model for the simulation of COVID-19 transmission. The transmission process is shown in Figure 2.
For the simulation study, the discrete-time stochastic compartment model for COVID-19 infection was constructed as:
The parameters are defined in Table 1. As the Japan authority has proposed that “During your stay in Japan, you will be expected to limit your activities to what is required to carry out your role” (6). Therefore, the range of activities of athletes outside the Olympic Village was not considered, and it was reasonable to set the total population number N = 11,000. Based on the historical data of Wuhan, China between January 11 and March 13, we used the Markov Chain Monte Carlo Method (MCMC) to estimate the recovery rate of symptomatic infected individuals. It turned out that the results of data-driven methods were consistent with the real-world results (25), where the 95% CI of γ1 and γ2 was (0.0085, 0.0085) and (0.0085, 0.0085), respectively. More details can be found in Supplementary Material A. In addition, we also initialized the other simulation parameters in Table 1 through literature review (11–18).
Since the Tokyo Olympic Games lasted 17 days, a total of t = 1, 2, ⋯ , 17 data points were simulated. Figures 3–16 show the simulated results. The goal for this simulation was to investigate the transmission of COVID-19 during the Tokyo Olympic Games. According to the general principles of the simulation study design (26), we implemented our simulation study following the details as presented in Table 2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/730611/fpubh-09-730611-HTML/image_m/fpubh-09-730611-g005.jpg)
Figure 5. Vaccinated but not quarantined. Assume pre-flight screening can identify 20% of asymptomatic infectors.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/730611/fpubh-09-730611-HTML/image_m/fpubh-09-730611-g006.jpg)
Figure 6. Vaccinated but not quarantined. Assume pre-flight screening can identify 40% of asymptomatic infectors.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/730611/fpubh-09-730611-HTML/image_m/fpubh-09-730611-g007.jpg)
Figure 7. Vaccination and 7-day compulsory quarantine. Assume pre-flight screening can identify 20% of asymptomatic infectors, 7-day compulsory quarantine has about an 18% risk of missed detection.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/730611/fpubh-09-730611-HTML/image_m/fpubh-09-730611-g008.jpg)
Figure 8. Vaccination and 7-day compulsory quarantine. Assume pre-flight screening can identify 40% of asymptomatic infectors, 7-day compulsory quarantine has about an 18% risk of missed detection.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/730611/fpubh-09-730611-HTML/image_m/fpubh-09-730611-g009.jpg)
Figure 9. Vaccination and 14-day compulsory quarantine. Assume pre-flight screening can identify 20% of asymptomatic infectors, 14-day compulsory quarantine has about a 4% risk of missed detection.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/730611/fpubh-09-730611-HTML/image_m/fpubh-09-730611-g010.jpg)
Figure 10. Vaccination and 14-day compulsory quarantine. Assume pre-flight screening can identify 40% of asymptomatic infectors, 14-day compulsory quarantine has about a 4% risk of missed detection.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/730611/fpubh-09-730611-HTML/image_m/fpubh-09-730611-g011.jpg)
Figure 11. Vaccinated but not quarantined. Assume pre-flight screening can identify 20% of asymptomatic infectors.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/730611/fpubh-09-730611-HTML/image_m/fpubh-09-730611-g012.jpg)
Figure 12. Vaccinated but not quarantined. Assume pre-flight screening can identify 40% of asymptomatic infectors.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/730611/fpubh-09-730611-HTML/image_m/fpubh-09-730611-g013.jpg)
Figure 13. Vaccination and 7-day compulsory quarantine. Assume pre-flight screening can identify 20% of asymptomatic infectors, 7-day compulsory quarantine has about an 18% risk of missed detection.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/730611/fpubh-09-730611-HTML/image_m/fpubh-09-730611-g014.jpg)
Figure 14. Vaccination and 7-day compulsory quarantine. Assume pre-flight screening can identify 40% of asymptomatic infectors, 7-day compulsory quarantine has about an 18% risk of missed detection.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/730611/fpubh-09-730611-HTML/image_m/fpubh-09-730611-g015.jpg)
Figure 15. Vaccination and 14-day compulsory quarantine. Assume pre-flight screening can identify 20% of asymptomatic infectors, 14-day compulsory quarantine has about a 4% risk of missed detection.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/730611/fpubh-09-730611-HTML/image_m/fpubh-09-730611-g016.jpg)
Figure 16. Vaccination and 14-day compulsory quarantine. Assume pre-flight screening can identify 40% of asymptomatic infectors, 14-day compulsory quarantine has about a 4% risk of missed detection.
The current regulations required all participants to take two COVID-19 tests before their flight to Japan, and the number of contacts should be strictly controlled after entry (7, 27, 28). To make the experiment more realistic, we defined the number of susceptible individuals exposed by infected [r1(t)] and exposed [r2(t)] as piecewise functions of time t, where t is the number of days after entry. The function forms were given as follows.
Sensitivity Analysis
In order to verify the robustness of the SEIARH model, we analyzed the sensitivity of the model. The parameter β1 and β2 were increased or decreased by 5% at the same time. These two parameters directly affect the final number of infectors. The cumulative infectors and daily new infectors were observed to determine whether the model was effective.
Assessing the Risks of COVID-19 Under Different Prevention Measures
Based on the results of the SEIARH model, we calculated the risk assessment index of COVID-19: peak hour of onset and secondary infectors on the 17th day. Since the Tokyo Olympic Games lasted for 17 days, secondary infectors on the 17th day were used as indicators. Secondary infectors on the 17th day refers to the number of daily new infectors at the end of the Tokyo Olympic Games; hour of onset is the time t corresponding to the maximum number of daily new infectors.
Results
Transmission Simulation Without Intervention
In this case, we ignored the activities of athletes outside the Olympic Village, and then we used the parameters of the SEIARH model in Table 1 to simulate the transmission within the Olympic Village in Japan without intervention, as shown in Figure 3. The number of daily new secondary infectors reached the peak on the 12th day, which was 1,683. The cumulative infectors were in an S-shaped curve.
The red-shadowed region in Figure 3 is the result of sensitivity analysis. We found that the cumulative infectors and daily new infectors changed little, indicating that the model was robust.
Transmission Simulation Under Current Interventions of the Japan Olympic Committee (JOC)
According to the current prevention and control measures of JOC, the parameters of the average number of infected class contact with susceptible class (r1), average number of infected class contact with exposed class (r2), initial asymptomatic infected population [A(0)] were adjusted, and the specific settings are shown in Table 3.
As can be seen from Figure 4, based on the measures currently proposed by JOC, it was assumed that the screening before take-off can reduce 20 and 40% of asymptomatic infectors, respectively, but the prevention and control effect was not ideal. The total number of secondary infectors still reached the level of no intervention measures, but the duration increased; the peak value of daily new secondary infectors decreased, reaching 1,159 on the 19th day and 1,141 on the 20th day, respectively. Secondary infectors on the 17th day reached 857.6 and 729.7, respectively. Compared with no intervention measures, the situation of current interventions of JOC improved.
The shadowed region in Figure 4 is the result of sensitivity analysis. It can be seen that cumulative infectors and daily new infectors were almost unchanged, indicating that the model was robust. More importantly, through the simulations of the current prevention and control measures, it could be found that it was not enough to rely solely on the existing measures. Other measures (such as vaccination) must be supplemented to effectively reduce the risk of the epidemic during the Tokyo Olympic Games. We will discuss it further in the next section.
Transmission Simulation Under Other Interventions
According to the simulation of intervention measures set in the above sections, we found that the current prevention and control effect was not ideal. But for various infectious diseases, vaccines are the most effective way to eradicate the transmission of infectious diseases. Therefore, based on the intervention measures proposed by JOC, we considered the transmission of athletes after vaccination. The transmission simulation was carried out by adding prevention and control scenarios such as “vaccinated but not quarantined,” “vaccination and 7-day compulsory quarantine,” and “vaccination and 14-day compulsory quarantine.” The protective efficacy of any vaccine cannot reach 100%, and the protective efficacy of the COVID-19 vaccine ranges from 50.7 to 95.0% (29–39), as shown in Table 4. Due to the limitations of clinical trials, vaccine development, and the fact that any vaccines have been used for <1 year, the protective efficacy of the COVID-19 vaccine is still under continuous observation. Therefore, we assumed that the protective efficacy of the vaccine is the median in Table 4, that was 78.1% (95% CI: 64.8–86.3%) (29). So the probability of transmission β1, β2 among athletes vaccinated fell by 21.9% (95% CI: 13.7–35.2%). Variants of COVID-19 can lead to a decrease in neutralizing activity of the vaccine-induced antibodies, especially the mutant strain B.1.351, but the antibodies still possess a certain or high neutralizing ability to the mutant strain (33, 40–70). At present, there is not enough evidence that viral variation has a significant impact on the protection of the new coronavirus vaccine. So, we did not consider the mutations of COVID-19 from different countries. The specific parameter adjustments are shown in Table 5.
It can be found from Figures 5–10 that the assumption of the protective efficacy of the vaccine was 78.1%. Whether isolated or not, screening before take-off can identify 20% or 40% of asymptomatic infectors, which significantly reduced the number of secondary infectors on the 17th day and delayed the peak hour of onset compared to no interventions and current interventions by the JOC. In the comparison of different vaccination rates of various measures, it can be found that when the vaccination rate of athletes reached 80%, the number of secondary infectors on the 17th day decreased significantly. When vaccination rates were 80, 90, and 95%, cumulative infectors still reached a relatively high value at the end, but the first 17 days were very low. It showed that the vaccine was effective, and that extensive vaccination would be useful for the prevention and control of COVID-19 in the Tokyo Olympic Games.
Comparing Figure 5 and Figure 6, Figure 7 and Figure 8, and Figure 9 and Figure 10, it can be seen that under the intervention measures of “vaccinated but not quarantined,” “vaccination and 7-day compulsory quarantine,” and “vaccination and 14-day compulsory quarantine,” the proportion of asymptomatic infected individuals identified by pre-flight screening had little effect on the number of secondarily infected individuals, and the outbreak time was basically the same under different vaccination proportions. However, for the intervention measures in Figures 5–10, as the vaccination rate increased, the secondary infectors on the 17th day and the peak hour of onset significantly decreased. When the vaccination rate of athletes reached 80%, the number of secondary infectors decreased significantly, and it was controlled within 1 when the vaccination rate of athletes reached 90%, indicating that the immune barrier could be established when the vaccination rate reached 80–90%.
Considering that athletes are generally physically fit, we took differences in population resistance into account. We assumed that the protective efficacy of the vaccine of athletes is 15% higher than the average person (78.1%), that is 89.8% (95% CI: 74.5–99.2%). Under this assumption, we developed new simulations. The optimal parameters are shown in Table 6.
From Figures 11–16, we know that when the vaccination rate of athletes reached 80%, the number of secondary infectors on the 17th day decreased significantly and was controlled within 1, indicating that the immune barrier could be established when the vaccination rate reached 80%.
Discussion
The strengths of our study were that the topic is timely, and we used an epidemic method to model the spread of COVID-19. We also discussed the sensitivity of parameters. And we discussed the relationship between vaccination services and spread of COVID-19, considering different protective efficacy of vaccines and vaccination rates. We took differences in population resistance into account as athletes are generally more physically fit. We tested the robustness of the SEIARH model and found that the model was stable. The model we constructed considered the time, environment, and behavior at the same time. The limitations were that we did not consider the transmission of COVID-19 from athletes to the audience, because the policy of JOC clearly limited the scope of activities of athletes and the distance from the audience (6, 7, 71). But we considered the latest results of protective efficacy of the vaccine to make the results more scientific and credible (28, 29). There was no specific number of participants from each country, so it was not possible to accurately estimate the initially infected individuals. The parameters of the model refer to other research based on real data, which may be different from the real situation in Japan, but we set the interval of parameters to reduce bias.
Through the simulation of no intervention measures and various intervention measures, we found that the prevention and control effect of the measures currently proposed by JOC was not ideal. The total number of secondary infectors could still reach the level of no intervention measures, but with a longer duration. However, the prevention and control effect of the intervention measures proposed by us was significantly better than that of no intervention measures or the epidemic response measures proposed by JOC. The simulated number of secondary infectors without intervention measures on the 17th day was 587.8. As shown in Tables 3, 5, 6, we simulated no interventions and different levels of intervention. It showed that these measures can delay the outbreak of COVID-19 to some extent, but the number of infectors was not well-controlled. If the protective efficacy of the vaccine was 78.1% and the vaccination rate reached 90%, the number of secondary infectors under all intervention measures was <1; if the protective efficacy of vaccine was 89.8% and the vaccination rate reached 80%, the number of secondary infectors under all intervention measures was <1, indicating that the epidemic situation was effectively controlled, and vaccination can effectively prevent COVID-19. Comparing Table 5 and Table 6, it was found that under the same intervention circumstance, if the protective efficacy of vaccine was 89.8%, the number of secondary infectors decreased significantly. But the peak hour of onset for each control measure in Table 5 kept being postponed, while it was postponed first and then advanced in Table 6. This was because when the protective efficacy of vaccine was 89.8%, it could not trigger the transmission of COVID-19 and the transmission period was short. Vaccines gave a protection ability to COVID-19 infectors. Studies have found that vaccine uptake was the “main driver” of the decline in COVID-19 infectors rather than lockdown (72). Although there are variant strains of COVID-19, the vaccines still possessed a neutralizing ability to the mutant strain (33, 40–70). Both the Oxford AstraZeneca and Pfizer BioNTech COVID-19 vaccines were effective in reducing the risk of SARS-CoV-2 infection and COVID-19 hospitalization in people with the Delta Variant (55). Because people who participate in the Tokyo Olympic Games will have closed contact (including athletes, staff, etc.), and the vaccine supply can also be guaranteed, we recommend that they should be vaccinated.
As of June 10, 2021, only 4.3% of Japan's population of about 126 million was fully vaccinated against COVID-19, and only 12.6% of the total population was partly vaccinated against COVID-19 (73). At the same time, the Japanese government is trying to complete the task of vaccinating one million times a day, so as to vaccinate all 36 million elderly people in the country before the Tokyo Olympic Games (74). So the prevention and control measures against COVID-19 in the Tokyo Olympic Games need to be further strengthened. In the face of infectious diseases such as COVID-19 with strong infectivity, in addition to the prevention and control measures listed in this study, the combination of multiple intervention measures will bring better results. For example, closed-loop management can be adopted during the Olympic Games, sites can be disinfected before the athletes arrive, and the staff can be trained on epidemic prevention and control.
Risk assessment includes risk identification, dose-response relationship, exposure assessment, etc. COVID-19 is mainly transmitted through direct transmission, aerosol transmission, and contact transmission. Therefore, we started from the potential risks, namely, intervening in the contact between people, then simulated the potential transmission risks, and finally formed a risk assessment for COVID-19. In this study, a simulation method of public health intervention prevention and control based on the dynamic model of COVID-19 was proposed, which can be extended to all kinds of large-scale activities or infectious disease prevention and control research. The intervention measures in this study were based on the epidemic prevention and control measures and ideas proposed by JOC and the IOC. At present, many countries can produce the COVID-19 vaccine, and IOC will also bear the cost of the vaccine, so it is feasible (8). Vaccination is not compulsory, so it does not involve human rights interventions, and it is a complex regulation. People who have been vaccinated still run the risk of infection, and there are variant strains, i.e., Alpha, Beta, Gamma, Delta, etc. (75). Therefore, we could consider reclassifying the population (e.g., adding vaccinated and unvaccinated, or/and according to COVID-19 virus subtype) and combining the latest data on the protective efficacy of vaccines in the future, with the purpose of exploring effective prevention and control measures for COVID-19.
COVID-19 has an incubation period, and there were studies showing that the median incubation period was 4 days (interquartile range, 2–7) (76, 77). So athletes may be infectious after returning, and countries need to take relevant prevention and control measures for athletes, such as nucleic acid detection, isolation, etc.
Conclusions
Based on the dynamic model, this paper simulated different prevention and control measures of the Tokyo Olympic Games, comparing the number of secondary infectors under different measures, and found that vaccination had the best prevention and control effect. When the protective efficacy of the vaccine was 78.1 or 89.8% and if the vaccination rate reached 90%, then the number of secondary infectors can be controlled within 1. Our study will contribute to the formulation of relevant measures by JOC and IOC. In summary, compared with the current public health interventions, mass vaccination will become a milestone in the control of COVID-19.
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
WZ consulted the literature, analyzed the data, wrote the programs, and was a major contributor in writing the manuscript. JF wrote part of the manuscript. CL collected the parameters of the model. HW made comments on the content of the article. YZ collected part of the data. LZ gave advice on setting parameters. XZ made constructive comments on the manuscript. TZ contributed significantly to analysis and manuscript preparation. All authors contributed to the article and approved the submitted version.
Funding
This research work was funded by Sichuan Science and Technology Program (Grant nos. 2020YFS0015, 2020YFS0091, and 2021YFS0001-LH), Health Commission of Sichuan province (Grant no. 20PJ092), National Natural Science Foundation of China (Grant nos. 81602935 and 82041033), Chongqing Science and Technology Program (Grant no. cstc2020jscx-cylhX0003), Central government funding items (Grant no. 2021zc02), and Liangshan Yi Autonomous Prefecture Center for Disease Control and Prevention (Grant no. H210322). The funders played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
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.
Acknowledgments
Many thanks are expressed to those who contributed to the prevention and control of COVID-19. The authors also want to thank Hongli Wan for the submission of the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2021.730611/full#supplementary-material
Supplementary Material A. Introduces a way to estimate parameters by data-driven methods.
Supplementary Material B. Shares our codes.
Abbreviations
COVID-19, coronavirus disease 2019; SEIARH, susceptible-exposed-symptomatic-asymptomatic-recovered-hospitalized; SEIR, susceptible-exposed-infectious-recovered; IOC, International Olympic Committee; JOC, Japan Olympic Committee.
References
1. IOC IPC. Tokyo 2020 Oorganising Committee and Tokyo Metropolitan Government Announce New Dates for the Olympic and Paralympic Games Tokyo 2020. (2020). Available online at: https://www.olympic.org/news/ioc-ipc-tokyo-2020-organising-committee-and-tokyo-metropolitan-government-announce-new-dates-for-the-olympic-and-paralympic-games-tokyo-2020 (accessed March 30, 2020).
2. World Health Organization. Coronavirus Disease (COVID-19) Weekly Epidemiological Update and Weekly Operational Update. (2021). Available online at: https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports (accessed April 10, 2020).
3. National Health Commission of the People's Republic of China. Coronavirus Statistics Reports of the World. (2021). Available online at: https://voice.baidu.com/act/newpneumonia/newpneumonia/?from=osari_aladin_banner#tab4 (accessed March 30, 2021).
4. IOC. What Goes Into Postponing the Olympic Games? (2020). Available online at: https://tokyo2020.org/en/news/what-goes-into-postponing-the-olympic-games (accessed March 30, 2021).
5. IOC. Statement on Overseas Spectators for the Olympic and Paralympic Games Tokyo 2020. (2020). Available online at: https://tokyo2020.org/en/news/statement-on-overseas-spectators-for-the-olympic-and-paralympic-games-tokyo-2020 (accessed March 30, 2021).
6. IOC. First Playbook Published Outlining Measures to Deliver Safe and Successful Games. (2020). Available online at: https://tokyo2020.org/en/news/first-playbook-published-outlining-measures-to-deliver-safe-and-successful-games (accessed March 30, 2021).
7. IOC IPC. Tokyo 2020, Tokyo Metropolitan Government and the Government of Japan. (2020). Available online at: https://olympics.com/tokyo-2020/en/news/joint-statement-by-the-ioc-ipc-tokyo-2020-tokyo-met-gov-and-the-gov-of-japan (accessed June 8, 2021).
8. IOC. Tokyo 2020 to Have ‘Reasonable' Crowd Size, IOC President Bach Says. (2020). Available online at: https://tokyo2020.org/en/news/tokyo-2020-to-have-reasonable-crowd-size-ioc-president-bach-says (accessed March 30, 2021).
9. Wikipedia. Summer Olympic Games. (2021). Available online at: https://en.wikipedia.org/wiki/Summer_Olympic_Games#Popularity_of_Olympic_sports (accessed March 30, 2021).
10. Wikipedia. 2016 Summer Olympics_wikipedia. (2021). Available online at: https://en.wikipedia.org/wiki/2016_Summer_Olympics#Number_of_athletes_by_National_Olympic_Committee (accessed March 30, 2021).
11. Yang ZF, Zeng ZQ, Wang K, Wong SS, Liang WH, Zanin M, et al. Modified SEIR and AI prediction of the epidemics trend of COVID-19 in China under public health interventions. J Thorac Dis. (2020) 12:165–74. doi: 10.21037/jtd.2020.02.64
12. Peak CM, Kahn R, Grad YH, Childs LM, Li R, Lipsitch M, Buckee CO. Individual quarantine versus active monitoring of contacts for the mitigation of COVID-19: a modelling study. Lancet Infect Dis. (2020) 20:1025–33. doi: 10.1016/S1473-3099(20)30361-3
13. Song QQ, Zhao H, Fang LQ, Liu Y, Zheng C, Zhang Y. Study on assessing early epidemiological parameters of COVID-19 epidemic in China. Chinese J Epidemiol. (2020) 4:461–5. doi: 10.3760/cma.j.cn112338-20200205-00069
14. Wei YY, Lu ZZ, Du ZC, Zhang ZJ, Zhao Y, Sheng SP, et al. Fitting and forecasting the trend of COVID-19 by SEIR+CAQ dynamic model. Chinese J Epidemiol. (2020) 4:470–5. doi: 10.3760/cma.j.cn112338-20200216-00106
15. Kinoshita R, Anzai A, Jung SM, Linton NM, Miyama T, Kobayashi T, et al. Containment, contact tracing and asymptomatic transmission of novel coronavirus disease (COVID-19): a modelling study. J Clin Med. (2020) 9:3125 doi: 10.3390/jcm9103125
16. 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
17. Tang B, Wang X, Li Q, Bragazzi NL, Tang SY, Xiao YN, et al. Estimation of the Transmission Risk of the COVID-19 and Its Implication for Public Health Interventions. J Clin Med. (2020) 9:462. doi: 10.2139/ssrn.3525558
18. Tang B, Xia F, Bragazzi NL, Wang X, He S, Sun XD, et al. Lessons drawn from China and South Korea for managing COVID-19 epidemic: insights from a comparative modeling study. Bull World Health Organ. (2020) 98:150. doi: 10.2471/BLT.20.257238
19. Lloyd AL, May RM. Epidemiology. How viruses spread among computers and people. Science. (2001) 292:1316–7. doi: 10.1126/science.1061076
20. Keeling MJ, Eames K. Networks and epidemic models[J]. J R Soc Interface. (2005) 2:295–307. doi: 10.1098/rsif.2005.0051
21. Hethcote HW. The mathematics of infectious diseases. SIAM Rev. (2000) 42:599–653. doi: 10.1137/S0036144500371907
22. Brauer F. Mathematical epidemiology: past, present, and future. Infect Dis Model. (2017) 2:113–27. doi: 10.1016/j.idm.2017.02.001
23. He SB, Peng YX, Sun KH. SEIR modeling of the COVID-19 and its dynamics. Nonlinear Dyn. (2020) 101:1667–80. doi: 10.1007/s11071-020-05743-y
24. Keeling MJ, Rohnai P. Modeling Infectious Diseases in Humans and Animals. Princeton University Press (2008). doi: 10.1515/9781400841035
25. He S, Tang SY, Rong L. A discrete stochastic model of the COVID-19 outbreak: forecast and control. Math Biosci Eng. (2020) 17:2792–804. doi: 10.3934/mbe.2020153
26. Morris TP, White IR, Crowther MJ. Using simulation studies to evaluate statistical methods. Stat Med. (2019) 38:2074–102. doi: 10.1002/sim.8086
27. IOC. Spectators, Athletes and Games Staff: Plans in Place to Ensure Successful Olympic and Paralympic Games Tokyo. (2020). Available online at: https://tokyo2020.org/en/news/plans-in-place-to-ensure-successful-tokyo-2020 (accessed March 30, 2021).
28. IOC. Fifth Coordination Meeting for COVID-19 Countermeasures at the Olympic and Paralympic Games Tokyo 2020. (2020). Available online at: https://tokyo2020.org/en/news/fifth-coordination-meeting-for-covid-19-countermeasures-at-tokyo-2020 (accessed March 30, 2021).
29. Al KN, Zhang YT, Xia SL, Yang YK, Al QMM, Abdulrazzaq N, et al. Effect of 2 inactivated SARS-CoV-2 vaccines on symptomatic COVID-19 infection in adults: a randomized clinical trial. JAMA. (2021) 326:35–45. doi: 10.1001/jama.2021.8565
30. Palacios R, Batista AP, Albuquerque CSN, Patio EG, Santos JDP, Mnica TRPC, et al. Efficacy and safety of a COVID-19 inactivated vaccine in healthcare professionals in Brazil: the PROFISCOV Study. Soc Sci Electron Publ. (2021). doi: 10.2139/ssrn.3822780
31. Thompson MG, Burgess JL, Naleway AL, Tyner H, Yoon SK, Meece J, et al. Prevention and attenuation of Covid-19 with the BNT162b2 and mRNA-1273 vaccines. N Engl J Med. (2021) 385:320–9. doi: 10.1056/NEJMoa2107058
32. Chung H, He S, Nasreen S, Sundaram ME, Buchan SA, Wilson SE, et al. Effectiveness of BNT162b2 and mRNA-1273 COVID-19 vaccines against symptomatic SARS-CoV-2 infection and severe COVID-19 outcomes in Ontario, Canada. BMJ. (2021) 374:n1943. doi: 10.1136/bmj.n1943
33. Sadoff J, Gray G, Vandebosch A, Cárdenas V, Shukarev G, Grinsztejn B, et al. Safety and efficacy of single-dose Ad26.COV2.S vaccine against Covid-19. N Engl J Med. (2021) 384:2187–201. doi: 10.1056/NEJMoa2101544
34. Voysey M, Clemens SAC, Madhi SA, Weckx LY, Folegatti PM, Aley PK, et al. Safety and efficacy of the ChAdOx1 nCoV-19 vaccine (AZD1222) against SARS-CoV-2: an interim analysis of four randomised controlled trials in Brazil, South Africa, and the UK. Lancet. (2021) 397:99–111. doi: 10.1016/S0140-6736(20)32661-1
35. Logunov DY, Dolzhikova IV, Shcheblyakov DV, Tukhvatulin AI, Zubkova OV, Dzharullaeva AS, et al. Safety and efficacy of an rAd26 and rAd5 vector-based heterologous prime-boost COVID-19 vaccine: an interim analysis of a randomised controlled phase 3 trial in Russia. Lancet. (2021) 397:671–81. doi: 10.1016/S0140-6736(21)00234-8
36. Turner JS, Kim W, Kalaidina E, Goss CW, Rauseo AM, Schmitz AJ, et al. SARS-CoV-2 infection induces long-lived bone marrow plasma cells in humans. Nature. (2021) 595:421–5. doi: 10.1038/s41586-021-03647-4
37. Pušnik J, Richter E, Schulte B, Dolscheid-Pommerich R, Bode C, Putensen C, et al. Memory B cells targeting SARS-CoV-2 spike protein and their dependence on CD4+ T cell help. Cell Rep. (2021) 35:109320. doi: 10.1016/j.celrep.2021.109320
38. Doria-Rose N, Suthar MS, Makowski M, O'Connell S, McDermott AB, Flach B, et al. Antibody persistence through 6 months after the second dose of mRNA-1273 vaccine for Covid-19. N Engl J Med. (2021) 384:2259–61. doi: 10.1056/NEJMc2103916
39. Dan JM, Mateus J, Kato Y, Hastie KM, Yu ED, Faliti CE, et al. Immunological memory to SARS-CoV-2 assessed for up to 8 months after infection. Science. (2021) 371:6529. doi: 10.1126/science.abf4063
40. Muik A, Wallisch AK, Sänger B, Swanson KA, Mühl J, Chen W, et al. Neutralization of SARS-CoV-2 lineage B.1.1.7 pseudovirus by BNT162b2 vaccine-elicited human sera. Science. (2021) 371:1152–3. doi: 10.1126/science.abg6105
41. Gu H, Chen Q, Yang G, He L, Fan H, Deng YQ, et al. Adaptation of SARS-CoV-2 in BALB/c mice for testing vaccine efficacy. Science. (2020) 369:1603–7. doi: 10.1126/science.abc4730
42. Xie X, Zou J, Fontes-Garfias CR, Xia HJ, Swanson KA, Cutler M, et al. Neutralization of N501Y mutant SARS-CoV-2 by BNT162b2 vaccine-elicited sera. bioRxiv [Preprint]. (2021). doi: 10.1101/2021.01.07.425740
43. Garcia-beltran WF, Lam EC, Denis KS, Nitido AD, Garcia ZH, Hauser BM, et al. Multiple SARS-CoV-2 variants escape neutralization by vaccine-induced humoral immunity. Cell. (2021) 184:2372–83.e2379. doi: 10.1016/j.cell.2021.03.013
44. Supasa P, Zhou D, Dejnirattisai W, Liu C, Mentzer AJ, Ginn HM, et al. Reduced neutralization of SARS-CoV-2 B.1.1.7 variant by convalescent and vaccine sera. Cell. (2021) 184:2201–11.e2207. doi: 10.1016/j.cell.2021.02.033
45. Ho D, Wang PF, Liu LH, Iketani S, Luo Y, Guo YC, et al. Increased resistance of SARS-CoV-2 variants B.1.351 and B.1.1.7 to antibody neutralization. Res Sq. (2021) 372:1413–8. doi: 10.21203/rs.3.rs-155394/v1
46. Wall EC, Wu M, Harvey R, Kelly G, Warchal S, Sawyer C, et al. Neutralising antibody activity against SARS-CoV-2 VOCs B.1.617.2 and B.1.351 by BNT162b2 vaccination. Lancet. (2021) 397:2331–3. doi: 10.1016/S0140-6736(21)01290-3
47. Dejnirattisai W, Zhou D, Supasa P, Liu C, Mentzer AJ, Ginn HM, et al. Antibody evasion by the P.1 strain of SARS-CoV-2. Cell. (2021) 184:2939–54.e2939. doi: 10.1016/j.cell.2021.03.055
48. Anichini G, Terrosi C, Gori SG, Gandolfo C, Franchi F, Cusi MG, et al. Neutralizing antibody response of vaccinees to SARS-CoV-2 variants. Vaccines. (2021) 9:517. doi: 10.3390/vaccines9050517
49. Zhou D, Dejnirattisai W, Supasa P, Liu C, Mentzer AJ, Ginn HM, et al. Evidence of escape of SARS-CoV-2 variant B.1.351 from natural and vaccine-induced sera. Cell. (2021) 184:2348–61.e2346. doi: 10.1016/j.cell.2021.02.037
50. Wang P, Casner RG, Nair MS, Wang M, Yu J, Cerutti G, et al. Increased resistance of SARS-CoV-2 variant P.1 to antibody neutralization. Cell Host Microbe. (2021) 29:747–51. doi: 10.1016/j.chom.2021.04.007
51. Abu RLJ, Chemaitelly H, Butt AA. Effectiveness of the BNT162b2 Covid-19 vaccine against the B.1.1.7 and B.1.351 variants. N Engl J Med. (2021) 385:187–9. doi: 10.1056/NEJMc2104974
52. Bernal JL, Andrews N, Gower C, Gallagher E, Simmons R, Thelwall S, et al. Effectiveness of COVID-19 vaccines against the B.1.617.2 variant. N Engl J Med. (2021) 385:585–94. doi: 10.1056/NEJMoa2108891
53. Lopez BJ, Andrews N, Gower C, Robertson C, Stowe J, Tessier E, et al. Effectiveness of the Pfizer-BioNTech and Oxford-AstraZeneca vaccines on covid-19 related symptoms, hospital admissions, and mortality in older adults in England: test negative case-control study. BMJ. (2021) 373:n1088. doi: 10.1136/bmj.n1088
54. Vasileiou E, Simpson CR, Robertson C, Shi T, Kerr S, Agrawal U, et al. Effectiveness of first dose of COVID-19 vaccines against hospital admissions in Scotland: national prospective cohort study of 5.4 million people. SSRN Electron J. (2021). doi: 10.2139/ssrn.3789264
55. Sheikh A, Mcmenamin J, Taylor B, Robertson C, Public Health Scotland and the EAVE II Collaborators. SARS-CoV-2 Delta VOC in Scotland: demographics, risk of hospital admission, and vaccine effectiveness. Lancet. (2021) 397:2461–2. doi: 10.1016/S0140-6736(21)01358-1
56. Wu K, Werner AP, Moliva JI, Koch M, Choi A, Stewart-Jones GBE, et al. mRNA-1273 vaccine induces neutralizing antibodies against spike mutants from global SARS-CoV-2 variants. bioRxiv [Preprint]. (2021). doi: 10.1101/2021.01.25.427948
57. Choi A, Koch M, Wu K, Dixon G, Oestreicher J, Legault H, et al. Serum neutralizing activity of mRNA-1273 against SARS-CoV-2 Variants. bioRxiv [Preprint]. (2021). doi: 10.1101/2021.06.28.449914
58. Emary KRW, Golubchik T, Aley PK, Ariani CV, Angus B, Bibi S, et al. Efficacy of ChAdOx1 nCoV-19 (AZD1222) vaccine against SARS-CoV-2 variant of concern 202012/01 (B.1.1.7): an exploratory analysis of a randomised controlled trial. Lancet. (2021) 397:1351–62.
59. Madhi SA, Baillie V, Cutland CL, Voysey M, Koen AL, Fairlie L, et al. Efficacy of the ChAdOx1 nCoV-19 Covid-19 vaccine against the B.1.351 variant. N Engl J Med. (2021) 384:1885–98. doi: 10.1056/NEJMoa2102214
60. Planas D, Vever D, BaidaliukA, Staropoli I, Guivel-Benhassine F, Rajah MM, et al. Reduced sensitivity of infectious SARS-CoV-2 variant B.1.617.2 to monoclonal antibodies and sera from convalescent and vaccinated individuals. bioRxiv [Preprint]. (2021). doi: 10.1101/2021.05.26.445838
61. Katherine RWE, Tanya G, Parvinder KA, Cristina VA, Brian A, Sagida B, et al. Efficacy of ChAdOx1 nCoV-19 (AZD1222) vaccine against SARS-CoV-2 variant of concern 202012/01 (B.1.1.7): an exploratory analysis of a randomised controlled trial. Lancet. (2021) 397:1351–62. doi: 10.1016/S0140-6736(21)00628-0
62. Hossain MK, Hassanzadeganroudsari M, Apostolopoulos V. The emergence of new strains of SARS-CoV-2. What does it mean for COVID-19 vaccines? Expert Rev Vaccines. (2021) 20:635–8. doi: 10.1080/14760584.2021.1915140
63. Alter G, Yu J, Liu J, Chandrashekar A, Borducchi EN, Tostanoski LH, et al. Immunogenicity of Ad26.COV2.S vaccine against SARS-CoV-2 variants in humans. Nature. (2021) 596:268–72. doi: 10.1038/s41586-021-03681-2
64. Yu J, Tostanoski LH, Mercado NB, McMahan K, Liu J, Jacob-Dolan C, et al. Protective efficacy of Ad26.COV2.S against SARS-CoV-2 B.1.351 in macaques. Nature. (2021) 596:423–7. doi: 10.1038/s41586-021-03732-8
65. AbdoolK SS, De OT. New SARS-CoV-2 variants - clinical, public health, and vaccine implications. N Engl J Med. (2021) 384:1866–8. doi: 10.1056/NEJMc2100362
66. Toback S, Heath PT, Galiza EP, Baxter DN, Boffito M, Smith K, et al. Efficacy of the NVX-CoV2373 Covid-19 vaccine against the B.1.1.7 Varian. medRxiv [Preprint]. (2021). doi: 10.1101/2021.05.13.21256639
67. Shinde V, Bhikha S, Hoosain Z, Archary M, Bhorat Q, Fairlie L, et al. Efficacy of NVX-CoV2373 Covid-19 vaccine against the B.1.351 variant. N Engl J Med. (2021) 384:1899–909. doi: 10.1056/NEJMoa2103055
68. Wang GL, Wang ZY, Duan LJ, Meng QC, Jiang MD, Cao J, et al. Susceptibility of circulating SARS-CoV-2 variants to neutralization. N Engl J Med. (2021) 384:2354–56. doi: 10.1056/NEJMc2103022
69. Huang BY, Dai LP, Wang H, Hu ZY, Yang XM, Tan WJ, et al. Neutralization of SARS-CoV-2 VOC 501Y.V2 by human antisera elicited by both inactivated BBIBP-CorV and recombinant dimeric RBD ZF2001 vaccines. bioRxiv [Preprint]. (2021). doi: 10.1101/2021.02.01.429069
70. Hitchings MDT, Ranzani OT, Torres MSS, de Oliveira SB, Almiron M, Said R, et al. Effectiveness of CoronaVac among healthcare workers in the setting of high SARS-CoV-2 Gamma variant transmission in Manaus, Brazil: a test-negative case-control study. Lancet Reg Health Am. (2021) 1:100025. doi: 10.1101/2021.04.07.21255081
71. IOC. Safe and Secure Olympic and Paralympic Games Tokyo 2020: Version 2 of Tokyo 2020 Playbooks Released. (2020). Available online at: https://olympics.com/tokyo-2020/en/news/safe-and-secure-oly-and-para-games-tokyo-2020-v2-of-tokyo-2020-playbooks (accessed June 6, 2021).
72. Haas EJ, Angulo FJ, McLaughlin JM, Anis E, Singer SR, Khan F, et al. Impact and effectiveness of mRNA BNT162b2 vaccine against SARS-CoV-2 infections and COVID-19 cases, hospitalisations, and deaths following a nationwide vaccination campaign in Israel: an observational study using national surveillance data. Lancet. (2021) 397:1819–29. doi: 10.1016/S0140-6736(21)00947-8
73. Mathieu E, Ritchie H, Ortiz-Ospina E, Roser M, Hasell J, Appel C, et al. A global database of COVID-19 vaccinations. Nat Hum Behav. (2021) 5:956–9. doi: 10.1101/2021.03.22.21254100
74. Netease times weekly. No Tokyo Olympics? Kan Yiwei Said That If the Epidemic Is Serious, It Will Not Be “Forced”. (2021). Available online at: https://www.163.com/dy/article/GC82UVHV0519APGA.html?clickfrom=w_yw (accessed June 12, 2021).
75. World Health Organization. WHO Announces Simple, Easy-To-Say Labels for SARS-CoV-2 Variants of Interest and Concern. (2021). Available online at: https://www.who.int/news/item/31-05-2021-who-announces-simple-easy-to-say-labels-for-sars-cov-2-variants-of-interest-and-concern (accessed August 22, 2021).
76. 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
Keywords: COVID-19, SEIARH model, interventions, Tokyo Olympic Games, risk assessment
Citation: Zhu W, Feng J, Li C, Wang H, Zhong Y, Zhou L, Zhang X and Zhang T (2021) COVID-19 Risk Assessment for the Tokyo Olympic Games. Front. Public Health 9:730611. doi: 10.3389/fpubh.2021.730611
Received: 25 June 2021; Accepted: 14 September 2021;
Published: 25 October 2021.
Edited by:
Reza Lashgari, Shahid Beheshti University, IranReviewed by:
Haijun Zhang, Harbin Institute of Technology, ChinaDr. Varun Tyagi, Eurofins Agroscience Services Pvt. Ltd., India
Khondoker Nazmoon Nabi, Bangladesh University of Engineering and Technology, Bangladesh
Copyright © 2021 Zhu, Feng, Li, Wang, Zhong, Zhou, Zhang and Zhang. 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: Tao Zhang, c3RhdHpoYW5ndGFvQHNjdS5lZHUuY24=