Skip to main content

Identification and analysis of lipid metabolism-related genes in allergic rhinitis

A Correction to this article was published on 04 August 2023

This article has been updated



Studies have shown that the lipid metabolism mediator leukotriene and prostaglandins are associated with the pathogenesis of allergic rhinitis (AR). The aim of this study was to identify key lipid metabolism-related genes (LMRGs) related to the diagnosis and treatment of AR.

Materials and methods

AR-related expression datasets (GSE75011, GSE46171) were downloaded through the Gene Expression Omnibus (GEO) database. First, weighted gene co-expression network analysis (WGCNA) was used to get AR-related genes (ARRGs). Next, between control and AR groups in GSE75011, differentially expressed genes (DEGs) were screened, and DEGs were intersected with LMRGs to obtain lipid metabolism-related differentially expressed genes (LMR DEGs). Protein-protein interaction (PPI) networks were constructed for these LMR DEGs. Hub genes were then identified through stress, radiality, closeness and edge percolated component (EPC) analysis and intersected with the ARRGs to obtain candidate genes. Biomarkers with diagnostic value were screened via receiver operating characteristic (ROC) curves. Differential immune cells screened between control and AR groups were then assessed for correlation with the diagnostic genes, and clinical correlation analysis and enrichment analysis were performed. Finally, real-time fluorescence quantitative polymerase chain reaction (RT-qPCR) was made on blood samples from control and AR patients to validate these identified diagnostic genes.


73 LMR DEGs were obtained, which were involved in biological processes such as metabolism of lipids and lipid biosynthetic processes. 66 ARRGs and 22 hub genes were intersected to obtain four candidate genes. Three diagnostic genes (LPCAT1, SGPP1, SMARCD3) with diagnostic value were screened according to the AUC > 0.7, with markedly variant between control and AR groups. In addition, two immune cells, regulatory T cells (Treg) and T follicular helper cells (TFH), were marked variations between control and AR groups, and SMARCD3 was significantly associated with TFH. Moreover, SMARCD3 was relevant to immune-related pathways, and correlated significantly with clinical characteristics (age and sex). Finally, RT-qPCR results indicated that changes in the expression of LPCAT1 and SMARCD3 between control and AR groups were consistent with the GSE75011 and GSE46171.


LPCAT1, SGPP1 and SMARCD3 might be used as biomarkers for AR.


Allergic rhinitis (AR) is a noninfectious inflammatory disease mediated by IgE that affects approximately 10%~20% of the global population [1]. Exposure to inhaled allergens in susceptible individuals is a frequent precipitating factor for AR. The most common clinical symptoms include paroxysmal sneezing, nasal obstruction, rhinorrhea, and nasal itching, sometimes in association with conjunctivitis, such as eye itching and tearing. Persistent severe rhinitis may predispose patients toward asthma [2]. Moreover, AR patients often have decreased learning and work efficiency, impaired sleep and quality of life, and even psychological disorders such as depression, leading to a huge economic burden on society. For treatment, nasal corticosteroids, antihistamines, and leukotriene receptor antagonists are currently the most recommended drugs [3]. However, their long-term drug use leads to a range of side effects, including epistaxis and drowsiness. Furthermore, the sustained poor efficacy of available drugs causes recurrent illnesses. Thus, it is crucial to find effective therapeutic targets for AR treatment.

Lipids are composed of fats and lipoids and play an important role in different organelles as a second messenger for intracellular signaling [4]. Lipid metabolism refers to the digestion, synthesis, and disassembly of lipids, with the help of various enzymes related to the processing of substances necessary for the body to ensure normal physiologic function. Previous studies have shown that lipid metabolism-related genes (LMRGs) are associated with several systemic diseases. For instance, Li et al. found that LMRGs in circulation have good predictive value for early diagnosis of intervertebral disc degeneration (IDD) [5]. LMRGs are also involved in lung cancer development and might serve as biomarkers for lung cancer [6]. The lipid compound prostaglandins show immunological activity in allergic airway diseases, and meanwhile some prostaglandin D2 (PGD2) receptor antagonists (Ramatroban, AMG 853, etc.) have been shown potentially beneficial effects on allergic inflammation [7]. Leukotriene is a well-recognized lipid inflammatory mediator in allergic diseases, and leukotriene receptor antagonists are one of the major medications for AR [8]. In addition, AR patients have a high level of apolipoprotein in nasal mucus, which may be involved in lipid metabolism and have immunomodulatory properties [9]. As a bioactive lipid, platelet-activating factor (PAF) has the chemotactic property that amplifies mucosal inflammation and causes increased vascular permeability, promoting rhinorrhea and mucus secretion [10]. Peroxisome proliferator-activated receptor gamma (PPAR-γ) may promote activation of type 2 immune response and affect target gene expression by regulating lipid metabolism in allergic diseases [11]. The metabolite of eicosapentaenoic acid (EPA), 15-hydroxyeicosapentaenoic acid (15-HEPE) may reduce allergic rhinitis symptoms via intranasal injection by interacting with PPAR-γ and inhibition of mast cell degranulation [12]. Nevertheless, how LMRGs contribute to AR pathogenesis has remained unclear.

Rapid advances in ‘genomics’ and ‘omics’ yield vast amounts of data and gene-by-gene analysis is insufficient to meet the demands of biological cognition. Compared with laboratory experiments, bioinformatics technology has a great advantage in collating, analyzing, and visualizing large amount of biological information quickly and accurately by computer, with a good reproducibility [13]. Biomarkers and therapeutic targets of many diseases as well as tumors were identified by bioinformatics [14,15,16]. In this study, AR-related public datasets (GSE75011, GSE46171) were downloaded from Gene Expression Omnibus (GEO) database, and weighted gene co-expression network analysis (WGCNA) was performed to gain AR-related genes (ARRGs). Following the lipid metabolism-related differentially expressed genes (LMR DEGs) were screened, the protein-protein interaction (PPI) networks were constructed for hub genes. Enrichment analysis and receiver operating characteristic (ROC) curve evaluation were subsequently performed to identify diagnostic genes, and single-set gene set enrichment analysis (ssGSEA) was used to analyze the relationship between diagnostic genes and immune cells. Simultaneously, the predictive value of the nomogram based on the diagnostic genes was created and assessed for clinical decision making. This study was aimed to identify potential LMRGs with diagnostic value for AR, providing a potential treatment of AR patients.

Materials and methods

Collection of LMRGs

According to the previous literature [17], LMRGs were gained from the Reactome database and Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The LMRG1 within the “Metabolism of lipids” pathway (R-HSA-556,833) located in the Reactome database were extracted and overlapped with the LMRG2 in the 16 lipid metabolism-related gene sets from the KEGG database to obtain 750 LMRGs for subsequently analysis [5].

Microarray datasets acquisition in AR-related public datasets

Considering the important significance of skin or blood allergy testing in AR diagnosis, the gene expression profiles of GSE75011 from Th2-enriched CD4 + T cells in blood and GSE46171 from nasal epithelial cells within AR samples were selected and downloaded through GEO database. The Illumina HiSeq 2500 sequencing dataset GSE75011, with the largest number of samples was conducted as the training set, containing 15 control and 25 AR blood samples [18]. Besides, GSE46171 dataset containing three control and six AR samples of nasal mucosa was used as an external validation set to enrich disease information and avoid bias among different researches [19].

Identification of ARRGs

To gain ARRGs in GSE75011, WGCNA was performed [20]. First, the samples were clustered to remove outliers. Thereafter, the determination of soft threshold (β) was performed. Modules were segmented via dynamic tree cutting based on optimal β. Correlations were analyzed between modules and AR. The genes of the highest relevance module with |gene significance (GS)| > 0.3, |module membership (MM)| > 0.6, and P < 0.05 were defined as ARRGs [21, 22].

Screening and functional analysis of LMR DEGs

First, the read count data matrix in GSE75011 was transformed into log-counts per million (log-cpm) value to fit a linear model through voom approach in “limma” R package (version 3.48.3) for differentially expressed analysis. The mRNA expression levels between control and AR groups in the GSE75011 dataset were contrasted with the screening threshold of P < 0.05 via the “limma” R package (version 3.48.3) [23]. DEGs and LMRGs were taken to intersect to get LMR DEGs. Subsequently, enrichment analysis for LMR DEGs via Metascape database (P < 0.05) [24]. In addition, the online database WebGestalt was used to study the Disease Ontology (DO) function of LMR DEGs [25].

Creation of the PPI networks of LMR DEGs and screening of hub genes

The PPI network of LMR DEGs were created via Search Tool for the Retrieval of Interacting Genes (STRING) with confidence score more than 0.4 and visualized by Cytoscape. Subsequently, hub genes were further obtained by intersecting the top 30 genes which were calculated and ranked on the basis of among stress, radiality, closeness, and edge percolated component (EPC) within CytoHubba algorithm.

Screening of diagnostic genes

First, candidate genes were obtained by intersecting hub genes with ARRGs, and enrichment analysis were applied on them. Second, ROC curves of the candidate genes were mapped via “pROC” R package (version 1.18.0) in the GSE75011 and GSE46171 [26]. Candidate genes with area under the curve (AUC) ≥ 0.7 were regarded diagnostic genes. A nomogram was constructed with diagnostic genes for clinical utilize, and a calibration curve of the nomogram was drawn to verify its validity. The diagnostic worth of age, sex, and time point was assessed in GSE46171 via ROC curves.

Immune analysis

The ssGSEA algorithm was utilized to assess infiltrating richness of immune cells between control and AR groups in the training set. Differences of the control and AR groups were compared by the Wilcox test. In addition, immune cells relevant to the diagnostic genes were detected via Spearman algorithm between diagnostic genes and differential immune cells.

Analysis of clinical correlation

Relevance between diagnostic genes and clinical characteristics (age, sex, time point) was analyzed using Pearson correlation analysis in the “corrplot” R package [27].

Gene Set Enrichment Analysis (GSEA) of diagnostic genes

On the basis of the median value of the diagnostic gene expressions, the samples of GSE75011 were grouped into high and low expression groups. All genes in two expression groups were performed GSEA with |normalized enrichment score (NES)| > 1, nominal (NOM) P value < 0.05, and q < 0.25 [28].

Expression profiles of diagnostic genes in external validation datasets

To further demonstrate the reliability of the results above, expression levels of the diagnostic genes between control and AR samples were compared in the GSE75011 and GSE46171 datasets for external validation.

Patients and tissue preparation

Ten AR patients and ten patients without AR or significant underlying disease were selected from people visiting to Shanghai Changzheng Hospital. There were no marked variation in sex and age between the groups (Table 1). Blood samples were acquired from these patients with informed consent and carried out real-time fluorescence quantitative polymerase chain reaction (RT-qPCR). The Medical Ethics Committee of Shanghai Changzheng Hospital endorsed this study.

Table 1 Basic information of the patients


Firstly, we conducted the total RNA extraction utilizing TRIzol (Ambion, Austin, USA). Then, reverse transcription of total RNA to cDNA was made via First-strand-cDNA-synthesis-kit (Servicebio, Wuhan, China). RT-qPCR was made utilizing the 2xUniversal Blue SYBR Green qPCR Master Mix (Servicebio, Wuhan, China). Specific experimental steps were carried out on the basis of instructions. The primer sequences were showcased in Additional file 1. Internal reference gene was GAPDH. The 2−ΔΔCt method was utilized to calculate the expression of diagnostic genes [29]. Levels of expression of diagnostic genes between the control and AR groups were compared by the T test.

Statistical analysis

Statistical analysis was carried out through GraphPad Prism 5 and R software (version 4.2.0). P < 0.05 represented a significant difference. Differences between groups were analyzed via the Wilcoxon test.


Acquisition of ARRGs

To identify ARRGs, WGCNA was performed with the GSE75011 dataset. Sample clustering analysis showed no outliers in the dataset (Fig. 1A). The β was 4 (Fig. 1B), and each gene module contained a minimum of 100 genes. Three modules were eventually identified, each with a unique color (Fig. 1C-D). The blue module correlated markedly with AR (cor = -0.35, P = 0.03) (Fig. 1D). Finally, 66 ARRGs were gained and utilized for further analysis (Fig. 1E).

Fig. 1
figure 1

Gene co-expression network of AR. (A) Sample clustering analysis of the GSE75011 dataset. (B) Analysis of network topology for various soft-thresholding powers showed that the soft threshold was 4. (C) Clustering dendrograms and modules identified by WGCNA, where the minimum gene number was 100. (D) Module-trait relationships demonstrated that the blue module correlated markedly with AR (cor = -0.35, P = 0.03) in three modules, and each square contains the corresponding correlation and P value. (E) Correlation scatterplot of 66 ARRGs with |GS| > 0.3, |MM| > 0.6, and P < 0.05 in the blue module

Acquisition and functional enrichment of LMR DEGs

The 25 samples were standardized for the GSE75011 dataset and are presented as box plots in Fig. 2A-B. The volcano plot and heatmap show 1621 DEGs between the AR and control groups, including 810 upregulated and 811 downregulated genes (Fig. 2C-D). A total of 73 LMR DEGs (Additional file 2) were obtained by Venn analysis with LMRGs (750 genes) and DEGs (1621 genes), with a significant difference detected based on a heatmap (Fig. 2E-F). Enrichment analysis of the 73 LMR DEGs by Metascape showed a total of 334 functional pathways (Fig. 2G-H) to be related to the LMR DEGs, such as metabolism of lipids, lipid biosynthetic process, and sterol regulatory element-binding protein (SREBP) signaling. DO enrichment results showed that the LMR DEGs are significantly associated with 10 diseases, namely, xanthomatosis, increased serum pyruvate, decreased high-density lipoprotein, hypoalphalipoproteinemia, myoglobinuria, insulin-resistant diabetes, neonatal death, thin skin, myalgia and cardiomegaly (Fig. 2I).

Fig. 2
figure 2

Acquisition and analysis of LMR DEGs. (A-B) The 25 samples of the GSE75011 dataset were standardized and presented as box plots. (C) Volcano plot of 1621 DEGs between the AR and control groups (P < 0.05). (D) Heatmap of 1621 DEGs between the AR and control groups (P < 0.05). (E) Venn diagram of 73 LMR DEGs by overlapping LMRGs and DEGs. (F) Heatmap of the 73 LMR DEGs. (G-H) Enrichment analysis of the 73 LMR DEGs in Metascape was performed, and the associated interaction network is shown. (I) DO function analysis of LMR DEGs using the WebGestalt database

Acquisition of hub genes

The PPI network was created for LMR DEGs. As illustrated in Fig. 3A-B, SREBF1 interacts with multiple proteins, such as LPIN1, GPAM, and MED1. To identify the most important genes, the 22 genes common to the four algorithms were used as hub genes (Fig. 3C), and a PPI network of hub genes was created (Fig. 3D). The results showed that GPAM interacts with seven genes, namely, PPARG, NFYA, SREBF1, ACSL3, LPIN1, HMGCS1 and AACS.

Fig. 3
figure 3

Acquisition of hub genes. (A) The protein-protein interaction (PPI) network for 73 LMR DEGs through STRING. (B) The degree of connectivity of each gene in the PPI network demonstrated that SREBF1 interacted with multiple proteins. (C) Venn diagram of hub genes common to four algorithms. (D) PPI network of hub genes

Acquisition of diagnostic genes

Four candidate genes associated with diagnosis of AR were obtained by 66 ARRGs with 22 hub genes taking intersections: LPCAT1, SREBF1, SMARCD3, and SGPP1 (Fig. 4A). The four candidate genes were involved in 133 Gene Ontology (GO) items, including 114 GO biological process (BP), nine GO cellular component (CC) and 10 GO molecular function (MF), such as retina development in camera-type eye, npBAF complex, and transcription coregulator binding (Fig. 4B).

The diagnostic value of four candidate genes was assessed via ROC curve in GSE75011 and GSE46171. The AUC values for the three genes (LPCAT1, SMARCD3, and SGPP1) were greater than 0.7 in both datasets, suggesting that the three genes have diagnostic value for AR (Fig. 4C-D). The AUCs for age, sex and time point were 0.4492, 0.7166, and 0.4893, respectively, in GSE46171, revealing that sex might be a diagnostic factor for AR (Fig. 4E).

Fig. 4
figure 4

Acquisition of diagnostic genes. (A) Venn diagram of four candidate genes by taking the intersections of ARRGs and hub genes. (B) Gene Ontology (GO) enrichment analysis of the four candidate genes. (C-D) AUC curves of the three genes for diagnostic prediction in the GSE75011 and GSE46171 datasets. (E) AUCs of age, sex, and time point in GSE46171

Finally, the nomograms were created containing the three diagnostic genes in GSE75011 and GSE46171 (Fig. 5A-B), and the AUC values in both datasets were above 0.6 (Fig. 5C-D). The results demonstrated that the nomogram has good prediction ability for AR.

Fig. 5
figure 5

(A-B) The nomogram drawn based on the three diagnostic genes for the diagnostic efficacy of AR in GSE75011 and GSE46171. (C-D) Calibration curves of the nomogram in GSE75011 and GSE46171

Immuno-infiltration analysis in control and AR groups

Analysis of the percentage of immune cells by ssGESA in all samples showed the highest for T cells (Fig. 6A). Differences in infiltrating immune cells between the control and AR groups were illustrated by a violin plot (Fig. 6B). The results suggested that infiltration of regulatory T cells (Tregs) and T follicular helper cells (TFHs) was markedly lower in the AR group. There was significant relevance between SMARCD3 and TFHs. However, neither LPCAT1 nor SGPP1 correlated with differential immune cells (Tregs and TFHs); therefore, SMARCD3 was selected for further analysis (Fig. 6C-E).

Fig. 6
figure 6

Immuno-infiltration analysis in AR. (A) Analysis of the percentage of immune cells in all samples of GSE75011 using ssGESA. (B) Differences in immune cell infiltration between control and AR groups (Wilcox. Test). **P < 0.01. Lollipop chart demonstrating the correlation between diagnostic genes and immune cells, including (C) LPCAT, (D) SGPP1, and (E) SMARCD3

Correlation analysis of clinical features, enrichment analysis and infiltration analysis of SMARCD3

Pearson correlation analysis demonstrated that SMARCD3 was significantly associated with clinical characteristics (age and sex) (Fig. 7A-C). Then, GSEA for SMARCD3 was performed, revealing 256 GO enrichment (Additional file 3) and 33 KEGG pathways (Additional file 4) (Fig. 7D-E). Overall, SMARCD3 was involved in immune-related pathways, for instance, the B-cell receptor signaling pathway and T-cell receptor signaling pathway. Four immune cells displayed marked variations between the high and low expression groups, namely, macrophages, T helper cells, central memory T cell (Tcm), and TFH cells, reflecting the strong relevance between SMARCD3 and the immune microenvironment (Fig. 7F).

Fig. 7
figure 7

Comprehensive analysis of SMARCD3 in AR. (A-C) Pearson correlation analysis of SMARCD3 clinical characteristics (age, sex and time point). GSEA of SMARCD3 showed that 256 GO enrichment pathways (D) and 33 KEGG pathways (E) were enriched in AR. (F) Differences in immune cell infiltration between high and low SMARCD3 expression groups (Wilcox test). *P < 0.05; **P < 0.01

mRNA levels of diagnostic genes

The significant differences in expression of SGPP1, LPCAT1 and SMARCD3 between control and AR in GSE75011 and GSE46171 were clearly observed via visualized data (Fig. 8A-B). Moreover, the changes of the three genes expression were consistent in blood and nasal mucosal tissues, suggesting that these three genes are of high diagnostic value.

Fig. 8
figure 8

Expression levels of three diagnostic genes in online datasets. (A) GSE75011. (B) GSE46171. *P < 0.05; **P < 0.01

To verify the three diagnostic genes expression, we collected blood samples to assess mRNA expression levels of these genes via RT-qPCR. The expression trends of LPCAT1 and SMARCD3 were consistent with public databases, and the expression was lower in AR group (Fig. 9A-B). However, SGPP1 exhibited the opposite trend compared to the results of public database, possibly due to different experimental designs or analysis methods (Fig. 9C).

Fig. 9
figure 9

Validation of diagnostic genes expression by RT-qPCR. (A) SMARCD3. (B) LPCAT1. (C) SGPP1. *P < 0.05; ***P < 0.001


AR is an airway allergic disease with a high incidence, affecting billions of people in the world. Nevertheless, the effect of current therapies for AR is unsatisfactory due to its complex pathogenesis. LMRGs are involved in the maintenance of systemic physiology and play an important role in diverse diseases, especially in malignant tumors. Moreover, lipid-related inflammatory mediators such as prostaglandins and leukotrienes have been implicated in AR pathogenesis. To our knowledge, this is the first study to identify and analyze LMRGs in AR.

