Abstract
Background The pandemic of coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has rapidly emerged to seriously threaten public health. We aimed to investigate whether white blood cell traits have potential causal effects on severe COVID-19 using Mendelian randomization (MR).
Methods To evaluate the causal associations between various white blood cell traits and severe COVID-19, we conducted a two-sample MR analysis with summary statistics from recent large genome-wide association studies.
Results Our MR results indicated potential causal associations of white blood cell count, myeloid white blood cell count, and granulocyte count with severe COVID-19, with odds ratios (OR) of 0.84 (95% CI: 0.72-0.98), 0.81 (95% CI: 0.70-0.94), and 0.84 (95% CI: 0.71-0.99), respectively. Increasing eosinophil percentage of white blood cells was associated with a higher risk of severe COVID-19 (OR: 1.22, 95% CI: 1.03-1.45).
Conclusions Our results suggest the potential causal effects of lower white blood cell count, lower myeloid white blood cell count, lower granulocyte count, and higher eosinophil percentage of white blood cells on an increased risk of severe COVID-19.
- COVID-19
- white blood cells
- eosinophil
- Mendelian randomization
Background
Coronavirus disease 2019 (COVID-19) is caused by infection of a novel virus called Severe Acute Respiratory Syndrome Coronavirus 2 (SARS⍰CoV⍰2) [1]. SARS⍰CoV⍰2 has rapidly spread, causing damage or even death [2-4]. Besides age and gender, some pre-existing conditions are also well-known to be associated with an increased risk of severe COVID-19, such as cardiovascular disease, diabetes, chronic respiratory disease, hypertension, and cancers [4-7]. Recently, genetic studies, including two genome-wide association studies (GWAS), have identified multiple genetic loci to be associated with the susceptibility and severity of COVID-19 [8-10]. However, causal risk factors for severe COVID-19 remain unclear. Identifying host factors predisposing individuals to severe COVID-19 is urgently needed to improve primary prevention and to develop treatment strategies.
Elevated white blood cell and neutrophil counts, and depleted lymphocyte count have been repeatedly observed in COVID-19 patients with severe outcomes, and neutrophil-to-lymphocyte ratio have been proposed as a prognostic biomarker [1, 11-14]. However, findings from recent observational studies are inconsistent and the exact roles of white blood cells and its various subtypes in severe COVID-19 remain elusive [15-17]. Most existing studies measured blood cell counts in patients with confirmed infection of SARS-CoV-2 and as a result, the hematological indices could have been modified by immune responses [18]. It is unknown if blood cell counts before infection are associated with the risk of developing severe COVID-19. Even if white blood cells are measured before infection, they are influenced by many exogenous and endogenous factors (e.g., age, gender, disease status, and medications), which will confound the observational associations [19, 20]. No previous research has been able to interrogate the causal role of white blood cells in severe COVID-19.
While traditional observational associations often suffer from confounding, reverse causality, and various biases, a complementary approach, Mendelian randomization (MR), utilizes genetic variants as instrumental variables to approximate the lifetime status of exposure and enable causal inference in observational data. The random allocation of the allele at conception and the natural directional effects of genetic variants empower MR estimates to be less plagued by confounding and reverse causality [21, 22]. In this study, we aimed to test the causal effects of 19 white blood cell traits on severe COVID-19 by performing a two-sample MR analysis.
Methods
Study design and data source
A two-sample MR study was conducted to examine the causal effects of 19 white blood cell traits on severe COVID-19. These blood cell traits were chosen based on a GWAS meta-analysis in 173,480 European ancestry individuals across three cohorts, and their summary statistics were available in the IEU OpenGWAS database [23, 24]. Genetic instruments for each white blood cell trait were chosen based on the following criteria: 1) p < 8.31×10−9 for association with the exposure; and 2) linkage disequilibrium (LD) clumping based on r2 > 0.001. The data source and details for 19 white blood cell traits are available in Supplementary Table 1. The instrument-outcome effects were retrieved from the largest GWAS meta-analysis of COVID-19 to date, by the COVID-19 Host Genetics Initiative (HGI, release 3, accessed on July 2, 2020) [10]. We used the summary statistics based on the comparison of hospitalized COVID-19 patients (N = 3,199) with the general population (N = 897,488).
Statistical analysis
The causal effect of a white blood cell trait on severe COVID-19 was evaluated using the inverse variance-weighted (IVW) method with a multiplicative random-effects model [22, 25, 26]. Horizontal pleiotropy occurs when SNPs exert a direct effect on the severe COVD-19 through pathways other than the hypothesized exposure. To evaluate the possible presence of horizontal pleiotropy, we calculated Cochran’s Q statistic for heterogeneity and conducted the intercept test associated with the MR-Egger method. Additional sensitivity analyses were performed with MR-Egger [22, 26], weighted median (WM) [27], and Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) test [28]. The MR-Egger estimates allowed directional or unbalanced horizontal pleiotropic effects. The weighted median method provides robust causal estimates even when up to 50 % SNPs are invalid genetic instruments [27]. The MR-PRESSO test was utilized to correct for the presence of specific horizontal pleiotropic outlier variants via detected outlier removal [28]. All MR analyses were conducted in R with the TwoSampleMR package [23].
Resources
The COVID-19 Host Genetics Initiative: https://www.covid19hg.org/
The IEU OpenGWAS database: https://gwas.mrcieu.ac.uk/
Results
The counts of white blood cell, myeloid white blood cell, and granulocyte are negatively associated with severe COVID-19
By applying a two-sample MR approach, we first investigated the causal associations of the counts of white blood cell and its subpopulations with severe COVID-19. A relatively large number of independent SNPs, ranging from 79 for basophil count to 185 for monocyte count, were selected as genetic instruments for each blood cell count (Supplementary Tables 2-11). Based on the IVW MR estimates under a multiplicative random-effects model, we identified potentially causal, negative associations of white blood cell count (OR = 0.84, CI: 0.72-0.98, p = 0.031), basophil count (OR = 0.75, CI: 0.58-0.96, p = 0.023), myeloid white blood cell count (OR = 0.81, CI: 0.70-0.94, p = 0.0070), and granulocyte count (OR = 0.84, CI: 0.71-0.99, p = 0.040) with severe COVID-19 (Fig. 1, Table 1). A suggestive negative association was also found for sum neutrophil eosinophil counts (OR = 0.85, CI: 0.73-1.00, p = 0.051). No evidence of heterogeneity in causal estimates was found by the Cochran Q statistic, and no evidence of horizontal pleiotropy was reported by the MR-Egger intercept test, except for basophil count (p = 0.036). Causal estimates from MR-Egger and WM MR revealed broadly concordant effect directions, although they are mostly not statistically significant, probably due to the reduced power of these two approaches [26] (Supplementary Fig. 1, Supplementary Table 21). MR-PRESSO analysis did not identify any outlier SNPs and yielded significant causal estimates for white blood cell count (p = 0.033), myeloid white blood cell count (p = 0.0078), and granulocyte count (p = 0.041) (Supplementary Table 21). Taken together, we demonstrated that white blood cell count, myeloid white blood cell count, and granulocyte count had consistent, negative effects on the risk of severe COVID-19.
The percentage of eosinophil in white blood cell is positively associated with severe COVID-19
We further investigated the causal associations of the percentages of specific white blood cells with severe COVID-19. As genetic instruments, we used 158 SNPs for eosinophil percentage, 63 SNPs for basophil percentage, 135 SNPs for the neutrophil percentage, 191 SNPs for monocyte percentage, and 135 SNPs for lymphocyte percentage (Supplementary Tables 12-16).
Genetically predicted higher eosinophil percentage of white blood cells (OR = 1.22, CI: 1.03-1.45, p = 0.023) was associated with increased risk of severe COVID-19 using IVW with the random-effect model (Figure 1). The WM MR method and the MR-PRESSO analysis both revealed a consistent, positive effect (OR = 1.41, CI: 1.07-1.87, p = 0.015; and OR = 1.22, CI: 1.03-1.45, p = 0.024; respectively. Supplementary Table 21). No pleiotropy or outlier SNPs were identified in the MR-Egger test and MR-PRESSO analysis. Cochran Q statistics indicated no heterogeneity among the genetic instruments (Supplementary Table 21).
The percentage of eosinophil in granulocyte is positively associated with severe COVID-19
We further focused on the granulocyte and evaluated if it specific compositions are associated with severe COVID-19. As genetic instruments, there were 168 SNPs for granulocyte percentage of myeloid white blood cells, 150 SNPs for eosinophil percentage of granulocytes, 59 SNPs for basophil percentage of granulocytes, and 141 SNPs for the neutrophil percentage of granulocytes (Supplementary Tables 17-20). Our results found suggestive evidence for a risk-increasing effect of eosinophil percentage of granulocytes on severe COVID-19 (OR = 1.18, CI: 1.00-1.39, p = 0.053). The WM MR method and the MR-PRESSO analysis both revealed a consistent positive effect (OR = 1.35, CI: 1.03-1.75, p = 0.029; and OR = 1.18, CI: 1.00-1.39, p = 0.055; respectively). No evidence of heterogeneity and pleiotropy was found. No statistically significant causal associations were identified for other specific granulocyte percentages (Figure 1 and Supplementary Table 21).
Discussion
To our knowledge, this is the first MR study evaluating the causal roles of white blood cell traits in severe COVID-19 risk. Overall, our results suggest potential causal protective effects of increasing white blood cell count, myeloid white blood cell count, and granulocyte count on severe COVID-19. Our novel findings also include that genetically predicted higher eosinophil percentage in white blood cells or in granulocytes was associated with a higher risk of severe COVID-19.
Previous observational studies have frequently pointed out increased white blood cell count in the severe COVID-19 patients, when compared to healthy controls or mild COVID-19 patients [1, 11, 12, 29]. However, there are also reports that normal or decreased white blood cell count is more common in COVID-19 patients when compared to the reference range or healty controls [10, 15-17, 30-32]. Our MR analysis showed that lower white blood cell count, myeloid white blood cell count, and granulocyte count may play a causal role in increasing the risk of severe COVID-19. The mechanism by which they contribute to severe COVID-19 remains unclear. Immune system disorders have been suspected of playing roles in severe COVID-19 risk [33, 34]. The complete elucidation of the potential mechanism warrants further investigation.
Persistent eosinopenia after admission was associated with COVID-19 risk in previous retrospective and prospective observational studies [12, 15, 31, 32, 35-37]. There are also reports that eosinophil remains stable in severe COVID-19 patients [13, 29, 38, 39]. A multiparametric flow cytometry analysis found that high eosinophil count was associated with an increased risk of severe COVID-19 [13]. Another observational study suggested that increasing eosinophil count might be an indicator of COVID-19 improvement [40]. Our MR analysis supported that high eosinophil percentage of white blood cells may be causal in increasing the risk of severe COVID-19 [40, 41], calling for future mechanistic studies.
Lymphopenia, as a response to viral infection, has been frequently associated with severe COVID-19 risk [1, 2, 13, 29, 41-51]. However, we did not detect a causal effect of lymphocyte count on severe COVID-19. This discrepancy may reflect reverse causality in retrospective and prospective observational studies, with depleted lymphocyte count as a result of immune response to SARS⍰CoV⍰2 infection [52]. Due to the limited number of SNPs associated with severe COVID-19, performing reverse MR analyses between severe COVID-19 and white blood cell traits is challenging. Using linkage disequilibrium (LD) assessment, only one SNP at locus 3p21.31 was retained in the previously published genome-wide significant SNPs [10]. We recently showed that SNPs at this locus are associated with multiple blood cell traits, suggesting they may have pleiotropic effects and are not suitable to be used as genetic instruments [53]. As more COVID-19-associated SNPs are identified in the future, reverse MR analysis will be valuable to understand the effect of COVID-19 on blood cell traits.
One assumption of MR is that the instrumental variable influences severe COVID-19 risk only through its effect on a specific white blood cell trait. The Cochran Q statistic did not reveal heterogeneity among our genetic instruments, and the MR-Egger intercept test also indicated no presence of pleiotropic effects except those for basophil count. Our MR study was performed with strong instrumental variables and adequate statistical power, and we have conducted extensive sensitivity analyses. Still, we emphasize that our results should be interpreted with caution and future studies are needed to elucidate the mechanistic roles of white blood cell traits in severe COVID-19.
Conclusions
Our results suggest that lower white blood cell count, lower myeloid white blood cell count, lower granulocyte count, and higher eosinophil percentage of white blood cells are causally associated with an increased risk of severe COVID-19.
Data Availability
All data generated or analyzed during this study are included in this published article and this supplementary information files.
Declarations
Ethics approval and consent to participate
Ethical approval from the North West Multi-Centre Research Ethics Committee and written informed consent from all participants were obtained.
Consent for publication
Not applicable.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and this supplementary information files.
Competing interests
The authors declare that they have no competing interests.
Funding
KY is supported by the University of Georgia Research Foundation. Funding sources had no involvement in the conception, design, analysis, or presentation of this work.
Authors’ contributions
YS, JZ, and KY conceived the study. YS performed data analysis and prepared visualizations. YS, JZ, and KY interpreted the results. YS and KY wrote the first draft of the manuscript. All authors critically read, revised, and approved the final version of the manuscript.
Supplemental Data
Supplementary Table 1. Summary of GWAS information for white blood cell phenotypes
Supplementary Table 2. 167 SNPs used IVs in MR analyses testing effect of white blood cell count on severe COVID-19 risk.
Supplementary Table 3. 169 SNPs used IVs in MR analyses testing effect of eosinophil count on severe COVID-19 risk.
Supplementary Table 4. 79 SNPs used IVs in MR analyses testing effect of basophil count on severe COVID-19 risk.
Supplementary Table 5. 140 SNPs used IVs in MR analyses testing effect of neutrophil count on severe COVID-19 risk.
Supplementary Table 6. 185 SNPs used IVs in MR analyses testing effect of monocyte count on severe COVID-19 risk.
Supplementary Table 7. 160 SNPs used IVs in MR analyses testing effect of lymphocyte count on severe COVID-19 risk.
Supplementary Table 8. 146 SNPs used IVs in MR analyses testing effect of myeloid white blood cell count on severe COVID-19 risk.
Supplementary Table 9. 144 SNPs used IVs in MR analyses testing effect of granulocyte count on severe COVID-19 risk.
Supplementary Table 10. 139 SNPs used IVs in MR analyses testing effect of sum neutrophil eosinophil counts on severe COVID-19 risk.
Supplementary Table 11. 142 SNPs used IVs in MR analyses testing effect of sum basophil neutrophil counts on severe COVID-19 risk.
Supplementary Table 12. 158 SNPs used IVs in MR analyses testing effect of eosinophil percentage of white blood cells on severe COVID-19 risk.
Supplementary Table 13. 63 SNPs used IVs in MR analyses testing effect of basophil percentage of white blood cells on severe COVID-19 risk.
Supplementary Table 14. 135 SNPs used IVs in MR analyses testing effect of neutrophil percentage of white blood cells on severe COVID-19 risk.
Supplementary Table 15. 191 SNPs used IVs in MR analyses testing effect of monocyte percentage of white blood cells on severe COVID-19 risk.
Supplementary Table 16. 135 SNPs used IVs in MR analyses testing effect of lymphocyte percentage of white blood cells on severe COVID-19 risk.
Supplementary Table 17. 168 SNPs used IVs in MR analyses testing effect of granulocyte percentage of myeloid white blood cells on severe COVID-19 risk.
Supplementary Table 18. 150 SNPs used IVs in MR analyses testing effect of eosinophil percentage of granulocytes on severe COVID-19 risk.
Supplementary Table 19. 59 SNPs used IVs in MR analyses testing effect of basophil percentage of granulocytes on severe COVID-19 risk.
Supplementary Table 20. 141 SNPs used IVs in MR analyses testing effect of neutrophil percentage of granulocytes on severe COVID-19 risk.
Supplementary Table 21. MR estimates from each method of assessing the causal effects of white blood cells on severe COVID-19 risk.
Supplementary Figure 1. Scatter plots of the genetic associations of white blood cell count (A), or myeloid white blood cell count (B), or granulocyte count (C), or basophil count (D) or eosinophil percentage in white blood cells (E) associated SNPs against the genetic associations of severe COVID-19. The slopes of each line represent the causal association using different MR methods. The light blue line represents the random-effects inverse variance-weighted estimate, the dark blue line represents the MR-Egger estimate, and the green line represents the weighted median estimate.
Acknowledgments
We would like to express our gratitude to all other Ye lab members for stimulating discussions.
List of abbreviations
- COVID-19
- coronavirus disease 2019
- SARS-CoV-2
- severe acute respiratory syndrome coronavirus 2
- MR
- Mendelian randomization
- GWAS
- genome-wide association studies
- SNPs
- single nucleotide polymorphisms
- LD
- linkage disequilibrium
- HGI
- Host Genetics Initiative
- IEU
- Integrative Epidemiology Unit
- IVW
- inverse variance-weighted
- WM
- weighted median
- MR-PRESSO
- Mendelian randomization pleiotropy residual sum and outlier