Skip to main content

Analysis of immune cell components and immune-related gene expression profiles in peripheral blood of patients with type 1 diabetes mellitus



Type 1 diabetes mellitus (T1DM) is a chronic autoimmune disease caused by severe loss of pancreatic β cells. Immune cells are key mediators of β cell destruction. This study attempted to investigate the role of immune cells and immune-related genes in the occurrence and development of T1DM.


The raw gene expression profile of the samples from 12 T1DM patients and 10 normal controls was obtained from Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified by Limma package in R. The least absolute shrinkage and selection operator (LASSO)—support vector machines (SVM) were used to screen the hub genes. CIBERSORT algorithm was used to identify the different immune cells in distribution between T1DM and normal samples. Correlation of the hub genes and immune cells was analyzed by Spearman, and gene-GO-BP and gene-pathway interaction networks were constructed by Cytoscape plug-in ClueGO. Receiver operating characteristic (ROC) curves were used to assess diagnostic value of genes in T1DM.


The 50 immune-related DEGs were obtained between the T1DM and normal samples. Then, the 50 immune-related DEGs were further screened to obtain the 5 hub genes. CIBERSORT analysis revealed that the distribution of plasma cells, resting mast cells, resting NK cells and neutrophils had significant difference between T1DM and normal samples. Natural cytotoxicity triggering receptor 3 (NCR3) was significantly related to the activated NK cells, M0 macrophages, monocytes, resting NK cells, and resting memory CD4+ T cells. Moreover, tumor necrosis factor (TNF) was significantly associated with naive B cell and naive CD4+ T cell. NCR3 [Area under curve (AUC) = 0.918] possessed a higher accuracy than TNF (AUC = 0.763) in diagnosis of T1DM.


The immune-related genes (NCR3 and TNF) and immune cells (NK cells) may play a vital regulatory role in the occurrence and development of T1DM, which possibly provide new ideas and potential targets for the immunotherapy of diabetes mellitus (DM).


