key: cord-0960088-cbz5c2iu authors: Shah, N. H.; Suthar, A.; Satia, M. H.; Shah, Y.; Shukla, N.; Shukla, J.; Shukla, D. title: Modeling the Impact of Nationwide BCG Vaccine Recommendations on COVID-19 Transmission, Severity, and Mortality date: 2020-05-14 journal: nan DOI: 10.1101/2020.05.10.20097121 sha: 24555c3155ec0f920b66f47e590f87bc5eedb854 doc_id: 960088 cord_uid: cbz5c2iu COVID-19 was declared as pandemic on 11th March 2020 by WHO. There are apparent dissimilarities in incidence and mortality of COVID-19 cases in different parts of world. Developing countries in Asia and Africa with fragile health system have shown lower incidence and mortality compared to developed countries with superior health system in Europe and America. Most countries in Asia and Africa have national BCG vaccination program while Europe and America do not have such program, or have ceased it. At present, there is no known therapy to treat COVID-19 disease. There is no vaccine available currently to prevent COVID- 19 disease. As mathematical modelling is ideal for predicting the rate of disease transmission as well as evaluating efficacy of possible public health prevention measures, we have created a mathematical model with seven compartments to understand nationwide BCG vaccine recommendation on COVID-19 transmission, severity and mortality. We have computed two basic reproduction number, one at vaccine free equilibrium point and other at non-vaccine free equilibrium point and carried out local stability, sensitivity and numerical analysis. Our result showed that individuals with BCG vaccinations have lower risk of getting COVID-19 infection, shorter hospital stays, and increased rate of recovery. Furthermore, countries with long-standing universal BCG vaccination policies have reduced incidence, mortality, and severity of COVID-19. Further research will focus on exploring the immediate benefits of vaccination to healthcare workers and patients as well as benefits of BCG re-vaccination. In December of 2019, a novel strain of coronavirus was found in Wuhan, China and identified as Coronavirus Disease 2019 (COVID-19) [1] . Transmission of COVID-19, similar to numerous airborne respiratory viruses including tuberculosis and influenza, occurs through direct contact with an infected person through respiratory droplets when a person coughs or sneezes or indirect contact through interaction with contaminated surfaces with respiratory droplets from infected person [2] . As of April 18, 2020, the disease has spread to 210 countries with 2,324,731 total confirmed new cases reported and 160,434 deaths [3] . Since currently there is no definitive treatment or an effective vaccine for disease control, there is a crucial need to find measures that cure or reduce morbidity due to COVID-19. Figure 1 : Geographical distribution of COVID-19 cases worldwide as of 18 th April 2020 [4] . There are prominent dissimilarities in how COVID-19 is affecting different countries. Although developing countries like India, Philippines, Sri Lanka reported their initial case in January, they have yet to experience extensive community spread [5] . In contrast, developed countries like Italy and United Kingdom have experienced high mortality and widespread infection despite strong curtailing of social interactions. Studies show a disproportionately smaller number of cases have been reported from disadvantaged and low-income countries [6] . It is puzzling that countries with more fragile health systems report less incidence of COVID-19. Figure 1 shows a figure from the European Centre for Disease Prevention and Control which indicates that much of the incidence of COVID-19 is in the developing countries such as United Kingdom and United States [4] . A possible explanation for lower number of cases detected in these developing countries with extensive travel and trade links to China may stem from the induced heterologous protective activity of the Bacillus Calmette-Guerin (BCG) immunization [5] . The country used to recommend BCG vaccination for everyone, but currently does not. C: The country never had universal BCG vaccination programs [7] . Adapted from BCG World Atlas. Vaccines provide protection from particular pathogens by inducing effector mechanisms directed towards that pathogen. Certain live attenuated vaccines like the BCG vaccine increase immunity against not only a particular pathogen, but also numerous unrelated pathogens that present with similar acute respiratory tract infections [8] . Previous studies have established that BCG vaccination confers non-specific protection via the induction of innate immune memory [9] . BCG is used widely across the world as a vaccine for Tuberculosis, with many developing nations having a universal BCG vaccination policy for new-born [10] . Studies have indicated that countries with uniform BCG vaccination policy show significantly lower COVID-19 cases and deaths per million people compared to countries where BCG vaccination policy was ceased or was never in place [5] . Figure-2 from the BCG world atlas shows that numerous developing countries currently employ a universal BCG vaccination program, while many developed countries have either ceased their BCG vaccination recommendation or never recommended it [7] . Comparison of Figures 1 and 2 also indicates that the largest number of COVID-19 cases are in countries that currently have no BCG vaccination recommendation. Furthermore, the year that universal BCG vaccination was established within a country has a significant correlation with the mortality rate, indicating that earlier implementation of the policy leads to protection of a larger fraction of the elderly population. For instance, although Iran currently employs a universal BCG vaccination policy, it was implemented as recently as 1984. Consequently, it reports an elevated mortality with 19.7 deaths per million inhabitants [6] . In contrast, Japan started its universal BCG policy in 1947 and has around 100 times less deaths per million people, with 0.28 deaths. Similarly, Brazil started universal vaccination in 1920 and also has an even lower mortality rate of 0.0573 deaths per million inhabitants [6] . Exploring this relationship between BCG vaccination and COVID-19 further is crucial in establishing a potential short-or long-term prevention measure that is effective in reducing incidence and mortality of COVID-19. While previous studies have explored a possible link . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted May 14, 2020. . https://doi.org/10.1101/2020.05. 10.20097121 doi: medRxiv preprint between the BCG vaccine and incidence of COVID-19, they have not modelled the impact of implementing a nationwide BCG immunization policy [6, 7] . Mathematical modelling is ideal for predicting the rate of disease transmission as well as evaluating efficacy of possible public health prevention measures [11] . This study aims to create a model to evaluate the effect of nationwide BCG vaccine recommendations on incidence, severity, and mortality rates of COVID-19. Furthermore, our study aims to model the protentional implications of an immediate implementation of the universal BCG vaccination policy. A novel corona virus pandemic is occurring around the world, both in countries with and without the BCG vaccine. To analyse the importance of this vaccine against the fight of COVID-19, we have developed a mathematical model. This model has seven various compartments. First compartment is the class of exposed individuals to COVID-19 noted as From these individuals, some are vaccinated with BCG vaccine noted as and those who have not taken it are noted as . Now, with the certain rate these individuals transfer into Infected stage noted as and some can become critically infected noted as compartment. The individuals suffering from COVID-19 that need to be hospitalised are noted as . The infected and hospitalised individual that recover are noted as . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 14, 2020. then we get new system as follow , , , , , , : CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 14, 2020. (2) The dynamical behaviour of system (1) is equivalent to the system (2). Hence, every solution of system (2) will remain in the region . The above system (2) has three equilibrium point by setting the equation as zero. 1. Diseases-free equilibrium point where 2. Vaccine free equilibrium point where 3. Non-vaccine free equilibrium point , ,0, , , , CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 14, 2020. The diseases-free equilibrium point is not valid according to current scenario and therefore we will discuss only two equilibrium points. In this section, basic reproduction number is derived using next generation matrix method [12] . . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 14, 2020. . https://doi.org/10.1101/2020.05.10.20097121 doi: medRxiv preprint These Jacobian matrices are used to derive whose spectral radius is our desired basic reproduction number . Here, we have computed two basic reproduction numbers; one at vaccine free equilibrium point and other at non-vaccine free equilibrium point. The expression of basic reproduction number for non-vaccinated individuals i.e. at vaccinefree equilibrium point is denoted by given by Whereas, the expression of basic reproduction number for vaccinated individuals i.e. at nonvaccine-free equilibrium point is denoted by given by Moreover, the Basic reproduction for the system is . In this section, we will be discussing the local stability for vaccine and non-vaccine free equilibrium point utilizing the concept of eigenvalues. To find the eigen values, the Jacobian matrix of the system (2) is calculated as CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 14, 2020. Here, sensitivity analysis is carried out to see which parameter is conferring a positive effect on the model. The sensitivity analysis is conducted on the basic reproduction number. It is to be done using Christoffel formula which is where is the model parameter. CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 14, 2020. . https://doi.org/10.1101/2020.05.10.20097121 doi: medRxiv preprint In this section, results of transmission of COVID-19 outbreak is simulated numerically, helping validate model results. Figure 4 shows the change in the behaviour of individuals. This suggests that exposed individuals quickly progress into infected ones. Infected individuals increase in first 1.5 weeks and then start to decrease whereas critically infected individuals stay in lower numbers. Hospitalisation increases rapidly for at least 3 weeks and the it flattens. Recovery rate increases but at lower rate. . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 14, 2020. . https://doi.org/10.1101/2020.05. 10.20097121 doi: medRxiv preprint To get precise results, four clusters of figure 5 is magnified. Our data indicates that individuals should go into isolation earlier to prevent exposure to COVID-19. BCG vaccinated individual may get exposed in 0.2 week while, at the same time, non-vaccinated individual may get critically infected. Moreover, recovered individual also may become exposed and get reinfected after 2 weeks. Figure 6 : Cyclic behaviour of recovered individuals Figure 6 indicates the cyclic behaviour of recovered individuals. The first graph shows vaccinated individuals while the other displays non-vaccinated individuals. The graphs show that the cycle of vaccinated individuals is larger than non-vaccinated individuals. Furthermore, vaccinated individuals have greater time to reinfection when compared to that of nonvaccinated individuals. . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 14, 2020. . https://doi.org/10.1101/2020.05. 10.20097121 doi: medRxiv preprint Figure 7 : Behaviour of exposed individuals Figure 7 compares the behaviour of exposed individuals with vaccinated individuals, nonvaccinated individuals, and infected individuals, respectively. The rate of exposure is consecutive for non-vaccinated individuals when compared to vaccinated individuals. Furthermore, the transmission from exposed to infected is quite quick i.e. almost every individual who is exposed will get infected. is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 14, 2020. . https://doi.org/10.1101/2020.05. 10.20097121 doi: medRxiv preprint Figure 9 shows the behaviour of infected individuals towards critical, hospitalised and recovered individuals, respectively. If individuals are vaccinated early, the chances of becoming critically infection are low. However, this also indicates that infected individuals must receive treatment in hospitals before the infection progresses to critical stage. This is further elaborated by the last graph, which indicates that the chances of recovering by staying at home are negligible. Figure 10 : Behaviour of hospitalised individuals Figure 10 shows the behaviour of hospitalised individuals and their recovery. This figure further emphasizes that individuals recover only after proper hospitalization is received. Our model indicates that a universal recommendation of BCG vaccination is beneficial in reducing the incidence rate of COVID-19. Moreover, the vaccination reduces the severity of the disease by reducing the hospitalization time for infected individuals. It also contributed to a higher rate of recovery for vaccinated individuals when compared to non-vaccinated individuals. The risk of reinfection is also reduced in vaccinated populations. Finally, the overall disease specific mortality of infected individuals is also reduced in vaccinated populations. Consequently, we validate the effectiveness of universal BCG vaccination policy in increasing herd immunity against COVID-19. Considering that the BCG vaccination is established in preventing acute respiratory tract infections even in the elderly [13] , vulnerable populations may be immunized with BCG vaccines until a COIVD-19 specific vaccine is developed. This strategy would be especially beneficial in the new-borns, elderly, and frontline essential workers. This study validates the call for clinical trials of BCG vaccine against . Given that our model shows a significant reduction in incidence of COVID-19 in countries with universal BCG vaccination, we recommend all countries to conduct more clinical trials to establish and consider policy regarding BCG vaccinations for its citizens. Our study may be subject to a few limitations. The study does not take into account confounding variables such as the potential differences in geographical and biological factors such as temperature, humidity, life expectancy, average income, social-cultural norms, ethnical genetic background, and mitigation between developing and developed countries. Due to these confounding variables, we cannot establish causality. The epidemiological data comparing vaccinated countries and non-vaccinated countries must be monitored throughout this . CC-BY-NC-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted May 14, 2020. . https://doi.org/10.1101/2020.05. 10.20097121 doi: medRxiv preprint pandemic to establish a causal impact in reducing incidence and mortality. Overall, our study serves to establish that individuals with BCG vaccinations have lower risk of getting COVID-19 infection, shorter hospital stays, and increased rate of recovery. Furthermore, countries with long-standing universal BCG vaccination policies have reduced incidence, mortality, and severity of COVID-19. Further research will focus on exploring the immediate benefits of vaccination to healthcare workers and patients as well as benefits of BCG re-vaccination. The results of this study will serve to inform and encourage further research and, ultimately, the creation of policy regarding universal BCG vaccination. The authors declare that we have no competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Our conclusion is based on the data derived from our above mathematical calculations which are available to readers. New SARS-like virus in China triggers alarm First case of 2019 novel coronavirus in the United States World Population Clock: Worldometer Situation update worldwide, as of Is Global BCG Vaccination Coverage Relevant To The Progression Of SARS-CoV-2 Pandemic? Correlation between universal BCG vaccination policy and reduced morbidity and mortality for COVID-19: an epidemiological study Aaron The BCG world atlas: A database of global BCG vaccination policies and practices Long-lasting effects of bcg vaccination on both heterologous th1/th17 responses and innate trained immunity Trained immunity: A memory for innate host defense Non-specific effects of BCG vaccine on viral infections Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study The construction of next-generation matrices for compartmental epidemic models BCG Vaccination Protects against Experimental Viral Infection in Humans through the Induction of