Next Article in Journal
What Happened to People with Non-Communicable Diseases during COVID-19: Implications of H-EDRM Policies
Next Article in Special Issue
A Look Behind the Scenes at COVID-19: National Strategies of Infection Control and Their Impact on Mortality
Previous Article in Journal
Burnout Syndrome in Police Officers and Its Relationship with Physical and Leisure Activities
Previous Article in Special Issue
Experiential Learning Program to Strengthen Self-Reflection and Critical Thinking in Freshmen Nursing Students during COVID-19: A Quasi-Experimental Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

BCG Vaccination and Mortality of COVID-19 across 173 Countries: An Ecological Study

1
Division of Molecular Epidemiology, The Jikei University School of Medicine, Tokyo 105-8461, Japan
2
Advanced Therapies Innovation Department, Siemens Healthcare K.K., Tokyo 141-8644, Japan
3
Hitachi, Ltd. Research & Development Group, Tokyo 185-8601, Japan
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(15), 5589; https://doi.org/10.3390/ijerph17155589
Submission received: 4 July 2020 / Revised: 24 July 2020 / Accepted: 31 July 2020 / Published: 3 August 2020
(This article belongs to the Special Issue COVID-19: Prevention, Diagnosis, Therapy and Follow Up)

Abstract

:
Ecological studies have suggested fewer COVID-19 morbidities and mortalities in Bacillus Calmette–Guérin (BCG)-vaccinated countries than BCG-non-vaccinated countries. However, these studies obtained data during the early phase of the pandemic and did not adjust for potential confounders, including PCR-test numbers per population (PCR-tests). Currently—more than four months after declaration of the pandemic—the BCG-hypothesis needs reexamining. An ecological study was conducted by obtaining data of 61 factors in 173 countries, including BCG vaccine coverage (%), using morbidity and mortality as outcomes, obtained from open resources. ‘Urban population (%)’ and ‘insufficient physical activity (%)’ in each country was positively associated with morbidity, but not mortality, after adjustment for PCR-tests. On the other hand, recent BCG vaccine coverage (%) was negatively associated with mortality, but not morbidity, even with adjustment for percentage of the population ≥ 60 years of age, morbidity, PCR-tests and other factors. The results of this study generated a hypothesis that a national BCG vaccination program seems to be associated with reduced mortality of COVID-19, although this needs to be further examined and proved by randomized clinical trials.

1. Introduction

Currently, more than four months since declaration of the coronavirus disease 2019 (COVID-19) as a pandemic by the World Health Organization (WHO) on March 11th, 2020, more than 14 million people have been infected with the virus and more than half a million have died worldwide. Marked differences in COVID-19 mortalities have been observed in different countries. For example, the mortality per million population is till now several tens of times or even higher in Western countries, e.g., Belgium (845), the United Kingdom (UK, 664), the United States of America (USA, 426) and Germany (109), than in Asian countries, e.g., India (19), Japan (8) and China (3), as of 17 July 2020. This is quite the opposite of what was reported during the 1918–20 influenza pandemic, the so called Spanish flu, in which the population mortality was over 30-fold higher, with excess death rates observed in low-income countries, such as India, than in high-income countries, such as those in the West [1]. Higher morbidities and mortalities may be explained by easy access to diagnostic polymerase chain reaction tests (PCR-tests) for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in these Western countries. In contrast, they may be underestimated in middle- and low-income countries due to low capacities of PCR-testing and poor access to health care. Alternatively, there is growing concern that racial and ethnic minority, as well as socioeconomic and biological factors may influence the high morbidity and mortality. A retrospective study of integrated health care systems suggested that Black race compared with White race, increasing age, male sex, having a greater burden of comorbidities, type of public insurance, residence in a low-income area and obesity were associated with increased odds of hospital admission [2]. On the other hand risk of in-hospital mortality increased only in aged patients and was not associated with the Black race, sex, comorbidities, obesity and other factors after multivariate adjustment [2]; this phenomenon was also confirmed in other races, i.e., Asians and Hispanics, compared with the White race [3,4]. Moreover, since COVID-19 is an infectious disease that spreads mainly through the droplet route by close contact in dense human societies, metropolitan areas, such as New York City in the USA [5] and Lombardy in Italy [6], Paris in France [7], Sao Paulo in Brazil [8], and so on, have tended to be regional epicenters. However, associations between population dynamics, e.g., population size, density, migrants and urbanization and the morbidity/mortality of COVID-19 have not yet been well examined.
In addition to the factors listed above, several scientists proposed a hypothesis [9,10,11,12] that Bacillus Calmette–Guérin (BCG) vaccination has preventive effects not only against tuberculosis, i.e., the target disease of the vaccine, but also other non-specific infectious diseases, i.e., off-target diseases such as COVID-19, through epigenetic mechanisms [13,14,15]. In fact, ecological studies have suggested that both COVID-19 morbidities and mortalities were less in BCG-vaccinated countries than in BCG-non-vaccinated countries [16,17,18]. However, because these studies only obtained data during the early phase of the pandemic, by which time the disease load had escalated in Western countries, but not yet escalated in low- and middle-income countries where BCG is given at birth and did not adjust for any potential confounders, including PCR-test number per million population, the BCG-hypothesis needs to be reexamined now, more than four months after declaration of the pandemic, since the number of PCR-confirmed COVID-19 cases are still growing in many countries. We therefore aimed to explore whether recent BCG vaccine coverage is associated with COVID-19 morbidity and/or mortality rates, using linear regression models to explore associations between the two continuous random variables adjusted for a variety of potential confounders, such as median age and body mass index (BMI) in individual countries through this ecological study.

2. Materials and Methods

2.1. Ethics Statement

Institutional review board approval for this work was not sought due to the use of publicly available data obtained from open resources.

2.2. Study Design

This ecological study compared population data of each country.

2.3. Data Resources

All data were obtained from open resources: Information on total number of cases, total number of deaths and number of PCR-tests performed per million population were obtained from ‘Coronavirus Update’ [19], data regarding population dynamics were obtained from ‘World Population Clock’ [20], socioeconomic covariates were obtained from ‘the United Nations database’ [21], BCG vaccine coverage data were obtained from ‘The BCG World Atlas’ [22] and ‘The global health observatory’ on the WHO homepage [18] and other health related data were from the same WHO site [23]. Definitions of each of the covariates are described in the Supplement. Only countries that had data of both total deaths and BCG vaccine coverage were included for analyses in this study.

2.4. Outcomes and PCR-Tests

The outcomes evaluated in this study were COVID-19 morbidity, i.e., total COVID-19 cases per million population and COVID-19 mortality, i.e., total COVID-19-related deaths per million population, in each country, obtained from ‘Worldometer COVID-19 Data’ on July 17th, 2020 [19]. Data on the total number of PCR-tests performed per million population were simultaneously obtained from the same ‘Worldometer COVID-19 Data’ website. PCR-test positivity was simply calculated as the total COVID-19 case number divided by the total number of PCR-tests performed in each country.

2.5. BCG Vaccine Coverage

Recent BCG vaccine coverage was defined as the percentage of the vaccinated population among one-year-olds in each country (World Health Data Platform, the World Health Observatory, BCG-Immunization coverage estimates) [23]. For countries that have already stopped a national BCG vaccine immunization program [22], which includes a significant number of countries, their BCG vaccine coverage rate was counted as zero percent in this study.

2.6. Sixty-One Covariates

The covariates evaluated in this study, the definitions of which are described in the Supplement (Table S1), included: (1) Population (n); (2) yearly change in population (%); (3) net change in population (n); (4) population density (n/km2); (5) land area (km2); (6) net number of migrants (n); (7) fertility rate (n); (8) median age (years); (9) urban population percentage (%); (10) world share (%); (11) population between age 0 to 14 years (%); (12) population ≥ 60 years of age (%); (13) population ≥ 70 years of age (%); (14) gross domestic product (GDP) (million US dollars); (15) GDP per capita (US dollars); (16) total unemployment rate (%); (17) male unemployment rate (%); (18) female unemployment rate (%); (19) total labor force participation rate (%); (20) male labor force participation rate (%); (21) female labor force participation rate (%); (22) annual incidence of tuberculosis per 100,000 population (n); (23) international health regulation (IHR) score (%); (24) universal health coverage (UHC) index; (25) hospital beds (n) per 10,000 population; (26) medical doctors (n) per 10,000 population; (27) nursing midwifery (n) per 10,000 population; (28) licensed qualified anesthesiologists who usually cover management of intensive care units actively working (n) per 10,000 population; (29) total expenditure on health as a percentage of GDP (%); (30) population with household expenditures on health greater than 10% of total household expenditure or income (%); (31) prevalence of high blood pressure (systolic ≥ 140 or diastolic ≥ 90 mmHg) (%); (32) prevalence of elevated fasting blood glucose levels (≥7.0 mmol/L or on anti-diabetic medication); (33) prevalence of elevated total cholesterol levels (≥5.0 mmol/L) (%); (34) mean BMI (body weight [kg]/height2 [m2]); (35) prevalence of obesity among adults, BMI ≥ 30 kg/m2 (%); (36) prevalence of ‘overweight’ people among adults, BMI ≥ 25 kg/m2 (%); (37) alcohol drinking, total per capita (≥15 years of age) consumption (in liters of pure alcohol over a calendar year); (38) prevalence of smoking any tobacco product among males aged ≥ 15 years (%); (39) prevalence of smoking any tobacco product among females aged ≥ 15 years (%); (40) prevalence of insufficient physical activity among adults aged ≥ 18 years (%); (41) estimated population-based prevalence of depression (%); (42) neonatal mortality rate (n per 1000 live births); (43) infantile mortality rate (n per 1000 live births); (44) under-five mortality rate (probability of dying by the age of 5 years per 1000 live births); (45) mortality rate for 5–14-year-olds (probability of dying per 1000 children aged 5–14 years); (46) adult mortality rate (probability of dying between 15 and 60 years of age per 1000 population); (47) probability of dying between age 30 and exact age 70 years from any of the following causes: cardiovascular disease, cancer, diabetes or chronic respiratory disease; (48) life expectancy at birth (years); (49) life expectancy at age 60 years (years); (50) healthy life expectancy (HALE) at birth (years); (51) HALE at age 60 years (years); (52) death due to chronic obstructive pulmonary disease (%); (53) death due to ischemic heart disease (%); (54) death due to lower respiratory infections (%); (55) death due to stroke (%); (56) death due to tracheal, bronchial and lung cancers (%); (57) total of (52) to (56) as ambient and household air pollution-attributable death rate (n per 100,000 population); and (58) annual mean concentration of particulate matter less than 2.5 microns in diameter (PM2.5) [µg/m3] in urban areas; and (59) coverage rate with the first dose of a measles-containing-vaccine (MCV1) among one-year-olds (%) as well as (60) recent BCG coverage and (61) PCR-tests number.

2.7. Statistics

Linear Regression Models

For preprocessed data, outcomes, i.e., morbidity and mortality per million population were transformed to the common logarithm (log10) to adjust for normality of the distribution, which was verified by means of kurtosis tests. When the number of total deaths was zero, these were changed to 0.01 per million population, because zero cannot be transformed to the common logarithm. Variance inflation factor (VIF) was used to detect the presence of multicollinearity. Only one variable among biologically similar variables, e.g., ‘median age (years)’, ‘≥ 60 years of age (%)’ and ‘≥ 70 years of age’, was selected to maximize adjusted R2 in the multi-linear regression models to avoid the influence of collinearity. If the variance inflation factor (VIF) of certain covariates was more than 10, then the covariates were avoided in multivariate analyses because of a collinearity issue. Multiple linear regression models were used to screen potential risk or preventive factors associated with morbidity by adjusting for PCR-test numbers transformed to the common logarithm (log10) and those associated with mortality were screened by adjusting for PCR-test numbers and morbidity per million population. Considering type I error due to a multiple testing, the significance level of alfa was set as p < 0.001. Then, all the screened factors were assessed in a multi-linear regression model to determine significant factors with p < 0.05 as the cutoff point. Each model was evaluated by adjusted R2 as a coefficient of determination. Pearson’s correlation coefficient for variables with normal distributions or Spearman’s rank correlation for variables with non-normal distributions, represented as rho, was used to quantify the strengths of associations between morbidity, mortality and significant factors determined by the final models, as: absolute value of rho ≥ 0.5: very strong; rho ≥ 0.4: strong; 0.4 > rho ≥ 0.2: moderate; and rho < 0.2: weak associations. Data were analyzed using Stata version 14.0 software (StataCorp LP, College Station, TX, USA).

3. Results

3.1. Variability of COVID-19 Morbidity and Mortality across 173 Countries

A total of 173 countries that had data of both total COVID-19 deaths and BCG vaccine coverage were included for analyses in this study. The 20 countries with highest and lowest COVID-19 morbidities (Table 1) and mortalities (Table 2), as well as their PCR-test rate per million population are shown below. Marked differences in morbidities and mortalities were observed among these countries, ranging from 1 (Papua New Guinea) to 37,566 (Qatar) and from 0 (Vietnam, etc.) to 845 (Belgium), respectively. Six and thirteen countries that do not have a BCG national vaccine program at present (indicated with bold and mark “∫”) were included in the 20 countries with the highest COVID-19 morbidity and mortality, respectively.
Histograms of morbidities and mortalities were drawn as normal density plots (Figure 1). Although the histograms of morbidity and mortality were skewed to the right, they followed a normal distribution by transformation with the common logarithm (log10).

