Systematic review with meta-analysis of the epidemiological evidence in Europe, Israel, America and Australasia on smoking and COVID-19

Abstract

BACKGROUND

Previous meta-analyses related smoking to death or severe infection from coronavirus disease 2019 (COVID-19) in hospitalized patients, but considered only a few studies, did not adjust for demographics and comorbidities, and inadequately defined smoking.

AIM

To review and meta-analyse epidemiological evidence on smoking and COVID-19, considering a range of endpoints, populations and smoking definitions and the effect of adjustment.

METHODS

Studies were identified from publications in English up to 30 September, 2020 involving at least 100 individuals, carried out in Europe, Israel, America or Australasia, not restricted to those with specific other diseases, and providing information relating smoking to various COVID-related endpoints. Meta-analyses were carried out for combinations of population and endpoint, with variation studied by smoking definition, adjustment level and other factors.

RESULTS

From 96 publications, 74 studies were identified, 37 in the United States, 10 in the United Kingdom, with up to four in the other countries. Three involved over a million individuals, and 37 involved less than a thousand. Adjusted results for smoking were available in 42 studies, with adjustment not considered in 20 studies. Results were considered by endpoint. No significant effect of smoking on COVID-19 positivity was seen in the general population, but there was a reduced risk in those tested. Best-adjusted estimates for current (vs never) smoking were 0.87 (95% confidence interval: 0.52-1.47) in the general population and 0.52 (0.43-0.64) in those tested. For those hospitalized due to COVID-19, unadjusted rates were significantly increased in current smokers (1.20, 1.01-1.42) and ever smokers (1.64, 1.41-1.91), but those adjusted for comorbidities showed no increase for current (0.82, 0.52-1.30) or ever smokers (1.00, 0.76-1.32). There was little evidence to suggest that smoking was associated with intensive care admission. For those hospitalized with COVID-19, best-adjusted estimates were 0.88 (0.72-1.08) for current smokers and 1.10 (0.99-1.22) for ever smokers. In those hospitalized with COVID-19, smoking was not significantly related to subsequent mechanical ventilation, with best-adjusted estimates of 1.12 (0.60-2.09) for current smokers and 1.05 (0.88-1.25) for ever smokers. For those hospitalized with severe COVID-19, best-adjusted estimates were 0.74 (0.49-1.12) for current smokers and 1.15 (0.87-1.51) for ever smokers; few estimates were adjusted for comorbidities. While smoking was associated with increased mortality in unadjusted analyses, the association disappeared after adjustment for comorbidities. For example, in those hospitalized with COVID-19, the unadjusted estimate for ever smokers of 1.59 (1.37-1.83) reduced to 1.07 (0.82-1.38) when adjusted for comorbidities. Studies on those with severe COVID-19 showed that smoking tended to be associated with worsening of the disease. However, no estimate was adjusted, even for demographics. Estimates did not clearly vary by location or study size, and there was too little evidence to usefully study variations by age, amount smoked or years quit.

CONCLUSION

The increased COVID-19 death rate in smokers seen in unadjusted analyses disappears following adjustment for demographics and comorbidities. Among those tested, smoking is associated with lower COVID-19 infection rates.

Key Words: Smoking; COVID-19; Meta-analyses; Review; Europe; America

Core Tip: Detailed analyses of 74 studies related smoking to being tested for coronavirus disease 2019 (COVID-19), having COVID-19, or suffering death or severe disease due to COVID-19. Various smoking indices were studied, as were the effects of adjusting for other factors. Although many studies provided limited unreliable results, consistent evidence showed that of those tested, smokers were less likely to have COVID-19. Among those positive for or hospitalized with COVID-19, there was a clear association between smoking and COVID-19 death and severity in unadjusted analyses, which disappeared following adjustment for comorbidities and demographics. Any adverse effects in smokers appear to derive from their poorer prior health status.

INTRODUCTION

In a previous project commenting on publications on smoking and coronavirus disease 2019 (COVID-19), we considered over 100 papers published up to the end of September 2020. Among these were various meta-analyses falling into two groups.

Eight publications[1-8] considered smoking prevalence in hospitalized patients, generally agreeing it was substantially less than expected from national statistics. This evidence does not necessarily show smoking protects against acquiring COVID-19. Smoking may be markedly under-reported in studies based on medical records. Also, among those with COVID-19, smokers might be less likely than non-smokers to be hospitalized. These meta-analyses ignored relevant information from studies of those tested for COVID-19, or of the general population, where smoking habits were collected pre-pandemic, as well as more recent studies of hospitalized patients.

The other meta-analyses[3,6,9-28] concerned hospitalized patients, relating smoking to severity, progression or death from COVID-19, mainly from studies in China. While these generally reported positive associations which were often statistically significant, many meta-analysis estimates were unadjusted even for age, with comorbidities present pre-pandemic rarely considered. These meta-analyses also varied on the index of smoking used, which was not always clearly defined.

Here we describe meta-analyses aimed at avoiding the limitations of the early meta-analyses by considering more studies, not limiting attention to hospitalized patients, and paying particular attention to the definition of smoking and the effects of adjustment, as well as the reliability of the smoking data.

To limit the scope of the study and provide timely results various restrictions to the studies were made, as described in the methods section. Notably studies in China and other parts of Asia, except for Israel, were excluded. A recent review[29] classified few studies from Asia as “good” or “fair”, most having much missing data and/or not reliably distinguishing current, former, ever and never smoking status. We were also aware of a large study in Israel[30] without these weaknesses.

Due to the more comprehensive data included, and better quality of the studies considered, the meta-analyses we describe should provide much better insight into the relationships of smoking to various COVID-related outcomes than do the meta-analyses referred to above.

MATERIALS AND METHODS

Full details of the methods used are given in Supplementary material 1 and are summarized below.

Study inclusion and exclusion criteria

Pre-defined criteria stipulated for practical reasons that studies should be described in English and were detected in searches up to September 2020. They should also be conducted in Europe, Israel, America, or Australasia, as discussed above. America here includes all the countries in South and Central America, as well as the United States and Canada. Studies restricted to individuals with specific other diseases were excluded as being less generalizable. Studies of less than 100 individuals were excluded as they provided inadequate power to detect reliable results, and as there were adequate numbers of larger studies. Studies should provide information relating smoking to the probability of one or more of the following relevant endpoints: being tested for COVID-19, having confirmed COVID-19, having self-reported COVID-19, being hospitalized with COVID-19, requiring mechanical ventilation for COVID-19, requiring intensive care for COVID-19, having severe/progressive COVID-19, or dying from COVID-19 or from any cause. The studies may concern various at-risk populations, including the general population, those tested for COVID-19, those positive for COVID-19, those hospitalized with COVID-19, or those with severe COVID-19.

Literature searches

As part of our earlier project, we carried out a first PubMed search on April 7, 2020, and then carried out further daily searches up to September 30, 2020.

Publications identified in our study as being of initial interest were then examined to identify studies satisfying our inclusion/exclusion criteria, and relevant meta-analyses. The meta-analyses were then examined for additional relevant studies meeting our inclusion criteria. Further studies were then sought from meta-analyses identified in further searches, from a further more detailed look at our original searches, and from examining reference lists of publications identified as relevant.

Multiple publications from the same study

To avoid double-counting results from the same study reported in multiple publications, the relevant publications were examined to identify publications from the same study.

Data recorded

Data from each publication were entered onto a study database and a linked effect estimate database. The study database recorded information on: publications considered; study title; study location; sexes, ages and races considered; study dates; study type; nature of population studied; sample size; definition of severe COVID-19 (if relevant); method of COVID-19 diagnosis; whether adjusted effect estimates were available; the confounding variables studied; and whether adjusted results on smoking were reported and, if not, why not (e.g., smoking not significant in the adjusted model). It also recorded information on the smoking index for which results were available (e.g., current vs never smoking), the source of the smoking data (e.g., medical records), the extent of missing data, the percentage of smokers in the population studied, and whether dose-response data were available, as well as details of the endpoints and at-risk populations studied.

