Introduction

The Coronavirus Disease 2019 (COVID-19) pandemic, caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), has imposed a huge burden on global health and the economy [1,2,3,4]. According to the World Health Organization, there have been 767,518,723 cumulative COVID-19 cases and 6,947,192 cumulative COVID-19 deaths worldwide as of July 3rd, 2023 (https://covid19.who.int/). A collection of risk and protective factors have been identified which correlate to susceptibility, hospitalization, and severity of COVID-19 [5,6,7,8,9]. Among these, the correlations between COVID-19 and allergic diseases have been repeatedly reported [10,11,12,13]. Allergic asthma (AA) is commonly defined as asthma associated with sensitization to aeroallergens [14]. Approximately 11% of the general population suffer from AA [15]. Previous studies reveal that AA decreases the risk of hospitalization and death in COVID-19 patients compared with non-AA [16, 17]. Patients with AA have a lower probability of SARS-CoV-2 test positivity than those with non-AA [16]. However, the associated mechanism between these two diseases remains largely unknown.

Immune cells play crucial roles in regulating AA attacks and defending against COVID-19 infections [18, 19]. As the most efficient antigen-presenting cells, dendritic cells (DCs) link innate and adaptive immunity [20]. There are five main subtypes of DCs in human lung tissue: conventional DC1 (cDC1), conventional DC2 (cDC2), plasmacytoid DC (pDC), monocyte-derived DC (Mo-DC), and Pre-DC [21]. Among these, cDC2 is derived from Pre-DCs, accounting for the most prominent subset of DCs [22]. The molecular markers of cDC2 include CD1C, ITGAM, CLEC10A, and others [21, 22]. In AA lungs, the number of cDC2s is significantly up-regulated [23]. cDC2s can promote T helper (Th) cell differentiation and regulate humoral immunity in asthma [19, 24]. After SARS-CoV-2 infection, cDC2s are activated and accumulate in patients lungs [25]. The activated cDC2s can induce a positive immune response in COVID-19 [20]. Furthermore, in some severe infections, the number and maturation of DCs are reduced, indicating DC dysfunction or failure [20]. Although some studies explored the immune landscape in AA and COVID-19, to our knowledge, the shared immune genetic architectures have not been investigated.

Herein, we used 18 GWAS datasets to examine the causal links between allergic diseases and COVID-19 by employing bi-directional two-sample Mendelian randomization (MR) analysis and considering blood parameters as risk factors. Among these, AA and white cell count stood out. We further performed machine learning using 3 bulk transcriptome datasets to screen the shared discriminative immunocytes in AA and COVID-19, and identified cDC2s. Moreover, the quantity, activity, intercellular communication, and differentiation of cDC2s were investigated in 2 single-cell transcriptome datasets, and the antigen presentation process of cDC2s was exhibited using molecular docking analysis. In summary, this study proposes for the first time, that cDC2s link the genetic causal association between AA and COVID-19.

Materials and methods

Ethics approval

An overview of the study design is presented in Fig. 1. As this study utilized pre-existing publications and public databases, no further ethical approval or consent was necessary.

Fig. 1
figure 1

Study design overview. Genomic and transcriptomic data of COVID-19 and allergic diseases were collected from public platforms. The causal correlations between allergy and COVID-19 were detected using the Mendelian randomization (MR) analysis, and blood parameters were used as mediating variables. Based on MR results, the study further applied bulk-tissue transcriptome, single-cell transcriptome, and molecular docking analyses to explore the potential cellular and molecular mechanism that bridges the allergic asthma (AA) and COVID-19

Data source

GWAS data for COVID-19, including 3 phenotypes (susceptibility, severity, and hospitalization), were extracted from the COVID-19 Host Genetics (hg) GWAS meta-analyses (Round 7, April 2022) (https://www.covid19hg.org/results/r7/). GWAS data for allergy, consisting of 8 phenotypes, and for blood, comprising 7 parameters, were obtained from the FinnGen database (https://www.finngen.fi/en/access_results) and GWAS Catalog/ IEU Open GWAS Project (https://www.ebi.ac.uk/gwas, http://gwas.mrcieu.ac.uk), separately. All GWAS was restricted to individuals of European ancestry. Moreover, transcriptome data of blood and airway mucosa in COVID-19 and AA cohorts were accessed from Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/). A total of 630 samples from 5 datasets (including 3 bulk and 2 single-cell datasets) were collected. Bulk-tissue datasets included blood samples from AA/COVID-19 patients and healthy controls (HCs). Single-cell datasets incorporated airway mucosa samples from AA patients and AA controls (AC) with or without allergen (Ag) challenge, as well as airway mucosa samples from severe (S) COVID-19 patients and HCs. For detailed information see Additional file 2:: Table S1. Additionally, 3 protein structures were downloaded from the Protein Data Bank (PDB, https://www1.rcsb.org/), including the COVID-19 spike glycolprotein (6VXX), TLR1 (6NIH), and OX40L-OX40 (TNFSF4-TNFRSF4) complex (2HEV).

MR analysis

The bi-directional two-sample MR study followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement [26]. The study design fulfilled 3 crucial assumptions [26]: (1) instrumental variables should be strongly correlated to exposures; (2) instrumental variables should be independent of underlying confounders; (3) instrumental variables influence the outcome only through the risk factors (Fig. 1). For forward analysis, allergic diseases were designated as exposure factors, while COVID-19 was set as the outcome factor. Reverse analysis was conducted for the causal relationships identified by the forward analysis, where COVID-19 (hospitalization) was designated as the exposure factor and allergic disease (AA) was set as the outcome factor. Next, instrumental variable SNPs were selected based on the following criteria: (1) high exposure correlation: P < 5 × 10–6; (2) high independence: linkage disequilibrium (LD) r2 < 0.001; (3) high statistical strength: F > 10; (4) no direct outcome correlation: P > 5 × 10–6; (5) no palindromic and incompatible structure: harmonization analysis [27]. The harmonization analysis between exposure and outcome data was conducted using the function ‘harmonise_data’ of the R package ‘TwoSampleMR’. MR-Steiger filtering [28] was applied to exclude SNPs which implied reverse causal direction. Moreover, since some phenotypes had a limited number of SNPs to meet the criteria, we relaxed their correlation threshold to P < 5 × 10–5. [29] Furthermore, 4 MR analytical algorithms were applied to evaluate the causal effects between exposures and outcomes. Among these, inverse-variance weighted (IVW) functioned as the primary method, while the other three, namely MR Egger, Weighted Median (WM), and MR PRESSO, served as supplementary tools [30]. Statistical significance was determined, only when the P-value of IVW was < 0.05 and the OR directions of the other three algorithms were consistent with that of IVW.

Sensitivity analyses were performed, which included tests for pleiotropy and heterogeneity and leave-one-out analysis. Directional pleiotropy was measured by the Egger-intercept test [31] (no pleiotropy: P > 0.05). Heterogeneity was assessed through the Cochran Q test (no heterogeneity: P > 0.05) and funnel pot. Leave-one-out analysis was conducted to determine whether the causal correlation was driven by any single SNP (no driven: results are all on one side). Furthermore, to explore the mechanisms genetically linking allergy and COVID-19, we evaluated some underlying mediators, namely seven blood parameters that have been previously reported to exhibit causal correlations to COVID-19 [32,33,34,35]. Allergy phenotypes were treated as exposures, and blood parameters were designated as outcomes. Their MR analyses were conducted based on the methods mentioned above.

Bulk transcriptome analysis

The immunocyte proportion of blood in AA and COVID-19 cohorts was estimated using 4 independent algorithms: CIBERSORT [36], MCPCounter [37], ssGSEA [38], and xCell [39]. The discriminating abilities of immunocytes between AA/COVID-19 patients and HCs were evaluated and ranked through three machine learning algorithms: Boruta, Random Forest, and Support Vector Machine Recursive Feature Elimination (SVM RFE) [40, 41]. Next, we intersected the immunocytes that exhibited high discriminating abilities across three algorithms. The overlapping cells between the AA and COVID-19 were identified as target cells.

Single-cell transcriptome analysis

The data were initially normalized using the functions ‘NormalizeData’ and ‘ScaleData’, and its dimensionality was reduced by the functions ‘FindVariableFeatures’ and ‘RunPCA’ (R package ‘Seurat) [42]. Batch effects were corrected with the function ‘RunHarmony’ (R package ‘harmony’) [43]. Cells were clustered and visualized through the functions ‘FindClusters’ and ‘UMAP’, separately (R package ‘Seurat) [42]. Cell lineage identity was annotated according to the previous studies [23, 44]. DCs were further divided into different subtypes (pDC, cDC1, cDC2, Mo DC, and Pre DC) using sets of marker genes [21]. Furthermore, the degree of cell activation was assessed by quantifying energy metabolism using the R package 'scMetabolism' [45]. Intercellular crosstalk was calculated by the R package ‘CellChat’ [46]. Cell differentiation was evaluated through pseudotime analysis with the R package ‘monocle3′ [47] and stemness analysis with the R package ‘cytoTRACE’ [48]. Additionally, the corresponding genes and functional consequences of instrumental variable SNPs for allergy were obtained from the dbSNP website (https://www.ncbi.nlm.nih.gov/snp/). For these genes, functional enrichment was conducted using the Metascape website (https://metascape.org/), and their role in cell differentiation was assessed.

Molecular docking

Protein docking was conducted using the ZDOCK module of Discovery Studio software (version 2019, available at https://discover.3ds.com/discovery-studio-visualizer-download). The docking structures with the top ZDOCK score were selected. Their 3D interaction diagram was visualized using the PyMOL software (version 1.8, available at http://www.pymol.org/pymol), and the 2D interaction figure was generated using the LigPlot + software (version 2.2.8, available at https://www.ebi.ac.uk/thornton-srv/software/LIGPLOT/).

Statistical analysis

Statistical analysis was conducted using the R software (version 4.1.2, available at https://www.r-project.org/). The Shapiro-Wilks test was applied for the normality test. The two-tailed t-test (normal distribution) and Mann–Whitney test (non-normal distribution) were utilized to compare the continuous variables between the two groups, separately. Additionally, MR analysis was conducted using the R package ‘TwoSampleMR’. Forest plots were produced through the R package ‘forestplot’. Bubble, bar, and violin plots were generated by the R package ‘ggplot2’. Venn diagrams were drawn through the R package ‘VennDiagram’. NS indicates no statistical significance, while ‘*’ indicates statistical significance, namely P < 0.05.

Results

AA exhibits a negative causal correlation to COVID-19

Instrumental variable SNPs were screened for different allergy phenotypes and COVID-19 hospitalization, which has been detailed in Additional file 2: Table S2. The number of SNPs varied from 5 to 76, and their summary data are shown in Additional file 2: Tables S3-S11. Among the allergy phenotypes, IVW analysis showed that AA decreased the risk of COVID-19 hospitalization (OR = 0.968, 95% CI 0.940—0.997, P = 0.031). The result directions of the other 3 MR methods were consistent (Fig. 2A, B, and Additional file 2: Table S12). Additionally, no reverse causal relationship was detected from COVID-19 hospitalization to AA (Additional file 2: Table S14). We also failed to observe causal associations between the other allergy phenotypes and COVID-19 risk (Fig. 2A and Additional file 2: Table S12). To evaluate the robustness and reliability of the above results, we conducted multiple types of sensitivity analyses. MR-Egger intercept tests indicated that no horizontal pleiotropy existed (all P > 0.05) (Additional file 2: Tables S13 and S14). However, heterogeneity was observed in Cochran's Q test among some allergy and COVID-19 phenotypes, like the association from AA to COVID-19 (Q = 52.293, P = 0.030) (Fig. 2C and Additional file 2: Tables S13 and S14). Despite the presence of heterogeneity, the MR estimates in the current study were not invalidated, as the random-effect IVW method was employed [49, 50]. Moreover, the existence of heterogeneity did not introduce any pleiotropic bias to MR estimates. Furthermore, the leave-one-out analysis demonstrated that no single SNP drove the results, indicating that none of the estimates were violated (Additional file 1: Fig. S1A).

Fig. 2
figure 2

AA exhibits a negative causal correlation to COVID-19 hospitalization. A Bubble plot showing the causal correlation between allergic diseases and COVID-19 phenotypes. B Scatter plots displaying the causal correlation between AA and COVID-19 phenotypes. C Funnel plots assessing the heterogeneity of the correlation between AA and COVID-19 phenotypes

White cell counts bridge the correlation between AA and COVID-19

Based on logical reasoning, the correlation between AA and COVID-19 hospitalization may not be direct. To explore the underlying mediators linking AA and COVID-19 hospitalization, we evaluated the effect of AA on COVID-19-related blood variables through MR analysis. AA exhibited a positive correlation to white cell count (IVW, β = 0.010, 95% CI 0.003—0.017, P = 0.008), which was consistent with the direction of results obtained from the other three MR methods (Fig. 3). However, no significant association was observed between AA and other blood variables (Fig. 3). Furthermore, the sensitivity analysis did not reveal any evidence of pleiotropy but reflected the presence of heterogeneity, which was accepted under the random-effect IVW method (Fig. 3). The leave-one-out analysis further excluded significant impact from a single SNP on the MR estimates (Additional file 1: Fig. S1B). These data suggest that white cell count bridged AA and COVID-19.

Fig. 3
figure 3

White blood cell counts bridge the correlation between AA and COVID-19. Blood parameters were used as mediating variables between AA and COVID-19. Forest plots show the causal correlation between AA and blood parameters. The causal correlation between these blood parameters and COVID-19 has been identified in previous studies

DC is the most discriminative white cell in AA and COVID-19

Because the range of white cells remains too large for the mediator variables between AA and COVID-19, we tried to narrow the scope. We used 4 independent algorithms to identify the particular type of white blood cell that bridged AA and COVID-19. Overall, compared to healthy subjects, AA/COVID-19 patients had significantly higher levels of granulocytes and monocytes/macrophages, as well as lower levels of B cells and T cells in the blood (Fig. 4A and Additional file 2: Table S15). The expression of DCs was up-regulated in the blood of AA patients and down-regulated in the blood of COVID-19 patients (Fig. 4A and Additional file 2: Table S15). Subsequently, three machine learning algorithms (Boruta, Random Forest, and SVM RFE) were separately applied to identify the sixteen most discriminative immunocytes between AA/COVID-19 patients and healthy subjects (Fig. 4B and Additional file 2: Table S16–S18). Through the intersection, we identified DCs as the target cells (Fig. 4C).

Fig. 4
figure 4

Dendritic cells (DCs) are the most discriminative immune cell in AA and COVID-19. A Heatmap showing the expression of immune cells in AA and COVID-19. B Three machine learning algorithms separately identify the top 16 discriminative immune cells in AA and COVID-19. C Venn plots showing the intersected immunocytes identified by three machine learning algorithms for distinguishing AA/COVID-19 and controls

cDC2 is upregulated/activated in AA exacerbation and downregulated/activated in severe COVID-19

To further uncover the correlation between the two diseases at the cellular subtype level, we performed the single-cell RNA sequencing (scRNA-seq) analysis on the airway mucosa of AA and COVID-19 cohorts. First, AA and COVID-19 datasets were integrated and corrected for batch effects (Additional file 1: Fig. S2A and B). Eight types of cells were annotated, including T cells, macrophages, mast cells, epithelial cells, B cells, NK cells, DCs, and neutrophils (Fig. 5A). Their cell markers were confirmed (Additional file 1: Fig. S2C). For AA, AAs with Ag challenge had the highest proportion of DCs, followed by ACs with Ag challenge and AAs/ACs without Ag (samples at baseline). For COVID-19, patients with severe infection exhibited a lower proportion of DCs compared to HCs (Fig. 5B). These two results validated the findings regarding DCs in the bulk sequencing analysis. Subsequently, we further divided the DCs into 5 cellular subtypes, including pDCs, cDC2s, Mo DCs, Pre DCs, and cDC1s (Fig. 5C). Their cell markers were confirmed (Fig. 5D). When comparing the expression of DC subtypes in different groups of AA/COVID-19, only cDC2s exhibited statistically significant disparities (P < 0.05). In terms of AA, AAs with Ag challenge exhibited the highest proportion of cDC2s, followed by ACs with Ag challenge and AAs/ACs without Ag. Regarding COVID-19, patients with severe infection demonstrated a significantly lower proportion of cDC2s compared to HCs (Fig. 5E). These results suggest that the level of cDC2s is up-regulated in AA exacerbation and down-regulated in the severe COVID-19.

Fig. 5
figure 5

Conventional dendritic 2 cells (cDC2s) are increased/activated in AA exacerbation and downregulated/activated in severe COVID-19. A UMAP plot of total cells from AA and COVID-19, with each cell color coded for cell type. B The proportion of DC in AA/COVID-19 and controls. C UMAP plot of DC populations from AA and COVID-19, with each cell color coded for cell subtype. D Bubble plot showing marker genes for cell subtype annotation. E The proportion of DC subtypes in AA/COVID-19 and controls. F The level of energy metabolism of cDC2 in AA/COVID-19 and controls. G Intercellular communication from cDC2 to other cells. H The number of interactions from cDC2 to T cells in AA/COVID-19 and controls. Bln samples at baseline, Ag samples receiving antigen, AC allergic non-asthmatic controls, HC healthy controls, S patients with severe COVID-19 infection

Next, we investigated whether there were differences in the activity of cDC2s among different groups of AA/COVID-19. We found that the levels of energy metabolism were increased in AAs in response to Ag and were higher in AAs in comparison with ACs (Fig. 5F). Additionally, severe COVID-19 exhibited higher levels of energy metabolism than HCs (Fig. 5F). cDC2s have been found to possess the ability to present antigens to T cells [22]. Thus, we examined the cell communication between cDC2s and other cells. In the context of antigen exposure, cDC2s from AAs sent more signals to T cells compared to that from ACs. Moreover, cDC2s from COVID-19 patients trended towards increased signal transmission to T cells in comparison to HCs, despite lacking statistical significance due to limited sample size (Fig. 5G, H). Collectively, these data suggest the presence of more activated cDC2s in exacerbation status of AA and COVID-19.

EXOC4, ZBTB46, TLR1 and TNFSF4 variations promote the cDC2 differentiation and cDC2 mediates SARS-CoV-2 antigen presentation in AA

Because the causal correlation between AA and COVID-19 had been identified in the aforementioned results, we further explored the shared genetic mechanism. Within 36 instrumental variable SNPs of AA, a total of 25 had corresponding genes, including TLR1, ZBTB46, EXOC4, and TNFSF4. Among these, 16 SNPs exhibited the functional consequences of transcriptional regulatory region variants, and 13 SNP corresponding genes were highly expressed in DCs (min.pct > 0.01, log2FC > 0.01, and P < 0.05) (Fig. 6A and Additional file 2: Tables S19, S20). Metascape enrichment analysis showed that these 13 genes focused on AA (childhood asthma), virus (response to virus), and immune cells (negative regulation of leukocyte activation) (Fig. 6B).

Fig. 6
figure 6

EXOC4, ZBTB46, TLR1, and TNFSF4 variations promote cDC2 differentiation and cDC2 mediates SARS-CoV-2 antigen presentation in AA. A The basic characteristics of instrumental variable SNPs of AA. B Metascape functional enrichment based on corresponding genes of instrumental variable SNPs of AA. C, D The differentiation potential of Pre-DCs in AA/COVID-19 and controls. E Pseudotime trajectories of Pre-DCs, cDC1s, and cDC2s in the UMAP plot. F The differentiation correlation levels of ZBTB46 and EXOC4 in AA and controls. G The expression alterations of ZBTB46 and EXOC4 in cDC2s as related to the increase in Pseudotime. H, I 3D protein docking showing the process of SARS-CoV-2 antigen presentation of cDC2s. J The expression level of TNFSF4 and TLR1 between AA and controls. S spike protein

Next, we investigated the differentiation and antigen presentation of DCs. CytoTRACE analysis predicted a higher developmental potential for Pre-DCs from AA Ag than those from AC Ag (Fig. 6C, D). Pre-DCs from COVID-19 patients exhibited decreased differentiation capacity compared to those from HCs (Fig. 6C, D). Monocle analysis drew the differentiation trajectory from Pre-DC to cDC1/cDC2 in AAs and ACs (Fig. 6E). Differentiation correlations of those genes were presented along the pseudotime trajectories (Additional file 2: Table S20). Among these, both ZBTB46 and EXOC4 exhibited dramatically higher correlations with differentiation in AAs than ACs (Fig. 6F and Additional file 1: Table S20). As pseudotime advanced, ZBTB46 was up-regulated in both AAs and ACs, while EXOC4 displayed down-regulation in AAs but up-regulation in ACs (Fig. 6G). Moreover, protein docking analysis presented the process of antigen presentation mediated by cDC2s, including the identification of SARS-CoV-2 (protein docking between TLR1 of cDC2 and spike protein on SARS-CoV-2) and presentation to T cells (protein docking between TNFSF4 of cDC1 and TNFRSF4 on T cells). Interactive forces and distances were displayed (Fig. 6H, I, and Additional file 1: Fig. S3). When antigen was exposed, both TLR1 and TNFSF4 exhibited higher expression in AAs than ACs, indicating a more pronounced activation of cDC2-mediated SARS-CoV-2 antigen presentation in AAs (Fig. 6J). Collectively, the above findings uncovered that EXOC4, ZBTB46, TLR1, and TNFSF4 variations promote cDC2 differentiation and cDC2 mediated SARS-CoV-2 antigen presentation in AA.

Discussion

AA has been suggested to be a protective factor for COVID-19, with the potential mechanism being largely unknown [5, 16, 17]. In this study, we found that AA exhibits negative genetic correlations to COVID-19 hospitalization, which is mediated by white cell count. Among them, cDC2s are the most discriminative immune cells between AA/COVID-19 and controls. When diseases are exacerbated, cDC2s exhibit upregulation, activation, and more differentiation in AA, while they show downregulation, activation, and poorer differentiation in COVID-19. ZBTB46, EXOC4, TLR1, and TNFSF4 gene mutations promoted cDC2 differentiation and SARS-CoV-2 antigen presentation in AA.

The causal correlation between asthma and COVID-19 has been previously investigated. Wang et al. reported that asthma was negatively associated with susceptibility to and severe respiratory symptoms of COVID-19 [13]. Baranova et al. indicated that asthma decreased the risk both for COVID-19 hospitalization and infection [8]. Similarly, our study revealed a negative causal correlation from AA to COVID-19 hospitalization. However, the correlation between AA and COVID-19 susceptibility/severity was not observed. Such disparities might be due to the use of different datasets, which may generate varying instrumental variables and yield different MR results. Moreover, hospitalization represents a particularly severe condition. Therefore, in subsequent transcriptome analyses, we focused on the exacerbation of AA and COVID-19.

Cellular immunity plays many crucial roles in AA and COVID-19 diseases [18, 19]. In AA, epithelial cell-derived cytokines drive the type 2 immune response, which promotes the infiltration of lung tissue with eosinophils, Type 2 T helper (Th2) cells, inflammatory DCs, and others. The accumulation of these cells leads to the acute response and chronic features of AA [19]. In COVID-19, the number of immunocytes, including DC, NK, B, CD4+ T, CD8+ T, and other cells are dramatically altered. Defects in cytotoxicity of NK and T cells lead to cytokine storm, causing tissue damage and organ dysfunction [51]. In our study, AA exhibited a positive causal correlation to white cell count, which has been previously reported to possess causal protective effects in COVID-19 [35]. Accordingly, we inferred that the white cell count bridged the correlation between AA and COVID-19. Through machine learning, we further identified DCs, as the most discriminative immunocytes between AA/COVID-19 and controls. As the most efficient antigen-presenting cells, DCs connect innate immunity and adaptive immunity. After penetrating the submucosa, Ag or SARS-CoV-2 are captured by DCs, undergo processing, and are subsequently presented to T cells [19, 20]. Previous studies have reported that AA increased the number of DCs in the peripheral blood [52]. However, the quantity and maturation of DCs decreases in severe COVID-19 [20]. Similarly, we found that the expression of DCs was up-regulated in AA blood and airway mucosa and down-regulated in those of COVID-19.

DCs have five subtypes, mainly including cDC1s, cDC2s, pDCs, Mo-DCs, and Pre-DCs [21]. Among these, cDC2s constitute a significant proportion of the DCs [22]. Our single-cell analysis found that only cDC2s were differentially expressed in the airway mucosa between AA/COVID-19 exacerbation and controls. The number of cDC2s was upregulated in AA, whereas it was downregulated in COVID-19. Consistent with our findings, previous studies have revealed that AA patients had a higher percentage of cDC2s in bronchial tissues, while hospitalized COVID-19 patients had a lower proportion of cCD2s in blood [25, 53]. cDC2s are derived from the differentiation of Pre-DCs [21]. Accordingly, we conducted stemness analysis, and found that the Pre-DCs of AA exacerbation possessed higher differentiation potential, whereas the Pre-DCs of COVID-19 exacerbation had lower differentiation potential. The stemness disparity of pre-DCs may, at least partially, account for the number difference of cDC2s in these two diseases. In addition, pseudotime analysis revealed that the expression of AA mutation genes ZBTB46 and EXOC4 were highly correlated to cDC2 differentiation. Among them, the role of ZBTB46 in DC differentiation has been previously reported. ZBTB46, a zinc finger transcription factor, is a marker for conventional DCs and committed DC precursors [54, 55]. Satpathy et al. found that ZBTB46 overexpression could promote cDC development [56]. Wang et al. revealed that ZBTB46 characterizations were potentially involved in DC maturation and activation [57]. In addition, EXOC4 encodes the protein belonging to the exocyst complex component family. Yi eti al. reported that EXOC2 knockdown inhibits COVID-19 infection. However, the study of EXOC4 in DC differentiation is lacking and deserves further investigation.

cDC2s respond to multiple danger signals ranging from nucleotides to polysaccharides, and can promote a variety of immune responses, especially the CD4 + T cell response and Th cell response [20, 22]. In our study, more activated cDC2s were found in both exacerbated AA and severe COVID-19. These cells sent more signals to T cells, promoting the production of adaptive immune responses. However, it is worth noting that the roles of cDC2s in AA and COVID-19 were mainly opposite. In AA, cDC2s responded to Ag, which contributed to disease occurrence, while in COVID-19, cDC2s resisted the virus invasion, which inhibited the disease progression. Furthermore, we investigated the mutation genes of AA, and found that the TLR1 protein of cDC2s (encoded by AA mutation gene TLR1) could strongly bind with the spike protein of SARS-CoV-2, and the TNFSF4 protein of cDC2s (encoded by AA mutation gene TNFSF4) could bind with its receptor TNFRSF4 on T cells. Both TLR1 and TNFSF4 were more highly expressed in cDC2s of AA compared to those of controls, suggesting an increased capacity of cDC2s for antigen presentation in AA. Among these, TLR1, belonging to the Toll-like receptor (TLR) family, plays a crucial role in pathogen recognition and activation of innate immunity [58]. Previous studies have confirmed the TLR1 gene mutation is associated with asthma [59]. The TLR family can mediate antiviral DC responses to resist SARS-CoV-2 infection [60]. Additionally, TNFSF4, belonging to the tumor necrosis factor (TNF) ligand family, functions in the interaction between antigen-presenting cells and T cells [61]. Previous studies have confirmed that SNP polymorphisms in TNFSF4 are correlated to asthma [62]. TNFSF4 is primarily expressed on DCs, and its receptor TNFRSF4 is preferentially expressed by T cells [63]. TNFSF4/TNFRSF4 signaling plays a crucial role in resisting the infection of bacteria and viruses [64].

According to the study results, we proposed a potential theory explaining why patients with AA are less susceptible to COVID-19 hospitalization. For COVID-19, SARS-CoV-2 may inhibit the differentiation of Pre-DCs, reduce the number of cDC2s, attenuate the connection between the innate immune and adaptive immune system, and as a result, increase its ability to infect the body. AA might enhance the capacity of Pre-DC differentiation through EXOC4 and ZBTB46 gene mutation, and therefore increase the number of cDC2s, especially the number of active cDC2s. These cells highly express TLR1 and TNFSF4, recognize more SARS-CoV-2, and present more signals to T cells, which induces a stronger adaptive immune response to resist SARS-CoV-2 infection.

Our study has some limitations. First, the MR analysis was performed among European ancestry participants, while bioinformatic analysis was conducted among American, Asian, and European ancestry participants. Race differentiation belongs to confounding factors, which may affect the accuracy of results. Second, apart from blood parameters, there was a lack of analysis for other risk factors that bridged AA and COVID-19. Third, this study is lacking the experimental validation.

Conclusions

In conclusion, this study suggests that cDC2s link the genetic causality from AA to COVID-19. Future strategies of COVID-19 prevention, intervention, and treatment for the population may be stratified according to AA and guided with dendritic cell (DC)-based therapies.