According to the latest statistics from the International Diabetes Federation (IDF) (, the global prevalence of diabetes is about 9.3% (463 million people) in 2019, and it is forecasted to rise to 10.2% (578 million people) by 2030 and to 10.9% (700 million people) by 2045. Type 1 diabetes mellitus (T1DM) caused by autoimmune reaction is a major subtype of diabetes, and is mostly prevalent in adolescents and childhood [1]. Compared with type 2 diabetes mellitus (T2DM), T1DM has higher morbidity, mortality and health care cost [2, 3]. It is widely accepted that T1DM is the result of the interaction of both genetic and environmental factors, but its exact molecular mechanism is still unclear. Over the past 30 years, insulin therapy, immunotherapy and some potential therapies such as cell therapy are the main treatments for diabetes. However, these treatments usually have some problems, for instance, insulin therapy can cause a series of complications, while immunotherapy may impact acquired immunity and the efficacy is short-term [4,5,6].

To understand the mechanisms lying in the pathogenesis of T1DM, a wide range of studies have been carried out and found that immune responses play an important role in T1DM, which need the coordianated efforts of multiple immune related genes and various immune cells [7,8,9,10,11]. For instance, the onset phase of T1DM is characterized by the perturbations of NK cells. Patients with long-standing T1DM showed reduced NK cell activity due to decreased mRNA expression of the cell surface markers NKp30/p46 [also known as natural cytotoxicity triggering receptor 3/1 (NCR3/1)], as well as IFN-γ and perforin [12]. Th1 cells and regulatory T cells (Tregs), two hypotypes of CD4+ T cells, have also been reported to promote the pathogenesis of T1DM by destroying β cells via secreting cytokines like IL-1 [13].

However, the roles of immune related genes and immune cells, as well as the associations between them in T1DM have not been fully investigated. Thus, the current study aims to analyze the immune cell components, immune related genes and their correlations in peripheral blood of T1DM patients via bioinformatic strategies. Our findings will provide novel insights into diagnosis and improve our knowledge of the immunotherapies for T1DM patients.

Materials and methods

Dataset acquisition

The raw gene expression profile dataset (GSE55098) was obtained from the National Center of Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database ( The samples for this data were the peripheral blood mononuclear cells (PBMCs). The dataset based on GPL570 platform included 12 T1DM patients and 10 normal controls. The 1639 immune-related genes were downloaded from the ImmPort database (Immunology Database and Analysis Portal database,

Identification of immune-related DEGs

The dataset between T1DM patients and normal controls was identified using Limma package in R. The P < 0.05 and |Log2 fold change (FC)|≥ 0.585 were set as the cut-off for the DEGs. The DEGs were exhibited by the volcano plot and the heatmap. The DEGs were overlapped with the 1639 immune-related genes, getting 50 immune-related DEGs.

Analyses of the GO and KEGG

The GO enrichment and the KEGG pathway analyses were conducted using the immune-related DEGs respectively by GOseq package and hypergeometric test in R, which was used to explore the potential biological processes (BP), cellular components (CC), molecular functions (MF) and identified significantly relevant signal pathways of the immune-related DEGs. The histogram and bubble chart were plotted by the ggplot2 package (version 1.0.2) in R.

Analyses of LASSO and SVM

By constructing a penalty function, LASSO can compress variable coefficients and make some regression coefficients to become 0, thereby achieving the purpose of variable selection [14]. To screen the gene signatures, the 50 immune-related DEGs were performed the LASSO regression analysis in the glmnet package in R. SVM is supervised machine learning techniques widely used in pattern recognition and classification problems, which have been used in medical applications to predict whether a new gene falls into one category or the other, thereby classifying the genes [15]. The SVM was also used to screen the gene signatures. The overlapping genes after LASSO and SVM analyses were used as the hub genes.

CIBERSORT analysis

CIBERSORT algorithm could quantify the abundance of specific cell types [16]. To compare the difference between T1DM samples and normal samples in immune cells, the CIBERSORT analysis ( was used to estimate the percentage of LM22 (22 immune cell types) in each sample. Moreover, the fraction of 22 immune cells was compared between T1DM samples and normal samples and the violin plot was drawn by vioplot package in R.

Enriched GO-BP and KEGG network by 50 immune-related DEGs

Correlation of the hub genes and immune cells was analyzed by Spearman and visualized by heatmap in T1DM and normal samples. To systematically explore potential functions between the hub genes, 50 immune-related DEGs were imported into the Cytoscape software v3.7.2 ( to construct the genes and pathways interaction network by ClueGO plug-in. ClueGO was used to decipher the functionally grouped GO and pathway annotation networks to understand their implication in three different classifications (BP, MF and CC), in addition to the KEGG signaling pathway. The relationship between the terms was calculated using κ statistics and the ClueGO network was built based on the similarity of their related genes. In the present study, the enrichment analysis of gene-BP and gene-pathway was statistically validated using the ClueGO plug-in. BPs/signaling pathways were functionally split into several groups with κ score ≥ 0.4. In network, the node represented a BP/pathway, and the edge between two nodes indicated that the two BPs/pathways shared common genes.

Statistical analysis

The WilcoxTest were used to compare the fraction of the immune cells between T1DM samples and normal samples in CIBERSORT analysis. ROC curves were used to assess the diagnostic value of genes in T1DM,and the higher the AUC value is, the stronger the diagnostic value will be. The P < 0.05 was considerable as the statistical significance.


Screening of DEGs between the T1DM samples and normal samples

The 216 DEGs between the T1DM samples and normal samples were screened out (Fig. 1a), including 92 up-regulated genes (Additional file 1: Table S1) and 124 down-regulated genes (Additional file 2: Table S2). The distribution of the DEGs was displayed by the volcano plot (Fig. 1b), and the expressions of the DEGs in each sample were shown in Fig. 1c.

Fig. 1

The amount of 216 DEGs between the T1DM samples and normal samples. a Histogram; b volcano plot; c heatmap

The raw identification and functional enrichment analysis of the immune-related DEGs

The 1639 immune-related genes were overlapped with the DEGs using the Venn diagram, getting 50 immune-related DEGs (Fig. 2a) that were used for subsequent analysis.

Fig. 2

Raw identification and functional enrichment analysis of the immune-related DEGs. a Venn diagram of overlapping genes; b GO terms; c KEGG pathway

To further investigate the functions of the immune-related DEGs, the GO and KEGG analyses were performed. The GO terms with the top 5 of gene enrichment were protein binding, membrane, extracellular region, extracellular space, and integral component of membrane (Fig. 2b and Additional file 3: Table S3), suggesting that the immune-related DEGs may be involved in the biological activities of the cell membrane. Meanwhile, the analysis of KEGG pathway indicated that these genes were mainly enriched in Type 1 diabetes mellitus pathway (Fig. 2c and Additional file 4: Table S4), further verifying that the immune-related DEGs played a vital role in the occurrence and development of T1DM.

Identification of the optimal immune-related biomarkers

To identify immune-related hub genes, the LASSO regression analysis for the 50 immune-related DEGs was performed to screen the gene signatures, getting 11 gene signatures (Fig. 3a, b). Besides, the SVM for the 50 immune-related DEGs was also used to screen gene signatures, getting 6 gene signatures (Fig. 3c). Subsequently, the 11 gene signatures identified by LASSO were overlapped with the 6 gene signatures identified by the SVM, and ultimately obtaining 5 hub genes [C–C motif chemokine receptor 3 (CCR3), major histocompatibility complex, class II, DQ alpha 1 (HLA-DQA1), NCR3, toll like receptor 3 (TLR3) and tumor necrosis factor (TNF), Fig. 3d], which were considered as the optimal immune-related biomarkers.

Fig. 3

Identification of the optimal immune-related biomarkers. a, b LASSO regression analysis; c SVM analysis; d venn diagram of overlapping genes

Distribution of immune cells between T1DM and normal samples

To further investigate the correlation of the hub genes related to immune with immune cells, the distribution of the 22 immune cells between T1DM and normal samples was analyzed by the CIBERSORT algorithm. Figure 4a showed the percentage of immune cells in each sample, indicating that Monocytes, resting NK cells, resting memory CD4+ T cells, and CD8+ T cells had a larger proportion (Additional file 5: Table S5). Removing the four types of immune cells that were not in the sample, the proportion of the other immune cells in each sample was displayed by a heatmap (Fig. 4b). According to statistical results, the distribution of plasma cells (P = 0.05), resting mast cells (P = 0.013) and neutrophils (P = 0.018) in T1DM samples was significantly increased, while the distribution of resting NK cells (P = 0.009) in T1DM samples was significantly reduced, compared with that in normal samples (Fig. 4c). The above results suggested that these differential immune cells might be involved in the immune regulation process of T1DM pathogenesis.

Fig. 4

Distribution of immune cells between T1DM and normal samples. a Percentage of immune cells in each sample; b heatmap; c vioplot

Correlation analysis of the immune cells and the hub genes

To further analyze the correlation of the hub genes and the immune cells, the Spearman correlation heatmaps were plotted respectively in T1DM and normal samples. As shown in Fig. 5a, b, the differential immune cells were marked with red. NCR3 was only significantly related to the activated NK cells (P < 0.05) in the normal samples, while was significantly associated with multiple immune cells including M0 macrophages, monocytes, resting NK cells, and resting memory CD4+ T cells (all P < 0.05) in the T1DM samples. Moreover, TNF had no significant correlations to all immune cells (all P > 0.05) in the normal samples, but was significantly related to naive B cell (P < 0.05) and naive CD4+ T cell (P < 0.01) in the T1DM samples. Additionally, there were significant differences between T1DM and normal samples that TLR3, HLA-DQA1 and CCR3 were only related to CD8+ T cells (P < 0.05), resting mast cells (P < 0.05) and Neutrophils (P < 0.05), respectively. These results suggested that NCR3 and TNF might play an important role in the regulation of immune cell-mediated T1DM progression.

Fig. 5

Spearman correlation of the hub genes and the immune cells. a normal samples; b T1DM samples

Construction of the hub genes and pathways network

To in-depth investigate functions of the hub genes, a network of the 50 immune-related DEGs and GO-BP interaction was constructed using ClueGO Plug-in of Cytoscape software, in which the hub genes were marked as red frame. As shown in Fig. 6a, as expected, the pathways that interacted with TNF were the most, mainly cellular response to interferon-gamma-related anti-infection pathways. Next was NCR3, which was mainly related to natural killer cell mediated immunity-related innate immune response pathways. Further, Fig. 6b also indicated that proportions of genes in cellular response to interferon-gamma and natural killer cell mediated immunity terms were the largest, 59.52% and 16.67% respectively. To more clearly show the interaction of TNF and NCR3 with pathways, the up-regulated genes and down-regulated genes were separated to construct networks which were shown in Fig. 6c, d, e–f. Similar results were obtained, TNF participated in antimicrobial humoral response-related anti-infection pathways, NCR3 was still involved in natural killer cell mediated immunity. Besides, the network of the 50 immune-related DEGs and KEGG interaction also indicated similar results (Additional file 6: Figure S1).

Fig. 6

The network of hub genes and GO-BP interaction. a All genes; b the pie chart of all genes; c up-regulated genes; d the pie chart of up-regulated genes; e down-regulated genes; f the pie chart of down-regulated genes

The above results suggested that TNF and NCR3 might regulate T1DM progression respectively by anti-infection pathways and natural killer cell mediated immunity. NK cells recognize and kill virus-infected cells in the absence of antibodies and major histocompatibility complex (MHC), allowing for a much faster immune reaction. The role of NK cells in both the innate and adaptive immune responses is becoming increasingly important in research using NK cell activity as a potential therapy [17].

Diagnostic value of TNF and NCR3

To explore the accuracy of the TNF and NCR3 as the diagnostic biomarkers for T1DM, the ROC curves were plotted, respectively. The AUC was 0.763 (TNF, Fig. 7a) and 0.918 (NCR3, Fig. 7b), suggesting that NCR3 possessed a higher accuracy than TNF in diagnosis of T1DM. And this might be related to NK cell mediated innate immune system which was often ignored in the design of novel immune-based therapies [18].

Fig. 7

The ROC curves. a TNF; b NCR3


Immune related genes and immune cells are crucial in the pathogenesis of T1DM. However, their roles and interactions in T1DM have yet to be clarified. In this study, we identified five candidate immune-related biomarkers for T1DM and four types of immune cells differentially distributed between T1DM patients and the normal controls. Furthermore, we explored the correlations between these immune-related biomarkers and immune cells.

Firstly, we obtained 50 immune-related genes involved in T1DM by intersecting 1639 immune-related genes with 216 DEGs. Functional analysis revealed that these genes mainly participate in T1DM pathway, suggesting that T1DM is closely related to immunity. Besides, we also found some pathways other than T1DM or immunity, such as positive regulation of cytosolic calcium ion concentration and G protein-coupled receptor signaling pathway in GO enrichment analysis. Calcium ions play an important role in the maintenance of cell membrane biopotential and nerve conduction. Washburn RL et al. noted that C-peptide, a cleavage product of insulin processing unproduced by patients with T1DM, can elevate calcium concentrations by binding G protein-coupled receptors and thus treat the complications of T1DM [19]. As well, the Jak-STAT signaling pathway was present in the results of the KEGG pathway analysis. Jak-STAT signaling pathway activation can regulate the transcription of downstream genes and cell proliferation, differentiation and apoptosis process. Gurzov et al. [20] observed that the highly conserved and efficient Jak-STAT signaling pathway is necessary for normal homeostasis, and when dysregulated, it leads to the development of obesity and diabetes.

To get more robust immune-related biomarkers in T1DM, we performed LASSO and SVM analyses and found five candidate immune-related biomarkers, including CCR3, HLA-DQA1, NCR3, TLR3 and TNF. CCR3 encodes a receptor for C–C type chemokines which belongs to family 1 of the G protein-coupled receptors [21]. Chemokine receptors play an important role in the extracellular infiltration of inflammatory factors into the inflamed tissue in T1DM insulitis [22]. According to the study of Lohmann et al. [23], there was no difference for the Th2-associated chemokine receptor CCR3 in CD3+ lymphocytes between newly diagnosed T1DM patients, long-existing T1DM patients and healthy individuals. Likewise, another study also showed that there was no significant difference in the percentages of CCR3 in recently activated circulating T cells (CD3+, HLA-DR+) between newly diagnosed T1DM patients, long standing T1DM patients and health individuals [24]. In this paper, CCR3 was differentially expressed in T1DM and normal samples, probably due to the different cell types examined in this study compared to the above studies. Thus, CCR3 may indirectly regulate the pathological process of T1DM by some mechanisms which require further research. HLA-DQA1 belongs to the HLA class II α chain paralogues. The higher expression of HLA-DQA1 in T1DM than in normal samples affected the encoding level of DQ2.5 and DQ8 molecules on the APC surface, which can promote the presentation of self-antigens to induce autoimmunity, and present gluten antigens to homologous CD4+ T cells, thereby promoting T cell activation and proliferation [25]. In our study, HLA-DQA1 was also highly expressed in T1DM patients, suggesting that it may be involved in the pathogenesis of T1DM by inducing autoimmunity and could potentially be a therapeutic target. The protein family encoded by NCR3 is an active receptor that conveys effective signals to NK cells to lyse harmful cells and produce inflammatory cytokines [26]. Rodacki et al. [12] reported that the expression of NCR3 is reduced in patients with long-standing T1DM, which is consistent with our results. TLR3 plays a fundamental role in pathogen recognition and innate immune activation by recognizing dsRNA derived from viral replication. In some genetically susceptible individuals, this defence system does not work properly and instead induces excessive progressive inflammation and prolonged cell death, leading to the development of T1DM [27], which may explain the lower expression of TLR3 in T1DM patients found in our study. TNF encodes a multifunctional pro-inflammatory cytokine secreted basically by macrophages. Qiao et al. [28] confirmed that serum TNF level in T1DM patients significantly elevated among all age, disease duration and ethnicity groups, which is similar to our result observed in PBMCs, suggesting its important role in T1DM process.

Furthermore, we analyzed the distribution of immune cells in T1DM patients and normal ones, and found that plasma cells, resting mast cells, neutrophils and resting NK cells were differentially distributed between the T1DM and the normal controls, indicating these immune cells are more important in the occurrence and development of T1DM. For plasma cells, Isabel et al. [29] found that they were increased in the thymus of non-obese diabetic (NOD) mice compared to the control ones. In another study, plasma cells were also increased in patients with diabetic nephropathy (DN). They found a positive correlation between the number of plasma cells and serum IgG levels in DN patients, probably due to antigenic stimulation promoting the activation and differentiation of naive B cells and memory B cells towards plasma cells, leading to IgG production in DN patients [30]. Our results also showed that plasma cells were increased in the samples of T1DM patients compared to the normal ones, suggesting that it may play a role in the antigen presentation process. Mast cells are a type of innate immune cells that express MHC molecules [31]. It has been demonstrated that mast cells negatively regulate T1DM and other autoimmune related diseases [32, 33]. Our findings revealed increased resting mast cell in T1DM, which also points out from a cell subtype perspective that mast cells may negatively regulate T1DM through participating in MHC expression. Neutrophils gelatinase-associated lipocalin (NGAL) has been shown to be elevated in T1DM patients and has potential to become a biomarker for DN [34, 35]. Studies suggest that neutrophil cells are involved in T1DM pathological damage through the formation of neutrophil extracellular traps (NETs) [36, 37]. Likewise, our study also indicated that neutrophils are highly expressed in T1DM. Rodacki et al. found that recently diagnosed patients with T1DM may have more, or more active, circulating NK cells, whereas patients with the long-standing T1DM had reduced cellular activity of NK cells [12]. Likewise, Lorini et al. also indicated that NK cell cytotoxic activity was reduced in patients with long-standing T1DM [38]. The reduced activity of NK cells in patients with long-standing T1DM suggests that the reduced activity of NK cells is a consequence of the disease rather than a cause. In our study, we showed that NK cell resting was less prevalent in T1DM from a cell subtype perspective, which may be due to the fact that samples in our study were from newly diagnosed T1DM patients rather than long-term T1DM patients.

In consideration that immunity requires the coordinated efforts of immune related genes and immune cells, we analyzed the correlations between five T1DM biomarkers and immune cells. We found a significantly positive correlation between TNF and naive B cells and a significantly negative correlation between TNF and naive CD4+ T cells in T1DM samples. Lee et al. reported that TNF promotes the proliferation of naive CD4+ T cells in NOD mice. It is probably due to the differences of the sources and treatment of the samples [39]. NCR3 was only negatively associated with activated NK cells in normal samples, while positively related to resting memory CD4+ T cells and negatively correlated with M0 macrophages, monocytes and resting NK cells in T1DM samples. It has been reported that NCR3 positively regulates NK cell activity in T1DM [12], which may explain the negative correlation between NCR3 and resting NK cells in T1DM observed in our study. Ward et al. [40] showed that a few HIV-infected CD4+ T cells expressed NCR3 ligands, leading to lysis of CD4+ T cells by NK cells. In our study, NCR3 and resting memory CD4+ T cells were positively correlated in patients with T1DM, by a potential mechanism that may be similar to the above study. The results of this study regarding the negative correlation of NCR3 with M0 macrophages and monocytes in T1DM patients may be a novel finding. Taken together, TNF may affect T1DM by regulating naive B cells and naive CD4+ T cells, while NCR3, an NK cell receptor, may affect T1DM primarily by delivering effective signals to NK cells and enable them to lyse target cells.

Interestingly, in the network of hub genes and GO-BP, we found most biological processes were mostly interacted with TNF and NCR3, and it is also a sign that they play an essential role in T1DM. Also, the network revealed the possible mechanisms of TNF and NCR3 in regulating T1DM, by anti-infection pathways and natural killer cell mediated immunity, respectively. Finally, given the importance of the above two genes, we assessed their diagnostic value by ROC and found a strong accuracy for the NCR3 (AUC = 0.918), which also provides a new biomarker for the diagnosis of T1DM.

Although we used multiple analyses to systematically investigate immune related genes, immune cells and their relationships in T1DM, this study still has some limitations due to the lack of experimental validation, single dataset and incomplete database information. Further optimization and fundamental experiments are required to reveal the detailed molecular mechanisms of these immune related genes in T1DM.


In this study, we identified CCR3, HLA-DQA1, NCR3, TLR3 and TNF as potential immune-related biomarkers in T1DM. For the first time, the current study revealed the associations between immune related genes and immune cells. Our findings improve the understanding of molecular mechanisms involved in T1DM and provide novel information of diagnose and therapy for T1DM patients.

Availability of data and materials

The datasets supporting the conclusions of this article are available in the Gene Expression Omnibus (GEO) database and the ImmPort database (Immunology Database and Analysis Portal database).



Area under curve


Biological processes


Cellular components


C–C motif chemokine receptor 3


Differentially expressed genes


Diabetes mellitus


Diabetic nephropathy


Fold change


Gene Expression Omnibus


Major histocompatibility complex, class II, DQ alpha 1


International Diabetes Federation

ImmPort database:

Immunology Database and Analysis Portal database


Least absolute shrinkage and selection operator


Molecular functions


National Center of Biotechnology Information


Natural cytotoxicity triggering receptor 3


Neutrophil extracellular traps


Neutrophils gelatinase-associated lipocalin


Non-obese diabetic


Peripheral blood mononuclear cells


Receiver operating characteristic


Support vector machines


Type 1 diabetes mellitus


Type 2 diabetes mellitus


Toll like receptor 3


Tumor necrosis factor


Regulatory T cells


  1. 1.

    Katsarou A, Gudbjornsdottir S, Rawshani A, Dabelea D, Bonifacio E, Anderson BJ, et al. Type 1 diabetes mellitus. Nat Rev Dis Primers. 2017;3:17016.

    Article  Google Scholar 

  2. 2.

    Ogurtsova K, da Rocha Fernandes JD, Huang Y, Linnenkamp U, Guariguata L, Cho NH, et al. IDF Diabetes Atlas: global estimates for the prevalence of diabetes for 2015 and 2040. Diabetes Res Clin Pract. 2017;128:40–50.

    CAS  Article  Google Scholar 

  3. 3.

    Rewers M. Challenges in diagnosing type 1 diabetes in different populations. Diabetes Metab J. 2012;36(2):90–7.

    Article  Google Scholar 

  4. 4.

    Kroger CJ, Clark M, Ke Q, Tisch RM. Therapies to suppress beta cell autoimmunity in type 1 diabetes. Front Immunol. 2018;9:1891.

    Article  Google Scholar 

  5. 5.

    DiMeglio LA, Evans-Molina C, Oram RA. Type 1 diabetes. Lancet. 2018;391(10138):2449–62.

    Article  Google Scholar 

  6. 6.

    Skyler JS. Hope vs hype: where are we in type 1 diabetes? Diabetologia. 2018;61(3):509–16.

    Article  Google Scholar 

  7. 7.

    Vallianou NG, Stratigou T, Geladari E, Tessier CM, Mantzoros CS, Dalamaga M. Diabetes type 1: can it be treated as an autoimmune disorder? Rev Endocr Metab Disord. 2021.

    Article  PubMed  Google Scholar 

  8. 8.

    Babon JA, DeNicola ME, Blodgett DM, Crevecoeur I, Buttrick TS, Maehr R, et al. Analysis of self-antigen specificity of islet-infiltrating T cells from human donors with type 1 diabetes. Nat Med. 2016;22(12):1482–7.

    CAS  Article  Google Scholar 

  9. 9.

    Michels AW, Landry LG, McDaniel KA, Yu L, Campbell-Thompson M, Kwok WW, et al. Islet-derived CD4 T cells targeting proinsulin in human autoimmune diabetes. Diabetes. 2017;66(3):722–34.

    CAS  Article  Google Scholar 

  10. 10.

    Smith MJ, Simmons KM, Cambier JC. B cells in type 1 diabetes mellitus and diabetic kidney disease. Nat Rev Nephrol. 2017;13(11):712–20.

    CAS  Article  Google Scholar 

  11. 11.

    Marek-Trzonkowska N, Mysliwiec M, Iwaszkiewicz-Grzes D, Gliwinski M, Derkowska I, Zalinska M, et al. Factors affecting long-term efficacy of T regulatory cell-based therapy in type 1 diabetes. J Transl Med. 2016;14(1):332.

    Article  Google Scholar 

  12. 12.

    Rodacki M, Svoren B, Butty V, Besse W, Laffel L, Benoist C, et al. Altered natural killer cells in type 1 diabetic patients. Diabetes. 2007;56(1):177–85.

    CAS  Article  Google Scholar 

  13. 13.

    Mandrup-Poulsen T, Pickersgill L, Donath MY. Blockade of interleukin 1 in type 1 diabetes mellitus. Nat Rev Endocrinol. 2010;6(3):158–66.

    CAS  Article  Google Scholar 

  14. 14.

    Vasquez MM, Hu C, Roe DJ, Chen Z, Halonen M, Guerra S. Least absolute shrinkage and selection operator type methods for the identification of serum biomarkers of overweight and obesity: simulation and application. BMC Med Res Methodol. 2016;16(1):154.

    Article  Google Scholar 

  15. 15.

    Wang ZL, Zhou ZG, Chen Y, Li XT, Sun YS. Support vector machines model of computed tomography for assessing lymph node metastasis in esophageal cancer with neoadjuvant chemotherapy. J Comput Assist Tomogr. 2017;41(3):455–60.

    Article  Google Scholar 

  16. 16.

    Felix Garza ZC, Lenz M, Liebmann J, Ertaylan G, Born M, Arts ICW, et al. Characterization of disease-specific cellular abundance profiles of chronic inflammatory skin conditions from deconvolution of biopsy samples. BMC Med Genomics. 2019;12(1):121.

    Article  Google Scholar 

  17. 17.

    Watzl C. How to trigger a killer: modulation of natural killer cell reactivity on many levels. Adv Immunol. 2014;124:137–70.

    Article  Google Scholar 

  18. 18.

    Fraker C, Bayer AL. The expanding role of natural killer cells in type 1 diabetes and immunotherapy. Curr Diab Rep. 2016;16(11):109.

    Article  Google Scholar 

  19. 19.

    Washburn RL, Mueller K, Kaur G, Moreno T, Moustaid-Moussa N, Ramalingam L, et al. C-peptide as a therapy for type 1 diabetes mellitus. Biomedicines. 2021;9(3):270.

    CAS  Article  Google Scholar 

  20. 20.

    Gurzov EN, Stanley WJ, Pappas EG, Thomas HE, Gough DJ. The JAK/STAT pathway in obesity and diabetes. FEBS J. 2016;283(16):3002–15.

    CAS  Article  Google Scholar 

  21. 21.

    Combadiere C, Ahuja SK, Murphy PM. Cloning and functional expression of a human eosinophil CC chemokine receptor. J Biol Chem. 1995;270(28):16491–4.

    CAS  Article  Google Scholar 

  22. 22.

    Bradley LM, Asensio VC, Schioetz LK, Harbertson J, Krahl T, Patstone G, et al. Islet-specific Th1, but not Th2, cells secrete multiple chemokines and promote rapid induction of autoimmune diabetes. J Immunol. 1999;162(5):2511–20.

    CAS  PubMed  Google Scholar 

  23. 23.

    Lohmann T, Laue S, Nietzschmann U, Kapellen TM, Lehmann I, Schroeder S, et al. Reduced expression of Th1-associated chemokine receptors on peripheral blood lymphocytes at diagnosis of type 1 diabetes. Diabetes. 2002;51(8):2474–80.

    CAS  Article  Google Scholar 

  24. 24.

    Baker C, Chang L, Elsegood KA, Bishop AJ, Gannon DH, Narendran P, et al. Activated T cell subsets in human type 1 diabetes: evidence for expansion of the DR+ CD30+ subpopulation in new-onset disease. Clin Exp Immunol. 2007;147(3):472–82.

    CAS  Article  Google Scholar 

  25. 25.

    Farina F, Picascia S, Pisapia L, Barba P, Vitale S, Franzese A, et al. HLA-DQA1 and HLA-DQB1 alleles, conferring susceptibility to celiac disease and type 1 diabetes, are more expressed than non-predisposing alleles and are coordinately regulated. Cells. 2019;8(7):751.

    CAS  Article  Google Scholar 

  26. 26.

    Hudspeth K, Silva-Santos B, Mavilio D. Natural cytotoxicity receptors: broader expression patterns and functions in innate and adaptive immune cells. Front Immunol. 2013;4:69.

    Article  Google Scholar 

  27. 27.

    Assmann TS, Brondani Lde A, Boucas AP, Canani LH, Crispim D. Toll-like receptor 3 (TLR3) and the development of type 1 diabetes mellitus. Arch Endocrinol Metab. 2015;59(1):4–12.

    Article  Google Scholar 

  28. 28.

    Qiao YC, Chen YL, Pan YH, Tian F, Xu Y, Zhang XX, et al. The change of serum tumor necrosis factor alpha in patients with type 1 diabetes mellitus: a systematic review and meta-analysis. PLoS ONE. 2017;12(4):e0176157.

    Article  Google Scholar 

  29. 29.

    Pinto AI, Smith J, Kissack MR, Hogg KG, Green EA. Thymic B cell-mediated attack of thymic stroma precedes type 1 diabetes development. Front Immunol. 2018;9:1281.

    Article  Google Scholar 

  30. 30.

    Zhang N, Tai J, Qu Z, Zhang Z, Zhao S, He J, et al. Increased CD4(+)CXCR5(+)T follicular helper cells in diabetic nephropathy. Autoimmunity. 2016;49(6):405–13.

    CAS  Article  Google Scholar 

  31. 31.

    Agier J, Brzezinska-Blaszczyk E, Witczak P, Kozlowska E, Zelechowska P. The impact of TLR7 agonist R848 treatment on mast cell phenotype and activity. Cell Immunol. 2021;359:104241.

    CAS  Article  Google Scholar 

  32. 32.

    Geoffrey R, Jia S, Kwitek AE, Woodliff J, Ghosh S, Lernmark A, et al. Evidence of a functional role for mast cells in the development of type 1 diabetes mellitus in the BioBreeding rat. J Immunol. 2006;177(10):7275–86.

    CAS  Article  Google Scholar 

  33. 33.

    Benoist C, Mathis D. Mast cells in autoimmune disease. Nature. 2002;420(6917):875–8.

    CAS  Article  Google Scholar 

  34. 34.

    Lacquaniti A, Donato V, Pintaudi B, Di Vieste G, Chirico V, Buemi A, et al. “Normoalbuminuric” diabetic nephropathy: tubular damage and NGAL. Acta Diabetol. 2013;50(6):935–42.

    CAS  Article  Google Scholar 

  35. 35.

    Kapoula GV, Kontou PI, Bagos PG. Diagnostic accuracy of neutrophil gelatinase-associated lipocalin for predicting early diabetic nephropathy in patients with type 1 and type 2 diabetes mellitus: a systematic review and meta-analysis. J Appl Lab Med. 2019;4(1):78–94.

    Article  Google Scholar 

  36. 36.

    Lang J, Wang X, Liu K, He D, Niu P, Cao R, et al. Oral delivery of staphylococcal nuclease by Lactococcus lactis prevents type 1 diabetes mellitus in NOD mice. Appl Microbiol Biotechnol. 2017;101(20):7653–62.

    CAS  Article  Google Scholar 

  37. 37.

    Njeim R, Azar WS, Fares AH, Azar ST, Kfoury Kassouf H, Eid AA. NETosis contributes to the pathogenesis of diabetes and its complications. J Mol Endocrinol. 2020;65(4):R65–76.

    CAS  Article  Google Scholar 

  38. 38.

    Lorini R, Moretta A, Valtorta A, d’Annunzio G, Cortona L, Vitali L, et al. Cytotoxic activity in children with insulin-dependent diabetes mellitus. Diabetes Res Clin Pract. 1994;23(1):37–42.

    CAS  Article  Google Scholar 

  39. 39.

    Lee LF, Xu B, Michie SA, Beilhack GF, Warganich T, Turley S, et al. The role of TNF-alpha in the pathogenesis of type 1 diabetes in the nonobese diabetic mouse: analysis of dendritic cell maturation. Proc Natl Acad Sci U S A. 2005;102(44):15995–6000.

    CAS  Article  Google Scholar 

  40. 40.

    Ward J, Bonaparte M, Sacks J, Guterman J, Fogli M, Mavilio D, et al. HIV modulates the expression of ligands important in triggering natural killer cell cytotoxic responses on infected primary T-cell blasts. Blood. 2007;110(4):1207–14.

    CAS  Article  Google Scholar 

Download references


We thank the researchers who gave their data for this analysis. And we also would like to acknowledge the GEO and ImmPort database.


The study was supported by National Natural Science Foundation of China (Grants No. 81370497 and 81500635) and the Science and Technology Department of Jilin Province (Grants No. 20190304067YY).

Author information




(I) Conception and design: JL, YL and JM; (II) Administrative support: PJ and JM; (III) Collection and assembly of data: JL and BW; (IV) Data analysis and interpretation: JL and YL; (V) Manuscript writing: all authors. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Ping Jiao or Jie Ma.

Ethics declarations

Ethics approval and consent to participate

The data of this study are from GEO and ImmPort database, and do not involve animal experiments and human specimens, no ethics-related issues.

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.

The 92 up-regulated genes

Additional file 2.

 The 124 down-regulated genes

Additional file 3.

 GO terms

Additional file 4.

 KEGG pathway

Additional file 5.

 Distribution of immune cells

Additional file 6.

 The network of hub genes and KEGG interaction

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

Verify currency and authenticity via CrossMark

Cite this article

Lin, J., Lu, Y., Wang, B. et al. Analysis of immune cell components and immune-related gene expression profiles in peripheral blood of patients with type 1 diabetes mellitus. J Transl Med 19, 319 (2021).

Download citation


  • Type 1 diabetes mellitus
  • Peripheral blood
  • Immune cells
  • Gene expression profiles