The effect estimate database included details of each individual effect estimate entered. Effect estimates were entered for every available combination of endpoints within a population, smoking index and level of adjustment (separated into unadjusted (U), adjusted for demographics only (D), adjusted also for comorbidities but not post-infection variables (C), and adjusted also for variables including post-infection responses to COVID-19 (P). Where available, effect estimates were entered by sex, age group and smoking dose (amount, duration, time since quitting). Other factors recorded included the publication the effect estimate was derived from; the population and endpoint considered; the smoking comparison; the type of effect estimate [odds ratio (OR), relative risk (RR) or hazard ratio (HR)]; the adjustment factors considered; the number of cases and at-risk subdivided by the smoking variable; the effect estimate and its lower and upper 95% confidence interval; and whether the estimate was given in the paper or was derived from the data presented. Derivation could be from the 2 × 2 table of numbers for unadjusted data, or using the method described by Hamling et al[31] to derive adjusted estimates for other smoking indices from estimates given in the publication (e.g., ever vs never estimates from those for current vs never and former vs never).

All data were entered by Hamling JS or Coombs KJ and checked by Lee PN, with any disagreements discussed and resolved.

Meta-analyses

For each studied combination of endpoint within a population (e.g., died while hospitalized with COVID-19), meta-analyses were carried out relating the endpoint to each of six indices of smoking; ever vs never, current vs non-current, current vs never, former vs never, a combined index most closely approximating to current vs never, and a combined index most closely approximating to ever vs never. The combined index for current vs never smoking includes, from each study reporting the endpoint/population combination, results in the following preference order (most to least preferred) – current vs never, current vs non-current, smoker (undefined) vs non-smoker, tobacco use vs none, ever vs never, and former vs never. The combined index for ever vs never smoking uses the preference - ever vs never, former vs never, smoker (undefined) vs non-smoker, tobacco use vs none, current vs never, and current vs non-current.

For each endpoint within a combination of population and smoking index, the results to be meta-analysed were selected using a first preference on level of adjustment and then a second preference on type of effect estimate. For adjustment, the order of preference (first to last) was adjustment for factors including comorbidities, adjustment for demographics only, unadjusted, and adjustment for factors measuring responses to COVID-19 infection (the lowest preference, as this may be a form of over-adjustment). For type of effect estimate, the preference order was HR, then RR, then OR.

Where the numbers of estimates permitted, the meta-analyses compare estimates by level of adjustment, type of estimate, sex, location and study size.

Statistical analysis

Fixed-effect and random-effects meta-analyses were conducted using the method of Fleiss and Gross[32] with heterogeneity quantified by H, the ratio of the heterogeneity chisquared to its degrees of freedom. H is directly related to the I² statistic[33] by the formula I² = 100(H−1)/H. For all meta-analyses, Egger’s test of publication bias[34] was included.

All analyses were carried out using RoeLee release 63, build 52, available from RoeLee Statistics Ltd (http://www.roelee.co.uk).

RESULTS

Literature searches

Figure 1 summarizes the results of the literature searches, with fuller details, including reasons for rejecting papers (Supplementary material 2). Overall, 98 publications met the selection criteria.

Figure 1 Flowchart of literature search. The flowchart shows, for each of the four stages of the process, the numbers of papers and meta-analyses considered, together with the numbers of papers meeting the selection criteria or rejected, with reasons for rejection shown.

Studies and study characteristics

Twelve references represented multiple publications from the same study. Allowing for this, data were entered for 76 separate studies. Subsequently, during data entry, it became apparent that two studies[35,36] provided no data for any of the endpoints considered, and only compared smoking prevalence with published data in the population at large. These studies were not further evaluated.

Supplementary material 3 summarizes the details for the 74 studies, including references. Studies were identified by the six character codes shown there.

Thirty-seven studies were from the United States, 10 from the United Kingdom (including seven restricted to England), four each from France, Israel, Italy, Mexico and Spain, two from Switzerland and one each from Australia, Brazil and Denmark. Also one study was conducted in the United Kingdom and Italy, and one in multiple European countries.

In 31 studies, the population considered was patients hospitalized with COVID-19, with a further four patients admitted to the intensive care unit (ICU). Also 19 studies considered those with a positive COVID-19 test, 11 included those who were tested for COVID-19 and seven included the general population. One study included hospital patients and non-COVID-19 controls, while another included those tested for COVID-19 as well as control groups not tested.

All studies included both sexes, and none selected individuals on race or ethnicity. One study was restricted to adults aged 47-87 years, with a further 41 restricted to adults, generally with a minimum age of 18 years, but sometimes having lower limits ranging from 15 years to 23 years. The remaining 32 studies did not refer to any age restriction. As shown in Supplementary material 3, the number of individuals with smoking data varied widely between the studies. The largest was OPENSA in England involving over 16 million individuals, with two other United Kingdom studies over a million, and three others over 100000. In contrast, 37 studies involved less than 1000 individuals.

It is also shown in Supplementary material 3 that some studies reported analyses for subsets of their populations, and for various endpoints. These endpoints included being tested for COVID-19, having confirmed or having self-reported COVID-19, being hospitalized with COVID-19, requiring intensive care or mechanical ventilation, having severe or progressive COVID-19, and mortality, either from COVID-19 or from all causes combined.

Thirty-six studies allowed calculation of separate effect estimates for current, former and ever smokers. Eight reported results for current smoking only, 19 for smoking history only, and eight for smoking undefined. One (TWIGG) reported results for tobacco use undefined, one (CHAND) had two source papers - one providing results for current smokers only, one for ever smokers only - and another (GUPTA) gave numbers of current or former smokers combined and adjusted results for current vs non-current smokers. The smoking data were mainly extracted from medical records, although in some cases the data came from a questionnaire or other sources, with no details given in six studies. Percentages with missing data on smoking were available in 35 studies, with 13 being over 20%.

Adjusted results for smoking were presented for 42 studies. In 12, only predictors other than smoking appeared in the adjusted model, either because smoking was not significant in univariate analyses, so was not considered in multivariate analyses, or because smoking dropped out of the multivariate modelling. In 20 studies, no adjusted results were presented.

In 54 studies, COVID-19 diagnosis was by reverse transcriptase-polymerase chain reaction or simply by polymerase chain reaction, but there were various alternatives, as indicated in Supplementary material 3.

Effect estimates available

Effect estimates were available for all but one study (KNIGHT), which stated only that “Ever cigarette smoking was predictive of death (P < 0.05)” without providing quantitative detail.

Overall there were 738 effect estimates: 548 ORs, 122 RRs and 68 HRs. The studies providing most estimates were GU (100), BIOBNK (86) and VETERA (52), with 18 other studies providing 10 or more estimates. Fourteen studies provided only one estimate. There were 153 estimates for current vs never smoking, 202 for current vs non-smoking, 142 for former vs never smoking, 223 for ever vs never smoking, 16 for smoker vs non-smoker not otherwise stated, and two for tobacco vs no-tobacco not otherwise stated. One study (HIPPIS) provided 15 effect estimates by amount smoked, and one (TOOLKI) provided four estimates by years quit. No other study provided dose-response data.

Of the 738 estimates, 432 (58.5%) were unadjusted, 110 (14.9%) adjusted for demographics only, 169 (22.9%) adjusted for variables including co-morbidities but not responses to COVID-19, and 27 (3.7%) adjusted for a list including responses to COVID-19.

One study (YANOVE) provided effect estimates subdivided jointly by sex and age group. Four others (BIOBNK, GDEM, HOPKIN, MIYARA) provided estimates subdivided by sex only, and three others (GDEM, SINAI, VETERA) provided estimates subdivided by age group only.

The at-risk population was the general population for 155 (21.0%) effect estimates, those tested for COVID-19 for 106 (14.4%), those positive for COVID-19 for 222 (30.1%), those hospitalized with COVID-19 for 188 (25.5%), and those in intensive care for 22 (3.0%). Other populations were considered for 45 (6.1%) effect estimates, with GU providing 36 of these, all based on populations representing a combination of a group considered above (e.g., tested for COVID-19) and unmatched controls. Other such populations included a mixture of those positive for COVID-19 and those untested (two estimates from BIOBNK), those hospitalized but including non-COVID cases (five from MEINI), those tested for COVID-19 but not hospitalized (one from ADORNI) and those hospitalized with COVID-19 but not in ICU (one from PELLAU).

The endpoints considered in the 738 effect estimates were tested for COVID-19 in 41 (5.6%), positive for COVID-19 in 189 (25.6%), hospitalized for COVID-19 in 128 (17.3%), in intensive care in 92 (12.5%), mechanical ventilation in 58 (7.9%), severe COVID-19 (defined variously) in 58 (7.9%), and died in 172 (23.3%).

The effect estimates concerned many different combinations of endpoint within population, the commonest being hospitalized among those positive for COVID-19 (98 effect estimates), positive within those tested (90), died among those hospitalized (82) and positive among the general population (79).

The results below are considered by endpoint and within endpoint by population. Meta-analysis results are shown in the main tables for endpoint/population combinations that have data from at least five studies, with the individual study data summarized in Supplementary material 4 for combinations with data from fewer than five studies. Fuller details, including the individual effect estimates meta-analysed, the extent of heterogeneity between the estimates, and results of tests for publication bias, are provided in Supplementary material 5.

Endpoint: Tested for COVID-19

Four studies provided effect estimates (Supplementary Table 1). GU compared tested individuals with unmatched controls, and the remaining studies (BIOBNK, ADORNI and HOPKIN) compared those tested and not tested. BIOBNK provided estimates from multiple publications, the results from two being shown in SupplementaryTable 1, with the others providing little extra information.

Table 1 Main results for endpoint: Tested positive for coronavirus disease 2019.

Selection of estimates1		Statistic2	Ever vs never	Current vs non-current	Current vs never	Former vs never	Closest to current	Closest to never
Population = general3
Best-adjusted	All	EE (95%CI)	1.04 (0.79-1.37)	0.76 (0.56-1.02)	0.87 (0.52-1.47)	1.09 (0.94-1.28)	0.78 (0.57-1.06)	0.89 (0.73-1.09)
		n	5	7	5	5	7	7
	C	EE (95%CI)	1.10 (0.70-1.72)	0.93 (0.43-2.01)	0.96 (0.42-2.17)	1.10 (0.86-1.40)	0.96 (0.42-2.17)	1.10 (0.70-1.72)
		n	3	3	3	3	3	3
	U	EE (95%CI)	0.71 (0.63-0.80)	0.56 (0.43-0.75)	0.45 (0.36-0.55)	0.92 (0.80-1.05)	0.56 (0.42-0.75)	0.65 (0.55-0.77)
		n	1	3	1	1	3	3
All estimates	U	EE (95%CI)	1.04 (0.83-1.29)	0.79 (0.51-1.22)	0.90 (0.46-1.78)	1.13 (0.97-1.31)	0.80 (0.53-1.22)	0.89 (0.74-1.08)
		n	5	7	5	5	7	7
Population = tested for COVID-194
Best-adjusted	All	EE (95%CI)	0.71 (0.60-0.84)	0.61 (0.52-0.71)	0.52 (0.43-0.64)	0.97 (0.86-1.09)	0.55 (0.47-0.65)	0.70 (0.63-0.77)
		n	10	13	8	8	17	17
	C	EE (95%CI)	0.77 (0.56-1.05)	0.60 (0.45-0.79)	0.53 (0.34-0.83)	0.90 (0.76-1.08)	0.58 (0.43-0.80)	0.77 (0.64-0.92)
		n	4	5	4	4	5	5
	U	EE (95%CI)	0.67 (0.52-0.86)	0.62 (0.51-0.75)	0.49 (0.42-0.58)	1.04 (0.85-1.28)	0.55 (0.45-0.67)	0.67 (0.58-0.78)
		n	6	8	4	4	11	11
All estimates	U	EE (95%CI)	0.72 (0.62-0.85)	0.61 (0.50-0.73)	0.53 (0.46-0.61)	0.97 (0.87-1.09)	0.55 (0.46-0.66)	0.71 (0.64-0.79)
		n	10	12	8	8	17	17

¹Best-adjusted = estimates for a study selected in preference order of adjustment (C = comorbidities, D = demographics only, U = unadjusted, P = includes post-coronavirus disease responses).

²EE: Effect estimate from random-effects meta-analysis.

³For population = general, the results are from six studies, one providing sex-specific estimates. Four provided odds ratio (OR)s, one hazard ratios and one relative risk (RR)s and ORs.

⁴For population = tested, the results are from 15 studies, two providing sex-specific estimates, with 13 providing only ORs and two RRs, one of these also providing ORs. COVID-19: Coronavirus disease 2019; CI: Confidence interval.

The results presented are somewhat conflicting, with the BIOBNK study reporting estimates above 1.00, generally significant, and tending to decrease with increasing adjustment, while ADORNI and HOPKIN, both only provided unadjusted results, with significant estimates below 1.00. GU shows a reduced probability of testing for current smoking and an increased probability of testing for former and ever smoking, with the estimates reducing with increasing adjustment.

Meta-analysis of this conflicting data was not attempted.

Endpoint: Positive for COVID-19

Table 1 presents the meta-analysis results from six studies where the at-risk population was the general population, and 15 where it was those tested for COVID-19. The “best-adjusted” results are those where, for each study, an effect estimate was selected in the order of preference C, D, U and P for level of adjustment. Results are shown for all the best-adjusted results and for those where C and U were the best-adjusted results available, which form the great majority of the best-adjusted results. The “all estimates U” results give estimates for the totality of unadjusted results, including those not included in the best-adjusted results being superseded by results for the same study with level of adjustment C or D. The results are consistent with no effect of smoking on positivity in the general population, but a reduced risk of positivity, particularly among current smokers, in those tested for COVID-19. No clear effects of adjustment were seen in either analysis.

Four other studies presented results for COVID-19 positivity based on other populations (Supplementary Table 2). In BIOBNK, COVID-19 positives were compared with the untested population, in ADORNI the population was those non-hospitalized, in MEINI the population was those hospitalized, whether or not from COVID-19, while GU compared tested individuals and unmatched controls. Although the results from BIOBNK suggested smokers were more likely to be positive, those from the other studies did not, the results for current smokers showing a negative association.

Table 2 Main results for endpoint: Hospitalized with coronavirus disease 2019.

Selection of estimates1		Statistic2	Ever vs never	Current vs non-current	Current vs never	Former vs never	Closest to current	Closest to ever
Population = positive for COVID-193
Best-adjusted	All	EE (95%CI)	1.34 (1.13-1.60)	1.00 (0.86-1.16)	1.05 (0.84-1.32)	1.31 (0.93-1.85)	1.31 (1.12-1.54)	1.38 (1.17-1.63)
		n	19	13	11	11	22	22
	C	EE (95%CI)	1.00 (0.76-1.32)	0.90 (0.76-1.07)	0.82 (0.52-1.30)	0.87 (0.67-1.12)	0.96 (0.79-1.18)	0.93 (0.80-1.08)
		n	6	6	4	4	8	8
	U	EE (95%CI)	1.53 (1.14-2.06)	1.03 (0.75-1.41)	1.13 (0.81-1.59)	1.77 (1.04-3.04)	1.44 (1.03-2.01)	1.71 (1.24-2.36)
		n	10	6	6	6	11	11
All estimates	U	EE (95%CI)	1.64 (1.41-1.91)	0.97 (0.71-1.32)	1.20 (1.01-1.42)	1.75 (1.35-2.27)	1.47 (1.21-1.77)	1.72 (1.46-2.02)
		n	18	13	11	11	21	21

¹Best-adjusted = estimates for a study selected in preference order of adjustment (C = comorbidities, D = demographics only, U = unadjusted, P = includes post-COVID responses).

²EE = effect estimate from random-effects meta-analysis.

³The results are from 21 studies, one study providing sex-specific estimates. Nineteen studies provide odds ratio (OR)s, one relative risk (RR)s and hazard ratios, and one ORs and RRs. COVID-19: Coronavirus disease 2019; CI: Confidence interval.

Endpoint: Hospitalized for COVID-19

Table 2 shows the meta-analysis results based on 19 studies of those positive for COVID-19. While the unadjusted estimates show increased hospitalization rates in former and ever smokers, those adjusted for comorbidities show no indication of an increase for any index of smoking.

Supplementary Table 3 shows additional results from three studies. BIOBNK provided results based on the general population, GDEM and GU provided results on those tested for COVID-19, and GU provided results based on those hospitalized and unmatched controls. The results are rather conflicting, with GU showing markedly lower hospitalization rates in current smokers, and the other studies increased rates. In BIOBNK, adjustment tended to reduce the associations, although they remained statistically significant. GDEM reported only results adjusted for demographics and GU reported only unadjusted results.

Table 3 Main results for endpoint: Admitted to intensive care unit.

Selection of estimates1		Statistic2	Ever vs never	Current vs non-current	Current vs never	Former vs never	Closest to current	Closest to ever
Population = positive for COVID-193
Best-adjusted	All	EE (95%CI)	1.57 (1.14-2.17)	0.87 (0.67-1.12)	0.94 (0.60-1.48)	1.55 (1.00-2.42)	1.13 (0.82-1.54)	1.40 (1.09-1.80)
		n	7	6	5	5	8	8
	C	EE (95%CI)	1.68 (1.15-2.46)	0.85 (0.77-0.94)	0.61 (0.11-3.36)	1.58 (0.97-2.57)	1.14 (0.54-2.41)	1.30 (0.74-2.31)
		n	2	2	1	1	3	3
	U	EE (95%CI)	1.42 (0.94-2.14)	0.94 (0.57-1.57)	0.98 (0.60-1.60)	1.55 (0.91-2.65)	0.98 (0.60-1.60)	1.42(0.94-2.14)
		n	4	4	4	4	4	4
All estimates	U	EE (95%CI)	1.52 (1.04-2.22)	1.03 (0.56-1.88)	0.92 (0.60-1.43)	1.68 (1.05-2.68)	1.11 (0.62-2.00)	1.61 (1.13-2.29)
		n	5	6	5	5	6	6
Population = hospitalized with COVID-194
Best-adjusted	All	EE (95%CI)	1.10 (0.99-1.22)	0.90 (0.72-1.12)	0.88 (0.72-1.08)	1.21 (0.99-1.47)	1.02 (0.87-1.21)	1.11 (0.98-1.25)
		n	9	9	6	6	14	14
	C	EE (95%CI)	1.01 (0.88-1.16)	0.80 (0.61-1.04)	0.81 (0.61-1.07)	1.07 (0.92-1.24)	0.81 (0.61-1.07)	1.01 (0.88-1.16)
		n	1	1	1	1	1	1
	U	EE (95%CI)	1.22 (1.04-1.43)	1.09 (0.93-1.27)	0.96 (0.70-1.31)	1.30 (1.00-1.70)	1.16 (1.05-1.28)	1.19 (1.08-1.30)
		n	8	7	5	5	12	12
All estimates	U	EE (95%CI)	1.20 (1.09-1.32)	1.00 (0.87-1.13)	0.89 (0.74-1.06)	1.28 (1.10-1.48)	1.05 (0.95-1.17)	1.14 (1.04-1.24)
		n	9	9	6	6	14	14

¹Best-adjusted = estimates for a study selected in preference order of adjustment (C = comorbidities, U = unadjusted, P = includes post coronavirus disease responses).

²EE = effect estimate from random-effects meta-analysis.

³For population = positive, the results are from eight studies. Six provided odds ratio (OR)s and two relative risk (RR)s;

⁴For population = hospitalized, the results are from 14 studies. Twelve provided ORs and two RRs. COVID-19: Coronavirus disease 2019; CI: Confidence interval.

Endpoint: Admitted to ICU

Table 3 shows the meta-analysis results based on estimates from eight studies of those positive for COVID-19 and 14 of those hospitalized with COVID-19. SupplementaryTable 4 shows the results from one general population study, and one study comparing intensive care patients with unmatched controls. Most estimates considered in Table 3 are unadjusted, not even for age, and show little evidence of an association between smoking and ICU admission. Exceptionally, the data in Supplementary Table 4 show reduced admission rates in current smokers, and increased admission rates in former smokers, tending to diminish and become marginally significant after adjustment in the HIPPIS study.

Table 4 Main results for endpoint: Mechanically ventilated.

Selection of estimates1		Statistic2	Ever vs never	Current vs non-current	Current vs never	Former vs never	Closest to current	Closest to ever
Population = hospitalized for COVID-193
Best-adjusted	All	EE (95%CI)	1.05 (0.88-1.25)	1.11 (0.72-1.71)	1.12 (0.60-2.09)	1.16 (0.87-1.55)	1.08 (0.86-1.36)	1.11 (0.92-1.35)
		n	12	8	5	5	15	15
	C	EE (95%CI)	1.04 (0.79-1.36)	1.94 (0.40-9.49)	2.77 (0.28-27.2)	2.43 (0.32-18.5)	0.99 (0.77-1.29)	1.04 (0.79-1.36)
		n	7	2	2	2	7	7
	U	EE (95%CI)	1.06 (0.94-1.21)	1.12 (0.60-2.09)	0.76 (0.59-0.98)	1.10 (0.96-1.27)	1.16 (0.80-1.67)	1.24 (0.99-1.54)
		n	5	5	3	3	7	7
All estimates	U	EE (95%CI)	1.07 (0.92-1.26)	1.10 (0.80-1.51)	1.09 (0.60-1.95)	1.12 (0.98-1.28)	1.12 (0.98-1.38)	1.17 (1.00-1.37)
		n	10	8	5	5	13	13

¹Best-adjusted = estimates for a study selected in preference order of adjustment (C = comorbidities, D = demographics only, U = unadjusted, P = includes post-coronavirus disease responses).

²EE = effect estimate from random-effects meta-analysis.

³Fourteen studies provided results, one providing estimates for two age groups. Thirteen studies provided odds ratios and one relative risks. COVID-19: Coronavirus disease 2019; CI: Confidence interval.

Endpoint: Mechanically ventilated

Fourteen studies provided results where the population involved patients hospitalized with COVID-19. While the best-adjusted effect estimates in Table 4 were greater than 1.0 for each smoking index, none were statistically significant at P < 0.05.

Supplementary Table 5 summarizes the results from three studies where the population was those tested for COVID-19. VETERA, which provided the most detailed results, did not demonstrate any clear association, with estimates for ever and for former smoking significantly increased when unadjusted, but close to 1.0 and non-significant when adjusted for comorbidities. In ESSVRD, a significant unadjusted increase again was non-significant after adjustment for comorbidities. Exceptionally, MUNOZP reported a very high OR for ever smoking after adjustment for comorbidities.

Table 5 Results for endpoint: Severe coronavirus disease 2019¹.

Selection of estimates2			Statistic3	Ever vs never	Current vs non-current	Current vs never	Former vs never	Closest to current	Closest to ever
Population = positive4
Best-adjusted	All		EE (95%CI)	1.39 (1.07-1.80)	1.04 (0.91-1.20)	1.27 (0.82-1.95)	1.20 (0.62-2.32)	1.10 (0.97-1.25)	1.10 (0.97-1.25)
			n	7	5	4	4	8	8
	C		EE (95%CI)	-	1.02 (0.88-1.18)	-	-	1.02 (0.88-1.18)	1.02 (0.88-1.18)
			n	0	1	0	0	1	1
	U		EE (95%CI)	1.30 (0.94-1.79)	1.26 (0.82-1.93)	1.27 (0.82-1.95)	1.20 (0.62-2.32)	1.35 (0.95-1.92)	1.30 (0.94-1.79)
			n	5	4	4	4	5	5
All estimates	U		EE (95%CI)	1.63 (1.27-2.10)	1.17 (1.09-1.25)	1.27 (0.82-1.95)	1.20 (0.62-2.32)	1.48 (1.13-1.94)	1.45 (1.12-1.86)
			n	7	5	4	4	8	8
Population = hospitalized5
Best-adjusted		All	EE (95%CI)	1.15 (0.87-1.51)	1.05 (0.93-1.18)	0.74 (0.49-1.12)	1.07 (0.87-1.32)	1.08 (0.94-1.25)	1.09 (0.98-1.22)
			n	7	5	2	2	10	10
		C	EE (95%CI)	1.35 (0.66-2.75)	0.76 (0.50-1.15)	0.77 (0.51-1.17)	1.06 (0.86-1.31)	1.23 (0.46-3.27)	1.35 (0.66-2.75)
			n	2	1	1	1	2	2
		U	EE (95%CI)	1.38 (1.04-1.83)	1.12 (1.04-1.21)	0.29 (0.03-2.64)	1.25 (0.56-2.76)	1.17 (0.99-1.37)	1.14 (1.06-1.23)
			n	4	3	1	1	6	6
All estimates		U	EE (95%CI)	1.40 (1.20-1.63)	1.14 (1.04-1.24)	0.89 (0.54-1.46)	1.48 (1.23-1.80)	1.17 (1.08-1.27)	1.21 (1.12-1.29)
			n	6	5	2	2	9	9

¹In most studies the endpoint was severe coronavirus disease 2019 (COVID-19) as defined by the author, but it was admitted to ICU or died in two studies, pneumonia in one, increased oxygen required in one and COVID-19 progressed in one.

²Best-adjusted = estimates for a study selected in preference order of adjustment (C = comorbidities, U = unadjusted, P = includes post-COVID responses).

³EE = effect estimate from random-effects meta-analysis.

⁴For population = positive, the results are from five studies, one study providing four estimates by age and sex. All estimates are odds ratio (OR)s.

⁵For population = hospitalized, the results are from 10 studies, nine providing ORs and one relative risks and hazard ratios. COVID-19: Coronavirus disease 2019; CI: Confidence interval.

Endpoint: Severe COVID-19

Table 5 shows the meta-analysis results from five studies on those positive for COVID-19 and 10 studies on those hospitalized with COVID-19. As shown in Table 5, definitions of severity varied by study. Few effect estimates were adjusted for comorbidities. The smoking indices were generally associated with a small increase in severity, but this was only significant at P<0.05 in one of the 12 best-adjusted meta-analysis estimates.

Endpoint: Died

Table 6 summarizes the results using data from 25 studies of those hospitalized with COVID-19, and eight studies of those positive for COVID-19 regardless of hospitalization. The estimates adjusted for comorbidities were virtually never statistically significant and usually close to 1.00, but the unadjusted estimates were nearly always elevated and often statistically significant. This was very clearly illustrated by the results for “closest to ever smoking” where about a two-fold increase was seen for the unadjusted results, with little or no increase seen for the comorbidity adjusted results, regardless of the population studied. It is also clear that higher unadjusted estimates were seen for former or ever smoking than for current smoking, perhaps because former smokers tend to be older than current or never smokers.

Table 6 Main results for endpoint: Died.

Selection of estimates1		Statistic2	Ever vs never	Current vs non-current	Current vs never	Former vs never	Closest to current	Closest to ever
Population = hospitalized with COVID-193
Best-adjusted	All	EE (95%CI)	1.42 (1.19-1.70)	0.99 (0.89-1.10)	1.12 (0.98-1.29)	1.61 (1.34-1.92)	1.27 (1.08-1.49)	1.42 (1.22-1.65)
		n	18	10	10	10	25	25
	C	EE (95%CI)	1.07 (0.82-1.38)	0.64 (0.12-3.35)	1.02 (0.70-1.49)	1.15 (1.02-1.29)	1.09 (0.84-1.40)	1.12 (0.91-1.39)
		n	8	2	3	3	10	10
	U	EE (95%CI)	1.78 (1.40-2.26)	0.95 (0.77-1.16)	1.18 (0.91-1.52)	1.92 (1.48-2.50)	1.52 (1.22-1.91)	1.79 (1.47-2.19)
		n	9	6	6	6	13	13
All estimates	U	EE (95%CI)	1.59 (1.37-1.83)	1.00 (0.88-1.13)	1.12 (0.97-1.29)	1.80 (1.48-2.17)	1.33 (1.18-1.51)	1.58 (1.39-1.78)
		n	15	10	9	9	20	20
Population = with COVID-194
Best-adjusted	All	EE (95%CI)	1.41 (0.93-2.14)	0.96 (0.92-1.01)	0.98 (0.72-1.32)	1.61 (0.79-3.29)	1.08 (0.82-1.42)	1.34 (0.96-1.86)
		n	7	6	5	5	8	8
	C	EE (95%CI)	0.99 (0.87-1.13)	0.96 (0.92-1.01)	0.85 (0.67-1.07)	1.04 (0.91-1.19)	0.96 (0.92-1.01)	0.97 (0.93-1.02)
		n	3	4	3	3	4	4
	U	EE (95%CI)	2.07 (1.12-3.82)	0.93 (0.60-1.46)	1.49 (0.94-2.36)	2.70 (1.11-6.58)	1.59 (1.40-1.80)	2.07 (1.12-3.82)
		n	3	2	2	2	3	3
All estimates	U	EE (95%CI)	2.04 (1.39-3.02)	1.06 (0.39-2.90)	1.07 (0.72-1.58)	2.58 (1.31-5.07)	1.52 (0.89-2.60)	2.22 (1.51-3.27)
		n	5	5	4	4	6	6

¹Best-adjusted = estimates for a study selected in preference order of adjustment [C = comorbidities, D = demographics only, U = unadjusted, P = includes post-coronavirus disease (COVID) responses.

²EE = effect estimate from random-effects meta-analysis.

³For population = hospitalized, the results are from twenty four studies, one providing age-specific estimates. Seventeen provided only odds ratio (OR)s, two only hazard ratio (HR)s, two only relative risk (RR)s, and three RRs and HRs, one of these also providing ORs.

⁴For population = with COVID-19, the results are from eight studies, six providing ORs, one RRs and one HRs and RRs. COVID-19: Coronavirus disease 2019; CI: Confidence interval.

Table 6 does not include results where the population studied was those admitted to the ICU, as these form a subset of those reported in Table 7; see the next section.

Table 7 Results for endpoint: Worsened or died¹.

Selection of estimates2	Statistic3	Ever vs never	Current vs non-current	Current vs never	Former vs never	Closest to current	Closest to ever
Unadjusted4	EE (95%CI)	1.34 (1.00-1.79)	1.31 (1.06-1.62)	1.21 (0.35-4.17)	2.01 (0.67-6.09)	1.27 (1.15-1.39)	1.36 (1.15-1.62)
	n	5	5	3	3	9	9

¹Seven studies considered death among those admitted to the intensive care unit (ICU), and one of these and one other considered death among those mechanically ventilated. One considered admission to ICU or death among those hospitalized with severe coronavirus disease 2019.

²All estimates were unadjusted.

³EE = effect estimate from random-effects meta-analysis.

⁴The data are from nine studies, five reporting odds ratio (OR)s, two relative risk (RR)s and two RRs and hazard ratios, one of these also reporting ORs. CI: Confidence interval.

Results were also available from three studies based on other populations (Supplementary Table 6). OPENSA provided estimates based on the general population, and these seem consistent with the pattern shown in Table 6. For former smoking, for example, an unadjusted estimate of 2.53 [95% confidence interval (CI): 2.43-2.63] reduced to 1.19 (1.14-1.24) after adjustment for comorbidities. Two other studies only reported unadjusted results. PELLAU found no increase in smokers (undefined) when comparing hospitalized patients with and without COVID-19, while GU reported an increased risk of death in former and ever smokers and a decreased risk in current smokers, whether the tested population was considered or whether decedents were compared to unmatched controls.

Not considered above was the KNIGHT study, which provided no effect estimates, merely stating that ever cigarette smoking predicted death from COVID-19.

Other endpoints

As shown in Table 7, nine studies reported results on the endpoint worsened or died, based on those with severe disease. Most estimates relate to death among those in the ICU or those requiring mechanical ventilation. There was a tendency for smoking to be positively associated with the endpoint. However, each estimate was unadjusted even for demographic variables.

Two studies reported results for endpoints worse than hospitalization among those tested for COVID-19 (see Supplementary Table 7). Based on estimates adjusted for demographics only, GDEM reported a significant increased risk of pneumonia in current smokers, but no increase in intensive care admissions or need for mechanical ventilation. Similar to the results shown in Supplementary Table 6 for death, and based on unadjusted estimates, GU reported an increased risk of ICU admission in former and ever smokers and a decreased risk in current smokers.

Consistency of results in subgroups

Table 8 compares best-adjusted effect estimates by level of adjustment, effect estimate type, location and study size separately for the indices of smoking closest to current smoking and closest to ever smoking, and for the six combinations of endpoint and population where data were available for at least 10 studies.

Table 8 Variations in best-adjusted effect estimates by level of adjustment, effect estimate type, location and study size.

Factor/level	Endpoint	CP	H	IC	MV	S	M
	Population	T	CP	H	H	H	H
	Estimates	17	22	14	15	10	25
Estimates for closest to current smoking
Adjustment	U	0.55 (0.45-0.67)	1.44 (1.03-2.01)	1.16 (1.05-1.28)	1.16 (0.80-1.67)	1.17 (0.99-1.37)	1.52 (1.22-1.91)
	D	0.23 (0.09-0.59)	1.52 (1.18-1.95)	0.73 (0.56-0.96)	0.85 (0.67-1.07)	1.02 (0.89-1.16)	1.01 (0.87-1.17)
	C	0.58 (0.43-0.80)	0.96 (0.79-1.18)	0.81 (0.61-1.07)	0.99 (0.77-1.29)	1.23 (0.46-3.27)	1.09 (0.84-1.40)
	PC	-	-	-	-	0.71 (0.41-1.23)	-
	P	NS	< 0.05	< 0.001	NS	NS	< 0.05
Estimate type	OR	0.54 (0.46-0.63)	1.26 (0.97-1.64)	1.07 (0.87-1.32)	1.13 (0.88-1.44)	1.09 (0.93-1.27)	1.35 (1.07-1.70)
	RR	0.64 (0.31-1.34)	1.25 (1.20-1.31)	0.90 (0.71-1.15)	0.73 (0.56-0.95)	-	1.11 (0.87-1.40)
	HR	-	1.10 (0.97-1.24)	-	-	1.02 (0.64-1.63)	1.15 (0.89-1.50)
	P	NS	NS	NS	< 0.05	NS	NS
Location	United States	0.45 (0.37-0.54)	1.27 (1.11-1.46)	1.00 (0.81-1.22)	1.02 (0.83-1.26)	1.06 (0.72-1.56)	1.21 (0.95-1.53)
	Other	0.64 (0.55-0.73)	1.11 (0.48-2.57)	1.08 (0.79-1.49)	1.05 (0.65-1.71)	1.10 (1.02-1.17)	1.34 (1.07-1.69)
	P	< 0.01	NS	NS	NS	NS	NS
Study size	< 500	-	0.99 (0.58-1.69)	1.25 (0.91-1.71)	1.35 (1.02-1.80)	1.02 (0.65-1.59)	1.29 (0.97-1.72)
	500-	0.40 (0.27-0.61)	1.42 (1.17-1.72)	1.07 (0.69-1.65)	0.92 (0.34-2.46)	1.62 (1.14-2.29)	1.57 (1.20-2.05)
	5000-	0.51 (0.39-0.66)	1.43 (0.97-2.10)	0.64 (0.29-1.42)	0.83 (0.70-0.98)	0.81 (0.56-1.19)	1.03 (0.83-1.28)
	50,000+	0.67 (0.51-0.87)	1.00 (0.85-1.18)	0.97 (0.72-1.31)	0.98 (0.61-1.59)	1.09 (0.99-1.19)	0.94 (0.82-1.07)
	P	NS	< 0.05	NS	< 0.05	NS	< 0.01
Estimates for closest to ever smoking
Adjustment	U	0.67 (0.58-0.78)	1.71 (1.24-2.36)	1.19 (1.08-1.30)	1.24 (0.99-1.54)	1.14 (1.06-1.23)	1.79 (1.47-2.19)
	D	0.23 (0.09-0.59)	1.41 (1.04-1.92)	0.73 (0.56-0.96)	0.85 (0.67-1.07)	1.02 (0.89-1.16)	1.11 (0.81-1.53)
	C	0.77 (0.64-0.92)	0.93 (0.80-1.08)	1.01 (0.88-1.16)	1.04 (0.79-1.36)	1.35 (0.66-2.75)	1.12 (0.91-1.39)
	PC	-	-	-	-	0.71 (0.41-1.23)	-
	P	< 0.05	< 0.001	< 0.01	NS	NS	< 0.01
Estimate type	OR	0.67 (0.59-0.76)	1.39 (1.08-1.79)	1.16 (0.99-1.37)	1.14 (0.91-1.42)	1.10 (0.98-1.24)	1.53 (1.19-1.96)
	RR	0.87 (0.59-1.27)	1.25 (1.20-1.31)	1.01 (0.89-1.15)	1.02 (0.88-1.18)	-	1.16 (1.05-1.29)
	HR	-	1.03 (0.95-1.11)	-	-	1.02 (0.64-1.63)	1.28 (1.14-1.44)
	P	NS	< 0.001	NS	NS	NS	NS
Location	United States	0.65 (0.50-0.84)	1.38 (1.14-1.67)	1.10 (0.99-1.23)	1.08 (0.90-1.30)	1.15 (0.87-1.51)	1.43 (1.11-1.84)
	Other	0.74 (0.66-0.82)	1.31 (0.64-2.71)	1.08 (0.79-1.49)	1.05 (0.65-1.71)	1.10 (1.02-1.17)	1.43 (1.19-1.73)
	P	NS	NS	NS	NS	NS	NS
Study size	< 500	-	1.21 (0.84-1.75)	1.25 (0.93-1.69)	1.29 (0.95-1.76)	1.06 (0.74-1.52)	1.36 (1.05-1.75)
	500-	0.55 (0.37-0.80)	1.32 (0.99-1.76)	1.04 (0.92-1.18)	1.24 (0.77-2.00)	1.62 (1.14-2.29)	1.70 (1.34-2.16)
	5000-	0.64 (0.50-0.81)	1.71 (1.18-2.49)	1.48 (1.10-2.00)	0.83 (0.70-0.98)	1.00 (0.83-1.22)	1.27 (0.96-1.68)
	50000+	0.81 (0.73-0.89)	0.97 (0.88-1.07)	0.97 (0.73-1.28)	1.10 (0.79-1.52)	1.09 (0.99-1.19)	0.94 (0.82-1.07)
	P	< 0.05	< 0.01	NS	< 0.05	NS	< 0.001

The table shows effect estimates with 95% confidence interval. Variation between levels of a factor is coded as P < 0.001, P < 0.01, P < 0.05 or NS (P ≥ 0.05). C: Adjusted for comorbidities; CP: Coronavirus disease 2019 positive; D: Adjusted for demographics only; H: Hospitalized; HR: Hazard ratio; IC: Intensive care; M: Mortality; MV: Mechanical ventilation; OR: Odds ratio; PC: Adjusted for variables including post-coronavirus disease 2019 responses; RR: Relative risk; S: Severe coronavirus disease 2019; T: Tested for coronavirus disease 2019; U: Unadjusted.

For level of adjustment, the results echo those summarized above, with adjustment for comorbidities eliminating unadjusted associations of both current and ever smoking with hospitalization within those COVID-19 positive, and with ICU admission and death within those hospitalized for COVID-19.

Significant variation by type of effect estimate was only seen in two of the 12 analyses, and where it was seen may reflect the fact that the unadjusted estimates were typically ORs.

There is no convincing evidence that effect estimates vary by location.

Large studies, involving 50000 or more individuals, showed no significant increases with smoking in any of the analyses shown. Again, higher effect estimates seen in smaller studies may reflect a greater tendency for such studies to report unadjusted results.

Although meta-analyses were attempted by sex, none of them included sex-specific results from more than two studies, and the results are not shown in Table 8. There were even fewer studies reporting results by age, or by sex and age jointly, so meta-analyses by these factors were not attempted.

Dose-response results

Only one study, HIPPIS, reported results by amount smoked in current smokers, and only one, TOOLKI, by quit duration in former smokers. Neither study showed a significant dose-response relationship for any of the endpoint/population combinations considered (results not shown).

DISCUSSION

Table 9 summarizes the results from the meta-analyses of the 10 endpoint/population combinations shown in Tables 1 to 6. Results for the best estimates are shown, where estimates adjusted for comorbidities (and not responses to COVID-19) were preferred, with those adjusted for demographics preferred to unadjusted estimates, and those adjusted for factors including responses to the infection being least preferred. Results are also presented for comorbidity-adjusted estimates and for unadjusted estimates. The results show no consistent evidence of publication bias.

Table 9 Summary of results from meta-analyses¹.

Endpoint	Population	Smoking index = closest to current smoking			Smoking index = closest to ever smoking
		Best estimate2	Comorbidity adjusted3	Unadjusted4	Best estimate2	Comorbidity adjusted3	Unadjusted4
Positive	General	0.78 (0.57-1.06)5	0.96 (0.42-2.17)	0.56 (0.42-0.75)5	0.89 (0.73-1.09)	1.10 (0.70-1.72)	0.65 (0.55-0.77)
Positive	Tested	0.55 (0.47-0.65)	0.58 (0.43-0.80)	0.55 (0.45-0.67)	0.70 (0.63-0.77)	0.77 (0.64-0.92)	0.67 (0.58-0.78)
Hospitalized	Positive	1.31 (1.12-1.54)	0.96 (0.79-1.18)	1.44 (1.03-2.01)	1.38 (1.17-1.63)	0.93 (0.80-1.08)	1.71 (1.24-2.36)
ICU admission	Positive	1.13 (0.82-1.54)	1.14 (0.54-2.41)	0.98 (0.60-1.60)	1.40 (1.09-1.80)5	1.30 (0.74-2.31)	1.42 (0.94-2.14)
ICU admission	Hospitalized	1.02 (0.87-1.21)	0.81 (0.61-1.07)	1.16 (1.05-1.28)	1.11 (0.98-1.25)	1.01 (0.88-1.16)	1.19 (1.08-1.30)
Mechanically ventilated	Hospitalized	1.08 (0.86-1.36)	0.99 (0.77-1.29)5	1.16 (0.80-1.67)	1.11 (0.92-1.35)	1.04 (0.79-1.36)5	1.24 (0.99-1.54)
Severe	Positive	1.10 (0.97-1.25)5	1.02 (0.88-1.18)	1.35 (0.95-1.92)	1.10 (0.97-1.25)5	1.02 (0.88-1.18)	1.30 (0.94-1.79)
Severe	Hospitalized	1.08 (0.94-1.25)	1.23 (0.46-3.27)	1.17 (0.99-1.37)	1.09 (0.98-1.22)	1.35 (0.66-2.75)	1.14 (1.06-1.23)
Died	Positive	1.08 (0.82-1.42)	0.96 (0.92-1.01)	1.59 (1.40-1.80)	1.34 (0.96-1.86)	0.97 (0.93-1.02)	2.07 (1.12-3.82)
Died	Hospitalized	1.27 (1.08-1.49)5	1.09 (0.84-1.40)	1.52 (1.22-1.91)	1.42 (1.22-1.65)	1.12 (0.91-1.39)	1.79 (1.47-2.19)

¹See also Tables 1 to 6 for further details. Data shown are random-effects estimates with 95% confidence interval.

²Best estimates based on, respectively, 7, 17, 22, 8, 14, 15, 8, 10, 8 and 25 effect estimates (total 134).

³Comorbidity adjusted estimates based on, respectively, 3, 5, 8, 3, 1, 7, 1, 2, 4 and 10 (total 44).

⁴Unadjusted estimates based on, respectively, 3, 11, 11, 4, 12, 7, 5, 6, 3 and 13 (total 75). The unadjusted results shown are based only on those included in the best estimates;

⁵Significant publication bias (P < 0.05). ICU: Intensive care unit.

The clearest result shows that, of those tested for COVID-19, smoking was associated with a reduced risk of positivity (Figures 2 and 3), with less clear evidence of a negative association between smoking and positivity seen among the general population.

Figure 2 Forest plot of the relationship of the index “closest to current smoking” to positivity for coronavirus disease 2019 among those tested. The individual study effect estimates and random-effects combined effect estimates are shown, both as numbers and as horizontal lines with the weight of the estimate indicated by a box of proportional size. “Best estimates” are selected in the preference order (first to last) comorbidity-adjusted, demographic-adjusted and unadjusted. Estimates from ESSVRD, GDEM and HOPKIN are for current vs non-smokers, from FONTAN and OREILL for smokers (undefined) vs non-smokers, and from MUNOZP and ROZENF for ever vs never smokers, with estimates from the other eight studies being for current vs never smokers.

Figure 3 Forest plot of the relationship of the index “closest to ever smoking” to positivity for coronavirus disease 2019 among those tested. This figure is laid out as described for Figure 2. Estimates from ESSVRD, GDEM and HOPKIN are for current vs non-smokers, from FONTAN and OREILL for smokers (undefined) vs non-smokers, with estimates from the other 10 studies being for ever vs never smokers.

In contrast, all the best estimates for the other eight endpoint/population combinations, each of which relate to adverse events in those positive for, or hospitalized with COVID-19, were greater than 1 (ranging from 1.02 to 1.42), with five of the 16 estimates statistically significant (at P < 0.05). However, there was a clear difference between estimates unadjusted for other risk factors, where nine of the 16 estimates were significant, and those adjusted for comorbidities where none were significant and six were below 1.0. This difference is strikingly seen for the most commonly considered endpoint/population combination - died among those hospitalized - where (Figures 4 and 5) unadjusted estimates of 1.52 (95%CI: 1.22-1.91) for the smoking index closest to current smoking and 1.79 (1.47-2.19) for that closest to ever smoking can be contrasted with comorbidity adjusted estimates of, respectively, 1.09 (0.84-1.40) and 1.12 (0.91-1.39).

Figure 4 Forest plot of the relationship of the index “closest to current smoking” to death among those hospitalized for coronavirus disease 2019. This figure is laid out as described for Figure 2. Estimates from CARRIL are for current vs non-smokers, from APEA, HGUGM, PELLAU, SOARES and TORRES for smokers (undefined) vs non-smokers, and from CHAND, FRIED, GIACOM, IMAM, MUNOZP, SHAH, SINAI, THONE and ZHAO for ever vs never smokers, with estimates from the other nine studies being for current vs never smokers.

Figure 5 Forest plot of the relationship of the index “closest to ever smoking” to death among those hospitalized for coronavirus disease 2019. This figure is laid out as described for Figure 2. Estimates from CARRIL are for current vs non-smokers, from APEA, HGUGM, PELLAU, SOARES and TORRES for smokers (undefined) vs non-smokers, and from PETRIL for former vs never smokers, with estimates from the other 17 studies being for ever vs never smokers.

A major limitation of the available data is that, of the 73 studies which provided effect estimates, adjusted results were available for only 42. In most epidemiological contexts, effect estimates adjusted for age and sex are a basic starting point for analysis, but this was not so here. Given the different age distribution of current, former and never smokers and the strong age relationship to severe COVID-19 and death, unadjusted estimates would seem likely to be biased, as would the analyses which adjusted for variables representing a response to the virus. The most useful analyses were based on estimates adjusted for demographics only, or those adjusted also for comorbidities. These answer different questions. Analyses adjusted for demographics and comorbidities attempted to answer the question “Is a smoker more at risk of COVID-19 related outcome (such as hospitalization, admission to intensive care, undergoing mechanical ventilation or death) than an otherwise equivalent never smoker of the same age, sex and other relevant demographics and health status pre-pandemic?” This is a valid question and is somewhat equivalent to that investigated in a cohort study where smoking, demographics and health status were recorded at baseline, and smoking was related to an outcome occurring during follow-up. In analyses adjusting for demographics only, any increased risk of COVID-19 related outcomes in smokers may be due to their poorer status of health pre-pandemic. While clear answers to both questions would be nice to have, it must be noted that there are very few studies providing effect estimates adjusted only for demographics. Thus, considering deaths in the hospitalized population, only three out of 25 studies with relevant data provided estimates adjusted for demographics only, and none provided comparable unadjusted, demographic adjusted and comorbidity adjusted effect estimates.

Another problem is that some studies provided ORs, some RRs and some HRs. Where results for the relevant 2 × 2 table on exposure × outcome were available, and the ORs and RRs were both estimable, we generally used the OR, using the RR only where the source paper had reported adjusted RRs or HRs. Although we could have used an alternative strategy, it is doubtful whether this would have materially affected our results, given the general consistency of the results by type of effect estimate.

Another possible concern is with the mortality data. Much relates to deaths occurring in patients hospitalized with COVID-19, and many publications implicitly assume all those deaths were due to the virus, when some might have been due to other causes. However, given this proportion is small, it seems probable that this would only result in a minor bias.

More concern relates to the quality of the smoking data. There are two issues here. One is the way smoking was recorded and defined, with eight studies reporting results for smoking undefined, and only 36 distinguishing current and former smoking. Also, only one study reported results by amount smoked, and only one reported results by duration of quit. Furthermore, papers generally did not provide details on the smoking questions asked, or when they were asked.

The other issue is that many studies derived their smoking data from medical records, known to be incomplete and inaccurate[37-39]. While many studies gave no information concerning missing data on smoking, 35 did so, and in 13 the proportion with missing data exceeded 20%, giving concern about the validity of their effect estimates. It would not surprise us to find that, in some studies that did not mention missing data, the “non-smokers” included some individuals never actually asked about their smoking.

Recent publications, particularly by Farsalinos et al[2-4,6-8], have observed that the prevalence of smoking seen in the studies of hospitalized patients was substantially less than reported in national statistics by a factor of four or so, and have suggested that smokers might be protected against getting COVID-19. While the mean percentage of current and former smokers in the studies of hospitalized patients that we considered (current 7.76%, SE 1.00%; ever 33.24%, SE 2.01%) was clearly less than in the studies of the general population (current 15.14%, SE 1.34%; ever 44.0%, SE 1.75%), the difference was only by a factor of 2 or 1.3 rather than about 4. While we also showed a reduced risk in smokers of COVID-19 positivity in those tested (Table 1), the reduction was again much less than the factor reported by Farsalinos et al[6]. Although it is possible that some of the difference between our results and those of Farsalinos et al[6] has arisen as we excluded Asian studies, while their results mainly came from Asia, the fact that two out of the three studies on the general population that we considered found no reduced risk of hospitalization in smokers (SupplementaryTable 3) suggests to us that the low prevalence of smoking seen in hospitalized patients may largely result from incompleteness of data in hospital records, although it would also be consistent with smokers with COVID-19 tending to be less likely to appear for testing or report to hospital.

As noted in the introduction, there were, by the end of September 2020 (the final date of the searches used to produce these results), quite a number of published meta-analyses which relate smoking to adverse events such as deaths or severity of COVID-19[3,6,9-28]. These meta-analyses had various limitations, including little attention to possible data inadequacy, limiting attention to studies of hospitalized patients and considering few, if any, studies conducted outside Asia. They also include paying scant attention to the need for adjusting effect estimates for other risk factors. Many of the studies and meta-analyses we considered demonstrated that dying from COVID-19, for example, was strongly related to various factors associated with smoking, including age, obesity, and a history of respiratory, cardiovascular and other diseases, and yet they attempt to draw conclusions for smoking from unadjusted analyses.

By now, we are aware of further meta-analyses that have been conducted, including those that concentrate on smoking[11,19,29,40,41] and those that consider smoking as one of a list of factors considered[16,20,21,42-55]. While some of the reviews considered far more studies than the earlier reviews, including those in non-Asian populations, the weaknesses seen generally persist. Thus, for example, a recent review of 109 studies[41] limited attention to hospitalized patients, considered only unadjusted effect estimates, hardly mentioned lack of adjustment as a possible weakness, and paid very limited attention to the possibility that smoking data from hospital records may be inadequate.

That meta-analysis[41] described the 109 studies included as being all of moderate or high quality and some of the other reviews also attempted to evaluate study quality. We did not attempt classification of study quality, but given so many studies used medical records as the source of smoking data and failed to present results adjusted even for basic demographics, we doubt very much that we would have considered more than a few studies to be of high or even moderate quality. This view aligns with that of a recent meta-analysis[29] which rejected 201 of 233 studies as being of poor quality, with only one of the studies considered in their meta-analyses considered to be of good quality and the rest classified as fair.

Our meta-analyses have various strengths, including giving careful attention to adjustment, considering many combinations of outcome and population at-risk, and including meta-analyses comparing estimates by factors such as location, study size and type of estimate considered. Weaknesses relate mainly to the poor underlying data quality, much from medical records, and many studies failing to provide adjusted estimates. Almost complete lack of data for males and females separately, and by age group is also a limitation, as is the very limited data on amount smoked or quit duration. Other possible limitations relate to the fact that, with the exception of studies from Israel, we did not consider results from other Asian countries; thus, our conclusions may not necessarily apply to all locations. We also did not consider studies involving less than 100 cases as their results would be less reliable, and studies of patients with specific diseases as their results would not be generalizable.

However, we feel that our meta-analysis provides a good insight into the relationship between smoking and a variety of endpoints relevant to COVID-19.

CONCLUSION

Based on data from 74 studies conducted in Europe, Israel, America and Australasia, many providing only limited results, there is evidence that, among those tested for COVID-19, smokers are less likely to be positive for the virus. There is also less clear evidence of reduced positivity in smokers in the general population. Among those who are positive for, or hospitalized with, COVID-19 there is a positive association between smoking and both death and severity of COVID-19. This association is most clearly seen for effect estimates unadjusted for other risk factors, and is not evident for estimates adjusted for comorbidities and demographic variables. This suggests that any apparent adverse effect of smoking is due to the poorer prior health status of smokers and that smokers and non-smokers with equivalent demographics and prior health status have a very similar risk of adverse events linked to COVID-19.

ARTICLE HIGHLIGHTS

Research background

Previous meta-analyses relating smoking to coronavirus disease 2019 (COVID-19) are limited by considering few studies, restricting attention to hospitalized patients, giving limited or no attention to the definition of smoking or the reliability of smoking as recorded, and failing to properly consider the effect of adjustment for demographics and comorbidities.

Research motivation

We wished to gain a detailed insight into the effect of smoking on a variety of endpoints in different populations.

Research objectives

To carry out a systematic review, based on epidemiological studies in Europe, Israel, America and Australasia on the relationship of smoking to being tested for COVID-19, being positive for COVID-19, being hospitalized with COVID-19, having severe disease or dying.

Research methods

Literature searches based on publications in English up to September 30, 2020 identified studies of at least 100 individuals, carried out in Europe, Israel, America and Australasia, and unrestricted to those with specific other diseases, and providing information relating smoking to various COVID-related endpoints. Fixed-effect and random-effects meta-analyses were conducted for combinations of index of smoking, endpoint, population and level of adjustment with heterogeneity studied by level of adjustment, study location, and other factors.

Research results

Data were available from 74 studies of highly variable size: 37 in the United States, 10 in the United Kingdom, and up to four elsewhere, with populations most commonly studied being those hospitalized with COVID-19, positive for COVID-19, tested for COVID-19 and the general population. Only 36 studies distinguished current and former smokers, and adjusted results for smoking were only given in 42 studies. Positivity for COVID-19 was reduced among smokers in those tested, but not in the general population. Apparent increases in risk in smokers of hospitalization for COVID-19 among those positive, and of death among those positive and among those hospitalized disappeared following adjustment for pre-existing comorbidities, and there was little evidence of any relationship of smoking with admission to intensive care, being mechanically ventilated or having severe COVID-19, even in the unadjusted results.

Research conclusions

There is some evidence that smoking is associated with a reduced risk of being COVID-19 positive. Any apparent adverse effects of smoking on hospitalization rates among those positive, and on death rates seem due to the poorer prior health status of smokers.

Research perspectives

Evidence from later studies could consolidate these conclusions, and help to explain why, among those tested for COVID-19, current smokers are less likely to be positive.

ACKNOWLEDGEMENTS

We thank Yvonne Cooper for typing various drafts of the paper and obtaining relevant references, and John Fry for comments at various stages.