Next Article in Journal
GBP5 Serves as a Potential Marker to Predict a Favorable Response in Triple-Negative Breast Cancer Patients Receiving a Taxane-Based Chemotherapy
Next Article in Special Issue
Alterations of Left Ventricular Function Persisting during Post-Acute COVID-19 in Subjects without Previously Diagnosed Cardiovascular Pathology
Previous Article in Journal
A Quantitative Framework for Measuring Personalized Medicine Integration into US Healthcare Delivery Organizations
Previous Article in Special Issue
Statin Therapy and the Risk of COVID-19: A Cohort Study of the National Health Insurance Service in South Korea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

White Blood Cells and Severe COVID-19: A Mendelian Randomization Study

1
Department of Genetics, Franklin College of Arts and Sciences, University of Georgia, Athens, GA 30602, USA
2
School of Public Health, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
3
Institute of Bioinformatics, University of Georgia, Athens, GA 30602, USA
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2021, 11(3), 195; https://doi.org/10.3390/jpm11030195
Submission received: 22 February 2021 / Revised: 1 March 2021 / Accepted: 10 March 2021 / Published: 12 March 2021
(This article belongs to the Special Issue COVID-19 Related Complications)

Abstract

:
Increasing evidence shows that white blood cells are associated with the risk of coronavirus disease 2019 (COVID-19), but the direction and causality of this association are not clear. To evaluate the causal associations between various white blood cell traits and the COVID-19 susceptibility and severity, we conducted two-sample bidirectional Mendelian Randomization (MR) analyses with summary statistics from the largest and most recent genome-wide association studies. Our MR results indicated causal protective effects of higher basophil count, basophil percentage of white blood cells, and myeloid white blood cell count on severe COVID-19, with odds ratios (OR) per standard deviation increment of 0.75 (95% CI: 0.60–0.95), 0.70 (95% CI: 0.54–0.92), and 0.85 (95% CI: 0.73–0.98), respectively. Neither COVID-19 severity nor susceptibility was associated with white blood cell traits in our reverse MR results. Genetically predicted high basophil count, basophil percentage of white blood cells, and myeloid white blood cell count are associated with a lower risk of developing severe COVID-19. Individuals with a lower genetic capacity for basophils are likely at risk, while enhancing the production of basophils may be an effective therapeutic strategy.

1. Introduction

Coronavirus disease 2019 (COVID-19) is caused by infection of a novel virus called Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) that has rapidly spread, causing irreversible damage or even death [1,2,3,4]. A growing number of studies have found that immune cells may play important roles in the COVID-19 severity and susceptibility [5,6,7,8,9]. Since prior findings were based on observational studies, the causal associations between white blood cells (WBC) and COVID-19 remain unclear. Identifying host factors predisposing individuals to COVID-19 susceptibility and severity is urgently needed to improve primary prevention and develop treatment strategies.
Many observational studies have reported risk factors for COVID-19 susceptibility and severity. 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,10,11,12]. Furthermore, elevated WBC and neutrophil counts, and depleted lymphocyte count have been repeatedly observed in COVID-19 patients with severe outcomes, and the neutrophil-to-lymphocyte ratio has been proposed as prognostic biomarkers [1,5,13,14,15,16]. However, findings from recent observational studies are inconsistent, and it is difficult to infer the causal roles of WBC in severe COVID-19 due to confounding and reverse causation [17,18,19]. 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 [20]. It is unknown if blood cell counts before infection are associated with the risk of developing severe COVID-19. Even if WBC 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 [21,22]. To date, no existing research has been able to interrogate the causal role of WBC in COVID-19 severity and susceptibility.
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 alleles at conception and the natural directional effects of genetic variants on phenotypes empower MR estimates to be less plagued by confounding and reverse causality [23,24]. In this study, we used two-sample MR analyses to examine potential causal associations between various WBC traits and the risk of COVID-19 severity and susceptibility. Recent genome-wide association studies (GWAS) with thousands of cases and millions of controls have identified scores of genetic loci associated with the susceptibility and severity of COVID-19 [25]. Two recent GWAS on WBC traits enrolled approximately half a million European-ancestry participants and dramatically increased the number of instrumental variables [26,27,28]. Building on these largest-to-date GWAS for both the exposure and outcome, we conducted two-sample bidirectional MR analyses to disentangle the causal effects of various WBC traits on COVID-19 severity and susceptibility, or vice versa.

2. Materials and Methods

2.1. Study Design and Data Sources

A two-sample bidirectional MR study was conducted to examine the causal effects of various WBC traits on COVID-19, and vice versa (Figure 1). GWAS summary statistics of WBC traits were extracted from publicly available data (https://www.ebi.ac.uk/gwas/downloads/summary-statistics, accessed on 19 January 2021) [29]. Summary statistics were obtained from the most recent version of GWAS from the COVID-19 Host Genetics Initiative (https://www.covid19hg.org/ accessed on 20 January 2021) [25]. Publicly available data sources had previously received appropriate ethics and institutional review board approvals, and thus further ethical approvals were not required.

2.2. Instrumental Variables for WBC and COVID-19

In our primary discovery analysis, instrumental variables for WBC were selected based on two recent GWAS as independent genetic variants at the study-specific genome-wide significance (p < 5 × 10−9) [26,27]. The first is a GWAS meta-analysis of 11 WBC traits, which included 408,112 participants of European ancestry and was released by Vuckovic et al. [26]. The other is a GWAS of 6 WBC traits, which included 562,132 European-ancestry participants and was derived by Chen et al. [27]. In the replication analysis, we selected SNPs that were associated with 21 white blood cell traits at the study-specific genome-wide significance (p < 8.31 × 10−9) in a GWAS meta-analysis comprising 173,480 individuals of European ancestry. This GWAS is called the Astle et al. study [28].
The instrumental variables for COVID-19 were retrieved at the genome-wide significance (p < 5 × 10−8) from the largest GWAS meta-analysis of COVID-19 to date, by the COVID-19 Host Genetics Initiative (HGI, round 5, released on 18 January 2021) [25]. For COVID-19 severity, we used two GWAS with European-ancestry participants. The first, called study A2, compared patients confirmed as “very severe respiratory” COVID-19 (N = 5105) with the general population samples (N = 1,383,241). The other, called study B2, compared hospitalized COVID-19 patients (N = 9986) with the general population (N = 1,877,672). For COVID-19 susceptibility, summary statistics were obtained from a GWAS in participants of European ancestry, including patients with laboratory-confirmed COVID-19 infection or physician confirmed COVID-19 (N = 36,590) and the population controls (N = 1,668,938). This is called study C2.
To ensure independence among instrumental variables for each exposure, we applied linkage disequilibrium (LD) clumping based on r2 > 0.001 and discarded variants within 1-Mb distance from other SNPs with a stronger association. Effect alleles of genetic instruments were harmonized across the GWAS of the exposure and the outcome, and those not present in the GWAS of the outcome were removed. F statistics were calculated to assess instrument strength and F ≥ 10 indicates strong instruments [30].

2.3. Statistical Analysis

The causal association of each exposure with each outcome was evaluated using the inverse variance-weighted (IVW) method with a multiplicative random-effects model [24,31,32]. Horizontal pleiotropy occurs when SNPs exert a direct effect on the severe COVID-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 [24,32], weighted median (WM) [33], weighted mode methods [34], and Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) test [35]. The MR-Egger estimates allow directional horizontal pleiotropic effects. The weighted median method provides robust causal estimates even when up to 50% of SNPs are invalid genetic instruments [33]. The weighted mode method requires that the largest number of instruments that identify the consistent causal effect derived from valid instruments even when the majority of instruments are invalid [34]. The MR-PRESSO test was utilized to correct for the presence of specific horizontal pleiotropic outlier variants via detected outlier removal [35]. The false discovery rate (FDR) approach was adopted to correct for multiple testing [36]. The nominal p-value < 0.05 but above the FDR significance threshold (q-value < 0.05) was considered as suggestive evidence. All MR analyses were conducted in R with the TwoSampleMR package (version 3.6.7) [37].

3. Results

3.1. Analysis Pipeline

We first aimed to investigate the causal effects of WBC traits on COVID-19 severity and susceptibility via the forward MR analysis. In our main analysis, we selected uncorrelated SNPs as instrumental variables for each WBC trait from two recent GWAS in 408,112 and 562,132 participants of European ancestry [26,27]. We then performed a replication analysis using genetic instruments from a previous GWAS meta-analysis that included 173,480 individuals of European ancestry [28]. On the other hand, to evaluate the causal effects of COVID-19 on WBC traits, we conducted a reverse MR analysis using COVID-19-associated SNPs as instrumental variables. All SNPs for the WBC traits (F statistic > 77.82) and COVID-19 (F statistic > 3458.84) are strong instruments. A schematic of the analysis pipeline is shown in Figure 1. Details about GWAS, from which we extracted summary-level data, are presented in Supplementary Table S1.

3.2. Forward MR Analysis of the Effects of WBC Traits on COVID-19

By applying a two-sample MR approach, we first investigated the causal associations of 12 WBC traits with COVID-19. A relatively large number of independent SNPs, ranging from 115 for basophil percentage of white blood cells to 469 for monocyte count, were selected as genetic instruments for each WBC trait (Supplementary Tables S2–S18).
Using very severe respiratory COVID-19 as the outcome (A2), we identified suggestive causal, negative associations of basophil count (OR = 0.75, CI: 0.60–0.95, p = 0.015), myeloid WBC count (OR = 0.85, CI: 0.73–0.98, p = 0.022), and basophil percentage of WBC (OR = 0.70, CI: 0.54–0.92, p = 0.011) based on the IVW MR estimates under a multiplicative random-effects model (Figure 2) (Supplementary Table S22). Causal estimates from MR-Egger, WM, and weighted mode methods revealed broadly concordant effect directions, although they are mostly not statistically significant, probably due to the reduced power of these two approaches [32]. MR-PRESSO analysis did not identify any outlier SNPs and yielded significant causal estimates for basophil count (p = 0.015), myeloid WBC count (p = 0.0022), and basophil percentage of WBC (p = 0.011). 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. In the replication analysis using genetic instruments from another GWAS meta-analysis of WBC traits (i.e., the Astle et al. study), we similarly observed the significant negative effects of basophil count and myeloid WBC count. Additionally, nominally significant negative associations were observed for total WBC count (p = 0.038), neutrophil count (p = 0.016), granulocyte count (p = 0.028), sum neutrophil eosinophil counts (p = 0.032), sum basophil neutrophil counts (p = 0.012) (Supplementary Table S23).
Using hospitalized COVID-19 as the outcome (B2), we also found negative associations of basophil count (OR = 0.83, CI: 0.71–0.97, p = 0.020), basophil percentage of WBC (OR = 0.78, CI: 0.65–0.93, p = 0.005), and myeloid WBC count (OR = 0.90, CI: 0.73–0.996, p = 0.041) (Figure 2) (Supplementary Table S22). The causal effect of basophil count (q = 0.014) and basophil percentage (q = 0.042) surpassed the FDR correction in this analysis. Moreover, genetically predicted lower total WBC count (OR = 0.89, CI: 0.80–0.98, p = 0.016) was associated with increased risk of hospitalized COVID-19 using IVW with the random-effect model (Figure 2) (Supplementary Table S22). MR-Egger, WM, weighted mode methods, and MR-PRESSO analysis revealed broadly consistent effect directions. No pleiotropy was identified in the MR-Egger test and a potential causal effect of WBC count was retained after removing outlier SNPs identified by MR-PRESSO analysis. In the replication analysis, basophil count, basophil percentage, and total WBC count were also negatively associated with hospitalized COVID-19 (Supplementary Table S23). Taken together these two outcomes of COVID-19 severity, very severe respiratory and hospitalized COVID-19 (A2 and B2), we demonstrated that basophil count, basophil percentage, myeloid WBC count, and total WBC count had consistent, negative effects on the risk of severe COVID-19.
We further investigated the causal associations of WBC traits with COVID-19 susceptibility (C2). None of the WBC traits was causally associated with COVID-19 susceptibility (Figure 2) (Supplementary Table S22). Similar findings were observed in the replication analysis (Supplementary Table S23). Overall, there was no evidence of associations between any WBC trait and COVID-19 susceptibility.

