Loading [MathJax]/jax/output/SVG/jax.js
Research article

Infection spread simulation technology in a mixed state of multi variant viruses

  • Received: 16 August 2021 Accepted: 01 November 2021 Published: 05 November 2021
  • ATLM (Apparent Time Lag Model) was extended to simulate the spread of infection in a mixed state of the variant virus and original wild type. It is applied to the 4th wave of infection spread in Tokyo, and (1) the 4th wave bottoms out near the end of the state of emergency, and the number of infected people increases again. (2) The rate of increase will be mainly by d strain (L452R) virus, while the increase by a strain (N501Y) virus will be suppressed. (3) It is anticipated that the infection will spread during the Olympic Games. (4) When variant viruses compete, the infection of highly infectious virus rises sharply while the infection by weakly infectious ones has converged. (5) It is effective as an infection control measure to find an infected person early and shorten the period from infection to quarantine by PCR test or antigen test as a measure other than the vaccine.

    Citation: Makoto Koizumi, Motoaki Utamura, Seiichi Kirikami. Infection spread simulation technology in a mixed state of multi variant viruses[J]. AIMS Public Health, 2022, 9(1): 17-25. doi: 10.3934/publichealth.2022002

    Related Papers:

    [1] Mario Coccia . Sources, diffusion and prediction in COVID-19 pandemic: lessons learned to face next health emergency. AIMS Public Health, 2023, 10(1): 145-168. doi: 10.3934/publichealth.2023012
    [2] Ali Moussaoui, El Hadi Zerga . Transmission dynamics of COVID-19 in Algeria: The impact of physical distancing and face masks. AIMS Public Health, 2020, 7(4): 816-827. doi: 10.3934/publichealth.2020063
    [3] Saina Abolmaali, Samira Shirzaei . A comparative study of SIR Model, Linear Regression, Logistic Function and ARIMA Model for forecasting COVID-19 cases. AIMS Public Health, 2021, 8(4): 598-613. doi: 10.3934/publichealth.2021048
    [4] Rafat Zreiq, Souad Kamel, Sahbi Boubaker, Asma A Al-Shammary, Fahad D Algahtani, Fares Alshammari . Generalized Richards model for predicting COVID-19 dynamics in Saudi Arabia based on particle swarm optimization Algorithm. AIMS Public Health, 2020, 7(4): 828-843. doi: 10.3934/publichealth.2020064
    [5] Robin M. Kowalski, Kenzie Hurley, Nicholas Deas, Sophie Finnell, Kelly Evans, Chelsea Robbins, Andrew Cook, Emily Radovic, Hailey Carroll, Lyndsey Brewer, Gabriela Mochizuki . Protection motivation unmasked: Applying protection motivation theory to skepticism toward COVID-19 mask and vaccine mandates. AIMS Public Health, 2022, 9(3): 506-520. doi: 10.3934/publichealth.2022035
    [6] Casey T. Harris, Kevin Fitzpatrick, Michael Niño, Priya Thelapurath, Grant Drawve . Examining disparities in the early adoption of Covid-19 personal mitigation across family structures. AIMS Public Health, 2022, 9(3): 589-605. doi: 10.3934/publichealth.2022041
    [7] Kottakkaran Sooppy Nisar, Muhammad Wajahat Anjum, Muhammad Asif Zahoor Raja, Muhammad Shoaib . Recurrent neural network for the dynamics of Zika virus spreading. AIMS Public Health, 2024, 11(2): 432-458. doi: 10.3934/publichealth.2024022
    [8] Srikanth Umakanthan, Anuradha Chauhan, Madan Mohan Gupta, Pradeep Kumar Sahu, Maryann M Bukelo, Vijay Kumar Chattu . COVID-19 pandemic containment in the Caribbean Region: A review of case-management and public health strategies. AIMS Public Health, 2021, 8(4): 665-681. doi: 10.3934/publichealth.2021053
    [9] Musyoka Kinyili, Justin B Munyakazi, Abdulaziz YA Mukhtar . Mathematical modeling and impact analysis of the use of COVID Alert SA app. AIMS Public Health, 2022, 9(1): 106-128. doi: 10.3934/publichealth.2022009
    [10] Chun Nok Lam, William Nicholas, Alejandro De La Torre, Yanpui Chan, Jennifer B. Unger, Neeraj Sood, Howard Hu . Factors associated with parents' willingness to vaccinate their children against COVID-19: The LA pandemic surveillance cohort study. AIMS Public Health, 2022, 9(3): 482-489. doi: 10.3934/publichealth.2022033
  • ATLM (Apparent Time Lag Model) was extended to simulate the spread of infection in a mixed state of the variant virus and original wild type. It is applied to the 4th wave of infection spread in Tokyo, and (1) the 4th wave bottoms out near the end of the state of emergency, and the number of infected people increases again. (2) The rate of increase will be mainly by d strain (L452R) virus, while the increase by a strain (N501Y) virus will be suppressed. (3) It is anticipated that the infection will spread during the Olympic Games. (4) When variant viruses compete, the infection of highly infectious virus rises sharply while the infection by weakly infectious ones has converged. (5) It is effective as an infection control measure to find an infected person early and shorten the period from infection to quarantine by PCR test or antigen test as a measure other than the vaccine.



    Mutation of SARS-CoV-2 has occurred in the infected areas of the United Kingdom, Brazil, and India over the past year. These variant viruses have also been brought to Japan. Variant viruses include those with weakened infectivity and those with increased infectivity. In the situation where several variant viruses coexist, it is considered that highly infectious virus trend to become the mainstream of the epidemic. In order to suppress the infection spreads, it is necessary to take measures depending on the virus. In constructing infection control measures, it is necessary to predict whether the infection will spread or shrink. To predict the infection spread, several calculation models were proposed.

    First, we take a look at the prediction methods that have been used to date. SIR [1] and SEIR [2] models are often used in early infection stage. Because, they are simple mathematical structure and require short calculation time. As the progressing of infection spreads, measures such as vaccination and lockdown will be taken. Kuniya et al. use SEIR to evaluate the effect of the SOE (State Of Emergency) in the second wave in Tokyo. They conducted a parameter survey with varying coefficients of the equation [3]. Britton et al. applied improved SEIR to spread infection under the case of non-uniform population structure [4]. Muñoz-Fernández et al. applied the modified SIR model to analyze the wave of COVID-19. They used nonconstant parameter [5]. Biala et al. improved SEIR model to calculate the spread of COVID-19 pandemic [6]. On the other hand, ABM (Agent Based Model) has been developed [7][11]. This technology is a probabilistic method, and unlike the deterministic methods such as SIR and SEIR, it is a method that assumes various behaviors of a person and calculates the infection probability, and it takes a considerable amount of calculation time. Many calculation methods as mentioned above are constructed assuming a single virus, and do not consider the mixture of variant viruses. The objective of our research is to develop a technology for predicting the spread of infection in a mixed state of variant viruses.

    Next, we briefly describe the progress of research so far. The above methods do not have a time delay from infection to quarantine. We considered that the time required until isolated from infected is the important role of contribution in expanding infection, therefore we developed ATLM (Apparent Time Lag Model) with a delay until isolation time [12]. This model currently has an extended version with vaccine and lockdown effects [13]. We have expanded it to handle variant viruses. The infectivity of variant viruses has already been reported [14]. We use these data to simulate the fourth wave of infection spread in Tokyo and investigate the availability of the method.

    The ATLM we have developed [12] uses the following equation, which takes into account the time delay from infection to quarantine and the time delay from infection to loss of infectivity. We denote the cumulative number of infected people by x as unknown, and daily infected people is by dx/dt.

    dxdt=α(t){x(t)(1ϵ)x(tT)ϵx(tS)}(1x(t)M(t))

    α(t)=α0ρ(t)

    M(t)=M0(1µ(t))

    where, T: time delay from infection to quarantine, µ(t): vaccination rate, α: infectivity, ϵ: ratio of asymptomatic persons, S is time delay from infection to the extinction of infectivity. M indicates the sensitive population. ρ(t) is the rate of decrease in infectivity due to the restriction of human flow such as a lockdown. α0 is original infectivity of the virus. Equation (3) represents the decrease in the sensitive population due to vaccination. Subscript 0 indicate initial value. Details of these equations were shown in previous paper [12].

    Number of quarantined persons in y(t) and number of infected people during infection isolation in a community z(t) can be calculated by Eqs. (4) and (5) respectively.

    y(t)=(1ϵ)(x(tT)x(tS))

    z(t)=x(t)(1ϵ)x(tT)

    To extend the above equations for handling variant viruses, the following assumptions are taken into account.

    (1) Infected people are infected with only one type of virus, and there is no simultaneous infection.

    (2) Patients who have been infected with one variant virus in the past are not infected with another variant virus.

    (3) The infection rates between viruses are independent of each other and do not interfere with each other.

    (4) The effect of the vaccine is the same for each variant virus.

    (5) Both delay times until the onset and the infectivity disappear are the same for each variant virus.

    Under the above assumptions, cumulative number of people infected by variant virus i is expressed by subscript i. the differential Eq. (1) is rewritten as follows.

    dxidt=αi(t){xi(t)(1ϵ)xi(tT)ϵxi(tS)}(1X(t)M(t))

    X(t)=ixi(t)

    where ai is infectivity of variant virus i. Equation (7) is a limitation induced from assumptions (1) and (2). That is, it is shown that variant viruses have the common sensitive population M. Number of quarantined persons in each virus yi(t) and number of infected people during infection isolation in a community zi(t) are also able to be calculated by following equations.

    yi(t)=(1ϵ)(xi(tT)xi(tS))

    zi(t)=xi(t)(1ϵ)xi(tT)

    The analytical solution of the differential Eq. (4) is unknown. Therefore, to solve the above equations, the 4th-order Runge-Kutta method was used for numerical integration. In Eq. (4), the numerical values of xi(t-T ) and xi(t-S) at time t have been already computed and there is no problem in accuracy. However, since these values are not calculated at the start of the calculation, precaution must be taken at the start time. Therefore, the initial value is given a sufficiently small value compared to M0. Next, it should be noted that, when X(t) is increased, too close to M(t). Especially when the vaccination rate becomes high, X(t) > M(t) may occur, in which case the solution oscillates and becomes unstable. To avoid numerical unstable, if X(t)/M(t) > 1, we set the right-hand side of Eq. (4) equal zero.

    Table 1.  Calculation conditions.
    α strain (N501Y) δ strain (L452R)
    Initial Infected People 4400 5.5
    Infectivity 0.09 0.121
    Sensitive Population 500000
    Start Point 2021/3/1
    Vaccination Start 2021/5/15
    Vaccination Rate 0.42%/day
    Start of SOE 2021/4/25
    End of SOE 2021/6/20
    Time Delay until Quarantine 14 days
    Quarantine Period 14 days
    Ratio of Asymptomatic Persons 0

    Note: SOE: State Of Emergency.

     | Show Table
    DownLoad: CSV

    Variant virus that is prevalent in Tokyo and is seen as almost by a strain virus (N501Y) [14]. Currently, the virus of interest is d strain (L452R) found in India. Infectivity of a strain prevalent in Tokyo is believed as 1.32 times that of the original wild type [14]. Infectivity of d strain is estimated at 1.78 times higher than original ones [15]. Thus, the ratio between two becomes about 1.35. We set the infectivity of each virus based on this ratio. It is confirmed that 9 people are infected at 2021/5/31 by d strain. Its sampling percentage of infected people is 10%, therefore, about 90 people are infected [15]. Table 1 shows calculation conditions. The initial value by d strain was determined so as to satisfy the above conditions. See appendix A. As shown in the table, the effect of the vaccine is incorporated. The SOE (State Of Emergency) was scheduled by Japanese Government and Tokyo Metropolitan Government.

    Figure 1 shows the pattern of the 4th wave in Tokyo. The daily change in the number of infected people is illustrated in (a). The origin of the horizontal axis is set at 2021/3/1. The epidemic peak is located at around fifty days from the calculation start (the beginning of May or the end of April), about 800 people infected persons have been calculated. This number is roughly equivalent to the actual 7-day average for the 4th wave. (b) shows the number of quarantined persons calculated by Eq. (6), including home medical treatment and hotel medical treatment. It is said that 80% of the infected persons are mild according to the WHO. Hence, we estimated that remaining 20% would be in the hospital. Therefore, at the peak more 11,000 people has been quarantined and about 2200 people is considered in the hospital. According to the data of Tokyo [15] at the time of the fourth wave peak about 2400 patients were hospitalized, then the results are consistent with the actual data. (c) displays number of infected people until isolated in a community calculated by Eq. (9). The higher this number, the higher the probability of having the next infected person. The average infectivity of the variant virus is shown in (d), and as the infection progresses, the average value becomes closer to the infectivity of strain. It shows that the strain is becoming dominant.

    Figure 1.  The pattern of 4th wave in Tokyo.

    The changing number of infected people by each virus are displayed in Figure 2. (a), (b) and (c) show the number of daily infected persons, quarantined persons, and infected persons in the community respectively. Solid line is the number of people infected by a, and broken line is due to d. Each figure shows that the patients by d has increased sharply after the end of the SOE. (d) in Figure 2 shows a ratio of patients by a and d. The patients by d increases from the SOE declared the end and becomes dominant after the point of 128 days (2021/7/6) from the calculation start. In addition, from these figures, it can be seen that the infection by d rises sharply at the stage when the infection of a has converged and bottomed out. As described above, this analysis also indicated that the one with stronger infectivity became dominant when the infection spreads.

    Figure 2.  Changing number of infected people by each virus.

    We have already examined the sensitivity about time delay T in the previous work [12] and reproductive number R is proportion to product of infectivity α and time delay T. Therefor less T gives suppress of spread of infection. The effect of vaccination was also considered in the recent paper [13]. Then, in this section, we examine the effect of infectivity difference under the coexistence of two variant viruses. Calculation conditions are shown in Table 2. Infectivity of α strain is set to constant and that of δ strain is changed from −8.2% to 4.1%. The results are displayed in Figure 3. In strong case, maximum infected people after SOE becomes about 1800, on the other hand, in weak case1 and weak case 2, maximum values become about 500 and 700 respectively. These results suggest that the ratio of infectivity between δ strain and α strain over 1.3 accelerates replacement of α by δ.

    Table 2.  Infectivity for sensitive analysis.
    α strain δ strain Difference (%) Ratio (δ/α)
    Nominal Case 0.09 0.121 0 1.35
    Weak Case 1 0.09 0.111 −8.2 1.23
    Weak Case 2 0.09 0.116 −4.1 1.29
    Strong Case 0.09 0.126 4.1 1.4

     | Show Table
    DownLoad: CSV
    Figure 3.  Change of infected people due to differences in infectivity.

    The problems of reinfection and breakthrough infection will be more important. To solve these problems, rate of decrease in antibody level or probability of reinfection and breakthrough infection must be taken into account. In the present calculation model has not yet adopted the methodology to solve above problems. Improvement of the model is future task.

    The above calculation results predict the peak of infection will become during the Olympic Games in Tokyo. Then, we thought three measures as follows.

    (1) In the case of continuing of the current measures (Case 1).

    (2) Measure 1: In the case of extending the SOE to 6/30 (Case 2).

    (3) Measure 2: To shorten the period from infection to isolation by PCR test or antigen test (Case 3).

    The infection status of each case is plotted together in Figure 4. The horizontal axis is the date from 2021/3/1. The broken line, dotted line and solid line indicate Case 1, Case 2, and Case 3 respectively. Solid line with blue dots shows the infection status in Tokyo on a 7-day average [15].

    Case 2 shows the transition of infected persons when the SOE is extended to 6/30. Peak of infected persons decreases by about 100, however the big improvement of the infection situation is not observed. Therefore, the extension until 6/30 has little effect on suppressing infection.

    Figure 4.  Infection status in case.

    The last measure is not effective unless as many people as possible participate. It is good if you know that you will be infected, but usually you do not know, so you need to have many people check it regularly. For that purpose, a negative certificate with a time limit (up to one week) should be issued and confirmed at restaurants and event venues. This will allow many people to be tested. In this way, it is possible to shorten the period from infection to isolation for two days or one day. In the present study, this period is set to 14 days due to the consistency of the data. We set the period to 12 days. Considering five days as a preparation period after the end of SOE, the implementation date was set to 6/25. Case 3 shows that the spread of infection after the peak of the 4th wave is suppressed to about 600. Searching for infected people early and shortening the infection period in this way is the most effective method other than vaccines as an infection control measure.

    We have extended the ATLM which has been developed to simulate the status of infection with various variant viruses. The developed model was applied to the 4th wave of Tokyo and the following results were obtained.

    (1) The fourth wave will bottom out near the end of the state of emergency, and the number of infected people will increase again.

    (2) The rate of increase will be mainly by δ strain, while the increase in α strain will be suppressed.

    (3) It is anticipated that the infection will spread during the Olympic Games.

    (4) When variant viruses compete, the infection of strongly infectious one rises sharply while the infection by weak infectious ones has converged.

    (5) The results of sensitivity analysis suggest that the ratio of infectivity between δ strain and α strain over 1.3 accelerates replacing speed of dominant virus in infection spread.

    (6) It is effective as an infection control measure to find an infected person early and shorten the period from infection to quarantine by PCR test or antigen test as a measure other than the vaccine.


    Acknowledgments



    MK is the former researcher of Hitachi Ltd., MU is the former professor of Tokyo Institute of Technology and SK is the former engineer of Hitachi Ltd. M. Koizumi developed the epidemiological model. M. Utamura verified the numerical results. S. Kirikami identified parameter values from data. All authors have read and agreed to the published version of manuscript.

    Data availability



    We used time-series data of COVID-19 for March 1 through June 10, 2021 in Tokyo [15].

    Conflict of interest



    The authors declare that they have no conflict of interest related to this report or the study it describes.

    [1] Kermack WO, Mckendrick AG (1927) A Contribution To The Mathematical Theory Of Epidemics. Proc Royal Soi London 115: 700-721.
    [2] Patel MD, Rosenstrom E, Ivy JS, et al. (2021) The Joint Impact of COVID-19 Vaccination and Non-Pharmaceutical Interventions on Infections, Hospitalizations, and Mortality: An Agent-Based Simulation. JAMA Netw Open 4: e2110782. doi: 10.1001/jamanetworkopen.2021.10782
    [3] Kuniya T, Inaba H (2020) Possible effects of mixed prevention strategy for COVID-19 epidemic: massive testing, quarantine and social distancing. AIMS Public Health 7: 490-503. doi: 10.3934/publichealth.2020040
    [4] Britton T, Ball F, Trapman P (2020) A mathematical model reveals the influence of population heterogeneity on herd immunity to SARS-CoV-2. Science 369: 846-849. doi: 10.1126/science.abc6810
    [5] Muñoz-Fernández GA, Seoane JM, Seoane-Sepúlveda JB (2021) A SIR-type model describing the successive waves of COVID-19. Chaos Solitons Fractals 144: 110682. doi: 10.1016/j.chaos.2021.110682
    [6] Biala TA, Khaliq AQM (2021) A fractional-order compartmental model for the spread of the COVID-19 pandemic. Commun Nonlinear Sci Numer Simul 98: 105764. doi: 10.1016/j.cnsns.2021.105764
    [7] Erik Cuevas (2020) An agent-based model to evaluate the COVID-19 transmission risks in facilities. Comput Biol Med 121: 103827.
    [8] Pageaud S, Ponthus N, Gauchon R, et al. (2021) Adapting French COVID-19 vaccination campaign duration to variant dissemination. medRxiv .
    [9] Chang E, Moselle KA (2021) Agent-Based Simulation of COVID-19 Vaccination Policies in CovidSIMVL. medRxiv .
    [10] Alagoz O, Sethi A, Patterson B, et al. (2021) The Impact of Vaccination to Control COVID-19 Burden in the United States: A Simulation Modeling Approach. medRxiv .
    [11] Truszkowska A, Behring B, Hasanyan J, et al. (2021) High-Resolution Agent-Based Modeling of COVID-19 Spreading in a Small Town. Adv Theor Simul 4: 2000277. doi: 10.1002/adts.202000277
    [12] Utamura M, Koizumi M, Kirikami S (2020) An Epidemiological Model Considering Isolation to Predict COVID-19 Trends in Tokyo, Japan: Numerical Analysis. JMIR Public Health Surveill 6: e23624. doi: 10.2196/23624
    [13] Utamura M, Koizumi M, Kirikami S (2021) A novel deterministic epidemic model considering mass vaccination and lockdown against COVID-19 spread in Israel: Numerical study. medRxiv .
    [14]  Characteristics and latest information of mutant viruses Available from: https://www3.nhk.or.jp/news/special/coronavirus/newvariant/.
    [15]  Updates on COVID-19 in Tokyo Available from: https://stopcovid19.metro.tokyo.lg.jp/en/.
  • publichealth-09-01-002-s001.pdf
  • This article has been cited by:

    1. Motoaki Utamura, Makoto Koizumi, Seiichi Kirikami, Novel deterministic epidemic model considering mass vaccination and lockdown against coronavirus disease 2019 spread in Israel: a numerical study, 2022, 7, 2396-8923, 10.1093/biomethods/bpac023
    2. Michio Murakami, Kenkichi Fujii, Wataru Naito, Masashi Kamo, Masaaki Kitajima, Tetsuo Yasutaka, Seiya Imoto, COVID-19 infection risk assessment and management at the Tokyo 2020 Olympic and Paralympic Games: A scoping review, 2023, 18760341, 10.1016/j.jiph.2023.03.025
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2465) PDF downloads(95) Cited by(2)

Figures and Tables

Figures(4)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog