Skip to main content
Download PDF

Decreased TLR7 expression was associated with airway eosinophilic inflammation and lung function in asthma: evidence from machine learning approaches and experimental validation

European Journal of Medical Research volume 29, Article number: 116 (2024) Cite this article

Abstract

Background

Asthma is a global public health concern. The underlying pathogenetic mechanisms of asthma were poorly understood. This study aims to explore potential biomarkers associated with asthma and analyze the pathological role of immune cell infiltration in the disease.

Methods

The gene expression profiles of induced sputum were obtained from Gene Expression Omnibus datasets (GSE76262 and GSE137268) and were combined for analysis. Toll-like receptor 7 (TLR7) was identified as the core gene by the intersection of two different machine learning algorithms, namely, least absolute shrinkage and selector operation (LASSO) regression and support vector machine-recursive feature elimination (SVM-RFE), and the top 10 core networks based on Cytohubba. CIBERSORT algorithm was used to analyze the difference of immune cell infiltration between asthma and healthy control groups. Finally, the expression level of TLR7 was validated in induced sputum samples of patients with asthma.

Results

A total of 320 differential expression genes between the asthma and healthy control groups were screened, including 184 upregulated genes and 136 downregulated genes. TLR7 was identified as the core gene after combining the results of LASSO regression, SVM-RFE algorithm, and top 10 hub genes. Significant differences were observed in the distribution of 13 out of 22 infiltrating immune cells in asthma. TLR7 was found to be closely related to the level of several infiltrating immune cells. TLR7 mRNA levels were downregulated in asthmatic patients compared with healthy controls (p = 0.0049). The area under the curve of TLR7 for the diagnosis of asthma was 0.7674 (95% CI 0.631–0.904, p = 0.006). Moreover, TLR7 mRNA levels were negatively correlated with exhaled nitric oxide fraction (r = − 0.3268, p = 0.0347) and the percentage of peripheral blood eosinophils (%) (r = − 0.3472, p = 0.041), and positively correlated with forced expiratory volume in the first second (FEV1) (% predicted) (r = 0.3960, p = 0.0071) and FEV1/forced vital capacity (r = 0.3213, p = 0.0314) in asthmatic patients.

Conclusions

Decreased TLR7 in the induced sputum of eosinophilic asthmatic patients was involved in immune cell infiltration and airway inflammation, which may serve as a new biomarker for the diagnosis of eosinophilic asthma.

Background

Asthma is a common chronic disease in which airways become inflamed and narrow, causing airflow obstruction [1,2,3]. Asthma is a heterogeneous clinical syndrome that affects more than 300 million people worldwide [4]. The common symptoms of asthma in the acute phase include wheezing, coughing, chest tightness, and shortness of breath [1, 2]. Asthma is a complex and heterogenous respiratory diseases. The underlying pathogenetic mechanisms of asthma were poorly understood [5].

Induced sputum has several desirable characteristics as a noninvasive marker of airway inflammation [6]. In patients with asthma, sputum induction is generally a well-tolerated and safe method, and sputum can be used to measure various soluble mediators, including eosinophilic-derived proteins, cytokines, and remodeling-related proteins [6,7,8,9,10,11]. Induced sputum may be used to discover inflammatory cell profiles in patients with asthma and other airway diseases, and these profiles may be related to the patient’s response to treatment [12]. The gene expression profile of induced sputum cells is altered in patients with asthma [13].

Microarray technology and integrated bioinformatics analysis have been used in recent years to identify novel genes associated with various diseases that may serve as biomarkers for diagnosis and prognosis [14, 15]. Bioinformatics analysis has also been performed to identify the underlying mechanisms and hub genes of asthma [16, 17]. Studies have also shown that immune cell infiltration plays an increasingly important role in the occurrence and development of various diseases [18,19,20,21]. Previous studies demonstrated that the Th1/Th2-mediated immune imbalance is the main mechanism of asthmatic airway inflammatory response, and various immune cells are involved in the pathogenesis of asthma [22].

CIBERSORT, a method for characterizing cell composition of complex tissues from their gene expression profiles, has been widely used to evaluate the relative content of 22 kinds of immune cells [23]. CIBERSORT method has also been applied to study the immune cell infiltration and candidate diagnostic markers in asthma. It has been reported by Yang et al. that autophagy-related genes are involved in the progression and prognosis of asthma and regulate the immune microenvironment [24]. Least absolute shrinkage and selector operation (LASSO) regression and support vector machine-recursive feature elimination (SVM-RFE) are two machine learning algorithms. LASSO is a dimension-reduction algorithm that can analyze high-dimensional data compared with regression analysis [25]. SVM-RFE is a machine learning algorithm used to identify the best variables through classification method [26]. The combination of LASSO and SVM-RFE has been applied in previous research to identify diagnostic markers [20, 27, 28].

In the present study, bioinformatics analysis and experimental validation were performed to investigate the change of immune cell infiltration in asthma, and screen the biomarker for the diagnosis and treatment of asthma. Two datasets from Gene Expression Omnibus (GEO) database were combined, and differential expression gene (DEG) analysis, machine learning algorithms and CIBERSORT were performed. Toll-like receptor 7 (TLR7), a candidate gene that was found to be closely associated with immune infiltration in asthma, was also validated in another GEO dataset and induced sputum samples of asthmatic patients.

Material and methods

Subjects