3.3. Reverse MR Analysis of the Effects of COVID-19 on WBC Traits

We further explored the causal effects of COVID-19 severity and susceptibility on each WBC trait in the reverse direction using COVID-19-associated SNPs as instrumental variables. We used nine SNPs (A2) and five SNPs (B2) for COVID-19 severity, and seven SNPs (C2) for COVID-19 susceptibility (Supplementary Table S19–S21). The reverse MR analysis showed no significant association of genetically predicted COVID-19 severity or susceptibility with any WBC trait (Figure 3) (Supplementary Table S24). The null findings were confirmed in the replication analysis that used summary statistics from the Astle et al. study of WBC traits (Supplementary Table S25).

4. Discussion

To our knowledge, this is the first MR study evaluating the causal roles of WBC traits in COVID-19 severity and susceptibility. Overall, our results suggest potential causal protective effects of higher basophil count, basophil percentage of WBC, and myeloid WBC count on the severity of COVID-19. No evidence was found for causal effects of WBC traits on the susceptibility of COVID-19, or for reverse causation of COVID-19 on WBC traits.
Previous observational studies have frequently shown increased WBC count in the severe COVID-19 patients, when compared to healthy controls or mild COVID-19 patients [1,5,13,38]. However, observational studies may be influenced by confounding and cannot always distinguish symptoms from causes. Using genetically predicted WBC traits, our MR analysis showed that lower basophil count, basophil percentage of white blood cells, and myeloid white blood cell count may play causal roles in increasing the risk of severe COVID-19. COVID-19 postmortems have revealed that the loss of basophils from the blood is partially explained by its transfer into myelomonocytic lung infiltrates [7,39]. Basophils play important roles in modulating innate immune responses to infection and the subsequent tissue repair [40,41]. Individuals with low basophil count and basophil percentage may not present a strong enough innate immune response to the SARS-CoV-2 infection. Moreover, dysregulated activation of basophils also plays a role in regulating coagulation that may be pertinent to frequent thrombotic complications in severe cases of COVID-19 [42,43]. Immune system disorders have been suspected of playing roles in severe COVID-19 risk [44,45]. Basophils are known to potentiate humoral immune responses by producing interleukins 6 and 4 (IL-6 and IL-4) [46]. A recent immune-monitoring study of COVID-19 patients found that the relative abundance of basophils increases noticeably from acute to recovery phases and is positively associated with the level of anti-CoV-2 immunoglobulin G (IgG). These observations indicate that the extent of basophil depletion may affect the efficacy of IgG responses to the SARS-CoV-2 infection. The mediating factors between basophils and IgG responses may include IL-6, IL-4, or other unexplored factors, and the same study also observed a negative association between IL-6 and anti-CoV-2 IgG [47]. The complete elucidation of the potential mechanism warrants further investigation.
Previous systematic reviews and meta-analyses of observational studies have shown WBC traits to be associated with COVID-19 severity and susceptibility [5,6,7,9]. However, the associations between WBC traits and COVID-19 may be susceptible to reverse causation. Our reverse MR study did not find evidence for the effects of COVID-19 on any WBC trait. This result further supports the direct effect of basophils on severe COVID-19, ruling out the confounding of reverse causation. Therefore, further mechanistic understanding of basophil will shed light on the etiology of severe COVID-19 and provide multiple targets of intervention for prevention and treatment.
Lymphopenia, as a response to viral infection, has been frequently associated with severe COVID-19 risk [1,2,38,48]. Persistent eosinopenia after admission was associated with COVID-19 risk in previous retrospective and prospective observational studies [5,17,49,50,51]. However, these associations of WBC traits are more common in COVID-19 patients when compared to the reference range or healthy controls [17,18,19,25,50,51]. This discrepancy may reflect reverse causality in retrospective and prospective observational studies, with depleted lymphocyte and eosinophil counts as a result of the immune response to SARS-CoV-2 infection [52]. In our MR analysis, we found no evidence for the causal association between lymphocyte or eosinophils and the risk of COVID-19 severity and susceptibility. We showed that neutrophils are slightly associated with severe COVID-19, consistent with our recent MR results [53]. Our reverse MR analysis indicates no causal effects of severe COVID-19 on lymphocytes.
One assumption of MR is that the instrumental variable influences severe COVID-19 risk only through its effect on the specific WBC 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 for basophils. Pleiotropic instrumental variables, which are SNPs associated with other phenotypes directly or indirectly associated with the outcome, were detected by using the MR-Egger intercept and MR-PRESSO tests. Our MR study was performed with strong instrumental variables and adequate statistical power, and we have conducted extensive sensitivity analyses. The exposure and outcome GWAS were drawn from participants of European ancestry. Our MR results were replicated in at least one replication analysis with different GWAS of WBC to ensure robustness and reduce false positives. To assess whether WBC traits are associated with COVID-19 and to evaluate the direction of association, we applied a two-sample bidirectional MR study using the most up-to-date GWAS of WBC traits and COVID-19 severity and susceptibility.
Our study also has several limitations. Although we applied multiple sensitivity analyses, it is impossible to fully rule out that some genetic variants might be pleiotropic. Another limitation is that some GWAS of WBC traits and the HGI GWAS of COVID-19 have overlapped samples, which may bias the causal effect. F statistics of instrumental variables were estimated to mitigate this issue, and strong instruments are less susceptible to this bias. Although not all significant results remained after correction for multiple testing, we also utilized other GWAS of WBC and COVID-19 to replicate the MR results. The replicated observations across multiple analyses reduce the chance of false positives. Lastly, we emphasize that our results should be interpreted with caution, and future studies are needed to elucidate the mechanical roles of WBC traits in severe COVID-19.