In this study, three key LMRGs most associated with AR, i.e., LPCAT1, SGPP1, and SMARCD3, were identified, all of which are protein-coding genes. As one of the lysophosphatidylcholine acyltransferase (LPCAT) family, the LPCAT1 protein is an enzyme essential for phosphatidylcholine metabolism and regulation of phosphatidylcholine composition [30]. LPCAT1 is also used in the prognosis of multiple tumors, such as breast cancer, colorectal cancer, and hepatocellular carcinoma [31,32,33]. Little is known about LPCAT1 in allergic diseases. One study reported that LPCAT1 downregulates eosinophilic inflammation in asthmatic mice [34]. In the current study, LPCAT1 was significantly lower in AR blood samples, consistent with published results, suggesting that it may be essential for AR pathogenesis. SGPP1 can catalyze degradation of S1P, who can regulate diverse biological processes, as a bioactive sphingolipid metabolite [35]. SGPP1 is considered to be closely related to several tumors, especially regarding chemoresistance and radioresistance [36, 37]. There are currently no reports about the function of SGPP1 in allergic diseases, and the results in the current study are the first to show significant downregulation of SGPP1 in both blood and nasal mucosa samples in AR patients; conversely, RT-qPCR using blood samples showed the opposite result, possibly due to different experimental designs or analysis methods. Thus, the effect of SGPP1 in AR is still unclear. SMARCD3 is a chromatin-remodeling factor and a member of the SWI/SNF family, which present helicase and ATPase activities and are crucial in the transcription process of certain genes. Its related pathways include the circadian clock and transcriptional activation of mitochondrial biogenesis [38]. SMARCD3 was found to be downregulated in AR patients in this study, but how it participates in disease processes remains to be explored.

Immuno-infiltration analysis refers to studying the composition and quantification of immune cells in diseases. In this study, TFHs were extremely significantly reduced in the AR group. TFHs are CD4+ T cells that specialize in helping B cells and are involved in a wide range of diseases. An increasing number of theories have concluded that the antigen-related IgE response depends on more TFHs than Th2 cells [39, 40]. There are few reports about TFH and SMARCD3. A microarray model system identified that the SMARCD3 gene is upregulated in T-cell acute lymphoblastic leukemia [41]. In this study, only SMARCD3 correlated with differential immune cells (Tregs and TFHs), and TFHs and SMARCD3 were downregulated simultaneously in AR patients. Hence, it is hypothesized that SMARCD3 participates in the differentiation of T cells.

GSEA was performed to further investigate the role of SMARCD3 in AR, the results of which showed significant enrichment in the adipocytokine signaling pathway, B-cell receptor signaling pathway, and chemokine signaling pathway, among others. The adipocytokine signaling pathway refers to a series of cascade events via autocrine or paracrine adipocytokines, such as leptin and adiponectin, by adipocytes in the body [42, 43]. This pathway is not only crucial for obesity, insulin resistance, and type 2 diabetes mellitus but also plays an important role in inflammation and allergic diseases. Dysregulation of pulmonary adipocytokine/insulin signaling caused by early-onset obesity has been proven to induce asthma-like disease in mice [44]. The leptin/osteopontin axis promotes Th2 inflammation and Th17 responses in AR through the NF-кB, MAPK, JNK pathway and β3 integrin [45, 46]. Signaling through the B-cell receptor (BCR) is crucial for antigen recognition and subsequent biological effects, including B-cell activation, proliferation, and differentiation, which ensure host defense [47]. One study demonstrated that the BCR signaling pathway was significantly enriched among differentially expressed vesicle miRNAs in AR patient nasal mucus, consistent with the findings in the current study and further elucidating the importance of the BCR signaling pathway in AR development [48]. Chemokines are small molecule-scale cytokines that recruit leukocyte subsets under steady-state and pathological conditions; signaling pathways are activated by their binding to receptors on the cell surface and are involved in chronic inflammatory and autoimmune diseases. Multiple studies have shown that knockdown of the chemokine receptor CCR3 reduces eosinophilic inflammation and the Th2 immune response in AR [49,50,51]. In summary, our findings are in accordance with all of the above studies.

Comparisons with other studies and contribution of the current work to existing knowledge

To the best of our knowledge, previous exploration of AR based on GSE75011 and GSE46171 has mainly targeted key genes differentially expressed between control and AR samples [18, 19, 52, 53]. In the current study, the biological significance of lipid metabolism in AR was first systematically explored at the genetic level through these same datasets. Moreover, immune infiltration analysis was conducted to investigate the immune cell targets associated with SMARCD3 to elucidate the underlying role of immune-related treatment targeting the SMARCD3 gene in exploration of AR development.

Study strengths and limitations

Three key LMRGs with high diagnostic value for AR were identified and analyzed for the first time based on bioinformatics analysis of AR-related expression datasets. However, the limitations of this study cannot be ignored. First, small sample sizes and small datasets of AR may have introduced bias. Second, the mechanisms of these genes in AR development have not been clearly elucidated. Further research is needed for the possibility of clinical use in the future.


In summary, this is the first bioinformatics analysis of LMRGs in AR, and three key genes (LPCAT1, SGPP1 and SMARCD3) with high diagnostic value for AR were identified, where the expression of two of these genes were confirmed by clinical validation and are considered potential treatment targets. Not only the predictive capability of these key LMRGs were confirmed more excellent than that of other clinical features both in the GSE75011 and GSE46171, but also the nomogram by converting the expression of the key LMRGs into a score was constructed to take all of which into consideration for clinical utilize. Considering the samples regarding the AR datasets used in the study were all from peripheral blood, we make the case that the nomogram could be considered as an effectively clinical diagnostic and practical device with the help of the peripheral blood tests and evaluation of the key LMRGs expressions. Besides, correlation of SMARCD3 expression and immune cell infiltration was helpful to reveal future research directions of immune-related treatment targeting the SMARCD3 gene in AR.