3.2. Associations between Morbidity, Mortality and PCR-Tests per Million Population

First, associations represented by rho and VIF between morbidity (n = 173), mortality (n = 173) and PCR-tests (n = 155) are shown (Figure 2A). These three variables were predicted to have very strong and positive associations with each other. However, VIFs were less than 2 among these three factors. Multicollinearity is considered to be present when the VIF is higher than 5 to 10 [24]. Thus, any variable with a VIF < 5.0 was considered for inclusion in multiple linear regression analyses. Considering the number of PCR-tests per million population may exhibit associations with morbidity, adjustment was performed for the PCR-test number in every analysis when screening for the risk factors of COVID-19 morbidity per million population (Figure 2B). Considering morbidities and the number of PCR-tests per million population may exhibit associations with mortality, adjustment was performed for the morbidity and PCR-test number in every analysis when screening for the risk factors of COVID-19 mortality per million population (Figure 2C).
Since fewer PCR-tests may underestimate morbidity and mortality, the association between mortality as the outcome and morbidity as the exposure was adjusted for number of PCR-tests performed (Table 3). In this multiple regression analysis, higher morbidity was associated with higher mortality, whereas more PCR-tests were associated with lower mortality.
Evaluation of the association between PCR-test positivity and mortality, shown as a scatter plot, indicated a very strong association between them (rho = 0.54) (Figure 3). Countries with higher PCR-positivity rates tended to have higher mortality rates. PCR-test positivity rates of countries where no deaths due to COVID-19 were observed were less than 3.5%. Minimum positivity and no deaths were observed in Vietnam.

3.3. Screening Factors Associated with Morbidity

Among the 59 covariates, plus BCG vaccine coverage and PCR-testing rate, i.e., a total of 61 factors, ‘urban population’ and ‘insufficient physical activity’ were significantly (p < 0.001) associated with morbidity after adjustment for PCR-test rate (Table 4). Next, these two significant factors were used in a multi-linear regression model to eliminate confounding (Table 5). As a result, ‘urban population’ (p = 0.02) and ‘insufficient physical activity’ (p = 0.01) remained significant factors associated with COVID-19 morbidity, even after adjustment for PCR-tests. The adjusted R2 was 0.5037.
The association between ’urban population’ and morbidity, demonstrated below as a scatter plot, showed a very strong association (rho = 0.55) (Figure 4).
The association between ‘insufficient physical activity’ and morbidity, demonstrated below as a scatter plot, also showed a very strong association (rho = 0.52) (Figure 5).
COVID-19-related morbidity rates per million population on July 17th were transformed to the common logarithm (log10) in the graph. Countries that had never had or that had stopped a national program of BCG vaccination are indicated in red, while countries that currently follow a national BCG vaccine program are indicated in black. Selected country’s name was shown using three-letter country codes.

3.4. Screening Factors Associated with Mortality

Among the 58 covariates evaluated, plus BCG vaccine coverage, adjusted for morbidity and PCR-tests, age-related factors, i.e., median age, ≥60 years of age (%) and ≥70 years of age (%), were significantly (p < 0.001) associated with mortality (Table 4). Since these three age-related factors had collinearity for mortality, ‘≥ 60 years of age’ was selected to maximize adjusted R2 of the multi-linear regression models. Moreover, ‘BCG vaccine coverage’, ‘Elevated total cholesterol levels’ and ‘Life expectancy at 60 years of age’, were also significant (p < 0.001) factors associated with mortality. Next, these significant factors were used in a multi-linear regression model to eliminate confounding (Table 6). As a result, ‘≥60 years of age’ (p < 0.001) and ‘BCG vaccine coverage’ (p = 0.002) remained significant factors associated with COVID-19 mortality, even after adjustment for morbidity and PCR-tests. The adjusted R2 was 0.8254.
Evaluation of the association between ’median age’ and mortality, demonstrated as a scatter plot, showed a very strong association between these two variables (rho = 0.54) (Figure 6). Countries with a larger population ≥ 60 years of age (%) showed a tendency toward a higher mortality rate.
Finally, evaluation of the association between ‘BCG vaccine coverage’ and mortality, demonstrated as a scatter plot, indicated a moderately negative association (rho = −0.29) (Figure 7). Countries with higher BCG vaccine coverage showed a tendency toward lower mortality. Additionally, COVID-19 mortality rates did not have significant associations with BCG strain, e.g., Tokyo 172. Moreover, there were no significant associations between mortality and either the year of stopping or introducing a national BCG vaccine program (data not shown).

4. Discussion