5. Conclusions

Our MR results suggest that there are inverse causal associations between basophil count, basophil percentage, and myeloid WBC count and the severity of COVID-19. No causal effects on the susceptibility of COVID-19 were found. There is no evidence of reverse causation of COVID-19 severity and susceptibility on WBC traits. Overall, this study provides valuable insights that individuals with a lower genetic capacity for basophils are likely at risk of the severe forms of COVID-19, and that enhancing the production of basophils may be a potential prevention and treatment strategy.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4426/11/3/195/s1, Table S1: Details of the datasets used in the bidirectional Mendelian randomization study. Table S2: 401 SNPs significantly associated with white blood cell count used IVs in forward MR analyses derived from Vuckovic D et al. Table S3: 367 SNPs significantly associated with eosinophil count used IVs in forward MR analyses derived from Vuckovic D et al. Table S4: 388 SNPs significantly associated with eosinophil count used IVs in forward MR analyses derived from Chen MH et al. Table S5: 146 SNPs significantly associated with basophil count used IVs in forward MR analyses derived from Vuckovic D et al. Table S6: 166 SNPs significantly associated with basophil count used IVs in forward MR analyses derived from Chen MH et al. Table S7: 350 SNPs significantly associated with neutrophil count used IVs in forward MR analyses derived from Vuckovic D et al. Table S8: 383 SNPs significantly associated with neutrophil count used IVs in forward MR analyses derived from Chen MH et al. Table S9: 421 SNPs significantly associated with monocyte count used IVs in forward MR analyses derived from Vuckovic D et al. Table S10: 469 SNPs significantly associated with monocyte count used IVs in forward MR analyses derived from Chen MH et al. Table S11: 393 SNPs significantly associated with lymphocyte count used IVs in forward MR analyses derived from Vuckovic D et al. Table S12: 456 SNPs significantly associated with lymphocyte count used IVs in forward MR analyses derived from Chen MH et al. Table S13: 464 SNPs significantly associated with myeloid white cell count used IVs in forward MR analyses derived from Chen MH et al. Table S14: 347 SNPs significantly associated with eosinophil percentage of white cells used IVs in forward MR analyses derived from Vuckovic D et al. Table S15: 115 SNPs significantly associated with basophil percentage of white cells used IVs in forward MR analyses derived from Vuckovic D et al. Table S16: 295 SNPs significantly associated with neutrophil percentage of white cells used IVs in forward MR analyses derived from Vuckovic D et al. Table S17: 371 SNPs significantly associated with monocyte percentage of white cells used IVs in forward MR analyses derived from Vuckovic D et al. Table S18: 307 SNPs significantly associated with lymphocyte percentage of white cells used IVs in forward MR analyses derived from Vuckovic D et al. Table S19: 9 SNPs significantly associated with COVID-19 severity (very severe respiratory confirmed COVID-19 vs. population) used IVs in reverse MR analyses derived from the COVID-19 Host Genetics Initiative. Table S20: 5 SNPs significantly associated with COVID-19 severity (hospitalized COVID-19 vs. population) used IVs in reverse MR analyses derived from the COVID-19 Host Genetics Initiative. Table S21: 7 SNPs significantly associated with COVID-19 susceptibility (COVID-19 vs. population) used IVs in reverse MR analyses derived from the COVID-19 Host Genetics Initiative. Table S22: Primary discovery forward MR analyses assessing the causal effects of white blood cells on COVID-19 severity and susceptibility. Table S23: Replicated forward MR analyses assessing the causal effects of white blood cells on COVID-19 severity and susceptibility. Table S24: Primary discovery reverse MR analyses assessing the causal effects of COVID-19 severity and susceptibility on white blood cells. Table S25: Replicated reverse MR analyses assessing the causal effects of COVID-19 severity and susceptibility on white blood cells.

Author Contributions

