Abstract Several studies reveal that allergic rhinitis (AR) is a significant risk factor of systemic lupus erythematosus (SLE). However, studies investigating the common pathogenesis linking AR and SLE are lacking. Our study aims to search for the shared biomarkers and mechanisms that may provide new therapeutic targets for preventing AR from developing SLE. [36]GSE50223 for AR and [37]GSE103760 for SLE were downloaded from the Gene Expression Omnibus (GEO) database to screen differentially expressed genes (DEGs). The Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed to explore the functions of shared DEGs. Hub genes were screened by cytoHubba (a plugin of Cytoscape) and validated in another two datasets. Gene set enrichment analysis (GSEA) and single-sample Gene set enrichment analysis (ssGSEA) algorithm were applied to understand the functions of hub gene. ENTPD1 was validated as a hub gene between AR and SLE. GSEA results revealed that ENTPD1 was associated with KRAS_SIGNALING_UP pathway in AR and related to HYPOXIA, TGF_BETA_SIGNALING and TNFA_SIGNALING_VIA_NFKB pathways in SLE. The expression of ENTPD1 was positively correlated with activated CD8 T cell in both diseases. Thus, ENTPD1 may be a novel therapeutic target for preventing AR from developing SLE. Keywords: ENTPD1, CD4^+T cells, Allergic rhinitis, Systemic lupus erythematosus, Bioinformatics Subject terms: Immunology, Biomarkers, Diseases, Rheumatology, Risk factors Introduction Systemic lupus erythematosus (SLE) is a systemic autoimmune disease that is characterized by various clinical manifestations ranging from mild skin lesion to disastrous organ damage (e.g., ischemic heart disease, lupus nephritis and pulmonary embolism)^[38]1. The prevalence of SLE greatly varies in worldwide regions, such as 84.8, 47, 37.6 and 3299.5 cases per 100,000 persons in San Francisco, France, Chian and Kenya, respectively^[39]2–[40]5. Moreover, SLE is the leading death reason among young women^[41]6. Thus, early detection, early diagnosis and effective treatment of SLE are warranted. Currently, the clinical treatments for SLE include corticosteroids, immunomodulators, cytotoxic-immunosuppressants, and biologics such as rituximab^[42]7. Nevertheless, there are some limitations to current therapies. For example, the consumption of B cells induced by rituximab may result in the absence of persistent response to rituximab. Thus, further investigations into new therapeutic strategies are necessitated. Allergic rhinitis (AR) is a chronic nasal disorder mainly mediated by immunoglobulin E (IgE), whose mechanism involves various kinds of immune cells^[43]8. The symptoms of AR including nasal itching, nasal obstruction, sneezing and rhinorrhea have bothered 10–40% of the population globally^[44]9. Superficially, AR and SLE seem to be no correlation, but accumulated pieces of research studies suggest that AR is a significant risk factor of SLE^[45]10–[46]13. A meta-analysis study reviewed 1 cohort and 7 case–control studies, and indicated that individuals with AR had a 1.36-fold risk of SLE compared with individuals without AR^[47]13. These intriguing findings inspired us to investigate the potential common mechanisms between AR and SLE in search of novel biological targets. CD4^+T cells are mainly composed of type 1 T helper (Th1) cells, Th 2 cells, Th9 cells, Th17 cells, T follicular helper (Tfh) cells, and regulatory T (Treg) cells^[48]14. They play a vital role in the mechanisms of both AR and SLE. The dysregulation of Th1/Th2 cells, Th17/Treg cells imbalance and cytokines produced by Th2 and Th17 cells lead to occurrence and development of SLE^[49]15. Furthermore, AR is an allergic reaction mainly driven by Th2 cells^[50]8. Recent studies indicate that Tfh cells and Th17 cells are also involved in the mechanism of AR^[51]16,[52]17. Considering the important function of CD4^+T cells in both diseases, we then focused on the shared pathogenesis of AR and SLE based on CD4^+T cells. Rapid development of microarray technology has enabled scientific researchers to obtain the expression profiles of thousands of genes in a short time, contributing to a better understanding of disease mechanisms at the gene level. In this study, we hypothesized that AR and SLE had a shared pathogenesis and tested the assumption by analyzing genes related to CD4^+T cells. Finally, we identified common pathways, shared hub gene, and potential chemicals targeting the hub gene, potentially providing a new insight into clinical therapies for both diseases. Materials and methods Collection of datasets The array datasets were obtained from the Gene Expression Omnibus (GEO) database. With the aim to comprehensively acquire the datasets for AR and SLE regarding CD4^+T cells, we first employed the search terms “((allergic rhinitis) and T cell) and Homo sapiens” and “((systemic lupus erythematosus) and T cell) and Homo sapiens” to search AR datasets and SLE datasets, respectively. Next, the following selection criteria were applied: (1) The datasets should be CD4^+T cells expression profiling. (2) The datasets should include cases and controls. (3) The stimulated allergen to CD4^+T cells in selected datasets of AR should be the same. Finally, two array datasets for AR ([53]GSE50223 and [54]GSE44960) and two array datasets ([55]GSE103760 and [56]GSE10325) were collected. In particular, both [57]GSE50223 and [58]GSE44960 were divided into four groups: patients challenged with diluent (PD) group, patients challenged with allergen (PA) group, healthy controls challenged with diluent (HD) group, and healthy controls challenged with allergen (HA) group. Additionally, [59]GSE103760 and [60]GSE10325 were composed of SLE patients and healthy controls. Identification of shared DEGs between AR and SLE The four datasets were first normalized by the function of “normalizeBetweenArrays” during the preprocessing procedure. Then, the “hgu133a.db” R package was used to convert gene probe ID into gene symbol name, and probes that do not correspond to the genes were removed. Next, the expression value of genes with several gene probes was averaged by the function of “aggregate”. Finally, the limma R package were utilized to search for the differentially expressed genes (DEGs) in [61]GSE50223 and [62]GSE103760. In order to identify DEGs more comprehensively, we set P value < 0.05 and FC (fold change) > 1.5 or < 2/3 to identify upregulated and downregulated DEGs, respectively^[63]18,[64]19. The R package “ggplot2” was applied for drawing the volcano plots of the two datasets. Then, the shared DEGs between AR and SLE were visualized with the use of the web tool Draw Venn Diagrams (https:// bioinformatics.psb.ugent.be/webtools/Venn/). GO and KEGG analyses Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis^[65]20–[66]22 were undertaken to better understand the functions of shared DEGs via the “clusterProfiler” R package. GO analysis contained three terms: biological process (BP), cellular component (CC), and molecular function (MF). These outcomes were presented as bubble plots through “ggplot2” R package. PPI network establishment and candidate hub genes selection The protein–protein interaction (PPI) network was established based on STRING database ([67]https://cn.string-db.org/) and visualized by Cytoscape software. The potential hub genes were selected by cytoHubba, a plugin of Cytoscape, which contains local-based scoring methods and global-based scoring methods for exploring key nodes of PPI. Due to the heterogeneity of topological features analyzed by different scoring methods, we utilized two local-based methods (MCC and Degree) and two global-based methods (EcCentricity and BottleNeck) in cytoHubba to search for candidate hub genes^[68]23. The top 5 genes ranked by four methods were chosen and these genes were intersected based on Draw Venn Diagrams tool (https:// bioinformatics.psb.ugent.be/webtools/Venn/). Validation and ROC curves of hub gene in two independent datasets The mRNA expressions of candidate hub genes above were verified in [69]GSE44960 for AR and [70]GSE10325 for SLE. The R package “limma” was performed to determine whether these genes were differentially expressed in datasets. Additionally, the Comparative Toxicogenomics Database (CTD, http:// ctdbase.org/) includes information concerning gene-disease, chemical-disease, and chemical-gene/protein interactions^[71]24. Data from CTD was utilized to validate the relationships between hub gene and diseases. In order to evaluate the diagnostic efficacy of hub gene, we drew receiver operating characteristic (ROC) curves to calculate the area under curve (AUC) and the cutoff value of [72]GSE44960 and [73]GSE10325. Study population The participants were recruited from the Changhai Hospital. These participants were divided into four groups: healthy controls (n = 8), AR without SLE (n = 9), SLE without AR (n = 8), and SLE with AR (n = 5). The diagnosis of AR was based on Allergic Rhinitis and its Impact on Asthma (ARIA) guidelines^[74]8, as follows: (1) having at least one symptom of sneezing, rhinorrhea, nasal itching, and nasal obstruction; (2) having the signs of pale turbinate and nasal mucosa, and nasal secreta. (3) having at least one positive allergen by serum immunoglobulin E test. The diagnosis of SLE was based on the guideline^[75]25, including evaluation of seven clinical domains (constitutional condition, haematological condition, neuropsychiatric condition, mucocutaneous condition, serosal condition, musculoskeletal condition and renal condition) and three immunological domains (antiphospholipid antibodies, low complements, anti-Smith and anti-dsDNA). Every domain had the corresponding score and the total score was added up by the scores of domains. The diagnosis of SLE was made if the antinuclear antibodies were positive and total score was ≥ 10. For the AR without SLE group, they met the following criteria: (1) being diagnosed as AR; (2) no autoimmune diseases and other atopic diseases except AR; (3) no mental and malignant disorders. For the SLE without AR group, the inclusion criteria were as follows: (1) being diagnosed as SLE; (2) no atopic diseases and other autoimmune diseases except SLE; (3) no mental and malignant diseases. For the SLE with AR group, they matched the following criteria: (1) being diagnosed as AR and SLE; (2) no other autoimmune diseases and other atopic diseases except SLE and AR; (3) no mental and malignant disorders. The inclusion criteria of healthy controls were as follows: (1) no atopic diseases; (2) no autoimmune diseases; (3) no mental and malignant diseases. The clinical features of participants were summarized in Table [76]S1. Validation of ENTPD1 expression by quantitative real-time PCR analysis and flow cytometry Blood samples were collected from patients and controls to verify the hub gene expression. There were four groups in our study: healthy controls (n = 8), AR without SLE (n = 9), SLE without AR (n = 8), and SLE with AR (n = 5). The CD4 ^+ T cells were purified from peripheral blood mononuclear cells (PBMCs) by using CytoSinct™ CD4 Nanobeads (GeneScript, L00863-1) and a magnetic-activated cell separation system. The total RNA was extracted from CD4 ^+ T cells according to the manufacturer’s instructions (Vazyme, R711). Next, the Evo M-MLV RT Kit (Accurate Biology, AG11705) was utilized to transcribe mRNA into cDNA. Finally, qRT-PCR analysis of cDNA by SYBR Green Premix Pro HS qPCR Kit (Accurate Biology, AG11701) was performed. ENTPD1 gene primers: 5′-AGCAGCTGAAATATGCTGGC-3′(forward), 5′-GAGACAGTATCTGCCGAAGTCC-3′ (reverse). GAPDH gene primers: 5′-CATGAGAAGTATGACAACAGCCT-3′ (forward), 5′-AGTCCTTCCACGATACCAAAGT-3′ (reverse). GAPDH was considered as an internal reference gene. The gene relative expression was calculated by the 2^−ΔΔCT method. Flow cytometry was used to determine the percentage of CD11b^-CD4^+CD39^+T cells. Firstly, PBMCs were isolated from fresh peripheral blood. Subsequently, after the incubation of FITC anti-human CD11b (Biolegend, 982,614), CD4 (Biolegend, 980,802) and FITC anti-human CD39 antibody (Biolegend, 328,205), PBMCs was analyzed by flow cytometry using a FACSCalibur flow cytometer (BD Biosciences). CD11b^-CD4^+ CD39^+T cells percentage was analyzed using FlowJo software. FITC anti-human CD11b was used to exclude the interference induced by myeloid cells. Function analyzed by GeneMANIA GeneMANIA is an online tool used for predicting interaction networks and functions of specific genes, which can provide information physical interactions, co-expression, colocalization, genetic interactions, pathway, and shared protein domains^[77]26. Hence, we applied GeneMANIA to reveal the function of ENTPD1. Gene set enrichment analysis and immune infiltration analysis Gene set enrichment analysis (GSEA) was utilized to investigate the pathways that ENTPD1 was involved. In our study, we first divided the PA in [78]GSE50223 into high- and low-expression groups based on ENTPD1 median expression. Similarly, SLE patients in [79]GSE103760 were divided into groups with high and low expression of ENTPD1. Secondly, the gene sets of “H: hallmark gene sets” were downloaded from the Molecular Signatures Database (http:// [80]www.broad.mit.edu/gsea/msigdb/index.jsp) as reference sets. Finally, we performed GSEA between high and low expression groups using the R package “clusterProfiler”. The GSEA results included enrichment score (ES), normalized enrichment score (NES), P value, and false discovery rate (FDR). Enriched pathways with |NES|> 1, P value < 0.05, and FDR < 0.05 were defined significant. Single-sample Gene set enrichment analysis (ssGSEA) algorithm, an extension method of GSEA, is used to evaluate the immune cell infiltration in each sample using gene expression profile^[81]27. In our study, the ssGSEA algorithm was used to assess the immune infiltration in [82]GSE50223 for AR and [83]GSE103760 for SLE, respectively. Spearman’s correlation analysis was performed to evaluate the associations between immune cells and ENTPD1 and to calculate the correlation indices. Identification of putative chemical compounds As previously mentioned, the CTD contains data regarding chemical-gene interactions. The chemical-ENTPD1 interactions was investigated in the CTD. Ethical approval and consent to participate This study was approved by Changhai Hospital Ethics Committee (CHEC2023-264). All participants signed the informed consents. All methods were carried out in accordance with relevant guidelines and regulations. Statistical analysis In this research study, statistical analyses were implemented in R software (version 4.2.2). The correlations between the expression of ENTPD1 and immune cell were analyzed by using the Spearman method. Continuous variables were presented as mean ± standard deviation. Unless otherwise stated, a P value < 0.05 was considered statistically significant. Results The shared DEGs between AR and SLE By using limma R package, we identified 174 DEGs between PA and HA groups in [84]GSE50223, including 123 downregulated DEGs and 51 upregulated DEGs (Fig. [85]1a). However, only 8 genes differentially expressed between PD and HD groups (Table [86]S2), which might because that DEGs were generated by CD4^+T cells mainly after exposure to the allergen. Therefore, DEGs between the PA and HA groups were used for further study. Meanwhile, as shown in Fig. [87]1b, 464 DEGs containing 9 downregulated and 455 upregulated genes were screened between SLE patients and healthy controls in [88]GSE103760. Subsequently, a comparative analysis was applied to confirm the shared DEGs between AR and SLE. Finally, a total of 11 shared genes (CYBB, ENTPD1, FCER1G, FGL2, HMMR, IFI30, IFNG, LY96, MCOLN2, RNASE6, TYMS) were determined for the further analyses (Fig. [89]1c). Figure 1. [90]Figure 1 [91]Open in a new tab Identification of shared DEGs between AR and SLE. (a) The volcano plot illustrates 174 DEGs between the AR patients and healthy controls. (b) The volcano plot illustrates 464 DEGs between the SLE patients and healthy controls. (c) The Venn diagram shows a total of 11 shared DEGs between AR and SLE. DEG differentially expressed gene, AR allergic rhinitis, SLE systemic lupus erythematosus. Functional enrichment analyses of shared DEGs To further understand the functional overview of the shared DEGs, GO analysis and KEGG pathway enrichment were utilized. The top 10 GO terms of BP, CC, and MF were depicted in Fig. [92]2a and Table [93]S3. In terms of BP, the shared DEGs were mainly enriched in positive regulation of cytokine production, tumor necrosis factor production, antigen processing and presentation, and T cell activation and differentiation. In terms of CC, the shared DEGs were mostly distributed in secretory granule membrane, ficolin-1-rich granule, tertiary granule, and tertiary granule membrane. Besides, these DEGs involved in MF terms of superoxide-generating NAD(P)H oxidase activity, oxidoreductase activity, IgG binding, Toll-like receptor binding, folic acid binding, sequence-specific mRNA binding, nucleoside diphosphate phosphatase activity, mRNA regulatory element binding translation repressor activity, and calcium-release channel activity. Furthermore, KEGG analysis showed that shared DEGs were highly associated with antigen processing and presentation, HIF-1 signaling pathway, natural killer cell mediated cytotoxicity, and necroptosis (Fig. [94]2b and Table [95]S4). Figure 2. [96]Figure 2 [97]Open in a new tab Functional enrichment analyses of the shared DEGs. (a) GO terms including BP, CC, and MF were used for functional enrichment clustering analysis on shared DEGs. (b) KEGG pathway enrichment analyses were performed on the shared DEGs^[98]20–[99]22. DEG differentially expressed gene, GO gene ontology, BP biological process, CC cellular component, MF molecular function, KEGG Kyoto Encyclopedia of Genes and Genomes. PPI network construction and selection of candidate hub genes Using the STRING database, we established a PPI network to investigate the potential associations of proteins encoded by shared DEGs. As shown in Fig. [100]3a, the PPI network includes 10 nodes and 13 edges, with a PPI enrichment P value of 9.99E-10. Next, the PPI network was imported into Cytoscape software and analyzed in cytoHubba. Finally, the intersection of top 5 genes ranked by four methods (MCC, Degree, EcCentricity and BottleNeck) confirmed 3 candidate hub genes: IFNG, TYMS, and ENTPD1 (Fig. [101]3b and Table [102]1). Additionally, the expression level of these three genes were higher in both AR and SLE patients than those in corresponding controls (Figure S1). Figure 3. [103]Figure 3 [104]Open in a new tab PPI network and identification of candidate hub genes. (a) PPI network of the shared DEGs constructed by STRING. (b) Identification of 3 candidates for hub genes by four ranked methods in Cytoscape plugin cytoHubba. PPI protein–protein interaction, DEG differentially expressed gene, STRING Search Tool for the Retrieval of Interacting Genes. Table 1. Identification of candidate hub genes in shared DEGs ranked by four methods of cytoHubba. Category Ranked methods in cytoHubba Gene symbol MCC Degree EcCentricity BottleNeck FCER1G FCER1G IFNG FGL2 CYBB CYBB ENTPD1 ENTPD1 IFNG IFNG FGL2 FCER1G ENTPD1 ENTPD1 TYMS IFNG TYMS TYMS CYBB TYMS [105]Open in a new tab Bold gene symbols were the top 5 shared hub genes ranked by four methods in cytoHubba. DEG differentially expressed gene; MCC Maximal clique centrality. Validation and efficacy evaluation of hub gene in two independent datasets The expression of candidate hub genes obtained above were analyzed in [106]GSE44960 for AR and [107]GSE10325 for SLE to validate our findings. Notably, the expression of ENTPD1 gene was significantly upregulated in both AR and SLE. The expression level of IFNG was higher in SLE; however, such expression difference was not observed in AR. In addition, no significance of TYMS expression was found in both AR and SLE (Fig. [108]4). Furthermore, ENTPD1 was found to be correlated with AR and SLE in CTD (Table [109]S5). Hence, ENTPD1 was confirmed as the hub gene of AR and SLE. Figure 4. [110]Figure 4 [111]Open in a new tab Validation of candidate hub genes. (a) The violin plots show the expression levels of candidate hub genes in [112]GSE44960 for AR. (b) The violin plots show the expression levels of candidate hub genes in [113]GSE10325 for SLE. AR allergic rhinitis, SLE systemic lupus erythematosus. ROC curves were drawn to assess the diagnostic efficacy of ENPD1. ROC curve analysis revealed that the AUC was 0.74 for ENTPD1 in prediction of AR, with an optimal cutoff point of 4.532 (sensitivity = 70% and specificity = 90%) (Fig. [114]S2a). Furthermore, the AUC was 0.76 for ENTPD1 in prediction of SLE, with an optimal cutoff point of 6.531 (sensitivity = 100% and specificity = 57%) (Fig. [115]S2b). These results suggested that ENTPD1 might be a promising marker for diagnosing AR and SLE. Validation of ENTPD1 expression by qRT-PCR and flow cytometry The clinical characteristics of individuals included in this study were summarized in Table [116]S1. As shown in Fig. [117]5a, the expression level of ENTPD1 is higher in AR, SLE and SLE with AR than that in healthy control. Additionally, the ENTPD1 expression was the highest in SLE with AR group among 4 groups. Consistently, the percentage of CD11b^-CD4^+ CD39^+T cells in PBMCs was significantly lower in healthy control than that in other 3 groups (healthy control, 0.93 ± 0.37%; AR, 2.36 ± 0.97%; SLE, 2.69 ± 0.86%; SLE with AR, 8.46 ± 1.87%). The SLE with AR group had the highest percentage of CD11b^-CD4^+ CD39^+T cells (Fig. [118]5b,c). Figure 5. [119]Figure 5 [120]Open in a new tab Validation of ENTPD1 expression. (a) The relative mRNA expression of ENTPD1 in four groups: healthy control, AR, SLE, and SLE with AR groups. (b) (c) The percentage of CD11b^-CD4^+CD39^+T cells in PBMCs of four groups: healthy control, AR, SLE, and SLE with AR groups. *P < 0.05, ***P < 0.001. AR allergic rhinitis, SLE systemic lupus erythematosus, PBMC peripheral blood mononuclear cell. Functional annotations of ENTPD1 analyzed by GeneMANIA After confirming ENTPD1 as the hub gene, we next investigated the potential functions of ENTPD1 via GeneMANIA (Fig. [121]6). The hub gene is located in the central circle, which is surrounded by predicted related genes. A sum of 20 targeted genes were predicted to construct the networks with physical interactions of 77.64%, co-expression of 8.01%, colocalization of 3.63%, genetic interactions of 2.87%, pathway of 1.88%, and shared protein domains of 0.6%. As shown in Fig. [122]6, the networks implicate these genes are mainly enriched in pyrimidine-containing compound catabolic, nucleoside diphosphate metabolic, organophosphate catabolic, nucleoside phosphate catabolic and nucleobase-containing small molecule biosynthetic process, and nucleoside-diphosphatase activity. Figure 6. [123]Figure 6 [124]Open in a new tab The gene–gene interaction network and functions for ENTPD1 based on the GeneMANIA database. GSEA results of ENTPD1 To uncover the effect of ENTPD1 on AR and SLE, GSEA was utilized to explore the significantly different pathways involved in these two diseases between the low and high ENTPD1 expression groups. GSEA result reveals that ENTPD1 is significantly associated with KRAS_SIGNALING_UP pathway in AR (Fig. [125]S3a). In addition, ENTPD1 is strongly related to HYPOXIA, TGF_BETA_SIGNALING, and TNFA_SIGNALING_VIA_NFKB pathways in SLE (Fig. [126]S3b). The P-values, ES, NES and FDR for [127]GSE50223 and [128]GSE103760 were listed in Table [129]S6. Immune infiltration analyses of ENTPD1 We obtained the associations between the expression of ENTPD1 and immune cell infiltration in AR and SLE via the ssGSEA algorithm and spearman’s correlation analysis. The Fig. [130]S4 shows the activated CD4 T cell and effector memory CD4 T cell are more highly expressed in AR and SLE patients compared with controls. As revealed by Fig. [131]7a, the expression of ENTPD1 is positively related to activated CD8 T cell and eosinophil, and is negatively linked with regulatory T cell, macrophage, gamma delta T cell and activated dendritic cell in AR. Moreover, immune infiltration analysis was also performed in SLE, which indicated that ENTPD1 expression was positively associated with type 2 T helper cell, gamma delta T cell, memory B cell, regulatory T cell, activated CD4 T cell, effector memory CD8 T cell, plasmacytoid dendritic cell, activated CD8 T cell, and activated B cell (Fig. [132]7b). Clearly, the expression of ENTPD1 was positively correlated with activated CD8 T cell in both AR and SLE. Figure 7. [133]Figure 7 [134]Open in a new tab Distribution of immune cell infiltration in AR and SLE. The lollipop plots reveal the relationships between ENTPD1 expression level and immune cell subtypes in AR (a) and SLE patients (b). AR allergic rhinitis, SLE systemic lupus erythematosus. Candidate chemical compounds targeting at ENTPD1 in the CTD Since ENTPD1 might play a key role in the pathogenesis of AR and SLE, we next investigated to search for the potential compounds that regulated the expression of ENTPD1. Table [135]S7 describes candidate chemical compounds which could regulate ENTPD1 mRNA expression. There were 10 chemical compounds decreasing ENTPD1 expression and 13 chemical compounds increasing ENTPD1 expression. Discussion SLE usually causes catastrophic organ damage and is the primary reason of young female^[136]6. However, current therapies have no benefit in decreasing the mortality^[137]28. Recent research has revealed that the prevalence of SLE was significantly higher in AR patients compared with the general population^[138]13. Effective AR therapy may prevent or delay the onset of SLE. Hence, identifying the common pathogenesis between AR and SLE is vital for facilitating treatment and prevention strategies. CD4^+T cells are involved in not only the occurrence of AR but also the onset of SLE. Nevertheless, no study has focused on exploring the common pathogenesis of AR and SLE based on gene expression profiling of CD4^+T cells. In our study, we identified shared DEGs and pathways by analyzing transcriptomic data of CD4^+T cells in both AR and SLE. Further analyses indicated that ENTPD1 played a crucial role in AR and SLE, and was a potential biomarker of these two diseases. In this work, we found 11 shared DEGs between AR and SLE based on the GEO database, and performed GO and KEGG analyses on these shared DEGs. GO analysis of the shared DEGs showed that these genes were mainly enriched in positive regulation of cytokine production, tumor necrosis factor production, antigen processing and presentation, and T cell activation and differentiation in terms of BP. These results suggested that both AR and SLE were closely associated with immune response and inflammation, which were consistent with previous studies^[139]15,[140]17. Notably, in terms of MF, the shared DEGs were related to Toll-like receptor binding. Several studies have suggested that Toll-like receptor played a crucial role in the onset of AR and SLE by mediating the production of inflammatory cytokine^[141]29 or the loss of B cell tolerance to autoantigens^[142]30. A small number of shared DEGs involved in GO analysis might lead to unstable results, but the results might give us some idea about the functions of these genes. KEGG analysis also revealed that shared DEGs were involved in pathways regarding immune responses. Additionally, KEGG analysis indicated that shared DEGs were associated with HIF-1 signaling pathway and necroptosis, implying that these pathways might be crucial to the common pathogenesis of AR and SLE. Using bioinformatics analysis, we identified that ENTPD1 was the shared hub gene between AR and SLE. ENTPD1, ectonucleoside triphosphate diphosphohydrolase 1, encodes an ectoenzyme CD39 which hydrolyzes extracellular ATP (eATP) to AMP. CD39 is able to convert the eATP-induced pro-inflammatory responses into anti-inflammatory status in immune actions^[143]31. Notably, ATP level elevated in SLE patients, activating inflammasome and releasing cytokines related with disease consequently. As previous studies described, CD39 regulated the immune and inflammatory response in SLE via catalyzing ATP into adenosine, a potent immune regulator^[144]32,[145]33. And Becker et al. reported that the expression of CD39 increased in SLE compared with healthy controls and suggested that the increased CD39 expression could help control inflammation as a compensatory mechanism^[146]33. Moreover, Treg-mediated suppression was absolutely absent in SLE patients with CD39 defect compared with controls and SLE without CD39 defect^[147]34. Further experimental studies showed that deficiency of CD39 led to clinically overt pulmonary hemorrhage, potentiated splenomegaly and autoantibody production, and accelerated the release of neutrophil extracellular traps by neutrophils in SLE model compared with corresponding controls^[148]35. In lupus nephritis mice, Zheng et al. revealed that CD39 expression in iTreg cells was crucial to induce immune tolerance, and modulated the proliferation and differentiation of T cells, thus inhibiting inflammation. Besides, CD39 was effective in treating lupus nephritis and the inhibition of CD39 expression led to the failure of amelioration in lupus nephritis^[149]36,[150]37. Based on these, it could be concluded that ENTPD1/CD39 mediated its beneficial effects by suppressing inflammation in SLE. Previous study showed that genetic variant in ENTPD1 was related with an increased risk of AR^[151]38. The Treg cells had an effect on suppressing immune responses by increasing CD39 expression after intranasal allergen immunotherapy^[152]39, which indirectly suggested the anti-inflammatory function of CD39. Although few studies have conducted in-depth exploration regarding the role of ENTPD1/CD39 in AR, the role of ENTPD1/CD39 in allergic asthma was relatively widely explored. CD39 was found to be associated with reduced airway inflammation and hyperreactivity in allergic asthma^[153]40. Overexpression of ENTPD1 could also mitigate airway eosinophilia, mucin deposition, and Th2 cytokines production^[154]41, while ENTPD1 knockout in allergic asthma mice aggravated airway inflammation through increasing the cytokines level and eosinophil recruitment^[155]42. Likewise, allergic asthma model administered with apyrase, which has ectoenzymatic activity identical to CD39, was observed to has attenuated peribronchial eosinophilic infiltration and decreased levels of Th2 cytokines in the bronchoalveolar lavage fluid^[156]43. These results indicated that ENTPD1/CD39 might participate in the pathogenesis of allergic asthma. Moreover, Lu et al. reported that Treg cells transfected by ENTPD1 could secrete higher levels of IL10 and TGF-β, which are essential in modulating immune system. Furthermore, up-regulation of ENTPD1 enhanced the immunosuppressive function and FOXP3 expression of Treg cells, contributing to the maintenance of immune homeostasis^[157]44,[158]45. It is well known that Treg cells play a leading role in the mechanism of AR^[159]17. Given that the concept of “one-airway-one-disease”^[160]46, we speculated that ENTPD1/CD39 might be involved in the mechanism of AR via regulating the levels of mucin and Th2 cytokines, and eosinophil recruitment. In summary, the above studies provided compelling evidence for the crucial role played by ENTPD1 in AR and SLE. As described previously, GeneMANIA can predict interaction networks of ENTPD1^[161]26. GeneMANIA results showed that 20 genes were related to ENTPD1, which might provide a new insight into investigating how ENTPD1 participate in the mechanism of AR and SLE. Although CD4^+T cells are pivotal in AR and SLE, the role of CD8 ^+ T cells in these two diseases could not be underestimated. CD8 ^+ T cell, a subset of T lymphocytes, was reported to reduce allergen-induced airway eosinophilic infiltration by secreting IFN-γ^[162]47. Moreover, CD8^+T cells might mitigate inflammation in AR condition^[163]48. Li et al. found that the frequency of CD8^+T cells in peripheral blood was higher in SLE patients than in controls, and suggested that CD8^+T cells might be of importance in maintaining peripheral tolerance and suppressing self-reactive^[164]49. In addition, CD8^+T cells were observed to expand and could suppress pathogenic CD4^+T cells in multiple sclerosis, an autoimmune disease^[165]50. According to our result that ENTPD1 was positively associated with activated CD8^+T cells in both AR and SLE, it can be deduced that ENTPD1 might prevent the development of AR and SLE by activating CD8^+T cells. Notably, the expression of ENTPD1 was not correlated with the activated CD4^+T cells. The contradictory result might due to the relatively small size of study population. Further immune infiltration analysis with a larger sample scale is necessary. Additionally, the immune cell infiltration patterns of these two diseases were not similar. Previous study reported that other factors such as epigenetics and gut microbiota were involved in the mechanism of SLE^[166]51, which might contribute to the difference of immune patterns. The importance of KRAS pathway to the development of human tumors, especially pancreatic ductal adenocarcinomas is hard to overstate^[167]52. However, the association between the KRAS pathway and AR has not been investigated. In our work, GSEA result showed that ENTPD1 was involved in the KRAS_SIGNALING_UP pathway in AR, which might provide evidence for us to explore this association. In SLE, the low ENTPD1 expression group was correlated with TGF_BETA_SIGNALING and TNFA_SIGNALING_VIA_NFKB pathways, which have been recognized as key characteristics of the occurrence of SLE^[168]53,[169]54. Interestingly, recent studies showed that hypoxia-inducible factor-1α, a key protein of hypoxia pathway, participated in the pathogenesis of SLE by regulating microRNA, which supported our GSEA result that hypoxia pathway was involved in SLE. There are some limitations to our studies. Firstly, the shared pathways have not been verified by basic experiments. Thus, more experimental studies are warranted as direct evidence to demonstrate the exact role of these pathways on AR and SLE. Secondly, though we confirmed as many DEGs as possible in our research, some genes correlating with diseases might be excluded because they did not match the selection criteria. Further larger sample studies are warranted to identify DEGs comprehensively. Thirdly, we were unable to analyze the association between hub gene and clinical characteristics due to a lack of such data. Lastly, we could not obtain SLE dataset including male and female individuals from the GEO database, which might because women are 9 times more frequently to develop SLE than men^[170]55. Thus, the SLE datasets included in this study were consisted of female individuals. Notably, previous study has reported that the distribution and expression of CD4^+T cells was not associated with sex^[171]56. Hence, sex difference is unlikely to significantly affect the reliability of our findings. Conclusions In conclusion, we identified ENTPD1 as a shared hub gene in AR and SLE by applying bioinformatic analysis based on the GEO database. The ENTPD1 expression increased in AR and SLE patients, and the diagnostic efficacy of ENTPD1 in two diseases was validated. Our study may provide a new therapeutic target for preventing AR from developing SLE. Supplementary Information [172]Supplementary Information.^ (754.9KB, docx) Acknowledgements