- Research
- Open access
- Published:
A systematic two-sample and bidirectional MR process highlights a unidirectional genetic causal effect of allergic diseases on COVID-19 infection/severity
Journal of Translational Medicine volume 22, Article number: 94 (2024)
Abstract
Background
Allergic diseases (ADs) such as asthma are presumed risk factors for COVID-19 infection. However, recent observational studies suggest that the assumed correlation contradicts each other. We therefore systematically investigated the genetic causal correlations between various ADs and COVID-19 infection/severity.
Methods
We performed a two-sample, bidirectional Mendelian randomization (MR) study for five types of ADs and the latest round of COVID-19 GWAS meta-analysis datasets (critically ill, hospitalized, and infection cases). We also further validated the significant causal correlations and elucidated the potential underlying molecular mechanisms.
Results
With the most suitable MR method, asthma consistently demonstrated causal protective effects on critically ill and hospitalized COVID-19 cases (OR < 0.93, p < 2.01 × 10–2), which were further confirmed by another validated GWAS dataset (OR < 0.92, p < 4.22 × 10–3). In addition, our MR analyses also observed significant causal correlations of food allergies such as shrimp allergy with the risk of COVID-19 infection/severity. However, we did not find any significant causal effect of COVID-19 phenotypes on the risk of ADs. Regarding the underlying molecular mechanisms, not only multiple immune-related cells such as CD4+ T, CD8+ T and the ratio of CD4+/CD8+ T cells showed significant causal effects on COVID-19 phenotypes and various ADs, the hematology traits including monocytes were also significantly correlated with them. Conversely, various ADs such as asthma and shrimp allergy may be causally correlated with COVID-19 infection/severity by affecting multiple hematological traits and immune-related cells.
Conclusions
Our systematic and bidirectional MR analyses suggest a unidirectional causal effect of various ADs, particularly of asthma on COVID-19 infection/severity, but the reverse is not true. The potential underlying molecular mechanisms of the causal effects call for more attention to clinical monitoring of hematological cells/traits and may be beneficial in developing effective therapeutic strategies for allergic patients following infection with COVID-19.
Graphical Abstract
Key messages
-
1.
We performed a two-sample, bidirectional Mendelian randomization study for five types of allergic diseases with COVID-19 infection/severity (critically ill, hospitalized, and infection cases), and also further elucidated the potential underlying molecular mechanisms.
-
2.
Our MR analyses consistently suggest a unidirectional protective causal effect of asthma on COVID-19 infection/severity, but the reverse is not true.
-
3.
Food allergies such as shrimp allergy showed significant causal correlations with the risk of COVID-19 phenotypes.
-
4.
Allergic diseases may be causally correlated with COVID-19 infection/severity by inducing abnormal fluctuations of multiple immune-related cells and hematological traits in peripheral blood.
-
5.
The findings call for more attention to clinical monitoring of hematological traits and may be beneficial in developing effective therapeutic strategies for allergic patients following infection with COVID-19.
Introduction
According to the Coronavirus Disease 2019 (COVID-19) Weekly Epidemiological Update (accessed on 15th May 2023) published by the World Health Organization, approximately 766 million confirmed cases and 6.9 million deaths have been reported globally. It is worth noting that more than one-quarter of COVID-19 cases suffered from at least one comorbid condition (such as diabetes, hypertension, cardiovascular disease, and respiratory diseases) [1, 2] and the disease was more likely to be severe in patients with chronic conditions [3]. Among these comorbid conditions, respiratory diseases attracted the most attention as they showed the most obvious phenotypes associated with COVID-19. Asthma is a common paroxysmal allergic condition partially induced by infection with respiratory viruses [4]. Atopic dermatitis/rhinitis and food allergies can also be caused by an abnormal adaptive immune response to a specific allergen [5, 6].
It is assumed that a close correlation exists between COVID-19 infection and allergic diseases (ADs) as both conditions result from a dysfunctional immune response [7]. However, findings from observational studies about the correlation between COVID-19 infection/severity and ADs contradict each other [6, 8,9,10] and are under active debate [11,12,13,14]. The controversy surrounding traditional observational studies may derive from multiple confounding factors, such as hypertension, diabetes, obesity and other comorbidities in COVID-19 patients [1, 2]. Accumulating evidence have implied that the predisposition to COVID-19 infection and its severity is closely associated with a patient’s age, body mass index (BMI), and smoking status [15,16,17]. Moreover, results of conventional observational studies are usually affected by reverse causality and, because of that, may not easily be validated via experimental approaches [18]. Although randomized controlled trials (RCTs) are the gold standard for inferring causal relationships, they have been difficult to implement due to high costs and required enrollment of large number of eligible subjects [19].
In recent years, Mendelian randomization (MR) has been used to investigate genetic variants in terms of single-nucleotide polymorphisms (SNPs) or copy number variations as instrumental variables to infer causal correlations between exposures and outcomes [20]. As genetic variants are determined by nature and usually cannot be altered by acquired environmental and behavioral factors, they are considered reliable and practical instrumental variables for inferring causal correlations, which may not be easily affected by issues inherent in observational studies and RCTs [21]. In addition, host genetics has been reported to play a vital role in the predisposition to infectious diseases [22]. Based on the MR approach, severe COVID-19 symptoms were reported to be causally linked to multiple risk factors, including BMI and smoking intensity [23,24,25]. Our group performed a phenome-wide MR study, which suggested that diabetes and related traits may increase pulmonary angiotensin-converting enzyme 2 (ACE2) gene expression, which may in turn increase the rate or severity of COVID-19 infection [26]. However, the potential causal correlation between ADs and COVID-19 infection/severity has been less studied.
In the present study, we aimed to answer the following two questions: (1) Are there causal correlations between various ADs and COVID-19 infection/severity? (2) If so, how do they affect each other? Here, we employed a process of two-sample, bidirectional MR analysis to systematically explore the potential genetic correlations between ADs and COVID-19 infection/severity. Upon obtaining significant correlations, we would further validate the results using another available genome-wide association study (GWAS) summary datasets. We also applied the process to investigate the molecular mechanism underlying the identified significant correlations according to the putative theory about the potential connections between ADs and COVID-19 phenotypes.
Methods
GWAS summary statistics
We collated all GWAS summary statistics for COVID-19, ADs, peripheral blood hematological traits, immune-related cell count data, as well as peripheral blood ACE2 protein expression from publicly available databases (Additional file 1: Table S1). Notably, all the GWAS summary data were based on predominantly European samples (> 88%), except for the two GWAS summary datasets of ADs caused by shrimp and peach. Moreover, all the summary data were originally subjected to proper correction regarding population stratification to avoid the negative effects of different genetic backgrounds.
COVID-19 phenotypes and ADs
We obtained the latest 7th round of COVID-19 Host Genome Initiative GWAS meta-analysis datasets (https://www.covid19hg.org/results/r7/) (Additional file 1: Table S1), which included 18,152 critically ill COVID-19 positive cases and 1,145,546 controls (A2), 44,986 hospitalized COVID-19 positive cases and 2,356,386 controls (B2), and 159,840 reported COVID-19 infection cases and 2,782,977 controls (C2). Critically ill COVID-19 positive cases were defined as laboratory-confirmed SARS-CoV-2 infections AND hospitalization for COVID-19 AND death or respiratory support, while hospitalized COVID-19 positive cases were defined as laboratory-confirmed SARS-CoV-2 infections AND hospitalization for COVID-19. The mildest phenotype (called the reported COVID-19 infection) was defined as laboratory-confirmed SARS-CoV-2 infection OR a COVID-19 diagnosis based on electronic health record/International Classification of Disease coding/physician confirmation OR self-reported COVID-19 via questionnaire. A control was defined as any other person apart from the patients.
Regarding ADs, we comprehensively searched all available databases to systematically investigate correlations between all available types of ADs and COVID-19. Eventually, we obtained seven sets of summary datasets for five types of ADs, including one board allergic disease (BAD, a combination of hay fever, asthma, and eczema), asthma, allergic dermatitis (ADE), shrimp allergy (SA) and peach allergy (PA) (Additional file 1: Table S1) [27,28,29,30]. Two of these summary datasets for asthma (Asthma2018) and ADE (ADE2021) were used to validate the primarily significant correlations [31, 32].
Peripheral blood hematological traits and count data of immune-related cells
We collected a set of 13 different peripheral blood hematological traits, including total blood hemoglobin (HB), red blood cell count, mean corpuscular volume (MCV), platelet count (PLT), white blood cell count, hematocrit (HT), mean cell hemoglobin (MCH), mean corpuscular hemoglobin concentration, and neutrophil (NEUT), monocyte (MONO), eosinophil (EOS), basophil (BASO), and lymphocyte (LYMPH) counts [33] (Additional file 1: Table S1). We also obtained six GWAS summary datasets of count data of immune-related cells, comprising CD3+ T cells, CD4+ T cells, CD8+ T cells, CD19+ B cells, CD56+ natural killer cells, and the ratio of CD4+:CD8+ T cells (Additional file 1: Table S1) [34].
Peripheral blood ACE2 protein expression
Two sets of GWAS summary statistics of soluble ACE2 protein expression level in peripheral blood were obtained for this study. One set was obtained from the KORA (Cooperative Health Research in the Augsburg Region) study [35], and the other from a systematic meta-analysis of 14 European cohorts [36] (Additional file 1: Table S1).
Study design
The primary MR analyses contained three successive parts. In the first part, we treated each COVID-19 phenotype (A2, B2, and C2) as an exposure and each AD as an outcome to evaluate whether the COVID-19 outbreak would genetically affect the incidence rates of ADs (Fig. 1A, C). Conversely, we considered each type of AD as an exposure and COVID-19 phenotypes as outcomes to identify genetic risk factors for COVID-19 infection/severity (Fig. 1B, C). In the second part, we employed two available independent GWAS datasets (Asthma2018 and ADE2021) to further validate those significant causal correlations identified in the first part (Fig. 1D). Finally, we investigated the molecular mechanisms underlying those significant correlations based on well-known theories of pathogenesis for COVID-19 or ADs (Fig. 1E). In brief, we treated COVID-19 infection/severity or ADs as exposures and peripheral blood hematological traits and count data of immune-related cells as outcomes to evaluate whether COVID-19 infection/severity or ADs played a vital role in affecting hematological traits or immune-related cells. Conversely, we also assessed if there are significant effects of hematological traits or cell counts on COVID-19 infection/severity and ADs. In addition, we also considered various allergic diseases as exposures and peripheral blood ACE2 protein expression as an outcome to explore whether allergic diseases would affect COVID-19 infection/severity by altering ACE2 protein expression in peripheral blood.
Two-sample MR analysis
First, all the GWAS summary statistics went through multiple pre-processing steps, including transferring genomic coordinates, adding missing but necessary parameters and removing those strong genetically correlated SNPs (r2 < 0.1, distance = 1000 kb) [37] (see details in Additional file 1: Methods). Independent SNPs with an original p-value below 5 × 10–8 would be treated as valid instrumental variables. In addition, we also employed MR robust adjusted profile score (RAPS) method, which took into account multiple weak instruments simultaneously using a robust procedure [38].
All two-sample MR analyses were performed using TwoSampleMR (version 0.4.26) package in the R programming language (version 4.2.0). Unlike common applications, we applied the most suitable MR approach in different scenarios. We first employed the MR Egger regression method to determine if there was a horizontal pleiotropy for the genetic instruments and the MR inverse-variance weighted (IVW) method to evaluate if there was an obvious heterogeneity among multiple genetic instruments. Then, the most suitable MR approach was determined in four different scenarios. The fixed-effect IVW method was applied when there was neither horizontal pleiotropy nor significant heterogeneity (scenario 1, Additional file 1: S1). The random-effect IVW method was employed when there was obvious heterogeneity but no horizontal pleiotropy (Additional file 1: S2). The MR Egger method or weighted median was applied when there was an obvious horizontal pleiotropy (Additional file 1: S3). The Wald ratio method was employed when there was only one valid genetic instrument in MR analysis (Additional file 1: S4). With that, only one hypothesis was tested for each pair of exposure and outcome, from which a p-value below 0.05 was considered indicative of significant genetic correlations. In terms of effect size in causal correlation when the exposure was a binary trait, the odds ratio (OR) could be roughly considered to represent the likelihood that the outcome would occur when an individual is exposed to one specific condition. In addition, we employed the Steiger directionality test to evaluate if the causal directions between the hypothesized exposures and outcomes are true [39].
Results
Selection of genetic instruments
Suitable independent genetic instruments were first selected from all the exposure datasets according to the pre-set criteria (r2 < 0.1, distance = 1000 kb, p < 5 × 10–8). In MR analyses that treated ADs as exposures, independent genetic instruments were selected from seven GWAS datasets of ADs (BAD, n = 144; Asthma2020, n = 231; ADE2015, n = 20; PA, n = 5; SA, n = 19; ADE2021, n = 52; and Asthma2018, n = 36) (Tables 1 and 2). For COVID-19 infection and severity, we also obtained 73, 81, and 57 independent instruments for critically ill cases (A2), hospitalized patients (B2), and infection cases (C2), respectively (Additional file 1: Table S2). The detailed number of independent SNPs used in the subsequent MR analyses is also shown in Tables 1 and 2, in which five ADs were treated as exposures and COVID-19 infection/severity were considered as outcomes in each set of analyses.
Genetic causal correlations between five types of ADs and COVID-19 infection/severity
In this set of MR analyses, we treated each type of AD as an exposure and COVID-19 infection/severity as an outcome to explore whether there was a significant causal effect of any type of AD on COVID-19 infection/severity. By adopting the most suitable MR method, our MR analyses found that asthma was a causal protective factor for critically ill and hospitalized COVID-19, respectively (OR = 0.90, p = 2.00 × 10–4; OR = 0.93, p = 2.01 × 10–2) (Table 1, Fig. 2A, B). Meanwhile, the protective effects of asthma on COVID-19 phenotypes were further confirmed in another GWAS dataset (Asthma2018) with the most suitable IVW (random-effects) MR method, in which asthma was also identified as a protective factor for critically ill and hospitalized COVID-19 (OR = 0.88, p = 4.22 × 10–3; OR = 0.92, p = 5.61 × 10–3; IVW(random-effects)) (Table 2, Fig. 2C, D). Besides that, the significant causal effects were also suggested by the other two MR methods, including the IVW (fixed-effects) and RAPS methods (p < 6.44 × 10–8; Table 2). Regarding the broad allergic disease (BAD), our MR analysis found a significant protective effect of it on critically ill COVID-19 (OR = 0.93, p = 3.26 × 10–2) (Table 1, Fig. 2E). As for food allergies, our MR analyses observed significantly causal correlations of PA with COVID-19 infection (OR = 1.01, p = 2.66 × 10–2) and SA with all the three COVID-19 phenotypes (A2: OR = 1.04, p = 6.73 × 10–3; B2: OR = 1.03, p = 1.61 × 10–3; C2: OR = 1.02, p = 2.43 × 10–3) (Table 1, Fig. 2F, G, and I), implying that food allergies may act as risk factors for both COVID-19 infection and severity.
In contrast to the foregoing set, we treated COVID-19 infection/severity as exposures and each type of AD as an outcome to evaluate if the prevalence of ADs would be affected by COVID-19 infection/severity. With the most suitable MR method, our MR analyses only observed a marginal negative correlation between COVID-19 infection (C2) and PA (Additional file 1: Table S2). However, the acting direction of the marginal correlation was found to be false by the directionality test. Except that, we did not find any other significant causal effect of COVID-19 infection/severity on ADs (Additional file 1: Table S2).
Genetic causal correlations between COVID-19 infection/severity and peripheral blood immune-related cells and hematological traits
We aimed to investigate the molecular mechanism underlying the causal correlations observed between COVID-19 phenotypes and various ADs. In this set of MR analyses, we first treated 13 types of hematological traits and 6 kinds of immune-related cell counts as exposures and three COVID-19 phenotypes as outcomes to explore whether the immune-related traits/cells in peripheral blood would significantly affect COVID-19 infection/severity. The most suitable MR methods observed a causal risk effect of monocyte counts on critically ill COVID-19 (OR = 1.14, p = 9.39 × 10–3) and COVID-19 infection (OR = 1.05, p = 3.63 × 10–3), respectively (Table 3) (Fig. 3A, B). In addition, the ratio of CD4+ and CD8+ T cells and the CD4+ T cell counts were found to be risk factors for all COVID-19 phenotypes (CD4+ T/8+ T: OR = 1.02–1.14, p = 7.21 × 10–9–4.63 × 10–3, Fig. 3C–E; CD4+ T: OR = 1.07–1.31, p = 3.17 × 10–4–3.42 × 10–3), while CD8+ T cell counts had a protective causal effect on critically ill (OR = 0.84, p = 3.83 × 10–3) and hospitalized COVID-19 cases (OR = 0.91, p = 2.17 × 10–2), respectively (Table 3) (Fig. 3F, G). Besides that, the CD56+ natural killer cell counts showed a causal risk effect on hospitalized COVID-19 cases (OR = 1.12, p = 2.83 × 10–3). Notably, all this set of MR analyses were indicated to be true in the directionality tests (Table 3).
Conversely, we treated three kinds of COVID-19 phenotypes as exposures and the 19 immune-related traits/cells as outcomes and found some significant correlations (Additional file 1: Table S3). However, most of the directionality tests in MR analyses were found to be false, except for the mean cell hemoglobin (MCH). The critically ill COVID-19 cases showed a marginal significant association with an increased MCH level in peripheral blood (OR = 1.09, p = 4.79 × 10–2) (Additional file 1: Table S3) (Fig. 3H).
Genetic causal correlations between five ADs and peripheral blood hematological traits, immune-related cell counts and ACE2 protein expression
In this set of MR analyses, we first treated 13 types of hematological traits and 6 kinds of immune-related cell counts as exposures and five types of ADs as outcomes to explore if the immune-related traits in peripheral blood would affect the occurrence of various ADs. With the most suitable MR methods, we observed that two kinds of immune-related cells had protective causal effects on asthma (CD4+ T: OR = 0.78, p = 1.73 × 10–13; CD4+ T/CD8+ T: OR = 0.92, p = 2.42 × 10–02) (Table 4 and Fig. 4A). Besides that, the lymphocyte cell counts also showed a protective effect on asthma (OR = 0.68, p = 5.23 × 10–14). Regarding the broad allergic disease (BAD), CD4+ T cells and monocyte cells had significant causal effects on it (p = 3.90 × 10–4 and p = 2.46 × 10–3) (Fig. 4B). Besides that, the CD4+ T cell showed a risk causal effect on PA (OR = 3.92, p = 2.74 × 10–4) and the ratio of CD4+ T/CD8+ T also had a risk effect on atopic dermatitis (OR = 2.33, p = 4.70 × 10–2) (Fig. 4C). The directionality tests confirmed that all the directions of MR analyses were true. Conversely, we considered five types of ADs as exposures and found significant causal correlations of asthma with total blood hemoglobin and monocytes, respectively (p = 3.66 × 10–2 and p = 5.46 × 10–3) (Additional file 1: Table S3) (Fig. 4D, E). In addition, shrimp allergy had a significant causal effect on CD3+ T cells (p = 4.75 × 10–2, Fig. 4F) and multiple hematological traits, including lymphocytes, monocytes, neutrophils and white blood cells (p = 6.30 × 10–6–1.85 × 10–2) (Additional file 1: Table S3). Although our MR analyses observed other significant correlations in this set of MR analyses, the directionality tests identified the directions of these correlations were false.
As for ACE2 protein expression, we treated five types of ADs as exposures and peripheral blood ACE2 protein expression level from two GWAS datasets as outcomes to investigate if causal effects of ADs on COVID-19 phenotypes were mediated by the ACE2 expression level. Employing 220 strong genetic instruments, our MR analyses found a marginal correlation between asthma and ACE2 protein expression level (OR = 1.02, p = 5.74 × 10–2) (Additional file 1: Table S4). In addition, we also observed a marginal negative correlation of shrimp allergy with ACE2 protein expression (OR = 0.87, p = 8.93 × 10–2) (Additional file 1: Table S4). But beyond that, the MR analyses did not find any other significant/marginal correlation between various ADs and ACE2 protein expression level.
Discussion
In the present study, our systematic two-sample and bidirectional MR analysis process identified unidirectional causal effects of various ADs, particularly for a protective causal effect of asthma on COVID-19 infection/severity, but the reverse is not true. As for the underlying molecular mechanisms, the unidirectional causal effects may be mediated by abnormal fluctuations of immune-related cells and hematology traits in peripheral blood, such as CD4+ T, CD8+ T and monocyte cells. The change of ACE2 protein expression level might also play a non-negligible role on the causal correlations. The two-sample MR approach is a relatively well-established technique for evaluating causal correlations for diseases/traits of interest using GWAS summary statistics [40].
Here we utilized MR approach to show unidirectional causal effects of various ADs on COVID-19 phenotypes, including a protective causal effect of asthma on COVID-19, but not vice versa. These unidirectional causal effects may be mediated by fluctuations in amounts of immune-related cells in the peripheral blood, including CD4+ and CD8+ T cells as well as monocytes, and by other hematological traits. The circulating levels of ACE2 protein might also contribute to observed causal correlations.
Strong protective causal effects of asthma on COVID-19 infection/severity
Using the most suitable MR method, our MR analyses suggested that asthma was a causal protective factor for critically ill and hospitalized COVID-19 cases (OR < 0.93), which was further validated by another independent asthma GWAS dataset with multiple MR methods (OR, 0.87–0.98). The validation MR analysis confirmed a solid, negative causal correlation of asthma with severe course of COVID-19 (OR, 0.87–0.92), suggesting asthma patients may experience less severe symptoms following infection with COVID-19. Although it was commonly taken for granted that asthma patients showed increased susceptibility to COVID-19, or develop a disease of greater severity [41], the results of conventional observational studies either do not support this this presumption or contradict each other [6, 8, 10, 42]. Two previous unidirectional MR analyses also reported conclusions that were opposite to each other, possibly due to use of different GWAS datasets [43, 44], and, therefore, varying sets of independent genetic instruments. Moreover, lacking validation step in an independent dataset would to some extent decrease the reliability of outcome. Our MR analyses did not only include a validation step but also utilized the largest GWAS dataset to date for the analysis [27], which contained the largest number of independent genetic instruments to infer the correlation and generated highly promising results of protective causal effects.
Immune-related cells mediated protective causal effects of asthma on COVID-19 infection/severity
Our MR analyses implied that the protective causal effects of asthma on COVID-19 infection/severity were mediated by abnormal change of multiple immune-related cells, mainly including the CD4+ T cells, CD8+ T cells, the ratio of CD4+ T/CD8+ T and the monocyte cells in peripheral blood. The ratio of CD4+ T/CD8+ T cells in the peripheral blood of healthy adults and mice is approximately 2:1, and breaking the balance usually indicates diseases relating to immunodeficiency induced by an impaired immune system [45,46,47,48]. Recently Huang et al. performed a systematic meta-analysis and found that COVID-19 patients showed significantly declined count of B cells, natural killer cells and total lymphocytes [49]. Notably, both the CD4+ and CD8+ T cells had a far greater extent of decrease. Moreover, a study of 60 COVID-19 patients suggested that the decreased CD8+ T cells and the increased ratio of CD4+/CD8+ T cells were associated with poor treatment efficacy [50], which provided strong evidence to our findings about the positive genetic causal effect of an increased ratio and the inverse causal effect of a decreased CD8+ T cell count. Patients with asthma also showed an increased ratio and a decreased CD8+ T cell count in peripheral blood compared to healthy controls [51]. After 7-day antigen stimulation, the ratio would have a significant increase in atopic dermatitis syndrome patient’s peripheral blood mononuclear cells (PBMCs) as compared with healthy control’s cells [52].
Regarding the monocyte cells in peripheral blood, Ren et al. applied single-cell RNA sequencing to 284 samples from 196 COVID-19 patients and reported that the percentage of CD14+ monocytes in PBMCs was significantly elevated in COVID-19 patients, especially in COVID-19 hospitalized and critically ill patients [53]. Meanwhile, subtypes of macrophages and monocytes had the highest cytokine and inflammatory scores in the severe COVID-19 samples, implying that these cell subtypes might be the major sources driving the inflammatory storm in lung tissue. The main manifestation of COVID-19 severity was severe acute respiratory distress syndrome, which was mainly caused by pneumonia, sepsis, or pulmonary aspiration induced by a cytokine storm. During a cytokine storm, immune and non-immune cells would release a large number of pro-inflammatory cytokines, including IL-6, tumor necrosis factor-α (TNF-α), and chemokines, which would cause substantial damage to the host immune response and further induce various ADs such as allergic asthma [54, 55]. These findings are consistent with our present results about a positive causal correlation between monocytes and COVID-19 cases, particularly in those severe patients.
The protective effects were also partially explained by the MR analyses of asthma and ACE2 protein expression using two different MR methods, which showed that asthma was marginally correlated with the protein expression of soluble ACE2 in peripheral blood. The soluble ACE2 has a protective effect against virus-induced lung injury by increasing the amount of angiotensin 1–7, which produces local vasodilation, thereby reducing blood pressure [56, 57]. In contrast, the surface bound ACE2 is the main receptor for COVID-19 virus to invade into cells that increasing the incidence rate of COVID-19 infection [36, 58]. Several previous studies have suggested that allergic sensitization in asthma patients is associated with a low surface bound ACE2 expression in sputum cells in the upper and lower respiratory tracts [4, 59]. In addition, the ACE2 gene is less active in asthma patients and probably limits viral entry into the respiratory epithelium [60]. Moreover, it was reported that epithelial surface bound ACE2 expression is inversely associated with the levels of Th2 cytokines (IL-4, IL-5, and IL-13) [61] that could alleviate the viral-induced release of interferons and downregulate the cytokine storm typical of increasing COVID-19 severity [62].
Potential causal effects of BAD and food allergies on COVID-19 infection/severity
In addition to asthma, the present MR analyses also found a potential protective causal effect of BAD on critically ill COVID-19 (OR = 0.93). It is worth noting that the BAD is a combination of asthma, hay fever and atopic dermatitis. Because of the negative causal correlations between atopic dermatitis and COVID-19 infection/severity in our MR analyses, the causal effect is more likely driven by the protective causal effect of asthma. The exact effect and interaction roles are worthy to be further investigated if GWAS summary statistics for the three types of individual allergic diseases are publicly available. Besides that, we identified risk causal effects of food allergies on COVID-19 phenotypes (OR > 1.01), particularly for shrimp allergy on all three COVID-19 phenotypes (OR > 1.02). Although the findings could not be further validated due to a lack of another available GWAS dataset, our MR analyses partially explained the molecular mechanism underlying the causal effects. The shrimp allergy not only showed significant effects on multiple immune-related cells and hematology traits but also have a marginal association with ACE2 protein expression level. The close relationships between food allergies and immune-related cells obtained supports from one previous study [63]. In cord blood, infants who developed food allergies showed a higher ratio of monocyte/CD4+ T cells, which would secrete higher amounts of inflammatory cytokines, such as IL-1β, IL-6, and TNF-α that further suppressed IL-2 expression by CD4+ T cells. With the decrease of IL-2, multiple inflammatory cytokines further decreased the number of activated natural regulatory T cells and generated an IL-4 expression non-classical Th2 phenotype, which subsequently induced a series of allergic reactions. Our MR analyses for the first time found potentially causal correlations between BAD and food allergies with COVID-19 phenotypes, calling more attention to other types of ADs.
Limitations
Several limitations in the present study should be elucidated. First, we could not obtain another GWAS dataset for the two kinds of food allergies, which weakened the reliability of causal correlations between food allergies and COVID-19 phenotypes. Although GWAS analyses of food allergies were performed in eastern Asian populations, the analyses comprised multiple steps to control for population stratification to minimize its impact on the final results as much as possible. In addition, our MR analyses identified a series of significant causal effects of immune-related cells on ADs and COVID-19 phenotypes, but a part of analyses only depended on one valid genetic variant, which is likely to decrease the reliability of results and the acting direction. However, the Wald ratio method was designed to estimate genetic correlation of diseases/traits of interest and could ensure the reliability of outcome in this scenario [64]. Finally, although we obtained two GWAS datasets of ACE2 protein expression that could directly reflect the effective level of the protein that exerts a biological function in peripheral blood tissue, it only indirectly reflected ACE2 protein expression in lung or respiratory tract tissue. Notably, our MR analyses about the potential correlations of asthma on ACE2 protein expression were supported by previous observational studies that conducted in the respiratory tract [4, 59].
Conclusions
Our systematic two-sample and bidirectional MR analyses consistently indicated the existence of a unidirectional protective causal effect of asthma on COVID-19 infection/severity, which was confirmed by another GWAS dataset and could be further explained by the underlying molecular mechanisms. In brief, asthma symptoms were causally associated with abnormal fluctuations of immune-related cells and hematological traits in peripheral blood. Our findings also implied potential causal effects of food allergies on COVID-19 phenotypes, which require further testing in observational or fundamental studies. The findings of our MR analyses suggest that in addition to focusing on the widely studied asthma for infection/severity to develop effective preventive measures for asthma patients, we should also pay attention to the other allergic diseases, which may also induced by the dysfunction of host’s immune responses.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its Additional files.
Abbreviations
- ADs:
-
Allergic diseases
- BAD:
-
Board allergic disease
- ADE:
-
Allergic dermatitis
- SA:
-
Shrimp allergy
- PA:
-
Peach allergy
- A2:
-
COVID-19 critical ill cases
- B2:
-
COVID-19 hospitalized patients
- C2:
-
COVID-19 infection cases
- MR:
-
Mendelian randomization
- GWAS:
-
Genome-wide association study
- OR:
-
Odds ratio
- RCT:
-
Randomized controlled trial
- SNP:
-
Single nucleotide polymorphism
- ACE2:
-
Angiotensin-converting enzyme 2
- HB:
-
Hemoglobin
- MCV:
-
Mean corpuscular volume
- PLT:
-
Platelet
- HT:
-
Hematocrit
- MCH:
-
Mean cell hemoglobin
- NEUT:
-
Neutrophil
- MONO:
-
Monocyte
- EOS:
-
Eosinophil
- BASO:
-
Basophil
- LYMPH:
-
Lymphocyte
- RAPS:
-
Robust adjusted profile score
- IVW:
-
Inverse-variance weighted
References
Huang C, Wang Y, Li X, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020;395(10223):497–506.
Guan WJ, Ni ZY, Hu Y, et al. Clinical characteristics of coronavirus disease 2019 in China. N Engl J Med. 2020;382(18):1708–20.
Epidemiology Working Group for NCIP Epidemic Response, Chinese Center for Disease Control and Prevention. The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19) in China. Zhonghua Liu Xing Bing Xue Za Zhi. 2020;41(2):145–51.
Wang R, Bikov A, Fowler SJ. Treating asthma in the COVID-19 pandemic. Thorax. 2020;75(10):822–3.
Tordesillas L, Berin MC, Sampson HA. Immunology of food allergy. Immunity. 2017;47(1):32–50.
Yang JM, Koh HY, Moon SY, et al. Allergic disorders and susceptibility to and severity of COVID-19: a nationwide cohort study. J Allergy Clin Immunol. 2020;146(4):790–8.
Shahgolzari M, Yavari A, Arjeini Y, et al. Immunopathology and immunopathogenesis of COVID-19, what we know and what we should learn. Gene Rep. 2021;25: 101417.
Zhang JJ, Dong X, Cao YY, et al. Clinical characteristics of 140 patients infected with SARS-CoV-2 in Wuhan, China. Allergy. 2020;75(7):1730–41.
Goyal P, Choi JJ, Pinheiro LC, et al. Clinical characteristics of Covid-19 in New York City. N Engl J Med. 2020;382(24):2372–4.
Williamson EJ, Walker AJ, Bhaskaran K, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature. 2020;584(7821):430–6.
Johnston SL. Asthma and COVID-19: is asthma a risk factor for severe outcomes? Allergy. 2020;75(7):1543–5.
Kow CS, Capstick T, Hasan SS. Are severe asthma patients at higher risk of developing severe outcomes from COVID-19? Allergy. 2021;76(3):959–60.
Morais-Almeida M, Barbosa MT, Sousa CS, Aguiar R, Bousquet J. Update on asthma prevalence in severe COVID-19 patients. Allergy. 2021;76(3):953–4.
Song J, Zeng M, Wang H, et al. Distinct effects of asthma and COPD comorbidity on disease expression and outcome in patients with COVID-19. Allergy. 2021;76(2):483–96.
Lighter J, Phillips M, Hochman S, et al. Obesity in patients younger than 60 years is a risk factor for COVID-19 hospital admission. Clin Infect Dis. 2020;71(15):896–7.
Vardavas CI, Nikitara K. COVID-19 and smoking: a systematic review of the evidence. Tob Induc Dis. 2020;18:20.
Grasselli G, Zangrillo A, Zanella A, et al. Baseline characteristics and outcomes of 1591 patients infected with SARS-CoV-2 admitted to ICUs of the Lombardy Region, Italy. JAMA. 2020;323(16):1574–81.
Smith GD, Ebrahim S. Data dredging, bias, or confounding. BMJ. 2002;325(7378):1437–8.
Jones DS, Podolsky SH. The history and fate of the gold standard. Lancet. 2015;385(9977):1502–3.
Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG, Consortium E-I. Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30(7):543–52.
Lawlor DA, Harbord RM, Sterne JA, Timpson N, Davey SG. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. 2008;27(8):1133–63.
Chapman SJ, Hill AV. Human genetic susceptibility to infectious disease. Nat Rev Genet. 2012;13(3):175–88.
Covid-Host Genetics Initiative. Mapping the human genetic architecture of COVID-19. Nature. 2021;600(7889):472–7.
Rao S, Baranova A, Cao H, Chen J, Zhang X, Zhang F. Genetic mechanisms of COVID-19 and its association with smoking and alcohol consumption. Brief Bioinform. 2021;22(6): bbab284.
Zhang F, Baranova A. Smoking quantitatively increases risk for COVID-19. Eur Respir J. 2022;60(6):2101273.
Rao S, Lau A, So HC. 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(7):1416–26.
Olafsdottir TA, Theodors F, Bjarnadottir K, et al. Eighty-eight variants highlight the role of T cell regulation and airway remodeling in asthma pathogenesis. Nat Commun. 2020;11(1):393.
Ferreira MA, Vonk JM, Baurecht H, et al. Shared genetic origin of asthma, hay fever and eczema elucidates allergic disease biology. Nat Genet. 2017;49(12):1752–7.
Paternoster L, Standl M, Waage J, et al. Multi-ancestry genome-wide association study of 21,000 cases and 95,000 controls identifies new risk loci for atopic dermatitis. Nat Genet. 2015;47(12):1449–56.
Khor SS, Morino R, Nakazono K, et al. Genome-wide association study of self-reported food reactions in Japanese identifies shrimp and peach specific loci in the HLA-DR/DQ gene region. Sci Rep. 2018;8(1):1069.
Sliz E, Huilaja L, Pasanen A, et al. Uniting biobank resources reveals novel genetic pathways modulating susceptibility for atopic dermatitis. J Allergy Clin Immunol. 2022;149(3):1105-1112.e9.
Demenais F, Margaritte-Jeannin P, Barnes KC, et al. Multiancestry association study identifies new asthma risk loci that colocalize with immune-cell enhancer marks. Nat Genet. 2018;50(1):42–53.
Ferreira MA, Hottenga JJ, Warrington NM, et al. Sequence variants in three loci influence monocyte counts and erythrocyte volume. Am J Hum Genet. 2009;85(5):745–9.
Ferreira MA, Mangino M, Brumme CJ, et al. Quantitative trait loci for CD4:CD8 lymphocyte ratio are associated with risk of type 1 diabetes and HIV-1 immune control. Am J Hum Genet. 2010;86(1):88–92.
Suhre K, Arnold M, Bhagwat AM, et al. Connecting genetic risk to disease end points through the human blood plasma proteome. Nat Commun. 2017;8:14357.
Yang Z, Macdonald-Dunlop E, Chen J, et al. Genetic landscape of the ACE2 coronavirus receptor. Circulation. 2022;145(18):1398–411.
Purcell S, Neale B, Todd-Brown K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75.
Zhao Q, Wang J, Hemani G, Bowden J, Small DS. Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. Ann Statist. 2020;48(3):1742–69.
Hemani G, Tilling K, Davey SG. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 2017;13(11): e1007081.
Pierce BL, Burgess S. Efficient design for Mendelian randomization studies: subsample and 2-sample instrumental variable estimators. Am J Epidemiol. 2013;178(7):1177–84.
Corne JM, Marshall C, Smith S, et al. Frequency, severity, and duration of rhinovirus infections in asthmatic and non-asthmatic individuals: a longitudinal cohort study. Lancet. 2002;359(9309):831–4.
De Fabregues O, Vinas J, Palasi A, et al. Ammonium tetrathiomolybdate in the decoppering phase treatment of Wilson’s disease with neurological symptoms: a case series. Brain Behav. 2020;10(5): e01596.
Qiu LJ, Yin KJ, Pan GX, Ni J, Wang B. Non-causal effects of asthma on COVID-19 susceptibility and severity. Front Genet. 2022;12: 762697.
Baranova A, Cao H, Chen J, Zhang F. Causal association and shared genetics between asthma and COVID-19. Front Immunol. 2022;13: 705379.
Owen JA, Punt JA, Stranford SA. Kuby immunology. New York: W. H. Freeman and Company; 2013. p. p40.
McBride JA, Striker R. Imbalance in the game of T cells: what can the CD4/CD8 T-cell ratio tell us about HIV and health? PLoS Pathog. 2017;13(11): e1006624.
Aiello A, Farzaneh F, Candore G, et al. Immunosenescence and its hallmarks: how to oppose aging strategically? A review of potential options for therapeutic intervention. Front Immunol. 2019;10:2247.
Turner JE. Is immunosenescence influenced by our lifetime “dose” of exercise? Biogerontology. 2016;17(3):581–602.
Huang W, Berube J, McNamara M, et al. Lymphocyte subset counts in COVID-19 patients: a meta-analysis. Cytometry A. 2020;97(8):772–6.
Wang F, Nie J, Wang H, et al. Characteristics of peripheral lymphocyte subset alteration in COVID-19 pneumonia. J Infect Dis. 2020;221(11):1762–9.
Lee SY, Kim SJ, Kwon SS, et al. Distribution and cytokine production of CD4 and CD8 T-lymphocyte subsets in patients with acute asthma attacks. Ann Allergy Asthma Immunol. 2001;86(6):659–64.
Kosonen J, Luhtala M, Viander M, Kalimo K, Terho EO, Savolainen J. Candida albicans-specific lymphoproliferative and cytokine (IL-4 and IFN-gamma) responses in atopic eczema dermatitis syndrome. Evidence of CD4/CD8 and CD3/CD16+CD56 ratio elevations in vitro. Exp Dermatol. 2005;14(7):551–8.
Ren X, Wen W, Fan X, et al. COVID-19 immune features revealed by a large-scale single-cell transcriptome atlas. Cell. 2021;184(7):1895-1913.e19.
Sun X, Hou T, Cheung E, et al. Anti-inflammatory mechanisms of the novel cytokine interleukin-38 in allergic asthma. Cell Mol Immunol. 2020;17(6):631–46.
Lambrecht BN, Hammad H, Fahy JV. The cytokines of asthma. Immunity. 2019;50(4):975–91.
Keidar S, Kaplan M, Gamliel-Lazarovich A. ACE2 of the heart: from angiotensin I to angiotensin (1–7). Cardiovasc Res. 2007;73(3):463–9.
Chamsi-Pasha MA, Shao Z, Tang WH. Angiotensin-converting enzyme 2 as a therapeutic target for heart failure. Curr Heart Fail Rep. 2014;11(1):58–63.
Kuba K, Imai Y, Ohto-Nakanishi T, Penninger JM. Trilogy of ACE2: a peptidase in the renin-angiotensin system, a SARS receptor, and a partner for amino acid transporters. Pharmacol Ther. 2010;128(1):119–28.
Peters MC, Sajuthi S, Deford P, et al. COVID-19-related genes in sputum cells in asthma. Relationship to demographic features and corticosteroids. Am J Respir Crit Care Med. 2020;202(1):83–90.
Jackson DJ, Busse WW, Bacharier LB, et al. Association of respiratory allergy, asthma, and expression of the SARS-CoV-2 receptor ACE2. J Allergy Clin Immunol. 2020;146(1):203-206.e3.
Kimura H, Francisco D, Conway M, et al. Type 2 inflammation modulates ACE2 and TMPRSS2 in airway epithelial cells. J Allergy Clin Immunol. 2020;146(1):80-88.e8.
Ramasamy S, Subbian S. Critical determinants of cytokine storm and type I interferon response in COVID-19 pathogenesis. Clin Microbiol Rev. 2021;34(3):e00299-20.
Zhang Y, Collier F, Naselli G, et al. Cord blood monocyte-derived inflammatory cytokines suppress IL-2 and induce nonclassic “T(H)2-type” immunity associated with development of food allergy. Sci Transl Med. 2016;8(321):321ra8.
Burgess S, Small DS, Thompson SG. A review of instrumental variable estimators for Mendelian randomization. Stat Methods Med Res. 2017;26(5):2333–55.
Acknowledgements
We could like to express our sincere thanks for great suggestions from Prof. Hon-Cheong So and Prof. Mary Miu Yee Waye at the Chinese University of Hong Kong in writing and revision of this article.
Funding
This study was supported by Fujian Provincial Natural Science Foundation Youth Innovation Project (Grant No. 2021J05050), Fujian Province Joint Innovation Project (Grant No. 2021Y9030), Research start-up funds for high-level talents from Fujian Medical University (Grant No. XRCZX2021009), and National Natural Scientific Foundation of China (Grant Nos. 32072336 and 31871720). The funding sources had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and the decision to submit the manuscript for publication.
Author information
Authors and Affiliations
Contributions
STR and GML designed and supervised the study. XTC, SZ, and JWJ extracted the raw GWAS summary data from public databases. XTC and SZ contributed to the data analysis and interpretation. XTC, XYH, STR, and GML contributed to the drafting and revision of the article. The corresponding authors attest that all listed authors met the authorship criteria and that no others meeting the criteria were omitted. All listed authors read and approved the final version of the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors have no conflicts of interest to declare.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
Supplementary Methods: Pre-processing of GWAS summary statistics. Supplementary Results: Table S1. Basic characteristics of GWAS summary statistics for COVID-19, various types of allergic diseases, ACE2 protein expression, and peripheral blood hematological traits and immune-related cells. Table S2. Overall MR analyses with strong genetic instruments (p < 5E-08) for causal effects of COVID-19 infection/severity on five types of allergic diseases. Table S3. Overall MR analyses with strong genetic instruments (p < 5E-08) for causal effects of COVID-19 infection/severity or five types of ADs on peripheral blood hematological traits and immune-related cells. Table S4. MR analyses with strong genetic instruments (p < 5E-08) for causal effects of asthma and shrimp allergy on peripheral blood ACE2 protein expression level.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Chen, Xt., Zhi, S., Han, Xy. et al. A systematic two-sample and bidirectional MR process highlights a unidirectional genetic causal effect of allergic diseases on COVID-19 infection/severity. J Transl Med 22, 94 (2024). https://doi.org/10.1186/s12967-024-04887-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12967-024-04887-4