Conceptualization, Y.S. and K.Y.; methodology, Y.S. and K.Y.; formal analysis, Y.S.; writing—original draft preparation, Y.S.; writing—review and editing, K.Y.; visualization, Y.S.; supervision, J.Z. and K.Y.; funding acquisition, K.Y. All authors have read and agreed to the published version of the manuscript.

Funding

K.Y. is supported by the University of Georgia Research Foundation. Funding sources had no involvement in the conception, design, analysis, or presentation of this work.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available in Supplementary Files and upon request. Scripts for the MR analysis are freely available on GitHub via https://github.com/yitangsun/WBC_COVID19_MR.

Acknowledgments

We would like to express our gratitude to all other Ye lab members for stimulating discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, F.; Yu, T.; Du, R.; Fan, G.; Liu, Y.; Liu, Z.; Xiang, J.; Wang, Y.; Song, B.; Gu, X.; 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–1062. [Google Scholar] [CrossRef]
  2. Huang, C.; Wang, Y.; Li, X.; Ren, L.; Zhao, J.; Hu, Y.; Zhang, L.; Fan, G.; Xu, J.; Gu, X.; et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 2020, 395, 497–506. [Google Scholar] [CrossRef] [Green Version]
  3. Wenham, C.; Smith, J.; Morgan, R. COVID-19: The gendered impacts of the outbreak. Lancet 2020, 395, 846–848. [Google Scholar] [CrossRef] [Green Version]
  4. Guzik, T.J.; Mohiddin, S.A.; DiMarco, A.; Patel, V.; Savvatis, K.; Marelli-Berg, F.M.; Madhur, M.S.; Tomaszewski, M.; Maffia, P.; D’Acquisto, F.; et al. COVID-19 and the cardiovascular system: Implications for risk assessment, diagnosis, and treatment options. Cardiovasc. Res. 2020, 116, 1666–1687. [Google Scholar] [CrossRef] [PubMed]
  5. Qin, C.; Zhou, L.; Hu, Z.; Zhang, S.; Yang, S.; Tao, Y.; Xie, C.; Ma, K.; Shang, K.; Wang, W.; et al. Dysregulation of Immune Response in Patients with Coronavirus 2019 (COVID-19) in Wuhan, China. Clin. Infect. Dis. 2020, 71, 762–768. [Google Scholar] [CrossRef]
  6. Chen, R.; Sang, L.; Jiang, M.; Yang, Z.; Jia, N.; Fu, W.; Xie, J.; Guan, W.; Liang, W.; Ni, Z.; et al. Longitudinal hematologic and immunologic variations associated with the progression of COVID-19 patients in China. J. Allerg. Clin. Immunol. 2020, 146, 89–100. [Google Scholar] [CrossRef]
  7. Laing, A.G.; Lorenc, A.; Barrio, I.D.M.D.; Das, A.; Fish, M.; Monin, L.; Muñoz-Ruiz, M.; McKenzie, D.R.; Hayday, T.S.; Francos-Quijorna, I.; et al. A dynamic COVID-19 immune signature includes associations with poor prognosis. Nat. Med. 2020, 26, 1623–1635. [Google Scholar] [CrossRef]
  8. Vitte, J.; Diallo, A.B.; Boumaza, A.; Lopez, A.; Michel, M.; Allardet-Servent, J.; Mezouar, S.; Sereme, Y.; Busnel, J.-M.; Miloud, T.; et al. A Granulocytic Signature Identifies COVID-19 and Its Severity. J. Infect. Dis. 2020, 222, 1985–1996. [Google Scholar] [CrossRef]
  9. Sambataro, G.; Giuffrè, M.; Sambataro, D.; Palermo, A.; Vignigni, G.; Cesareo, R.; Crimi, N.; Torrisi, S.E.; Vancheri, C.; Malatino, L.; et al. The Model for Early COVID-19 Recognition (MECOR) Score: A Proof-of-Concept for a Simple and Low-Cost Tool to Recognize a Possible Viral Etiology in Community—Acquired Pneumonia Patients during COVID-19 Outbreak. Diagnosis 2020, 10, 619. [Google Scholar] [CrossRef] [PubMed]
  10. Teng, Y.; Xu, J.; Zhang, Y.; Liu, Z.; Zhang, S. Mendelian randomization in COVID-19: Applications for cardiovascular comorbidities and beyond. EBioMedicine 2020, 57, 102847. [Google Scholar] [CrossRef] [PubMed]
  11. Wu, Z.; McGoogan, J.M. Characteristics of and Important Lessons from the Coronavirus Disease 2019 (COVID-19) Outbreak in China: Summary of a Report of 72314 Cases from the Chinese Center for Disease Control and Prevention. JAMA 2020, 323, 1239–1242. [Google Scholar] [CrossRef] [PubMed]
  12. Rao, S.; Lau, A.; So, H.-C. Exploring Diseases/Traits and Blood Proteins Causally Related to Expression of ACE2, the Putative Receptor of SARS-CoV-2: A Mendelian Randomization Analysis Highlights Tentative Relevance of Diabetes-Related Traits. Diabetes Care 2020, 43, 1416–1426. [Google Scholar] [CrossRef] [PubMed]
  13. Shi, S.; Qin, M.; Shen, B.; Cai, Y.; Liu, T.; Yang, F.; Gong, W.; Liu, X.; Liang, J.; Zhao, Q.; et al. Association of Cardiac Injury with Mortality in Hospitalized Patients With COVID-19 in Wuhan, China. JAMA Cardiol. 2020, 5, 802. [Google Scholar] [CrossRef] [Green Version]
  14. Kuri-Cervantes, L.; Pampena, M.B.; Meng, W.; Rosenfeld, A.M.; Ittner, C.A.; Weisman, A.R.; Agyekum, R.S.; Mathew, D.; Baxter, A.E.; Vella, L.A.; et al. Comprehensive mapping of immune perturbations associated with severe COVID-19. Sci. Immunol. 2020, 5, eabd7114. [Google Scholar] [CrossRef]
  15. Liu, J.; Liu, Y.; Xiang, P.; Pu, L.; Xiong, H.; Li, C.; Zhang, M.; Tan, J.; Xu, Y.; Song, R.; et al. Neutrophil-to-lymphocyte ratio predicts critical illness patients with 2019 coronavirus disease in the early stage. J. Transl. Med. 2020, 18, 206. [Google Scholar] [CrossRef]
  16. Zhou, J.; Sun, Y.; Huang, W.; Ye, K. Altered blood cell traits underlie a major genetic locus of severe COVID-19. J. Gerontol. Ser. A Boil. Sci. Med. Sci. 2021, glab035. [Google Scholar]
  17. Xie, G.; Ding, F.; Han, L.; Yin, D.; Lu, H.; Zhang, M. The role of peripheral blood eosinophil counts in COVID-19 patients. Allergy 2021, 76, 471–482. [Google Scholar] [CrossRef]
  18. Liu, K.; Fang, Y.-Y.; Deng, Y.; Liu, W.; Wang, M.-F.; Ma, J.-P.; Xiao, W.; Wang, Y.-N.; Zhong, M.-H.; Li, C.-H.; et al. Clinical characteristics of novel coronavirus cases in tertiary hospitals in Hubei Province. Chin. Med. J. 2020, 133, 1025–1031. [Google Scholar] [CrossRef]
  19. Terpos, E.; Ntanasis-Stathopoulos, I.; Elalamy, I.; Kastritis, E.; Sergentanis, T.N.; Politou, M.; Psaltopoulou, T.; Gerotziafas, G.; Dimopoulos, M.A. Hematological findings and complications of COVID-19. Am. J. Hematol. 2020, 95, 834–847. [Google Scholar] [CrossRef] [Green Version]
  20. Siedlinski, M.; Jozefczuk, E.; Xu, X.; Teumer, A.; Evangelou, E.; Schnabel, R.B.; Welsh, P.; Maffia, P.; Erdmann, J.; Tomaszewski, M.; et al. White Blood Cells and Blood Pressure. Circulation 2020, 141, 1307–1317. [Google Scholar] [CrossRef] [PubMed]
  21. Zhang, X.; Tan, Y.; Ling, Y.; Lu, G.; Liu, F.; Yi, Z.; Jia, X.; Wu, M.; Shi, B.; Xu, S.; et al. Viral and host factors related to the clinical outcome of COVID-19. Nat. Cell Biol. 2020, 583, 437–440. [Google Scholar] [CrossRef] [PubMed]
  22. Zhao, K.; Li, R.; Wu, X.; Zhao, Y.; Wang, T.; Zheng, Z.; Zeng, S.; Ding, X.; Nie, H. Clinical features in 52 patients with COVID-19 who have increased leukocyte count: A retrospective analysis. Eur. J. Clin. Microbiol. Infect. Dis. 2020, 39, 2279–2287. [Google Scholar] [CrossRef] [PubMed]
  23. Davies, N.M.; Holmes, M.V.; Smith, G.D. Reading Mendelian randomisation studies: A guide, glossary, and checklist for clinicians. BMJ 2018, 362, k601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Burgess, S.; Smith, G.D.; Davies, N.M.; Dudbridge, F.; Gill, D.; Glymour, M.M.; Hartwig, F.P.; Holmes, M.V.; Minelli, C.; Relton, C.L.; et al. Guidelines for performing Mendelian randomization investigations. Wellcome Open Res. 2020, 4, 186. [Google Scholar] [CrossRef] [PubMed]
  25. The COVID-19 Host Genetics Initiative. A global initiative to elucidate the role of host genetic factors in susceptibility and severity of the SARS-CoV-2 virus pandemic. Eur. J. Hum. Genet. 2020, 28, 715–718. [Google Scholar] [CrossRef]
  26. Vuckovic, D.; Bao, E.L.; Akbari, P.; Lareau, C.A.; Mousas, A.; Jiang, T.; Chen, M.-H.; Raffield, L.M.; Tardaguila, M.; Huffman, J.E.; et al. The Polygenic and Monogenic Basis of Blood Traits and Diseases. Cell 2020, 182, 1214–1231. [Google Scholar] [CrossRef]
  27. Chen, M.-H.; Raffield, L.M.; Mousas, A.; Sakaue, S.; Huffman, J.E.; Moscati, A.; Trivedi, B.; Jiang, T.; Akbari, P.; Vuckovic, D.; et al. Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations. Cell 2020, 182, 1198–1213. [Google Scholar] [CrossRef]
  28. Astle, W.J.; Elding, H.; Jiang, T.; Allen, D.; Ruklisa, D.; Mann, A.L.; Mead, D.; Bouman, H.; Riveros-Mckay, F.; Kostadima, M.A.; et al. The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease. Cell 2016, 167, 1415–1429. [Google Scholar] [CrossRef] [Green Version]
  29. Buniello, A.; MacArthur, J.A.L.; Cerezo, M.; Harris, L.W.; Hayhurst, J.; Malangone, C.; McMahon, A.; Morales, J.; Mountjoy, E.; Sollis, E.; et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019, 47, D1005–D1012. [Google Scholar] [CrossRef] [Green Version]
  30. Burgess, S.; Thompson, S.G. CRP CHD Genetics Collaboration Avoiding bias from weak instruments in Mendelian randomization studies. Int. J. Epidemiol. 2011, 40, 755–764. [Google Scholar] [CrossRef] [Green Version]
  31. Burgess, S.; Butterworth, A.; Thompson, S.G. Mendelian Randomization Analysis with Multiple Genetic Variants Using Summarized Data. Genet. Epidemiol. 2013, 37, 658–665. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Bowden, J.K.; Smith, G.D.; Burgess, S. Mendelian randomization with invalid instruments: Effect estimation and bias detection through Egger regression. Int. J. Epidemiol. 2015, 44, 512–525. [Google Scholar] [CrossRef] [Green Version]
  33. Bowden, J.; Smith, G.D.; Haycock, P.C.; Burgess, S. Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet. Epidemiol. 2016, 40, 304–314. [Google Scholar] [CrossRef] [Green Version]
  34. Hartwig, F.P.; Smith, G.D.; Bowden, J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int. J. Epidemiol. 2017, 46, 1985–1998. [Google Scholar] [CrossRef] [Green Version]
  35. Verbanck, M.; Chen, C.-Y.; Neale, B.; Do, R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 2018, 50, 693–698. [Google Scholar] [CrossRef]
  36. Korthauer, K.; Kimes, P.K.; Duvallet, C.; Reyes, A.; Subramanian, A.; Teng, M.; Shukla, C.; Alm, E.J.; Hicks, S.C. A practical guide to methods controlling false discoveries in computational biology. Genome Biol. 2019, 20, 118. [Google Scholar] [CrossRef] [Green Version]
  37. Hemani, G.; Zheng, J.; Elsworth, B.; Wade, K.H.; Haberland, V.; Baird, D.; Laurin, C.; Burgess, S.; Bowden, J.; Langdon, R.; et al. The MR-Base platform supports systematic causal inference across the human phenome. eLife 2018, 7, 34408. [Google Scholar] [CrossRef] [PubMed]
  38. Zhang, J.-J.; Dong, X.; Cao, Y.-Y.; Yuan, Y.-D.; Yang, Y.-B.; Yan, Y.-Q.; Akdis, C.A.; Gao, Y.-D. Clinical characteristics of 140 patients infected with SARS-CoV-2 in Wuhan, China. Allergy 2020, 75, 1730–1741. [Google Scholar] [CrossRef] [PubMed]
  39. Carsana, L.; Sonzogni, A.; Nasr, A.; Rossi, R.S.; Pellegrinelli, A.; Zerbi, P.; Rech, R.; Colombo, R.; Antinori, S.; Corbellino, M.; et al. Pulmonary post-mortem findings in a series of COVID-19 cases from northern Italy: A two-centre descriptive study. Lancet Infect. Dis. 2020, 20, 1135–1140. [Google Scholar] [CrossRef]
  40. Crivellato, E.; Travan, L.; Ribatti, D. Mast Cells and Basophils: A Potential Link in Promoting Angiogenesis during Allergic Inflammation. Int. Arch. Allerg. Immunol. 2010, 151, 89–97. [Google Scholar] [CrossRef]
  41. Ye, Q.; Wang, B.; Mao, J. The pathogenesis and treatment of the “Cytokine Storm” in COVID-19. J. Infect. 2020, 80, 607–613. [Google Scholar] [CrossRef]
  42. Swystun, L.L.; Liaw, P.C. The role of leukocytes in thrombosis. Blood 2016, 128, 753–762. [Google Scholar] [CrossRef] [Green Version]
  43. Klok, F.A.; Kruip, M.J.H.A.; van der Meer, N.J.M.; Arbous, M.S.; Gommers, D.A.M.P.J.; Kant, K.M.; Kaptein, F.H.J.; van Paassen, J.; Stals, M.A.M.; Huisman, M.V.; et al. Incidence of thrombotic complications in critically ill ICU patients with COVID-19. Thromb. Res. 2020, 191, 145–147. [Google Scholar] [CrossRef] [PubMed]
  44. Jesenak, M.; Brndiarova, M.; Urbancikova, I.; Rennerova, Z.; Vojtkova, J.; Bobcakova, A.; Ostro, R.; Banovcin, P. Immune Parameters and COVID-19 Infection—Associations with Clinical Severity and Disease Prognosis. Front. Cell. Infect. Microbiol. 2020, 10, 364. [Google Scholar] [CrossRef]
  45. Azkur, A.K.; Akdis, M.; Azkur, D.; Sokolowska, M.; Van De Veen, W.; Brüggen, M.; O’Mahony, L.; Gao, Y.; Nadeau, K.; Akdis, C.A. Immune response to SARS-CoV-2 and mechanisms of immunopathological changes in COVID-19. Allergy 2020, 75, 1564–1581. [Google Scholar] [CrossRef]
  46. Denzel, A.; Maus, U.A.; Gomez, M.R.; Moll, C.; Niedermeier, M.; Winter, C.; Maus, R.; Hollingshead, S.K.; Briles, D.E.; Kunz-Schughart, A.L.; et al. Basophils enhance immunological memory responses. Nat. Immunol. 2008, 9, 733–742. [Google Scholar] [CrossRef] [PubMed]
  47. Rodriguez, L.; Pekkarinen, P.T.; Lakshmikanth, T.; Tan, Z.; Consiglio, C.R.; Pou, C.; Chen, Y.; Mugabo, C.H.; Nguyen, N.A.; Nowlan, K.; et al. Systems-Level Immunomonitoring from Acute to Recovery Phase of Severe COVID-19. Cell Rep. Med. 2020, 1, 100078. [Google Scholar] [CrossRef] [PubMed]
  48. Vabret, N.; Britton, G.J.; Gruber, C.; Hegde, S.; Kim, J.; Kuksin, M.; Levantovsky, R.; Malle, L.; Moreira, A.; Park, M.D.; et al. Immunology of COVID-19: Current State of the Science. Immunity 2020, 52, 910–941. [Google Scholar] [CrossRef] [PubMed]
  49. Wang, D.; Li, R.; Wang, J.; Jiang, Q.; Gao, C.; Yang, J.; Ge, L.; Hu, Q. Correlation analysis between disease severity and clinical and biochemical characteristics of 143 cases of COVID-19 in Wuhan, China: A descriptive study. BMC Infect. Dis. 2020, 20, 1–9. [Google Scholar] [CrossRef] [PubMed]
  50. Sun, D.-W.; Zhang, D.; Tian, R.-H.; Li, Y.; Wang, Y.-S.; Cao, J.; Tang, Y.; Zhang, N.; Zan, T.; Gao, L.; et al. The underlying changes and predicting role of peripheral blood inflammatory cells in severe COVID-19 patients: A sentinel? Clin. Chim. Acta 2020, 508, 122–129. [Google Scholar] [CrossRef]
  51. Li, Q.; Ding, X.; Xia, G.; Chen, H.-G.; Chen, F.; Geng, Z.; Xu, L.; Lei, S.; Pan, A.; Wang, L.; et al. Eosinopenia and elevated C-reactive protein facilitate triage of COVID-19 patients in fever clinic: A retrospective case-control study. EClinicalMedicine 2020, 23, 100375. [Google Scholar] [CrossRef] [PubMed]
  52. Sattar, N.; Preiss, D. Reverse Causality in Cardiovascular Epidemiological Research. Circulation 2017, 135, 2369–2372. [Google Scholar] [CrossRef] [PubMed]
  53. Sun, Y.; Zhou, J.; Ye, K. Prioritizing causal risk factors for severe COVID-19: An exhaustive mendelian randomization study. Res. Square 2021. [Google Scholar] [CrossRef]
