- Research
- Open access
- Published:
Correlation between immunity from BCG and the morbidity and mortality of COVID-19
Tropical Diseases, Travel Medicine and Vaccines volume 6, Article number: 17 (2020)
Abstract
Background
To investigate the association between parameters indicating immunity from BCG at country level (presence of BCG vaccination policy, BCG coverage, age-specific incidence of tuberculosis (TB)) and the morbidity and mortality of COVID-19.
Methods
Country-specific data for COVID-19 cases and deaths, demographic details, BCG coverage and policy, age-specific TB incidence and income level were obtained. The crude COVID-19 cases and deaths per 100,000 population were calculated and assessed against the parameters indicating immunity from BCG using linear regression analysis.
Results
Univariate analysis identified higher income level of a country to be significantly associated with COVID-19 cases (p < 0.0001) and deaths (p < 0.0001) but not with its case fatality rate. The association between COVID-19 and TB was strongest for TB incidence in patients > 65-years (Cases (rs = − 0.785,p = 0.0001)) and deaths (rs = − 0.647,p = 0.0001).
Multivariate analysis identified the higher income level of a country and not having a universal BCG vaccination policy to affect the COVID-19 cases. The deaths were inversely affected by the presence of BCG vaccination policy and coverage; and positively by the TB incidence in patients > 65-years.
Conclusion
Significant inverse correlations observed between cases and deaths of COVID-19 and BCG related parameters highlights immunity from BCG as a likely explanation for the variation in COVID-19 across countries.
Introduction
The COVID-19 outbreak started in December 2019 in Wuhan Province, China [1] and has by now spread to six continents. It is caused by SARS-CoV-2 [2]. It has overwhelmed the health system capacity in many countries. No specific treatment has been effective to date, while the thrust in controlling it is placed on public health interventions. Available literature suggests the mortality to be higher in patients of advanced age and comorbidities [1].
Bacillus Calmette–Guérin (BCG) is a live-attenuated vaccine primarily developed to protect against childhood meningitis and disseminated tuberculosis (TB) [3]. The World Health Organization (WHO) recommends BCG vaccination at birth in countries with both universal and selective vaccination policies [4]. Furthermore, multiple meta-analyses have confirmed the effectiveness of BCG against pulmonary and extra-pulmonary TB in adults [5,6,7]. It is however assumed that prior exposure to environmental mycobacteria would confer a partial resistance and thereby diminish the response to BCG vaccine. Thus, vaccination in older age groups is only recommended in a few selected scenarios, while there is no universal policy for BCG vaccination of immunocompromised or other high-risk individuals.
Several researchers have suggested an inverse relationship between immunity from BCG and the case and death rates from COVID-19 [8, 9]. Considering the high morbidity and mortality of COVID-19 in the elderly, a possible protective effect through BCG in this group would invariably be of clinical significance. Nevertheless, ecological studies conducted so far have failed to establish this relationship [10]. It could be that these studies have all used a single indicator to depict immunity from BCG, which is commonly the population coverage of BCG vaccination. This value does not reflect the immunisation status of different age groups. At the same time, the BCG immunisation status is also not routinely assessed in different age groups. Therefore, considering the protection that BCG offers against pulmonary and extra-pulmonary TB, we hypothesised that the incidence of TB in different age groups could be used as a surrogate marker to indicate poor immunity from BCG in those age groups.
BCG vaccination is believed to confer a non-specific increase in immunity [11], and acts via both innate and adaptive immune responses. The use of BCG as an immunotherapy in bladder cancer and melanoma [12] is an extension of this theory. We believe this increase in immunity may contribute to the differences in COVID-19 morbidity and mortality.
We further hypothesized that the income level of each country is likely to confound this relationship, as the morbidity and mortality from COVID-19 have been consistently higher in countries with a higher mean GDP. The majority of these countries had also abandoned universal BCG vaccination in the preceding decades.
Thus, the primary aim of this study was to investigate the population level associations of BCG coverage, incidence of TB in different age groups and income status of the countries with the morbidity and mortality of COVID-19.
Methods
Data sources
Country-specific data for COVID-19 cases and deaths were retrieved from the Johns Hopkins Coronavirus Resource Centre [13]. Data on population were obtained from the World Bank population data [14]. This included the total population, population of 15–64 years and > 65 years. BCG coverage data were obtained from the BCG Atlas [15] and the World Health Organization data on immunisation coverage [16]. Data on TB burden was obtained from the World Health Organization Global Health Observatory Data Repository [17]. The income level of the country was obtained from the World Bank population data and open data repository [18]. Data were retrieved on 11/04/2020.
To reduce the bias of the availability of testing, we calculated the crude COVID-19 mortality per 100,000 population using population data published by the World Bank.
Statistical analysis
Data were analysed using IBM SPSS statistics for Windows, version 26.0. Correlations were assessed using the Spearman correlation coefficient (rs). Tests applied to assess the significance between groups were Mann-Whitney test for scaled data, Chi-Square test for categorical data, and Kruskal-Wallis test for multiple groups of continuous data.
Linear regression analysis using backward LR method was performed to assess the relationship between COVID-19 and immunity from BCG adjusted for confounders. In the models, the morbidity or mortality were used as the dependent variable, while the independent variables were proxy indicators of the status of immunity from BCG of populations, namely, the presence or absence of a universal BCG vaccination policy, the WHO estimates of the BCG coverage of the total population, the BCG coverage in 1-year olds, TB incidence in 15–64 year age group and TB incidence in > 65-year age group. The income level of each country was also added to the model, as it could confound a true association between the dependent and independent variables. A probability value of < 0.05 was considered statistically significant.
Results
A total of 225 countries were eligible to be included in the analysis. Tuberculosis and population-level data were available for 214 countries. COVID-19 and BCG coverage data were available for 180 and 175 countries, respectively.
Distribution of COVID-19 cases and deaths by income level of countries
The COVID-19 cases were significantly different between countries of different income levels (Kruskal-Wallis test, p < 0.0001) (Fig. 1). All differences except that between low and lower-middle-income countries were statistically significant (Mann-Whitney test, p < 0.0005 for all). Similarly, COVID-19 deaths were significantly different between countries in different income levels (Kruskal-Wallis test, p < 0.0001) (Fig. 2). All differences except that between low and lower-middle-income countries were statistically significant (Mann-Whitney test, p < 0.005 for all). There was no statistically significant difference in the case fatality rate and the level of income (Fig. 3).
Relationship of COVID-19 cases and deaths with immunity from BCG at country level
The BCG coverage was inversely correlated with the incidence of TB (rs = − 0.164, p = 0.03). However, it showed no significant correlation with the cases or deaths per 100,000 from COVID-19 (p = 0.43 and 0.247, respectively). Instead, as shown in Table 1, the countrywide TB incidence had a strong inverse correlation with the COVID-19 cases per 100,000 population (rs = − 0.743, p = 0.0001) and COVID-19 deaths per 100,000 population (rs = − 0.611, p = 0.0001). Similar correlations were observed for both pulmonary (both bacteriologically confirmed and clinically confirmed) and extra-pulmonary TB. The association between COVID-19 and TB was strongest when comparing the TB incidence in patients over 65-years, in which both COVID-19 cases per 100,000 population (rs = − 0.785, p = 0.0001) and COVID-19 deaths per 100,000 population (rs = − 0.647, p = 0.0001) showed significant inverse correlations. The correlation between the younger age group (15–64 years) was also significant, but less strong (Cases (rs = − 0.742, p = 0.0001) and deaths (rs = − 0.621, p = 0.0001)).
We further assessed the relationship between the incidence of TB and COVID-19 cases and deaths stratified by countries with BCG coverage < 95% vs. ≥95% (Table 1). All correlations with COVID-19 deaths were higher in countries with BCG coverage < 95% than ≥95%. Similar observations were observed for COVID-19 cases except for three parameters, which were marginally higher in countries with higher BCG coverage. The Case Fatality Rate (CFR) showed a weak inverse correlation with TB deaths in countries (rs = − 0.208, p = 0.006).
When stratified according to income level, the BCG coverage showed a statistically significant moderate inverse correlation with both cases (rs = − 0.545, p < 0.0005) and deaths (rs = − 0.564, p < 0.0005) within the high-income group countries. No other income group showed a statistically significant correlation with BCG coverage. The correlation of COVID-19 cases and deaths with the incidence of TB in the population over 65-years progressively increased according to the income status of the country (Figs. 4 & 5).
Multivariate analysis
Regarding the linear regression performed for COVID-19 cases per 100,000 population, the model was statistically significant (F (3,62) = 27.028, p < 0.0005) and explained 54.6% of the variance observed. It identified a statistically significant effect on the caseload from the income status of a country; and a significant inverse effect from the presence of BCG vaccination policy. The WHO estimates of the BCG coverage of the population were not statistically significant (Table 2).
With regard to the linear regression performed for COVID-19 deaths per 100,000 population, the model was statistically significant (F(3,144) = 43.694, p < 0.0005) and explained 46.6% of the variance observed. It identified a statistically significant effect on deaths, from the incidence of TB in > 65-year old population, and an inverse effect from both the presence of BCG vaccination policy and WHO estimates of the BCG coverage of the population (Table 3).
The linear regression analysis using the case fatality rate as the dependent variable did not provide a statistically significant result.
Discussion
This study examines a previously observed association between the proxy measures of immunity from BCG of countries and the morbidity and mortality from COVID-19. We explored this association further by adopting the TB incidence in different age-groups as a surrogate of BCG coverage, as no age-specific BCG coverage data are currently available.
In summary, our findings indicate that the BCG coverage did not have a significant correlation with the cases or deaths from COVID-19. The countrywide TB incidence as well as pulmonary and extra-pulmonary TB incidence had strong inverse correlations with both cases and deaths. This association was strongest when comparing the TB incidence in patients over 65-years. Furthermore, all correlations with COVID-19 deaths were higher in countries with BCG coverage < 95% than ≥95%. Our multivariate analysis identified the presence of a universal BCG vaccination policy to inversely affect both cases and deaths from COVID-19.
The impact of BCG vaccination on the transmission of Mycobacterium tuberculosis is limited, as it does not prevent the reactivation of latent pulmonary infection, which is the principal source of bacillary spread in the community. There is significant variation in the effectiveness of BCG vaccination, and prior exposure to environmental mycobacteria is thought to influence this [19]. This variation may also explain the inconclusive evidence shown for the protection of BCG vaccination against COVID-19. In contrast, a reduced incidence of TB is a more direct indicator of the efficacy of BCG vaccination. This may explain our observations.
The presence of a universal BCG vaccination policy showed a stronger correlation than the BCG coverage, for both cases and deaths from COVID-19. In concurrence, in the study by Sala et al. [9], BCG coverage only accounted for 12.5% of the variance observed. The results of Miller et al. [20] suggest that factors increasing the fraction of the elderly population being vaccinated against BCG, like an earlier year of initiation of a universal BCG vaccination policy, had a good correlation with COVID-19 cases and deaths.
BCG vaccination is believed to confer a non-specific increase in immunity [11]. It acts via both innate and adaptive immune responses, and in the latter, via both T and B cells. This theory is applied when BCG is used as an immunotherapy in the management of melanoma and bladder cancer [12]. BCG induces long-lived memory B-cells [21]. BCG also induces and maintains tissue-resident CD4 T cells in the lungs [22]. COVID-19 is presumed to enter the body through the lungs [21]. These tissue-resident CD4 cells may therefore influence the entry into the body and therefore contribute to lower cases and deaths in countries with BCG vaccination. A clinical trial has already begun testing the utility of BCG for boosting immunity against COVID-19 [23].
Over 90% of the TB cases in the elderly are due to reactivation of primary infection [24]. Rarely, previously infected older persons who have eliminated the viable tubercle bacilli may revert to a “naïve” immunologic status and are at risk of reinfection [25]. These mechanisms may imply the loss of the initial protection from BCG vaccination. The better correlation of COVID-19 cases and deaths with the TB incidence in older ages groups than the younger groups also supports this. The changes occurring in the immune system with ageing may also contribute to this [26, 27]. The reduced incidence of TB could also represent an overall improvement in preventive care. The results we observed may be related to an unknown organism or antigen closely correlated to both TB and COVID-19, similar to the hygiene hypothesis [28].
The mortality from COVID-19 is reportedly higher in patients who are older and comorbid [1, 29]. In the study by Zhou et al. [1], the interquartile range (IQR) of the age of survivors and non-survivors was 45–58 and 63–76 years, respectively. This shows a clear demarcation with possibly no-overlap. Our analysis only identified the TB incidence in > 65-year olds to be significantly inversely correlated with deaths. The incidence in the total population and the 15–64-year group showed smaller correlations. We used the incidence of TB is a surrogate marker of the immunity offered by the BCG vaccination. Lesser deaths observed in populations where there is a lesser number of TB in the > 65-year olds support this hypothesis. There is no worldwide age-specific case and mortality data for COVID-19 and we were therefore unable to investigate this hypothesis further.
The relationship between the income and COVID-19 cases is paradoxical to observations during previous epidemics [30, 31]. High-income countries generally depend on larger budgetary allocations for secondary preventive strategies to contain an epidemic within the health system capacity. Their outcomes during the COVID-19 pandemic have remained worse than countries with lower income due to a multitude of reasons. The income level remained independently associated with the case number in our multivariate analysis and has been reported by other researchers as well [3, 8, 20, 32]. In the subgroup analysis, we observed that, for high-income countries, the correlation between cases and deaths and the degree of BCG coverage was higher. Most high-income countries do not have a universal BCG vaccination policy, but income level and BCG coverage were independent predictors in our multivariate analysis. One possible contributing factor to the higher cases and deaths in high-income countries could be the higher median age of the population in these countries [33]. A positive correlation between the median age of the population and COVID-19 morbidity and mortality has been reported [32].
We observed minimal correlation between the CFR and the variables studies. The CFR can be affected by the lack of testing to identify cases. There is a significant difference in the testing per-capita among countries, i.e. high-income countries performing more tests per million population [34]. There was also an inverse correlation between the percentage of the population tested and the positive rate [32]. This may have caused the apparent reduction in CFR compared to the deaths per 100,000 population in high- and upper-middle-income countries. The lag-between deaths and cases may also contribute to this [8], as the pandemic is still evolving. The CFR, therefore, may be unreliable in the COVID-19 pandemic.
Non-pharmacological methods like social distancing, quarantining and isolation and changes in population behaviour have been successful in reducing rates of transmission and hospitalization [35, 36]. There is, however, no reliable worldwide data available for these measures, and we have not included them in our hypothesis.
The main limitation of this analysis is the inability to link the exposures relate to BCG vaccination to the occurrence of COVID-19 in the same person (ecological fallacy) as the data considered are at the population level. Furthermore, there may be other factors between countries that are associated with the exposure, which may account for the differences, such as comorbidities [1] and other epidemiological characteristics of the vulnerable groups, which have not been considered in the model. Several countries who do not have a universal BCG vaccination at present have stopped the vaccination in the last 20-years. Their most recent estimated BCG coverage is still high. We have also not included other variables like temperature and humidity, which have already been proven to be associated with COVID-19 [37]. The use of both BCG coverage as well as surrogate indicators of BCG coverage (i.e. incidence of pulmonary and extra-pulmonary TB, TB incidence in different age groups) has enabled a more robust evaluation of the association between BCG vaccination and COVID-19 morbidity and mortality. We did not attempt to analyse a difference between males and females, as this data was only available for a few countries.
Conclusions
We saw a significant inverse relationship between Universal BCG vaccination policy, BCG coverage and the incidence of TB in > 65-year olds with COVID-19 cases and deaths. Although we cannot imply causality, this is a potentially useful intervention in a disease which presently has no definitive preventive or curative strategy.
Availability of data and materials
Available with corresponding author upon requesting.
References
Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet. 2020;395:1054–62.
Wu D, Wu T, Liu Q, Yang Z. The SARS-CoV-2 outbreak: what we know. Int J Infect Dis. 2020;94:44–8.
Calmette A. L'infection bacillaire et la tuberculose chez l'homme et chez les animaux: processus d'infection, et de defense, étude biologique et expérimentaleMasson, et cie; 1922.
World Health O. BCG vaccine: WHO position paper, February 2018 – recommendations. Vaccine. 2018;36:3408–10.
Mangtani P, Abubakar I, Ariti C, Beynon R, Pimpin L, Fine PE, et al. Protection by BCG vaccine against tuberculosis: a systematic review of randomized controlled trials. Clin Infect Dis. 2014;58:470–80.
Abubakar I, Pimpin L, Ariti C, Beynon R, Mangtani P, Sterne JA, et al. Systematic review and meta-analysis of the current evidence on the duration of protection by bacillus Calmette-Guerin vaccination against tuberculosis. Health Technol Assess. 2013;17:1–372 v-vi.
Roy A, Eisenhut M, Harris RJ, Rodrigues LC, Sridhar S, Habermann S, et al. Effect of BCG vaccination against Mycobacterium tuberculosis infection in children: systematic review and meta-analysis. BMJ (Clinical research ed). 2014;349:g4643.
Shet A, Ray D, Malavige N, Santosham M, Bar-Zeev N. Differential COVID-19-attributable mortality and BCG vaccine use in countries. medRxiv. 2020.
Sala G, Miyakawa T. Association of BCG vaccination policy with prevalence and mortality of COVID-19. medRxiv; 2020.
Cooper JA, Ryan R, Parsons N, et al. The use of electronic healthcare records for colorectal cancer screening referral decisions and risk prediction model development. BMC Gastroenterol. 2020;20:78. https://doi.org/10.1186/s12876-020-01206-1.
Covián C, Fernández-Fierro A, Retamal-Díaz A, Díaz FE, Vasquez AE, Lay MK, et al. BCG-induced cross-protection and development of trained immunity: implication for vaccine design. Front Immunol. 2019;10:2806.
Talat Iqbal N, Hussain R. Non-specific immunity of BCG vaccine: a perspective of BCG immunotherapy. Trials Vaccinol. 2014;3:143–9.
Coronavirus COVID-19 Global Cases by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University [https://coronavirus.jhu.edu/map.html]. Accessed 11 Apr 2020.
Population, total [https://data.worldbank.org/indicator/SP.POP.TOTL]. Accessed 11 Apr 2020.
Zwerling A, Behr MA, Verma A, Brewer TF, Menzies D, Pai MJP. The BCG world atlas: a database of global BCG vaccination policies and practices. PLoS Med. 2011;8:e1001012.
Baccille Calmette Guérin vaccine Reported estimates of BCG coverage [https://apps.who.int/immunization_monitoring/globalsummary/timeseries/tscoveragebcg.html]. Accessed 11 Apr 2020.
Global Health Observatory data repository [https://apps.who.int/gho/data/node.main.1315?lang=en]. Accessed 11 Apr 2020.
World Bank Population and Open Data repository [https://data.worldbank.org/income-level/low-income]. Accessed 11 Apr 2020.
Vaudry W. "to BCG or not to BCG, that is the question!". The challenge of BCG vaccination: why can't we get it right? Paediatr Child Health. 2003;8:141–4.
Miller A, Reandelar MJ, Fasciglione K, Roumenova V, Li Y, Otazu GH. Correlation between universal BCG vaccination policy and reduced morbidity and mortality for COVID-19: an epidemiological studymedrxiv; 2020.
Sahin AR, Erdogan A, Agaoglu PM, Dineri Y, Cakirci AY, Senel ME, et al. Novel coronavirus (COVID-19) outbreak: a review of the current literature. Eurasian J Med Oncol. 2019;2020(4):1–7.
Bull NC, Kaveh DA, Garcia-Pelayo MC, Stylianou E, McShane H, Hogarth PJ. Induction and maintenance of a phenotypically heterogeneous lung tissue-resident CD4+ T cell population following BCG immunisation. Vaccine. 2018;36:5625–35.
BCG Vaccination to Protect Healthcare Workers Against COVID-19 (BRACE) [https://clinicaltrials.gov/ct2/show/NCT04327206#moreinfo]. Accessed 11 Apr 2020.
Yoshikawa TT. Tuberculosis in aging adults. J Am Geriatr Soc. 1992;40:178–87.
Rajagopalan S, Yoshikawa TT. Tuberculosis in the elderly. Z Gerontol Geriatr. 2000;33:374–80.
Stervbo U, Meier S, Mälzer JN, Baron U, Bozzetti C, Jürchott K, et al. Effects of aging on human leukocytes (part I): immunophenotyping of innate immune cells. Age (Dordr). 2015;37:92.
Stervbo U, Bozzetti C, Baron U, Jürchott K, Meier S, Mälzer JN, et al. Effects of aging on human leukocytes (part II): immunophenotyping of adaptive immune B and T cell subsets. Age (Dordr). 2015;37:1–13.
Strachan DP. Hay fever, hygiene, and household size. Bmj. 1989;299:1259–60.
Ruan Q, Yang K, Wang W, Jiang L, Song J. Clinical predictors of mortality due to COVID-19 based on an analysis of data of 150 patients from Wuhan, China. Intensive Care Med. 2020;46:846–8.
Dawood FS, Iuliano AD, Reed C, Meltzer MI, Shay DK, Cheng P-Y, et al. Estimated global mortality associated with the first 12 months of 2009 pandemic influenza a H1N1 virus circulation: a modelling study. Lancet Infect Dis. 2012;12:687–95.
Murray CJL, Lopez AD, Chin B, Feehan D, Hill KH. 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–8.
Bhoj R, Singh RG. Are BCG vaccination, population density, median age and poverty important determinants of COVID-19 pandemic spread, morbidity and mortality? 2020.
Population Data by World Bank income group [https://apps.who.int/gho/data/view.main.POP2030]. Accessed 11 Apr 2020.
COVID-19 CORONAVIRUS PANDEMIC [https://www.worldometers.info/coronavirus/#countries]. Accessed 11 Apr 2020.
Cowling BJ, Ali ST, Ng TWY, Tsang TK, Li JCM, Fong MW, et al. Impact assessment of non-pharmaceutical interventions against coronavirus disease 2019 and influenza in Hong Kong: an observational study. Lancet Public Health. 2020;5:e279–88.
Prem K, Liu Y, Russell TW, Kucharski AJ, Eggo RM, Davies N, et al. The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study. Lancet Public Health. 2020;5:e261–70.
Sajadi MM, Habibzadeh P, Vintzileos A, Shokouhi S, Miralles-Wilhelm F, Amoroso A. Temperature, Humidity, and Latitude Analysis to Estimate Potential Spread and Seasonality of Coronavirus Disease 2019 (COVID-19). JAMA Netw Open. 2020;3(6):e2011834. https://doi.org/10.1001/jamanetworkopen.2020.11834.
Code availability
Not applicable.
Funding
None.
Author information
Authors and Affiliations
Contributions
All authors conceptualized the study. DW, NW and SK performed the data collection. All authors were involved in the data analysis and preparation of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not required.
Consent for publication
Not applicable.
Competing interests
None of the authors declare any conflicts of interest.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wickramasinghe, D., Wickramasinghe, N., Kamburugamuwa, S.A. et al. Correlation between immunity from BCG and the morbidity and mortality of COVID-19. Trop Dis Travel Med Vaccines 6, 17 (2020). https://doi.org/10.1186/s40794-020-00117-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40794-020-00117-z