Among the variety of parameters abstracted from open resources, ‘BCG vaccine coverage’ had a significant association with COVID-19 mortality, even after adjusting for morbidity, PCR-tests, age, universal health coverage, numbers of medical doctors, elevated total cholesterol and healthy life expectancy. On the other hand, BCG vaccination was not associated with COVID-19 morbidity. The main results of this study are consistent with a very recent article demonstrating that every 10% increase in the BCG index was associated with a 10% reduction in COVID-19 mortality [25]. Moreover, a retrospective cohort study suggested that BCG-vaccinated individuals were less likely to require hospital admission during the disease course [26]. In contrast to BCG, coverage of the measles vaccine, which is also considered to induce heterologous protection against infections through long-term boosting of innate immune responses [9], showed no association with the morbidity and mortality of COVID-19, which was also consistent with a very recent article showing a significant low risk of COVID-19 mortality in countries with higher BCG vaccine coverage, but not with measles vaccine coverage [27].
Moreover, SARS-CoV-2 is a single-stranded positive-sense RNA virus and the BCG vaccine has been shown in controlled trials to reduce the severity of infections by other viruses with such a structure [9]. For example, the BCG vaccine reduced yellow fever vaccine viremia by 71% in volunteers in the Netherlands [28]. However, some countries with a current national BCG vaccination policy have high mortality rates. Plausible reasons for this discrepancy may be: (1) low coverage of BCG vaccination in these countries (% coverage–mortality per million population), e.g., Ireland (18%—354); Portugal (32%—165); and Greece (50%—19): (2) late introduction of BCG vaccine program (year of introduction–mortality per million population), e.g., Iran, (1984–162): (3) oral delivery of the BCG vaccine, e.g., Brazil retained oral delivery of the vaccine until 1977: (4) UK administered the vaccine to older children (12 to 13 years of age) and (5) connection with endemic country by land, e.g., Mexico, Panama, Peru and Chile.
Higher morbidity, but fewer PCR-tests, were associated with higher mortality. Moreover, higher PCR-test positivity was associated with higher mortality, which was consistent with the report by Hisaka et al. [29]. Expanding application of PCR-tests not only to typical symptomatic cases, but also to mild or asymptomatic cases and to those who had close contact with patients, may decrease the PCR-test positivity rate. Thus, enhancing the capacity of PCR-testing may enable identification of cases, so that appropriate measures can be taken to prevent them spreading SARS-CoV-2 to others at home, at their work place or at events of mass gatherings.
In this study, the covariate of ‘urban population’ and ‘insufficient physical inactivity’ had a strong and positive association with morbidity, but not with mortality. People living in urban areas tend to have close contact with a greater number of people per day than those in rural areas, independent of age. A report on 4103 patients with COVID-19 in New York City found that obesity, which may strongly depend on the balance between physical activity and diet, was one of the clinical features leading to hospital admission [30]. On the other hand, older age was associated with higher mortality, which is consistent with previous articles [2,3,4]. From this study, the risk factors for morbidity seem to be different from those associated with mortality, suggesting that factors related to susceptibility may be different from those related to disease severity.
There are several limitations to this study. First, although we selected 61 covariates in this study, we did not evaluate range and timing of non-pharmaceutical interventions, e.g., school closures, workplace closures, cancellation of public events, restrictions on public gatherings, stay-at-home restrictions, restrictions on internal movement, international travel controls, etc., all of which would also have had significant effects on COVID-19-related morbidity and mortality. Therefore, the present study results are burdened with an extreme error. Second, the study design was ecological. Therefore, the outcome of this work should be considered highly limited, with a potential risk of high bias. Consequently, only the hypothesis that BCG vaccination mitigates COVID-19 mortality can be proposed here; cohort or case–control studies and randomized clinical trials, similar to the BCG–CORONA study [31], are required to test this hypothesis. Third, the COVID-19 pandemic is still ongoing, although we have confirmed the results using the latest data. However, the results may be different a few months from now. Fourth, it is clear that an extremely large number of covariates, i.e., 61, were selected for the limited number of 173 countries. This sample size could allow use of a maximum of 10 covariates in multiple regression analysis. Therefore, significant (p < 0.001) variables were at first screened after adjustment for PCR-tests and morbidity (Table 4). Then, multivariable linear regression using the screened variables were performed after adjustment for PCR-tests and morbidity (Table 5 for morbidity and Table 6 for mortality). Fifth, mortality in each country should be compared with excess deaths whether these two do not make a big difference. Sixth, the definition of COVID-19 cases may differ by country, e.g., including PCR confirmed, but asymptomatic cases and pneumonia cases with negative PCR-tests. Seventh, there were still residual confounders even after adjustment for PCR-tests.

5. Conclusions

Our results suggest the hypothesis that greater BCG vaccine coverage may reduce the risk of deaths due to COVID-19, which needs to be further studied by observational studies and confirmed by randomized clinical trials.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/17/15/5589/s1, Table S1: Definitions or descriptions of all covariates.

Author Contributions

Conceptualization, M.U.; methodology, M.U., K.O., Y.H.; validation, T.A., M.U.; formal analysis, M.U., T.A.; resources, T.A.; data curation, T.A.; writing—original draft preparation, M.U.; writing—review and editing, K.O., Y.H., T.A.; visualization, M.U. All authors have read and agreed to the published version of the manuscript. These four authors contributed equally to this work.

Funding

This research received no external funding.

Acknowledgments