Figure 1. Overview of two-sample MR study on the bidirectional associations between WBC and COVID-19. The forward MR analysis evaluating the causal effects of WBC on COVID-19 was indicated with orange arrows, while the reverse MR analysis evaluating the causal effects of COVID-19 on WBC was indicated with blue arrows. WBC, white blood cells; SNP, single nucleotide polymorphism.
Figure 1. Overview of two-sample MR study on the bidirectional associations between WBC and COVID-19. The forward MR analysis evaluating the causal effects of WBC on COVID-19 was indicated with orange arrows, while the reverse MR analysis evaluating the causal effects of COVID-19 on WBC was indicated with blue arrows. WBC, white blood cells; SNP, single nucleotide polymorphism.
Jpm 11 00195 g001
Figure 2. Forward MR estimates of the effects of WBC traits on COVID-19 severity and susceptibility. Odds ratios and 95% confidence intervals were derived using the IVW random-effects model. The p-value threshold (0.05) was used for the significant association. Detailed summary statistics could be found in Supplementary Table S22.
Figure 2. Forward MR estimates of the effects of WBC traits on COVID-19 severity and susceptibility. Odds ratios and 95% confidence intervals were derived using the IVW random-effects model. The p-value threshold (0.05) was used for the significant association. Detailed summary statistics could be found in Supplementary Table S22.
Jpm 11 00195 g002
Figure 3. Reverse MR estimates of the effects of COVID-19 severity and susceptibility on WBC traits. Odds ratios and 95% confidence intervals were derived using the IVW random-effects model. The p-value threshold (0.05) was used for the significant association. Detailed summary statistics could be found in Supplementary Table S24.
Figure 3. Reverse MR estimates of the effects of COVID-19 severity and susceptibility on WBC traits. Odds ratios and 95% confidence intervals were derived using the IVW random-effects model. The p-value threshold (0.05) was used for the significant association. Detailed summary statistics could be found in Supplementary Table S24.
Jpm 11 00195 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, Y.; Zhou, J.; Ye, K. White Blood Cells and Severe COVID-19: A Mendelian Randomization Study. J. Pers. Med. 2021, 11, 195. https://doi.org/10.3390/jpm11030195

AMA Style

Sun Y, Zhou J, Ye K. White Blood Cells and Severe COVID-19: A Mendelian Randomization Study. Journal of Personalized Medicine. 2021; 11(3):195. https://doi.org/10.3390/jpm11030195

Chicago/Turabian Style

Sun, Yitang, Jingqi Zhou, and Kaixiong Ye. 2021. "White Blood Cells and Severe COVID-19: A Mendelian Randomization Study" Journal of Personalized Medicine 11, no. 3: 195. https://doi.org/10.3390/jpm11030195

APA Style

Sun, Y., Zhou, J., & Ye, K. (2021). White Blood Cells and Severe COVID-19: A Mendelian Randomization Study. Journal of Personalized Medicine, 11(3), 195. https://doi.org/10.3390/jpm11030195

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

Article Metrics

Back to TopTop