We recruited 12 healthy controls and 36 newly diagnosed asthma patients with untreated asthma. The asthmatic patients included in this study and the control group were non-smokers, and the asthmatic patients were newly diagnosed and untreated. The asthmatic patients were from outpatients and were diagnosed with asthma by specialists. The characteristics of the subjects are summarized in Table 1. No significant differences were observed in terms of age, sex, and body mass index between the two groups. All subjects provided written informed consent. The study was approved by the Ethics Committee of the First Affiliated Hospital of Sun Yat-sen University (2021071).

Table 1 Characteristics of subjects
Full size table

Dataset acquisition and processing

The study design is shown in Additional file 1: Fig S1. The datasets GSE76262 and GSE137268 were downloaded from the GEO database (http://www.ncbi.nlm.nih.gov/geo). GSE76262 dataset, which is based on GPL13158 platform, included induced sputum samples from 118 asthmatic patients and 21 healthy controls. GSE137268 dataset, which is based on GPL6104 platform, included induced sputum samples from 54 asthmatic patients and 15 healthy controls. The series matrix files were annotated to the official gene symbols, and the two gene expression files were merged. The batch normalization was then conducted using combat method in “sva” R package. Finally, a merged file with 15,043 genes was prepared for the subsequent analysis.

Identification of DEGs and enrichment analysis

The “limma” R package was used to identify DEGs, and the |log2FC|> 0.5 and adjusted p value < 0.05 were filtered as statistically significant. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were then performed using “clusterProfiler” R package. Gene Set Enrichment Analysis (GSEA) was conducted to analyze the associated biological functions and pathways in asthma. Disease Ontology (DO) was also conducted using “DOSE” R package.

Identification of the core gene

First, two distinct machine learning algorithms, namely, least absolute shrinkage and selector operation (LASSO) regression and support vector machine-recursive feature elimination (SVM-RFE), were utilized in DEGs to screen the gene signatures. The LASSO is a regression analysis algorithm that uses regularization to improve the prediction accuracy. The LASSO analysis was undertaken using “glmnet” R package; the response type was set as binomial, and the alpha was set as 1. SVM is a supervised machine-learning technique widely utilized for both classification and regression. To avoid overfitting, an RFE algorithm was employed to select the optimal genes from the meta-data cohort. Therefore, to identify the set of genes with the highest discriminative power, SVM-RFE was applied to select the appropriate features. The SVM-RFE was performed using “e1071” and “caret” R package. Second, STRING database was used to construct the protein–protein interaction (PPI) network, and a core network was obtained through Cytoscape software and CytoHubba plugin. The top 10 hub genes were screened according to Degree algorithm. Finally, the results of LASSO regression, SVM-RFE algorithm, and hub genes were incorporated, and the overlapping gene (TLR7) was identified as the core gene.

Analysis of immune cell infiltration

The CIBERSORT algorithm was used to evaluate the percentage of 22 immune cell types in each sample. The fraction of 22 immune cells was compared between the asthma and healthy control groups, and the violin plot was drawn by “vioplot” R package. The correlation coefficient between immune cells was calculated using “corrplot” R package. Spearman correlation analysis was also performed to investigate the correlation of TLR7 and infiltrating immune cells.

Validation of TLR7 in a GEO dataset

The expression level of TLR7 in the merge dataset was visualized, and receiver operating characteristic (ROC) curve was applied to evaluate the diagnostic value of TLR7 for asthma. Furthermore, the GEO dataset GSE147878 with endobronchial biopsy samples from 60 asthmatic patients and 13 healthy controls was used to validate the expression level and diagnostic effectiveness of TLR7 in asthma.

Collection of induced sputum from subjects

A total of 48 subjects from First Affiliated Hospital of Sun Yat-sen University (Guangzhou, Guangdong, China) were enrolled in this study, including 12 healthy controls and 36 asthmatic patients. Patients with asthma met the diagnostic criteria for Global Asthma Initiative (GINA) guidelines [29] and were free of other respiratory diseases. People with normal lung function test results and no history of pulmonary disease, allergic disease, and autoimmune disease were included in the healthy control group. Sputum samples were collected from the participants. Participants were induced to cough by hypertonic saline. The above steps are completed by ultrasonic atomizer (Yuyue, Jiangsu, China). Sputum cell pellet was selected, weighed, and dissolved by adding 0.1% dithiothreitol (DTT) that is 4 times the weight. The pellet was then filtered through cell sieving [30,31,32]. After centrifugation, sputum cells were added with 1 ml TRIzol for subsequent RNA extraction. Additional clinical information was collected for each subject, including lung function, exhaled nitric oxide fraction (FeNO), and peripheral blood eosinophil percentage.

Quantitative real‐time polymerase chain reaction (qRT-PCR)

Total RNA was extracted from induced sputum cells using TRIzol reagent following the manufacturer’s instructions. Evo M-MLV RT Premix kit (AG, Hunan, China) was used for reverse transcription. The reaction conditions were 37 ℃ for 15 min and 85 ℃ for 5 s. Candidate gene expression was quantified using Biosystems Light Cycler 480 (Applied Biosystems, Massachusetts, USA) as standard procedure. The primers used were TLR7: forward, 5′- TCCTTGGGGCTAGATGGTTTC-3′, reverse, 5′- TCCACGATCACATGGTTCTTTG-3′ and GAPDH: forward, 5ʹ-ACCCAGAAGACTGTGGATGG-3ʹ, reverse, 5ʹ-TTCTAGACGGCAGGTCAGGT-3ʹ.

Statistical analysis

All data in this study were analyzed through GraphPad Prism 8. 0 (GraphPad, San Diego, California, USA). Normally distributed data were obtained through unpaired t-test and expressed as mean ± standard deviation. For non-normally distributed data, the results were obtained via a nonparametric test (i.e., Kruskal–Wallis test) and expressed as median (interquartile spacing). Fisher’s exact test was used to analyze classified data, and Spearman rank correlation was used for correlation analysis. ROC was generated to determine the diagnostic value of TLR7. P < 0.05 was considered statistically significant.

Results

Identification of DEGs and enrichment analysis

The inclusion criteria of the DEGs were |log2FC|> 0.5 and adjusted p value < 0.05. A total of 320 DEGs between the asthma and healthy control groups were screened, including 184 upregulated genes and 136 downregulated genes. The expressions of the DEGs in each sample are shown in the heatmap (Additional file 1: Fig S2A), and the distribution of the DEGs is illustrated through a volcano plot (Additional file 1: Fig S2B).

GO, KEGG, GSEA, and DO analyses were performed to further investigate the DEGs’ functions. GO enrichment analysis was conducted to analyze the gene function in terms of biological processes (BP), cellular component (CC), and molecular function (MF). The GO analysis results showed that BP is mainly enriched in the regulation of immune effector process, CC is mainly enriched in tertiary granule, and MF is mainly enriched in immune receptor activity (Additional file 1: Fig S3A). The KEGG pathway enrichment analysis demonstrated that the DEGs were mainly involved in cytokine–cytokine receptor interaction, tumor necrosis factor (TNF) signaling pathway, and nuclear factor (NF)-kappa B signaling pathway (Additional file 1: Fig S3B). As shown in Fig. 1A, GSEA also showed that the significantly enriched hallmark terms associated with asthma included chemokine signaling pathway, cytokine–cytokine receptor interaction, Janus kinase/signal transducer and activator of transcription (JAK/STAT) signaling pathway, mitogen-activated protein kinase (MAPK) signaling pathway, and neuroactive ligand receptor interaction. Furthermore, DO analysis revealed that the DEGs were mainly related to lung disease and obstructive lung disease (Fig. 1B).

Fig. 1
figure 1

GSEA and DO enrichment analyses. A Illustration of several important enrichment hallmark terms in asthma obtained through GSEA. B Top 20 terms in the DO enrichment analysis

Full size image

Identification of TLR7 as the core gene

To explore the biomarkers of asthma, two distinct machine learning algorithms, namely, the LASSO regression and SVM-RFE, were performed. The LASSO regression analysis identified 46 DEGs as signature genes in asthma (Fig. 2A). The SVM-RFE algorithm screened 28 DEGs as characteristic genes in asthma (Fig. 2B). In addition, a PPI network of DEGs was constructed using the STRING database. A core network was then obtained through Degree algorithm in the Cytohubba plugin, and 10 hub genes were identified (Fig. 2C). After combining the results of LASSO regression, SVM-RFE algorithm, and hub genes by Venn diagram, only an intersection gene was identified, i.e., TLR7 (Fig. 2D). TLR7 was thus identified as the core gene for the subsequent research.

Fig. 2
figure 2

Two machine learning algorithms and PPI network were performed for core gene selection. A LASSO regression analysis. B SVM-RFE algorithm. C Hub genes based on the Degree algorithm in the Cytohubba plugin. D Venn diagram showing the overlapping gene of LASSO regression, SVM-RFE algorithm, and hub genes

Full size image

Immune infiltration analyses

CIBERSORT algorithm was used to analyze the difference of immune cell infiltration between the asthma and healthy control groups in 22 subpopulations of immune cells. The total value of all immune cells in each sample was set at 100%, and the proportion of each immune cell in these samples is presented in Fig. 3A. The interaction between the immune cells was also analyzed. Average linkage clustering revealed that M1 macrophages and activated memory CD4 T cells have a significant positive correlation, whereas neutrophils and M0 macrophages are significantly negatively correlated (Fig. 3B). The violin plot showed marked differences in the distribution of 13 out of 22 immune cells (Fig. 3C). The fractions of naive CD4 T cells (p = 0.027), resting dendritic cells (p = 0.018), activated dendritic cells (p < 0.001), and eosinophils (p = 0.006) in the asthma group were remarkably higher compared with those of the healthy controls, while the fractions of memory B cells (p = 0.018), CD8 T cells (p = 0.015), activated memory CD4 T cells (p = 0.043), follicular helper T cells (p < 0.001), gamma delta T cells (p = 0.018), monocytes (p = 0.021), M0 macrophages (p = 0.006), M1 macrophages (p = 0.029), and M2 macrophages (p < 0.001) were lower in asthma. Taken together, these results suggest that the heterogeneity of infiltrating immune cells in asthma is evident and may play a role in the pathogenesis of asthma.

Fig. 3
figure 3

Landscape of immune infiltration between the asthma and healthy control groups. A The box plot diagram indicates the relative percentage of different types of immune cells in each sample. B The heatmap shows the correlation of infiltrating immune cells. C The violin plot shows the differences of immune infiltration between the asthma (red) and healthy control (blue) groups. D The lollipop chart presents the correlation of TLR7 and infiltrating immune cells on the basis of Spearman correlation analysis results. (p value < 0.05 indicated statistical significance)

Full size image

To further investigate the correlation of TLR7 and infiltrating immune cells, Spearman correlation was performed (Table 2) and plotted in a lollipop chart (Fig. 3D) and several scatter charts (Additional file 1: Fig S4). The results demonstrated that TLR7 was positively correlated with M2 macrophages (r = 0.59, p < 0.001), follicular helper T cells (r = 0.46, p < 0.001), memory B cells (r = 0.36, p < 0.001), CD8 T cells (r = 0.32, p < 0.001), M0 macrophages (r = 0.32, p < 0.001), M1 macrophages (r = 0.28, p < 0.001), resting memory CD4 T cells (r = 0.24, p < 0.001), and monocytes (r = 0.21, p < 0.01). Meanwhile, TLR7 was negatively correlated with activated dendritic cells (r = − 0.52, p < 0.001), naive CD4 T cells (r = − 0.40, p < 0.001), plasma cells (r = − 0.25, p < 0.001), eosinophils (r = − 0.22, p < 0.01), resting dendritic cells (r = − 0.20, p < 0.01), activated mast cells (r = − 0.19, p < 0.01), activated memory CD4 T cells (r = − 0.18, p < 0.05), and resting NK cells (r = − 0.17, p < 0.05). These results indicate that the core gene TLR7 is closely related to the level of immune cell infiltration and plays a crucial role in the immune microenvironment of asthma.

Table 2 Correlation of TLR7 and infiltrating immune cells
Full size table

Validation of TLR7 in a GEO dataset and the diagnostic value of TLR7 for asthma

In the merged dataset, the expression level of TLR7 in the asthma group significantly decreased compared with that of the healthy control group (p < 0.001, Fig. 4A). ROC curve analysis was conducted to evaluate the sensitivity and specificity of TLR7 for the diagnosis of asthma. As shown in Fig. 4B, the area under curve (AUC) value of TLR7 was 0.799 (95% CI 0.719–0.874). Moreover, the GSE147878 dataset was used to validate the expression and diagnostic effectiveness of TLR7 in asthma. Consistently, the TLR7 expression level in the asthma group of the GSE147878 dataset also significantly decreased (p < 0.01, Fig. 5A), and the AUC value of TLR7 was 0.783 (95% CI 0.645–0.897, Fig. 5B).

Fig. 4
figure 4

TLR7 expression level and its diagnostic value in asthma. A The expression level of TLR7 in the asthma (red) and healthy control (blue) groups in the merged dataset. (p value < 0.05 indicated statistical significance). B ROC curve analysis of TLR7 in the merged dataset

Full size image
Fig. 5
figure 5

Validation of the expression and diagnostic value of TLR7 in the GSE147878 dataset. A The expression level of TLR7 in the asthma (red) and healthy control (blue) groups in the GSE147878 dataset. (p value < 0.05 indicated statistical significance). B The ROC curve analysis of TLR7 in the GSE147878 dataset

Full size image

Validation of TLR7 mRNA expression in induced sputum cells of asthmatic patients

Detection of TLR7 mRNA levels via qRT-PCR showed that TLR7 mRNA levels were significantly downregulated in asthmatic patients compared with those in healthy controls (p = 0.0049, Fig. 6A). The AUC value was 0.7674 (95% CI 0.631–0.904, p = 0.006) (Fig. 6B). Our test results are thus consistent with those of the GEO dataset, and TLR7 has a satisfactory diagnostic ability for asthma.

Fig. 6
figure 6

Validation of the expression and diagnostic value of TLR7 in asthmatic patients. A TLR7 mRNA expression level in induced sputum cells of asthma. B ROC curve of TLR7 in induced sputum cells

Full size image

TLR7 mRNA expression is associated with airway eosinophilic inflammation and lung function

We investigated the correlation between TLR7 mRNA expression and clinical indicators such as FeNO, percentage of peripheral blood eosinophils (%), and lung function. The results showed that TLR7 mRNA expression was significantly negatively correlated with FeNO (r = − 0.3268, p = 0.0347) (Fig. 7A) and percentage of peripheral blood eosinophils (%) (r = − 0.3472, p = 0.041) (Fig. 7B), and positively correlated with forced expiratory volume in the first second (FEV1) (% predicted) (r = 0.3960, p = 0.0071) (Fig. 7C) and FEV1/forced vital capacity (FVC) (r = 0.3213, p = 0.0314) (Fig. 7D). These data suggest that TLR7 is involved in the pathogenesis of eosinophilic inflammation and bronchoconstriction in asthmatic patients.

Fig. 7
figure 7

Relationship between TLR7 mRNA expression level and clinical parameters. Relationship between TLR7 mRNA expression level in induced sputum cells and A FeNO, B percentage of peripheral blood eosinophils (%), C FEV1 (% predicted), and D FEV1/FVC (%)

Full size image

Discussion

Asthma is a common chronic disease [2]. Induced sputum may have some characteristics as a noninvasive marker of airway inflammation [12]. The gene expression profile of induced sputum cells is altered in patients with asthma [13]. In the current study, two datasets (i.e., GSE76262 and GSE137268), including induced sputum samples of 172 asthmatic patients and 36 healthy controls, were combined for analysis. The combat algorithm in “sva” R package was used to eliminate batch effect [33]. TLR7 was identified as the core gene through the intersection of two different machine learning algorithms (i.e., LASSO regression and SVM-RFE) and the top 10 core networks based on Cytohubba. The immune infiltration analysis results showed that TLR7 is closely related to the level of numerous infiltrating immune cells. Finally, the decreased TLR7 expression levels were validated in induced sputum samples of patients with asthma. The diagnostic value of TLR7 for eosinophilic asthma was evaluated, and its correlation with related clinical indicators was also analyzed.

In the present study, a total of 320 DEGs between the asthma and healthy control groups were obtained. GO and KEGG analyses revealed that DEGs between the asthma and healthy controls were primarily enriched in cytokine–cytokine receptor interaction and immune-related functions, such as immune effector process and immune receptor activity. GSEA is a threshold-free method that analyzes all genes on the basis of their differential expression rank, or other score, without prior gene filtering [34]. GSEA results coincided with the GO and KEGG results. Moreover, these DEGs were proven to be related to lung diseases, such as asthma, by DO analysis. Furthermore, two machine learning algorithms, the LASSO regression and SVM-RFE, were performed to identify the biomarkers of asthma. The combination of LASSO and SVM-RFE has been applied in previous research to identify diagnostic markers [20, 27, 28]. The traditional PPI network of DEGs was also constructed to identify hub genes. After combining the results of LASSO, SVM-RFE, and hub genes, decreased TLR7 was finally identified as the core gene of asthma.

Toll-like receptors (TLRs) play crucial roles in the recognition of invading pathogens and the immune system. The role of TLR signatures in asthma has been reported by Wu et al. that TLR2/TLR3/TLR4 pathway, MyD88-dependent/independent TLR pathway, positive regulation of TLR4 pathway and TLR binding signatures were correlated with asthma [35]. TLR7 is an endosomal receptor that recognizes microbial or self-antigen-derived single-stranded RNA ligands [36]. Currently, TLR7 has been reported to be involved in the pathogenesis of various immunological diseases [37,38,39,40,41,42,43]. Research reports that TLR7 agonists reduce Th2-mediated airway inflammation, airway hyperreactivity, and chronic airway remodeling in asthma [44]. Jha A and coworkers also achieved similar results [45]. TLR7 agonists can increase the expression of interferon and C-C motif chemokine ligand 13 (CCL13) in nasal mucosa of patients with asthma and allergic rhinitis [46]. Several research findings also revealed that TLR7 regulates RV1b-induced type I and type III interferon signaling pathways in allergic asthma [47]. TLR7 may confer predisposition to asthma and related atopic diseases [48]. A significant correlation was found between TLR7 single nucleotide polymorphism (SNP) and childhood asthma [49]. Furthermore, the expression of TLR7 in the airway of asthmatic mice was significantly decreased, and upregulation of TLR7 was found to inhibit the activation of NF-κB signaling pathway, reduce airway inflammation, inhibit the proliferation of airway smooth muscle cells (ASMCS), and promote apoptosis in asthmatic mice [50]. Recently, TLR7-nanoparticle adjuvants have been reported to improve the immune response to viral antigens [51]. TLR7 plays a key role in the pathogenesis of rosacea by activating the NFκB-mTORC1 axis [52]. Another study also showed that TLR7 expression is decreased in the lungs of patients with severe asthma [53]. The GSE147878 dataset confirmed that the TLR7 expression level in asthma is also significantly reduced and has good diagnostic value. The expression trend of our test result was consistent the GEO datasets, that is, TLR mRNA expression is significantly decreased in the induced sputum of asthmatic patients and has satisfactory diagnostic ability. TLR7 mRNA expression was significantly negatively correlated with FeNO and percentage of peripheral blood eosinophils (%) and positively correlated with FEV1 (% predicted) and FEV1/FVC. We thus inferred that TLR7 is involved in the pathogenesis of eosinophilic inflammation and bronchoconstriction in asthmatic patients.

In addition, immune infiltration analysis in this study demonstrated that the changes of infiltrating immune cells in asthma are evident. Significant differences were observed in the distribution of 13 out of 22 immune cells in asthma. The fractions of dendritic cells and eosinophils in the asthma group were remarkably higher, whereas the fractions of memory B cells, T cells, monocytes, and macrophages were lower compared with those of the healthy controls. Interestingly, TLR7 was also found to be closely related to the level of immune cell infiltration in the current study. Therefore, it could be concluded that TLR7 may play a critical role in asthma by regulating immune cells.

There are also inherent limitations in this study. First, the size of induced sputum samples was not sufficiently large. Further study should include more samples. Second, our sample size was small and we did not compare TLR7 protein levels across different asthma subtypes. Finally, the mechanism by which TLR7 affects eosinophilic asthma was not thoroughly studied. Therefore, further studies are warranted to confirm this mechanism as potential new therapeutic targets of eosinophilic asthma.

Conclusions

In conclusion, this study proved that decreased TLR7 in the induced sputum of eosinophilic asthmatic patients was involved in immune cell infiltration and airway inflammation, which may serve as a new biomarker for the diagnosis of eosinophilic asthma.

Availability of data and materials

The data that support the findings of this study are available in GEO database (http://www.ncbi.nlm.nih.gov/geo), reference number [GSE76262, GSE137268 and GSE147878].

Abbreviations

GEO:

Gene expression omnibus

DEG:

Differential expression gene

TLR7:

Toll-like receptor 7

GO:

Gene ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

GSEA:

Gene set enrichment analysis

DO:

Disease Ontology

LASSO:

Least absolute shrinkage and selector operation

SVM-RFE:

Support vector machine-recursive feature elimination

PPI:

Protein–protein interaction

ROC:

Receiver operating characteristic

FeNO:

Exhaled nitric oxide fraction

qRT-PCR:

Quantitative real‐time polymerase chain reaction

BP:

Biological processes

CC:

Cellular component

MF:

Molecular function

TNF:

Tumor necrosis factor

NF:

Nuclear factor

JAK/STAT:

Janus kinase/signal transducer and activator of transcription

MAPK:

Mitogen-activated protein kinase

AUC:

Area under curve

FEV1:

Forced expiratory volume in the first second

FVC:

Forced vital capacity

CCL13:

C-C motif chemokine ligand 13

SNP:

Single nucleotide polymorphism

References

  1. Beasley R, Braithwaite I, Semprini A, Kearns C, Weatherall M, Pavord ID. Optimal asthma control: time for a new target. Am J Respir Crit Care Med. 2020;201(12):1480–7.

    Article  PubMed  Google Scholar 

  2. Boulet LP, Reddel HK, Bateman E, Pedersen S, FitzGerald JM, O’Byrne PM. The global initiative for asthma (GINA): 25 years later. Eur Respir J. 2019;54(2):1900598.

    Article  PubMed  Google Scholar 

  3. Mo Y, Ye L, Cai H, Zhu G, Wang J, Zhu M, et al. SERPINB10 contributes to asthma by inhibiting the apoptosis of allergenic Th2 cells. Respir Res. 2021;22(1):178.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Stern J, Pier J, Litonjua AA. Asthma epidemiology and risk factors. Semin Immunopathol. 2020;42(1):5–15.

    Article  PubMed  Google Scholar 

  5. Ntontsi P, Photiades A, Zervas E, Xanthou G, Samitas K. Genetics and epigenetics in asthma. Int J Mol Sci. 2021;22(5):2412.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Davies AR, Hancox RJ. Induced sputum in asthma: diagnostic and therapeutic implications. Curr Opin Pulm Med. 2013;19(1):60–5.

    Article  PubMed  Google Scholar 

  7. Lemie`re C. The use of sputum eosinophils in the evaluation of occupational asthma. Curr Opin Allergy Clin Immunol. 2004;4(2):81–5.

    Article  Google Scholar 

  8. Hastie AT, Steele C, Dunaway CW, Moore WC, Rector BM, Ampleford E, et al. Complex association patterns for inflammatory mediators in induced sputum from subjects with asthma. Clin Exp Allergy. 2018;48(7):787–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Kim CK, Hagan JB. Sputum tests in the diagnosis and monitoring of asthma. Ann Allergy Asthma Immunol. 2004;93(2):112–22.

    Article  PubMed  Google Scholar 

  10. Mastalerz L, Celejewska-Wojcik N, Wojcik K, Gielicz A, Cmiel A, Ignacak M, et al. Induced sputum supernatant bioactive lipid mediators can identify subtypes of asthma. Clin Exp Allergy. 2015;45(12):1779–89.

    Article  CAS  PubMed  Google Scholar 

  11. Tanaka A, Sato H, Akimoto K, Matsunaga T, Sagara H. Spontaneous sputum discriminates inflammatory phenotypes in patients with asthma. Ann Allergy Asthma Immunol. 2021;126(1):54–60.

    Article  CAS  PubMed  Google Scholar 

  12. Bakakos P, Schleich F, Alchanatis M, Louis R. Induced sputum in asthma: from bench to bedside. Curr Med Chem. 2011;18(10):1415–22.

    Article  CAS  PubMed  Google Scholar 

  13. Kuo CS, Pavlidis S, Loza M, Baribaud F, Rowe A, Pandis I, et al. T-helper cell type 2 (Th2) and non-Th2 molecular phenotypes of asthma using sputum transcriptomics in U-BIOPRED. Eur Respir J. 2017;49(2):1602135.

    Article  PubMed  Google Scholar 

  14. Cao L, Chen Y, Zhang M, Xu DQ, Liu Y, Liu T, et al. Identification of hub genes and potential molecular mechanisms in gastric cancer by integrated bioinformatics analysis. PeerJ. 2018;6:e5180.

    Article  PubMed  PubMed Central  Google Scholar 

  15. He LL, Xu F, Zhan XQ, Chen ZH, Shen HH. Identification of critical genes associated with the development of asthma by co-expression modules construction. Mol Immunol. 2020;123:18–25.

    Article  CAS  PubMed  Google Scholar 

  16. Li M, Zhu W, Wang C, Zheng Y, Sun S, Fang Y, et al. Weighted gene co-expression network analysis to identify key modules and hub genes associated with paucigranulocytic asthma. BMC Pulm Med. 2021;21(1):343.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Zhang Z, Wang J, Chen O. Identification of biomarkers and pathogenesis in severe asthma by coexpression network analysis. BMC Med Genomics. 2021;14(1):51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Chen H, Pan Y, Jin X, Chen G. An immune cell infiltration-related gene signature predicts prognosis for bladder cancer. Sci Rep. 2021;11(1):16679.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  19. Li T, Wang T, Zhao X. Profiles of immune infiltration in abdominal aortic aneurysm and their associated marker genes: a gene expression-based study. Braz J Med Biol Res. 2021;54(11):e11372.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Liu Z, Li H, Pan S. Discovery and validation of key biomarkers based on immune infiltrates in Alzheimer’s disease. Front Genet. 2021;12:658323.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Xiu MX, Liu YM, Chen GY, Hu C, Kuang BH. Identifying hub genes, key pathways and immune cell infiltration characteristics in pediatric and adult ulcerative colitis by integrated bioinformatic analysis. Dig Dis Sci. 2021;66(9):3002–14.

    Article  CAS  PubMed  Google Scholar 

  22. Menzies-Gow A, Corren J, Bourdin A, Chupp G, Israel E, Wechsler ME, et al. Tezepelumab in adults and adolescents with severe, uncontrolled asthma. N Engl J Med. 2021;384(19):1800–9.

    Article  CAS  PubMed  Google Scholar 

  23. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Yue Xu, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Yang F, Kong J, Zong Y, Li Z, Lyu M, Li W, et al. Autophagy-related genes are involved in the progression and prognosis of asthma and regulate the immune microenvironment. Front Immunol. 2022;13:897835.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Tibshirani R. Regression shrinkage and selection via the lasso. J Roy Stat Soc: Ser B. 1996;58:267–88.

    MathSciNet  Google Scholar 

  26. Sanz H, Valim C, Vegas E, Oller JM, Reverter F. SVM-RFE: selection and visualization of the most relevant features through non-linear kernels. BMC Bioinformatics. 2018;19(1):432.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Zhao E, Xie H, Zhang Y. The combine of LASSO and SVM-RFE has also been applied in previous literatures to identify diagnostic markers. Front Cardiovasc Med. 2020;7:586871.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Zhang J, Rongguo Yu, Guo X, Zou Y, Chen S, Zhou K, et al. Identification of TYR, TYRP1, DCT and LARP7 as related biomarkers and immune infiltration characteristics of vitiligo via comprehensive strategies. Bioengineered. 2021;12(1):2214–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Reddel HK, Bacharier LB, Bateman ED, Brightling CE, Brusselle GG, Buhl R, et al. Global initiative for asthma strategy 2021: executive summary and rationale for key changes. Respirology. 2022;27(1):14–35.

    Article  PubMed  Google Scholar 

  30. Gershman NH, Wong HH, Liu JT, Mahlmeister MJ, Fahy JV. Comparison of two methods of collecting induced sputum in asthmatic subjects. Eur Respir J. 1996;9(12):2448–53.

    Article  CAS  PubMed  Google Scholar 

  31. Spanevello A, Beghé B, Bianchi A, Migliori GB, Ambrosetti M, Neri M, et al. Comparison of two methods of processing induced sputum: selected versus entire sputum. Am J Respir Crit Care Med. 1998;157(2):665–8.

    Article  CAS  PubMed  Google Scholar 

  32. Niessen NM, Gibson PG, Baines KJ, Barker D, Yang IA, Upham JW, et al. Sputum TNF markers are increased in neutrophilic and severe asthma and are reduced by azithromycin treatment. Allergy. 2021;76(7):2090–101.

    Article  CAS  PubMed  Google Scholar 

  33. Leek JT, Evan Johnson W, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28(6):882–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Reimand J, Isserlin R, Voisin V, Kucera M, Tannus-Lopes C, Rostamianfar A, et al. Pathway enrichment analysis and visualization of omics data using g:profiler, GSEA, cytoscape and enrichmentmap. Nat Protoc. 2019;14(2):482–517.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Xiaoyu Wu, Wang P, Zhang Y, Gao L, Zheng B, Youwen Xu, et al. Toll-like receptor characterization correlates with asthma and is predictive of diagnosis. DNA Cell Biol. 2020;39(7):1313–21.

    Article  Google Scholar 

  36. Tojo S, Zhang Z, Matsui H, Tahara M, Ikeguchi M, Kochi M, et al. Structural analysis reveals TLR7 dynamics underlying antagonism. Nat Commun. 2020;11(1):5204.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  37. Azar P, Mejia JE, Cenac C, Shaiykova A, Youness A, Laffont S, et al. TLR7 dosage polymorphism shapes interferogenesis and HIV-1 acute viremia in women. JCI Insight. 2020;5(12):e136047.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Chodisetti SB, Fike AJ, Domeier PP, Choi NM, Soni C, Rahman ZSM. TLR7 negatively regulates B10 cells predominantly in an IFNgamma signaling dependent manner. Front Immunol. 2020;11:1632.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Diebold SS, Kaisho T, Hemmi H, Akira S, Sousa CRe. Innate antiviral responses by means of TLR7-mediated recognition of single-stranded RNA. Science. 2012;303(5663):1529–31.

    Article  ADS  Google Scholar 

  40. Kobold S, Wiedemann G, Rothenfußer S, Endres S. Modes of action of TLR7 agonists in cancer therapy. Immunotherapy. 2014;6(10):1085–95.

    Article  CAS  PubMed  Google Scholar 

  41. Van Raemdonck K, Umar S, Palasiewicz K, Romay B, Volkov S, Arami S, et al. TLR7 endogenous ligands remodel glycolytic macrophages and trigger skin-to-joint crosstalk in psoriatic arthritis. Eur J Immunol. 2021;51(3):714–20.

    Article  PubMed  Google Scholar 

  42. Zheng N, Xie K, Ye H, Dong Y, Wang B, Luo N, et al. TLR7 in B cells promotes renal inflammation and Gd-IgA1 synthesis in IgA nephropathy. JCI Insight. 2020;5(14):e136965.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Celhar T, Lu HK, Benso L, Rakhilina L, Lee HY, Tripathi S, et al. TLR7 protein expression in mild and severe lupus-prone models is regulated in a leukocyte, genetic, and IRAK4 dependent manner. Front Immunol. 2019;10:1546.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Drake MG, Kaufman EH, Fryer AD, Jacoby DB. The therapeutic potential of toll-like receptor 7 stimulation in asthma. Inflamm Allergy Drug Targets. 2012;11(6):484–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Leaker BR, Singh D, Lindgren S, Almqvist G, Eriksson L, Young B, et al. Effects of the toll-like receptor 7 (TLR7) agonist, AZD8848, on allergen-induced responses in patients with mild asthma: a double-blind, randomised, parallel-group study. Respir Res. 2019;20(1):288.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Jha A, Thwaites RS, Tunstall T, Kon OM, Shattock RJ, Hansel TT, et al. Increased nasal mucosal interferon and CCL13 response to a TLR7/8 agonist in asthma and allergic rhinitis. J Allergy Clin Immunol. 2021;147(2):694–703.

    Article  CAS  PubMed  Google Scholar 

  47. Krug J, Kiefer A, Koelle J, Vuorinen T, Xepapadaki P, Stanic B, et al. TLR7/8 regulates type I and type III interferon signalling in rhinovirus 1b-induced allergic asthma. Eur Respir J. 2021;57(5):2001562.

    Article  CAS  PubMed  Google Scholar 

  48. Moller-Larsen S, Nyegaard M, Haagerup A, Vestbo J, Kruse TA, Borglum AD. Association analysis identifies TLR7 and TLR8 as novel risk genes in asthma and related disorders. Thorax. 2008;63(12):1064–9.

    Article  CAS  PubMed  Google Scholar 

  49. Shan L, Hou P, Kang X, Shang Y. Effects of single-nucleotide polymorphisms in the TLR7 and TLR9 genes of asthmatic children. Ann Clin Lab Sci. 2018;48(5):601–7.

    CAS  PubMed  Google Scholar 

  50. Song L, Luan B, Xu Q-R, Wang X-F. Effect of TLR7 gene expression mediating NF-κB signaling pathway on the pathogenesis of bronchial asthma in mice and the intervention role of IFN-γ. Eur Rev Med Pharmacol Sci. 2021;25(2):866–79.

    CAS  PubMed  Google Scholar 

  51. Yin Q, Luo W, Mallajosyula V, Bo Y, Guo J, Xie J, et al. A TLR7-nanoparticle adjuvant promotes a broad immune response against heterologous strains of influenza and SARS-CoV-2. Nat Mater. 2023. https://doi.org/10.1038/s41563-022-01464-2.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  52. Huang Y, Liu D, Chen M, Xu S, Peng Q, Zhu Y, et al. TLR7 promotes skin inflammation via activating NFκB-mTORC1 axis in rosacea. PeerJ. 2023. https://doi.org/10.7717/peerj.15976.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Shikhagaie MM, Andersson CK, Mori M, Kortekaas Krohn I, Bergqvist A, Dahl R, et al. Mapping of TLR5 and TLR7 in central and distal human airways and identification of reduced TLR expression in severe asthma. Clin Exp Allergy. 2014;44(2):184–96.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank the participants who volunteered in the study.

Funding

This study is not funded.

Author information

Authors and Affiliations

  1. Department of Geriatrics, The First Affiliated Hospital of Sun Yat-Sen University, Guangzhou, Guangdong, China

    Kemin Yan

  2. Department of Pulmonary and Critical Care Medicine, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China

    Yuxia Liang

  3. Department of Pulmonary and Critical Care Medicine, The First Affiliated Hospital of Sun Yat-Sen University, Guangzhou, Guangdong, China

    Yuxia Liang

Authors
  1. Kemin Yan
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yuxia Liang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

YL and KY conceived and designed the study. They also drafted and critically revised the manuscript. KY collected and analyzed the data for bioinformatic analysis. YL conducted the experiments, performed statistical analysis, and illustrated the results. All authors approved the manuscript’s final version.

Corresponding author

Correspondence to Yuxia Liang.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee of the First Affiliated Hospital of Sun Yat-sen University (2021071).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1: Figure S1.

Flow chart of the study design. Figure S2. Visualization of differentially expressed genes (DEGs). (A) Heatmap showed the expression of DEGs in each sample. (B) DEGs filtered by thresholds were presented in volcano map. Red dots represent upregulated genes and blue dots represent downregulated genes. Figure S3. Functional enrichment analysis of differentially expressed genes (DEGs). (A) GO analysis of DEGs. (B) KEGG analysis of DEGs. Figure S4. Scatter charts of the correlation of TLR7 and infiltrating immune cells. Spearman correlation of the correlation of TLR7 and infiltrating immune cells was performed. The results were presented in scatter charts.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yan, K., Liang, Y. Decreased TLR7 expression was associated with airway eosinophilic inflammation and lung function in asthma: evidence from machine learning approaches and experimental validation. Eur J Med Res 29, 116 (2024). https://doi.org/10.1186/s40001-023-01622-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40001-023-01622-5

Keywords

European Journal of Medical Research

ISSN: 2047-783X

Contact us