The authors thank the personnel who created the excellent open resources. We also thank Haruka Wada for data management at the Division of Molecular Epidemiology, Jikei University School of Medicine.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Murray, C.J.; Lopez, A.D.; Chin, B.; Feehan, D.; Hill, K.H. Estimation of potential global pandemic influenza mortality on the basis of vital registry data from the 1918-20 pandemic: A quantitative analysis. Lancet 2006, 368, 2211–2218. [Google Scholar] [CrossRef]
  2. Price-Haywood, E.G.; Burton, J.; Fort, D.; Seoane, L. Hospitalization and mortality among black patients and white patients with Covid-19. N. Engl. J. Med. 2020, 382, 2534–2543. [Google Scholar] [CrossRef] [PubMed]
  3. Patel, A.P.; Paranjpe, M.D.; Kathiresan, N.P.; Rivas, M.A.; Khera, A.V. Race, Socioeconomic deprivation, and hospitalization for COVID-19 in English participants of a National Biobank. medRxiv 2020. [Google Scholar] [CrossRef]
  4. Rentsch, C.T.; Kidwai-Khan, F.; Tate, J.P.; Park, L.S.; King, J.T.; Skanderson, M.; Hauser, R.G.; Schultze, A.; Jarvis, C.I.; Holodniy, M.; et al. Covid-19 by race and ethnicity: A national cohort study of 6 million United States veterans. medRxiv 2020. [Google Scholar] [CrossRef]
  5. Richardson, S.; Hirsch, J.S.; Narasimhan, M.; Crawford, J.M.; McGinn, T.; Davidson, K.W. And the Northwell COVID-19 Research Consortium. Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized with COVID-19 in the New York City area. JAMA 2020, 10, 2052–2059. [Google Scholar] [CrossRef]
  6. Grasselli, G.; Pesenti, A.; Cecconi, M. Critical Care utilization for the COVID-19 outbreak in Lombardy, Italy: Early experience and forecast during an emergency response. JAMA 2020, 323, 1545–1546. [Google Scholar] [CrossRef] [Green Version]
  7. Lapostolle, F.; Goix, L.; Vianu, I.; Chanzy, E.; De Stefano, C.; Gorlicki, J.; Petrovic, T.; Adnet, F. COVID-19 epidemic in the Seine-Saint-Denis Department of Greater Paris: One month and three waves for a tsunami. Eur. J. Emerg. Med. 2020, 27, 274–278. [Google Scholar] [CrossRef]
  8. Pereira, I.G.; Guerin, J.M.; Silva Júnior, A.G.; Garcia, G.S.; Piscitelli, P.; Miani, A.; Distante, C.; Gonçalves, L.M.G. Forecasting Covid-19 dynamics in Brazil: A data driven approach. Int. J. Environ. Res. Public Health. 2020, 17, 5115. [Google Scholar] [CrossRef]
  9. Netea, M.G.; Joosten, L.A.; Latz, E.; Mills, K.H.; Natoli, G.; Stunnenberg, H.G.; O’Neill, L.A.; Xavier, R.J. Trained immunity: A program of innate immune memory in health and disease. Science 2016, 352, 6284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Sharma, A.R.; Batra, G.; Kumar, M.; Mishra, A.; Singla, R.; Singh, A.; Singh, R.S.; Medhi, B. BCG as a game-changer to prevent the infection and severity of COVID-19 pandemic? Allergol. Immunopathol. 2020. [Google Scholar] [CrossRef] [PubMed]
  11. Koti, M.; Morales, A.; Graham, C.H.; Siemens, D.R. BCG vaccine and COVID-19: Implications for infection prophylaxis and cancer immunotherapy. J. Immunother. Cancer 2020, 8, e001119. [Google Scholar] [CrossRef] [PubMed]
  12. Ebina-Shibuya, R.; Horita, N.; Namkoong, H.; Kaneko, T. National policies for paediatric universal BCG vaccination were associated with decreased mortality due to COVID-19. Respirology 2020. [Google Scholar] [CrossRef] [PubMed]
  13. Curtis, N.; Sparrow, A.; Ghebreyesus, T.A.; Netea, M.G. Considering BCG Vaccination to reduce the impact of COVID-19. Lancet 2020, 395, 1545–1546. [Google Scholar] [CrossRef]
  14. Netea, M.G.; Giamarellos-Bourboulis, E.J.; Domínguez-Andrés, J.; Curtis, N.; van Crevel, R.; van de Veerdonk, F.L.; Bonten, M. Trained immunity: A tool for reducing susceptibility to and the severity of SARS-CoV-2 infection. Cell 2020, 181, 969–977. [Google Scholar] [CrossRef]
  15. O’Neill, L.A.J.; Netea, M.G. BCG-induced trained immunity: Can it offer protection against COVID-19? Nat. Rev. Immunol. 2020, 20, 335–337. [Google Scholar] [CrossRef]
  16. Ozdemir, C.; Kucuksezer, U.C.; Tamay, Z.U. Is BCG vaccination affecting the spread and severity of COVID-19? Allergy 2020. [Google Scholar] [CrossRef]
  17. Gursel, M.; Gursel, I. Is global BCG vaccination-induced trained immunity relevant to the progression of SARS-CoV-2 pandemic? Allergy 2020, 75, 1815–1819. [Google Scholar] [CrossRef]
  18. Miyasaka, M. Is BCG vaccination causally related to reduced COVID-19 mortality? EMBO Mol. Med. 2020, 12, e12661. [Google Scholar] [CrossRef]
  19. Coronavirus Update. Available online: https://www.worldometers.info/coronavirus/ (accessed on 17 July 2020).
  20. World Population Clock. Available online: https://www.worldometers.info/world-population/ (accessed on 22 May 2020).
  21. UN Data. A World of Information. Available online: https://data.un.org/ (accessed on 17 July 2020).
  22. The BCG World Atlas. Available online: http://www.bcgatlas.org/index.php (accessed on 20 May 2020).
  23. The Global Health Observatory. Available online: https://www.who.int/data/gho (accessed on 20 May 2020).
  24. Kim, J.H. Multicollinearity and misleading statistical results. Korean J. Anesthesiol. 2019, 72, 558–569. [Google Scholar] [CrossRef] [Green Version]
  25. Escobar, L.E.; Molina-Cruz, A.; Barillas-Mury, C. BCG vaccine protection from severe coronavirus disease 2019 (COVID-19). Proc. Natl. Acad. Sci. USA 2020, 117, 202008410. [Google Scholar] [CrossRef]
  26. Weng, C.H.; Saal, A.; Butt, W.W.; Bica, N.; Fisher, J.Q.; Tao, J.; Chan, P.A. Bacillus Calmette-Guérin vaccination and clinical characteristics and outcomes of COVID-19 in Rhode Island, United States: A cohort study. Epidemiol. Infect. 2020, 148, e140. [Google Scholar] [CrossRef] [PubMed]
  27. Klinger, D.; Blass, I.; Rappoport, N.; Linial, M. Significantly improved COVID-19 outcomes in countries with higher BCG vaccination coverage: A multivariable analysis. Vaccines 2020, 8, 378. [Google Scholar] [CrossRef] [PubMed]
  28. Arts, R.J.W.; Moorlag, S.J.C.F.M.; Novakovic, B.; Li, Y.; Wang, S.Y.; Oosting, M.; Kumar, V.; Xavier, R.J.; Wijmenga, C.; Joosten, L.A.B.; et al. BCG vaccination protects against experimental viral infection in humans through the induction of cytokines associated with trained immunity. Cell Host Microbe 2018, 23, 89–100. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Hisaka, A.; Yoshioka, H.; Hatakeyama, H.; Sato, H.; Onouchi, Y.; Anzai, N. Global comparison of changes in the number of test-positive cases and deaths by coronavirus infection (COVID-19) in the world. J. Clin. Med. 2020, 9, 1904. [Google Scholar] [CrossRef]
  30. Yates, T.; Razieh, C.; Zaccardi, F.; Davies, M.J.; Khunti, K. Obesity and risk of COVID-19: Analysis of UK biobank. Prim. Care Diabetes 2020. [Google Scholar] [CrossRef]
  31. Ten Doesschate, T.; Moorlag, S.J.C.F.M.; van der Vaart, T.W.; Taks, E.; Debisarun, P.; Ten Oever, J.; Bleeker-Rovers, C.P.; Verhagen, P.B.; Lalmohamed, A.; Ter Heine, R.; et al. BCG-CORONA study team: Two randomized controlled trials of Bacillus Calmette-Guérin vaccination to reduce absenteeism among health care workers and hospital admission by elderly persons during the COVID-19 pandemic: A structured summary of the study protocols for two randomised controlled trials. Trials 2020, 21, 555. [Google Scholar] [CrossRef]