Data Availability

GSE75011 and GSE46171 were downloaded from the GEO database (

Change history



allergic rhinitis


lipid metabolism-related genes


Gene Expression Omnibus


weighted gene co-expression network analysis


AR-related genes


differentially expressed genes


protein-protein interaction


edge percolated component


receiver operating characteristic


real-time fluorescence quantitative polymerase chain reaction


intervertebral disc degeneration


prostaglandin D2


platelet-activating factor


peroxisome proliferator-activated receptor gamma


eicosapentaenoic acid


15-hydroxyeicosapentaenoic acid


Kyoto Encyclopedia of Genes and Genomes


topological overlap matrix


gene significance


module membership


Gene Ontology


Disease Ontology


Search Tool for the Retrieval of Interacting Genes


area under the curve


single-set gene set enrichment analysis


gene set enrichment analysis


normalized enrichment score




glyceraldehyde-3-phosphate dehydrogenase


sterol regulatory element-binding protein


biological process


cellular component


molecular function


regulatory T cell


T follicular helper cell


central memory T cell


B-cell receptor


  1. Brozek JL, Bousquet J, Baena-Cagnani CE, Bonini S, Canonica GW, Casale TB, et al. Allergic Rhinitis and its impact on Asthma (ARIA) guidelines: 2010 revision. J Allergy Clin Immunol. 2010;126(3):466–76.

    Article  PubMed  Google Scholar 

  2. Wheatley LM, Togias A. Clinical practice. Allergic rhinitis. N Engl J Med. 2015;372(5):456–63.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Seidman MD, Gurgel RK, Lin SY, Schwartz SR, Baroody FM, Bonner JR, et al. Clinical practice guideline: allergic rhinitis. Otolaryngol Head Neck Surg. 2015;152(1 Suppl):1–43.

    Google Scholar 

  4. Olzmann JA, Carvalho P. Dynamics and functions of lipid droplets. Nat Rev Mol Cell Biol. 2019;20(3):137–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Li W, Ding Z, Zhang H, Shi Q, Wang D, Zhang S, et al. The roles of blood lipid-metabolism genes in Immune Infiltration could promote the development of IDD. Front Cell Dev Biol. 2022;10:844395.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Merino Salvador M, Gomez de Cedron M, Moreno Rubio J, Falagan Martinez S, Sanchez Martinez R, Casado E, et al. Lipid metabolism and lung cancer. Crit Rev Oncol Hematol. 2017;112:31–40.

    Article  PubMed  Google Scholar 

  7. Lee K, Lee SH, Kim TH. The Biology of Prostaglandins and Their Role as a Target for Allergic Airway Disease Therapy. Int J Mol Sci. 2020;21(5).

  8. Yokomizo T, Nakamura M, Shimizu T. Leukotriene receptors as potential therapeutic targets. J Clin Invest. 2018;128(7):2691–701.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Tomazic PV, Birner-Gruenberger R, Leitner A, Darnhofer B, Spoerk S, Lang-Loidolt D. Apolipoproteins have a potential role in nasal mucus of allergic rhinitis patients: a proteomic study. Laryngoscope. 2015;125(3):E91–6.

    Article  CAS  PubMed  Google Scholar 

  10. Schauberger E, Peinhaupt M, Cazares T, Lindsley AW. Lipid mediators of allergic disease: pathways, treatments, and emerging therapeutic targets. Curr Allergy Asthma Rep. 2016;16(7):48.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Stark JM, Coquet JM, Tibbitt CA. The role of PPAR-gamma in allergic disease. Curr Allergy Asthma Rep. 2021;21(11):45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Sawane K, Nagatake T, Hosomi K, Hirata SI, Adachi J, Abe Y et al. Dietary Omega-3 Fatty Acid Dampens Allergic Rhinitis via Eosinophilic Production of the Anti-Allergic Lipid Mediator 15-Hydroxyeicosapentaenoic Acid in Mice. Nutrients. 2019;11(12).

  13. Tao Z, Shi A, Li R, Wang Y, Wang X, Zhao J. Microarray bioinformatics in cancer- a review. J BUON. 2017;22(4):838–43.

    PubMed  Google Scholar 

  14. Wang T, Wu J, Wu Y, Cheng Y, Deng Y, Liao J, et al. A novel microRNA-based signature predicts prognosis among nasopharyngeal cancer patients. Exp Biol Med (Maywood). 2021;246(1):72–83.

    Article  CAS  PubMed  Google Scholar 

  15. Xiao Y, Xiao T, Ou W, Wu Z, Wu J, Tang J, et al. LncRNA SNHG16 as a potential biomarker and therapeutic target in human cancers. Biomark Res. 2020;8:41.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Ntontsi P, Photiades A, Zervas E, Xanthou G, Samitas K. Genetics and Epigenetics in Asthma. Int J Mol Sci. 2021;22(5).

  17. Li Y, Chen X, Zhang R, Chen M, Shen J, Wu J et al. Correlation analysis of lipid metabolism genes with the immune microenvironment in gastric cancer and the construction of a novel gene signature. Clin Transl Oncol. 2022.

  18. Ma Y, Liu C, Xi G, Guan Y, Tang Y, Zhang J, et al. Bioinformatic analysis and Cellular assays identify substance P influencing Th17/Treg differentiation via the MyD88 pathway as a potential contributor to the progression of Asthma and allergic Rhinitis. Dis Markers. 2022;2022:3843954.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Chang WC, Xu Y. [Integrated bioinformatics analysis of key genes in allergic rhinitis]. Zhonghua Er Bi Yan Hou Tou Jing Wai Ke Za Zhi. 2020;55(5):458–64.

    CAS  PubMed  Google Scholar 

  20. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Nicastro KR, Pearson GA, Ramos X, Pearson V, McQuaid CD, Zardi GI. Transcriptome wide analyses reveal intraspecific diversity in thermal stress responses of a dominant habitat-forming species. Sci Rep. 2023;13(1):5645.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Xi Y, Zhang T, Sun W, Liang R, Ganesh S, Chen H. GOLM1 and FAM49B: Potential Biomarkers in HNSCC Based on Bioinformatics and Immunohistochemical Analysis. Int J Mol Sci. 2022;23(23).

  23. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innov (Camb). 2021;2(3):100141.

    CAS  Google Scholar 

  25. Bello SM, Shimoyama M, Mitraka E, Laulederkind SJF, Smith CL, Eppig JT et al. Disease Ontology: improving and unifying disease annotations across species. Dis Model Mech. 2018;11(3).

  26. Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, et al. pROC: an open-source package for R and S + to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Hauben M, Hung E, Hsieh WY. An exploratory factor analysis of the spontaneous reporting of severe cutaneous adverse reactions. Ther Adv Drug Saf. 2017;8(1):4–16.

    Article  PubMed  Google Scholar 

  28. Ledley FD, Grenett HE, McGinnis-Shelnutt M, Woo SL. Retroviral-mediated gene transfer of human phenylalanine hydroxylase into NIH 3T3 and hepatoma cells. Proc Natl Acad Sci U S A. 1986;83(2):409–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  30. Harayama T, Eto M, Shindou H, Kita Y, Otsubo E, Hishikawa D, et al. Lysophospholipid acyltransferases mediate phosphatidylcholine diversification to achieve the physical properties required in vivo. Cell Metab. 2014;20(2):295–305.

    Article  CAS  PubMed  Google Scholar 

  31. Lebok P, von Hassel A, Meiners J, Hube-Magg C, Simon R, Hoflmayer D, et al. Up-regulation of lysophosphatidylcholine acyltransferase 1 (LPCAT1) is linked to poor prognosis in breast cancer. Aging. 2019;11(18):7796–804.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Mansilla F, da Costa KA, Wang S, Kruhoffer M, Lewin TM, Orntoft TF, et al. Lysophosphatidylcholine acyltransferase 1 (LPCAT1) overexpression in human colorectal cancer. J Mol Med (Berl). 2009;87(1):85–97.

    Article  CAS  PubMed  Google Scholar 

  33. Li L, Wang X, Ding Y, Hui N, Su B, Yang M. LPCAT1 acts as an independent prognostic biomarker correlated with immune infiltration in hepatocellular carcinoma. Eur J Med Res. 2022;27(1):216.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Cheng S, Chen H, Wang A, Xie M, Xie J, Osanai K, et al. Lentiviral vector-mediated delivery of lysophosphatidylcholine acyltransferase 1 attenuates airway inflammation in ovalbumin-induced allergic asthmatic mice. Asian Pac J Allergy Immunol. 2015;33(4):320–9.

    CAS  PubMed  Google Scholar 

  35. Le Stunff H, Giussani P, Maceyka M, Lepine S, Milstien S, Spiegel S. Recycling of sphingosine is regulated by the concerted actions of sphingosine-1-phosphate phosphohydrolase 1 and sphingosine kinase 2. J Biol Chem. 2007;282(47):34372–80.

    Article  PubMed  Google Scholar 

  36. Huang X, Taeb S, Jahangiri S, Emmenegger U, Tran E, Bruce J, et al. miRNA-95 mediates radioresistance in tumors by targeting the sphingolipid phosphatase SGPP1. Cancer Res. 2013;73(23):6972–86.

    Article  CAS  PubMed  Google Scholar 

  37. Zhang B, Gao S, Bao Z, Pan C, Tian Q, Tang Q. MicroRNA-656-3p inhibits colorectal cancer cell migration, invasion, and chemo-resistance by targeting sphingosine-1-phosphate phosphatase 1. Bioengineered. 2022;13(2):3810–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. He A, Kong SW, Ma Q, Pu WT. Co-occupancy by multiple cardiac transcription factors identifies transcriptional enhancers active in heart. Proc Natl Acad Sci U S A. 2011;108(14):5632–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Crotty ST, Follicular Helper. Cell Biology: a decade of Discovery and Diseases. Immunity. 2019;50(5):1132–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Kamekura R, Shigehara K, Miyajima S, Jitsukawa S, Kawata K, Yamashita K, et al. Alteration of circulating type 2 follicular helper T cells and regulatory B cells underlies the comorbid association of allergic rhinitis with bronchial asthma. Clin Immunol. 2015;158(2):204–11.

    Article  CAS  PubMed  Google Scholar 

  41. Hoffmann K, Dixon DN, Greene WK, Ford J, Taplin R, Kees UR. A microarray model system identifies potential new target genes of the proto-oncogene HOX11. Genes Chromosomes Cancer. 2004;41(4):309–20.

    Article  CAS  PubMed  Google Scholar 

  42. Galic S, Oakhill JS, Steinberg GR. Adipose tissue as an endocrine organ. Mol Cell Endocrinol. 2010;316(2):129–39.

    Article  CAS  PubMed  Google Scholar 

  43. Achari AE, Jain SK. Adiponectin, a Therapeutic Target for Obesity, Diabetes, and Endothelial Dysfunction. Int J Mol Sci. 2017;18(6).

  44. Dinger K, Kasper P, Hucklenbruch-Rother E, Vohlen C, Jobst E, Janoschek R, et al. Early-onset obesity dysregulates pulmonary adipocytokine/insulin signaling and induces asthma-like disease in mice. Sci Rep. 2016;6:24168.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Liu W, Zeng Q, Zhou L, Li Y, Chen Y, Luo R. Leptin/osteopontin axis contributes to enhanced T helper 17 type responses in allergic rhinitis. Pediatr Allergy Immunol. 2018;29(6):622–9.

    Article  PubMed  Google Scholar 

  46. Zeng Q, Luo X, Han M, Liu W, Li H. Leptin/Osteopontin Axis regulated type 2T helper cell response in allergic rhinitis with obesity. EBioMedicine. 2018;32:43–9.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Tanaka S, Baba Y. B cell receptor signaling. Adv Exp Med Biol. 2020;1254:23–36.

    Article  CAS  PubMed  Google Scholar 

  48. Wu G, Yang G, Zhang R, Xu G, Zhang L, Wen W, et al. Altered microRNA expression profiles of Extracellular vesicles in nasal mucus from patients with allergic Rhinitis. Allergy Asthma Immunol Res. 2015;7(5):449–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Fortin M, Wagner JG, Brault J, Harkema JR, Renzi PM, Paquet L, et al. Spatial and temporal expression of CCR3 and the common beta chain of the IL-3, IL-5 and GM-CSF receptor in the nasal epithelium and lymphoid tissues in a rat model of allergic rhinitis. Cytokine. 2010;52(3):194–202.

    Article  CAS  PubMed  Google Scholar 

  50. Dai M, Zhu X, Yu J, Yuan J, Zhu Y, Bao Y, et al. CCR3 gene knockout in bone marrow cells ameliorates combined allergic rhinitis and asthma syndrome (CARAS) by reducing airway inflammatory cell infiltration and Th2 cytokines expression in mice model. Int Immunopharmacol. 2022;104:108509.

    Article  CAS  PubMed  Google Scholar 

  51. Yuan J, Liu Y, Yu J, Dai M, Zhu Y, Bao Y, et al. Gene knockdown of CCR3 reduces eosinophilic inflammation and the Th2 immune response by inhibiting the PI3K/AKT pathway in allergic rhinitis mice. Sci Rep. 2022;12(1):5411.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Yan Z, Liu L, Jiao L, Wen X, Liu J, Wang N. Bioinformatics Analysis and Identification of underlying biomarkers potentially linking allergic Rhinitis and Asthma. Med Sci Monit. 2020;26:e924934.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Hao Y, Wang B, Zhao J, Wang P, Zhao Y, Wang X, et al. Identification of gene biomarkers with expression profiles in patients with allergic rhinitis. Allergy Asthma Clin Immunol. 2022;18(1):20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We thank the GEO database for providing high-quality data, as well as authors for making their valuable research results publicly available.


Not applicable.

Author information

Authors and Affiliations



QT and JW designed and conducted the study. YZ and TW collected and analyzed the data. QT and YZ collected the blood samples and carried out the RT-qPCR tests. QT and YD drafted the manuscript. HL and JW critically reviewed this manuscript. All authors contributed to and approved the final manuscript.

Corresponding authors

Correspondence to Huanhai Liu or Jian Wu.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Ethics approval and consent to participate

Informed consent was obtained from all participants, and all experimental protocols were approved by the Medical Ethics Committee of Shanghai Changzheng Hospital.

Consent for publication

Not applicable.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The original online version of this article was revised: the authors noticed that the equal contributor symbol to authors “Qilei Tao” and “Yajing Zhu” were missing.

Electronic supplementary material

Below is the link to the electronic supplementary material.


Supplementary Materials 1: Table S1. Primer sequences for real-time fluorescence quantitative polymerase chain reaction (RT-qPCR)

Supplementary Materials 2: Table S2. List of 73 lipid metabolism-related differentially expressed genes (LMR DEGs)


Supplementary Materials 3: Table S3. Results of the 256 Gene Ontology (GO) enriched terms using gene set enrichment analysis (GSEA)


Supplementary Materials 4: Table S4. Results of the 33 Kyoto Encyclopedia of Genes and Genomes (KEGG) enriched terms using GSEA

Supplementary Materials 5: iThenticate repor

Supplementary Materials 6: SNAS Editing Certificate

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tao, Q., Zhu, Y., Wang, T. et al. Identification and analysis of lipid metabolism-related genes in allergic rhinitis. Lipids Health Dis 22, 105 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Allergic rhinitis
  • Lipid metabolism
  • Biomarkers
  • Diagnostic model
  • Immune infiltration