Figure 1. Histograms of morbidities and mortalities drawn as a normal density plot. (A) Morbidity per million; (B) mortality per million; (C) log10 transformed morbidity; and (D) log10 transformed mortality.
Figure 1. Histograms of morbidities and mortalities drawn as a normal density plot. (A) Morbidity per million; (B) mortality per million; (C) log10 transformed morbidity; and (D) log10 transformed mortality.
Ijerph 17 05589 g001
Figure 2. (A) Associations between morbidity, mortality and PCR-tests. Either Pearson’s correlation coefficient or Spearman’s rank correlation was applied to calculate rho; (B) associations between morbidity and risk factors were adjusted for PCR-tests per million population (log10); (C) associations between mortality and risk factors were adjusted for morbidity (log10) and PCR-tests (log10) per million population.
Figure 2. (A) Associations between morbidity, mortality and PCR-tests. Either Pearson’s correlation coefficient or Spearman’s rank correlation was applied to calculate rho; (B) associations between morbidity and risk factors were adjusted for PCR-tests per million population (log10); (C) associations between mortality and risk factors were adjusted for morbidity (log10) and PCR-tests (log10) per million population.
Ijerph 17 05589 g002
Figure 3. Scatter plot showing the association between PCR-test positivity and mortality. Mortality per million population and PCR-test positivity rates (%) on 17 July 2020 were transformed to the common logarithm (log10) in the graph. Since the variable of ‘PCR-test positivity (log10)’ showed a normal distribution, Pearson’s correlation coefficient was applied to calculate rho, to quantify the strength of the association. Countries that never had or stopped a national BCG vaccine program are indicated in red, while countries with current national BCG vaccine programs are indicated in black. Selected country names are shown using three-letter country codes.
Figure 3. Scatter plot showing the association between PCR-test positivity and mortality. Mortality per million population and PCR-test positivity rates (%) on 17 July 2020 were transformed to the common logarithm (log10) in the graph. Since the variable of ‘PCR-test positivity (log10)’ showed a normal distribution, Pearson’s correlation coefficient was applied to calculate rho, to quantify the strength of the association. Countries that never had or stopped a national BCG vaccine program are indicated in red, while countries with current national BCG vaccine programs are indicated in black. Selected country names are shown using three-letter country codes.
Ijerph 17 05589 g003
Figure 4. Scatter plot showing the association between urban population and COVID-19 morbidity rates. COVID-19-related morbidity rates per million population on 17 July 2020 were transformed to the common logarithm (log10) in the graph. Since the variable of ‘urban population’ showed a non-normal distribution, Spearman’s rank correlation was applied to calculate rho, to quantify the strength of the association. Countries that had never had or that had stopped a national program of BCG vaccination are indicated in red, while countries that currently follow a national BCG vaccine program are indicated in black. Selected country names are shown using three-letter country codes.
Figure 4. Scatter plot showing the association between urban population and COVID-19 morbidity rates. COVID-19-related morbidity rates per million population on 17 July 2020 were transformed to the common logarithm (log10) in the graph. Since the variable of ‘urban population’ showed a non-normal distribution, Spearman’s rank correlation was applied to calculate rho, to quantify the strength of the association. Countries that had never had or that had stopped a national program of BCG vaccination are indicated in red, while countries that currently follow a national BCG vaccine program are indicated in black. Selected country names are shown using three-letter country codes.
Ijerph 17 05589 g004
Figure 5. Scatter plot showing the association between insufficient physical activity and COVID-19 morbidity rates.
Figure 5. Scatter plot showing the association between insufficient physical activity and COVID-19 morbidity rates.
Ijerph 17 05589 g005
Figure 6. Scatter plot showing the association between ≥60 years of age (%) of the population and mortality rate. Mortalities per million population on 17 July 2020 transformed to the common logarithm (log10) are presented in the graph. Since the variable of ‘percentage of population ≥ 60 years of age (%)’ showed a non-normal distribution, Spearman’s rank correlation was applied to calculate rho, to quantify the strength of the association. Countries that had never had or that had stopped their national BCG vaccine program are indicated in red, while countries that currently follow a national BCG vaccine program are indicated in black. Selected country names are shown using three-letter country codes.
Figure 6. Scatter plot showing the association between ≥60 years of age (%) of the population and mortality rate. Mortalities per million population on 17 July 2020 transformed to the common logarithm (log10) are presented in the graph. Since the variable of ‘percentage of population ≥ 60 years of age (%)’ showed a non-normal distribution, Spearman’s rank correlation was applied to calculate rho, to quantify the strength of the association. Countries that had never had or that had stopped their national BCG vaccine program are indicated in red, while countries that currently follow a national BCG vaccine program are indicated in black. Selected country names are shown using three-letter country codes.
Ijerph 17 05589 g006
Figure 7. Scatter plot showing the association between COVID-19-related mortality and BCG vaccine coverage. Mortalities per million population on 17 July 2020 are transformed to the common logarithm (log10) in the graph. Since the variable of ‘BCG vaccine coverage’ showed a non-normal distribution, Spearman’s rank correlation was applied to calculate rho, to quantify the strength of the association. Selected country names are shown using three-letter country codes.
Figure 7. Scatter plot showing the association between COVID-19-related mortality and BCG vaccine coverage. Mortalities per million population on 17 July 2020 are transformed to the common logarithm (log10) in the graph. Since the variable of ‘BCG vaccine coverage’ showed a non-normal distribution, Spearman’s rank correlation was applied to calculate rho, to quantify the strength of the association. Selected country names are shown using three-letter country codes.
Ijerph 17 05589 g007
Table 1. Twenty countries with the highest and lowest COVID-19 morbidity rates and their PCR-test rate per million population.
Table 1. Twenty countries with the highest and lowest COVID-19 morbidity rates and their PCR-test rate per million population.
The Highest 20 CountriesThe Lowest 20 Countries
RankCountryMorbidity *PCR-Test *RankCountry *Morbidity *PCR-Test *
1Qatar37,566153,380173Papua New Guinea1798
2Chile16,92770,696172Lao People’s Democratic Republic32990
3Kuwait13,496105,205171Vietnam42824
4Oman12,24450,512170Myanmar61770
5Panama11,40640,384169United Republic of Tanzania9
6Armenia11,32447,718168Cambodia102896
7Andorra 11,15648,531167Angola18304
8USA 11,118137,544166Timor-Leste181189
9Peru10,35560,747165Taiwan193319
10Brazil946423,098164Uganda235133
11Luxembourg 8438502,852163Burundi23563
12Singapore8053172,506162Syrian Arab Republic27
13Sweden 761067,495161Fiji294461
14Saudi Arabia698371,623160Gambia321477
15Belarus6945123,003159Mozambique441360
16Spain 6543128,893158Niger45276
17UAE5673436,262157Thailand468647
18South Africa546439,182156Burkina Faso50
19Belgium 5438121,891155Yemen524
20Maldives5234118,769154Chad54
* Numbers per million population, USA—United States of America, UAE—United Arab Emirates. Bold letters with mark “∫” mean no recent national BCG vaccination program for all.
Table 2. Twenty countries with the highest and lowest COVID-19 mortality rates and their PCR-test rates per million population.
Table 2. Twenty countries with the highest and lowest COVID-19 mortality rates and their PCR-test rates per million population.
The Highest 20 CountriesThe Lowest 20 Countries
RankCountryMortality *PCR-Test *RankCountry *Mortality *PCR-Test *
1Belgium 845121,891173Papua New Guinea0798
2Andorra 67348,531172Lao People’s Democratic Republic02990
3UK 664186,591171Vietnam02824
4Spain608128,893170Cambodia02896
5Italy 579100,954169Timor-Leste01189
6Sweden 55467,495168Uganda05133
7France 46239,868167Fiji04461
8USA 426137,544166Eritrea0
9Peru38260,747165Mongolia08688
10Chile38170,696164Bhutan051,107
11Brazil36123,098163Saint Lucia011,494
12Netherlands 35844,588162Greenland090,957
13Ireland 354107,707161Dominica011,820
14Ecuador 29511,023160Saint Vincent and the Grenadines010,932
15Mexico2865958159Saint Kitts and Nevis013,248
16Canada 23488,515158Seychelles0
17Switzerland 22781,825157Burundi0.08563
18Panama22740,384156Myanmar0.11770
19Armenia20547,718155Taiwan0.33319
20North Macedonia19238,771154Mozambique0.31360
* Numbers per million population, USA—United States of America, UK—United Kingdom, VCT—Saint Vincent and the Grenadines. Bold letters with mark “∫” mean no recent national BCG vaccination program for all.
Table 3. Association between mortality as the outcome and morbidity as the exposure, adjusted for number of PCR-tests performed.
Table 3. Association between mortality as the outcome and morbidity as the exposure, adjusted for number of PCR-tests performed.
VariableCoef.Std. Err.p-Value95% CIAdjusted R2
Morbidity per million population (log10)1.3950.083<0.0011.232 to 1.5580.7067
PCR-tests per million population (log10)−0.2410.0860.005−0.416 to −0.075
Coef.—Coefficient; Std. Err.—Standard Error; 95% CI—95% confidence interval.
Table 4. Associations between COVID-19 morbidity and mortality per million population and world data for each country.
Table 4. Associations between COVID-19 morbidity and mortality per million population and world data for each country.
VariableN * 3Median, IQRMin (Country)
Max (Country)
Morbidity * 1
Adjusted R2
Mortality * 2
Adjusted R2
Population, n1731.02 × 107
3.5 × 106–3.3 × 107
3.9 × 105 (Monaco)
1.4 × 109 (China)
0.37320.7120
Yearly change (%)1731.11
0.35–1.96
−1.35 (Lithuania)
3.84 (Niger)
0.38220.7172
Net change, n17389,516
8516–447,563
−383,840 (Japan)
1.40 × 107 (India)
0.37130.7073
Population density, n/km217383
34–167
0 (Greenland)
26,337 (Monaco)
0.37230.7048
Land area, km2173183,630
35,410–581,795
1 (Monaco)
1.6 × 107 (Russia)
0.37260.7079
Migrants, n (net)168−1590
−16,026–14,440.5
−653,249 (Venezuela))
954,806 (USA)
0.38850.7460
Fertility rate, n1682.27
1.71–3.65
1.11 (Korea)
6.95 (Niger)
0.37840.7516
Urban population (%)16860
43–78
13 (Papua New Guinea)
98 (Belgium)
0.4246 *0.7243
World share (%)1730.1
0.05–0.4
0 (Saint Lucia, etc.)
18.5 (China)
0.37320.7120
Age
Median age, years16830
21.5–39
15 (Niger)
48 (Japan)
0.37870.7811 *
0 to 14 years of age (%)16526.8
17.6–37.3
12.3 (Niger)
49.8 (Singapore)
0.37890.7574
≥ 60 years of age (%)1659.72
5.30–20.48
2.80 (UAE)
34.01 (Japan)
0.38020.7986 *
≥ 70 years of age (%)1623.41
2.03–9.17
0.53 (UAE)
18.49 (Japan)
0.38180.8027 *
Economy
GDP, million US dollars17040,729
11,135–223,780
393 (Sao Tome)
1.90 × 107 (USA)
0.37010.7095
GDP per capita, US dollars17051.4
1534–15,347
104 (Somalia)
165,421 (Monaco)
0.37220.7051
Total unemployment rate (%)1645.85
3.75–9.40
0.1 (Qatar)
28.2 (Lesotho)
0.39330.7436
Male unemployment rate (%)1645.5
3–8.5
0 (Qatar)
26 (Eswatini)
0.39050.7433
Female unemployment rate (%)1647
4–12
0 (Niger)
38 (Syria)
0.40780.7435
Total labor force participation rate (%)16462
55–68
38 (Yemen)
87 (Qatar)
0.39390.7597
Male labor force participation rate (%)16473
66.5–79
45 (Moldova)
95 (Qatar)
0.38090.7602
Female labor force participation rate (%)16453.5
45–60.5
6 (Yemen)
84 (Rwanda)
0.42160.7479
National BCG vaccine program
Annual incidence of tuberculosis, n/100,00016457
17–176
2 (UAE)
834 (South Africa)
0.36930.7053
Recent BCG vaccine coverage (%)17393
81–98
0 (Italy, etc.)
99 (Japan, etc.)
0.38140.7415 *
Global health observatory
Health policy
International Health Regulations score17064
48 to 82
17 (Central African Republic)
99 (Canada, Russia)
0.38170.7205
Universal health coverage index1660.41
0.36 to 0.47
0.17 (Saint Lucia)
0.70 (Peru)
0.38120.7644 *
Total expenditure on health as a percentage of gross domestic product1696.38
4.76–8.43
1.48 (Timor-Leste)
17.14 (USA)
0.39770.7283
Population with household expenditures on health greater than 10% of total household expenditure/income (%)1466.56
3.13–12.76
0.20 (Gambia)
54.2 (Sierra Leone)
0.39560.7906
Medical personnel
Hospital beds, n/10,000 population16421
11–40
1 (Mali)
134 (Japan)
0.37180.7409
Medical doctors, n/10,000 population16915.7
3.3–29.8
0.1 (Tanzania)
84.2 (Cuba)
0.38000.7351 *
Nursing and midwifery personnel, n/10,000 population16926.6
9.4–64.4
0.06 (Cameroon)
201.6 (Monaco)
0.38230.7203
Licensed qualified anesthesiologists actively working, n141250
18–1511
0 (Congo)
194,634 (China)
0.39040.6924
Health biomarkers
High blood pressure (SBP > 140 OR DBP > 90) (crude estimate) (%)16723
20–25
13 (Peru)
41 (Croatia)
0.40910.7199
Elevated fasting blood glucose (>7.0 mmol/L or on medication) (crude estimate) (%)1678.1
5.9–9.7
2.6 (Burundi)
16.6 (Fiji)
0.40690.7162
Elevated total cholesterol (≥5.0 mmol/L or 193 mg/dL)
(crude estimate) (%)
16737.6
25.2–52.5
14.8 (Niger)
69.7 (Denmark)
0.40700.7401 *
BMI
Mean BMI, kg/m216726.2
23.4–27.2
20.5 (Ethiopia)
30.0 (Saint Lucia)
0.42330.7102
BMI ≥ 30 kg/m2 (%)16719.9
7.5–25.7
2.1 (Vietnam)
37.3 (USA)
0.43930.7135
BMI ≥ 25 kg/m2 (%)16753.5
27.8–62.2
17.9 (Timor-Leste)
72.1 (Kuwait)
0.44540.7204
Alcohol drinking, total consumption per capita among persons aged ≥ 15 years, liters of pure alcohol over a calendar year1686.4
2.6–9.8
0 (Bangladesh)
15.2 (Moldova)
0.38340.7191
Prevalence of smoking any tobacco product among males aged ≥15 years (%)11932
23–44
9 (Ethiopia)
76 (Indonesia)
0.44100.8006
Prevalence of smoking any tobacco product among females aged ≥ 15 years (%)1198
3–19
0 (Niger)
40 (Serbia)
0.39320.8192
Prevalence of insufficient physical activity among adults aged ≥ 15 years (crude estimate) (%)14328.7
19.0–36.2
5.0 (Uganda)
65.3 (Kuwait)
0.4984 *0.7270
Estimated population-based prevalence of depression (%)1664.4
4.0–5.0
3.0 (Papua New Guinea)
6.3 (Ukraine)
0.40620.7486
Mortality according to age group
Neonatal mortality rate, n per 1000 live births1709.6
3.5–21.0
0.9 (Japan)
42.0 (Pakistan)
0.38110.7234
Infantile mortality rate, n per 1000 live births170
14.0
5.9–34.9
1.4 (Finland)
84.5 (Central African Republic)
0.38100.7507
Under-five mortality rate, probability of dying by age 5/1000 live births17014.4
7.0–44.2
1.7 (Finland)
121.5 (Somalia)
0.38130.7145
Mortality rate for 5–14-year-olds, n/1000 children aged 5–14 years1682.8
1.5–8.6
0.4 (Luxembourg)
37.3 (Niger)
0.37810.7117
Adult mortality rate, probability of dying between 15 and 60 years of age/1000 population167150
96–224
49 (Switzerland)
484 (Lesotho)
0.38120.7377
Probability of dying between age 30 and exact age 70 from cardiovascular disease, cancer, diabetes or chronic respiratory disease16618.4
14.7–22.6
7.8 (Korea)
30.6 (Yemen)
0.38860.7061
Life expectancy at birth, years
At birth16673.4
65.3–81.5
52.9 (Lesotho)
84.2 (Japan)
0.38060.7600
At age 60 years16619.6
17.2–22.1
13.4 (Sierra Leone)
26.4 (Japan)
0.38380.7630
Healthy life expectancy (HALE), years
At birth16665.3
57.5–68.1
44.9 (Central African Republic)
76.2 (Singapore)
0.37950.7599
At age 60 years16614.8
13.0–17.1
10.3 (Sierra Leone)
21.0 (Singapore)
0.38250.7652 *
Disease-specific mortality
Chronic obstructive pulmonary disease1666.13.9–8.30.6 (Qatar)
36.8 (India)
0.39750.7436
Ischemic heart disease16621.5
13.4–35.2
4.9 (Brunei Darussalam)
106.9 (Georgia)
0.37890.7431
Lower respiratory infections1668.2
3.7–31.6
0.22 (Finland)
130.3 (Chad)
0.37910.7455
Stroke16610.2
5.6–16.6
1.4 (Qatar)
49.3 (Georgia)
0.39630.7726
Tracheal, bronchial, lung cancers1661.8
0.8–4.2
0.2 (Niger)
20.5 (China)
0.39710.7541
Air pollution
Ambient and household air pollution attributable death rate /100,000 population16665.9
38.5–97.8
8.5 (Georgia)
184 (Brunei Darussalam)
0.38610.7427
Concentrations of fine particulate matter (PM2.5)17021.1
14.6–34.7
5.7 (New Zealand)
94.3 (Nepal)
0.39400.7210
Coverage rate with the first dose of a measles-containing-vaccine (MCV1) among one-year-olds (%)16993
84–97
30 (Equatorial Guinea)
99 (Mongolia and 23 other countries)
0.38510.7095
IQR—interquartile range; GDP—gross domestic product; BMI—body mass index—weight (kg)/height (m)2; * 1—adjusted for PCR tests per million population; * 2—adjusted for PCR tests and COVID-19 morbidity per million population; * 3—number of countries the data were able to abstract, *: p < 0.001. Bold letters mean statistically significant: p < 0.001.
Table 5. Factors associated with morbidity adjusted for PCR-test rates.
Table 5. Factors associated with morbidity adjusted for PCR-test rates.
VariableCoef.Std. Err.p-Value95% CIAdjusted R2
PCR-tests per million population (log10)0.5740.101<0.0010.374 to 0.7740.5037
Urban population (%)0.7640.3290.020.113 to 1.416
Insufficient physical activity (%)0.0150.0060.010.004 to 0.026
Coef.—Coefficient; Std. Err.—Standard Error; 95% CI—95% confidence interval.
Table 6. Factors associated with mortality adjusted for morbidity and PCR-test rates.
Table 6. Factors associated with mortality adjusted for morbidity and PCR-test rates.
VariableCoef.Std. Err.p-Value95% CIAdjusted R2
Morbidity per million population (log10)1.3420.067<0.0011.210 to 1.4750.8254
PCR-tests per million population (log10)−0.4850.099<0.001−0.680 to −0.290
Population ≥ 60 years of age (%)0.0300.0100.0030.011 to 0.050
BCG vaccine coverage−0.0040.0020.006−0.007 to −0.001
Universal health coverage1.4730.5490.0080.387 to 2.560
Medical doctors/10,000 population0.0030.0040.47−0.005 to 0.011
Elevated total cholesterol0.0020.0060.76−0.010 to 0.014
HALE at age 60−0.030.0290.27−0.090 to 0.025
HALE—healthy life expectancy, Coef.—Coefficient; Std. Err.—Standard Error; 95% CI—95% confidence interval.

Share and Cite

MDPI and ACS Style

Urashima, M.; Otani, K.; Hasegawa, Y.; Akutsu, T. BCG Vaccination and Mortality of COVID-19 across 173 Countries: An Ecological Study. Int. J. Environ. Res. Public Health 2020, 17, 5589. https://doi.org/10.3390/ijerph17155589

AMA Style

Urashima M, Otani K, Hasegawa Y, Akutsu T. BCG Vaccination and Mortality of COVID-19 across 173 Countries: An Ecological Study. International Journal of Environmental Research and Public Health. 2020; 17(15):5589. https://doi.org/10.3390/ijerph17155589

Chicago/Turabian Style

Urashima, Mitsuyoshi, Katharina Otani, Yasutaka Hasegawa, and Taisuke Akutsu. 2020. "BCG Vaccination and Mortality of COVID-19 across 173 Countries: An Ecological Study" International Journal of Environmental Research and Public Health 17, no. 15: 5589. https://doi.org/10.3390/ijerph17155589

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop