Immune cell infiltration characteristics and related core genes in lupus nephritis: results from bioinformatic analysis
BMC Immunology volume 20, Article number: 37 (2019)
Lupus nephritis (LN) is a common complication of systemic lupus erythematosus that presents a high risk of end-stage renal disease. In the present study, we used CIBERSORT and gene set enrichment analysis (GSEA) of gene expression profiles to identify immune cell infiltration characteristics and related core genes in LN.
Datasets from the Gene Expression Omnibus, GSE32591 and GSE113342, were downloaded for further analysis. The GSE32591 dataset, which included 32 LN glomerular biopsy tissues and 14 glomerular tissues from living donors, was analyzed by CIBERSORT. Different immune cell types in LN were analyzed by the Limma software. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis based on GSEA were performed by clusterProfiler software. Lists of core genes were derived from Spearman correlation between the most significant GO term and differentially expressed immune cell gene from CIBERSORT. GSE113342 was employed to validate the association between selected core genes and clinical manifestation. Five types of immune cells revealed important associations with LN, and monocytes emerged as having the most prominent differences. GO and KEGG analyses indicated that immune response pathways are significantly enriched in LN. The Spearman correlation indicated that 15 genes, including FCER1G, CLEC7A, MARCO, CLEC7A, PSMB9, and PSMB8, were closely related to clinical features.
This study is the first to identify immune cell infiltration with microarray data of glomeruli in LN by using CIBERSORT analysis and provides novel evidence and clues for further research of the molecular mechanisms of LN.
Systemic lupus erythematosus (SLE), one of the most complicated autoimmune diseases in the world, is caused by various endogenous antigens . Lupus nephritis (LN), a common and serious complication of SLE, is characterized by hematuria, proteinuria, and impaired glomerular filtration rate . The lack of understanding regarding the molecular mechanisms of LN hinders the development of specific targeted therapy for this progressive disease . Tracking the biological changes in LN at the genomic level is a worthwhile strategy . In recent years, gene sequencing technology combined with bioinformatic analysis has been conducted to identify genes relevant to diseases that might serve as prognostic biomarkers and be developed as therapeutic targets in the future . Bioinformatic analysis can process large amounts of samples within an extremely short time and provide valuable information about diseases, and several genes closely associated with SLE have been identified and driven research innovations in recent years [6,7,8]. However, few studies utilized bioinformatic analysis to characterize kidney tissue in the context of LN.
Many previous works found that immune cell infiltration is associated with treatment and clinical outcome in different types of cancer [9, 10]. Immune cells consisting of innate and adaptive immune populations, including dendritic cells, macrophages, neutrophils, T cells, and B cells, are associated with active and suppressive immune functions . However, given the functionally distinct cell types that comprise the immune response, assessing immune infiltration and determining whether differences in the composition of the immune infiltration can improve the development of novel immunotherapeutic drugs to target these cells is important. The CIBERSORT algorithm is an analytical tool whereby RNA-seq data can be used to assess the expression changes of immune cells and obtain the proportion of various types of immune cells from the samples. CIBERSORT offers 22 cell types encompassing monocytes, natural killer cells, B cells, T cells, eosinophils, macrophages, neutrophils, plasma cells, dendritic cells, and mast cells . It has been prevalently used to determine the immune cell landscapes in many malignant tumors such as breast cancer, hepatocellular carcinoma, and colorectal cancer [13,14,15]. In SLE pathogenesis, various immune cells have been widely evaluated and demonstrated to be harmful . Immune cell infiltration is also a hallmark of LN. Immune cells, such as monocytes, B cells, and T cells, are recruited to kidney tissue and produce cytokines and chemokines to cause tissue damage . However, the landscape of immune infiltration in LN has not been entirely revealed.
Although LN can affect all components of the kidney, the glomerulus is the most suitable tissue and is closely related to the pathogenesis and treatment of the disease . In our present study, the microarray data were downloaded from the Gene Expression Omnibus (GEO) database. By using CIBERSORT, we first investigated the difference in immune infiltration between LN kidney tissue and normal tissue in 22 subpopulations of immune cells. Gene set enrichment analysis (GSEA) was employed for functional enrichment analyses and to determine the most significant functional terms. A list of genes closely related to immune infiltration was screened out and validated against another dataset with clinical information from the GEO database. This study aimed to describe the characteristics of LN glomerular immune infiltration for the first time and to identify some key genes related to immune infiltration that affect clinical manifestation, so as to provide data resources for future research.
Bioinformatic analysis workflows and data description
Our workflows are shown in Fig. 1. We first investigated the difference of immune cell infiltration between normal glomerular tissues and LN glomerular tissues. Next, we discovered the most significant GO and KEGG functional term by GSEA. We screened out a list of genes closely related to immune infiltration and validated these genes against the clinical data. A total of 46 samples from GSE32591 were used in this study, including 32 LN glomerular biopsy tissues and 14 glomerular tissues from living donors. After data processing, the expression matrix of 30 LN glomerular samples and 6 normal control glomerular samples was obtained by screening the immune cell infiltration. GSE113342 contained 14 biopsy kidney tissues and 6 normal tissues.
Performance of CIBERSORT
Figure 2a shows the proportions of immune cells in 36 kidney tissues. Obviously, monocytes accounted for the majority of all infiltrating cells, especially in LN tissue. The differential expressional proportion of immune infiltration cells in the LN and control groups is shown in Fig. 2b. Five types of immune cells, namely, memory B cells, M0 macrophages, monocytes, activated NK cells, and follicular helper T (Tfh) cells, were differentially expressed. Monocytes, M0 macrophages, and activated NK cells were upregulated in LN tissue. The P-values of the five types of immune cells were 0.014, 0.04, 0.0001, 0.037, and 0.022, respectively. Among them, the increase in monocytes was the most significant. Memory B cells and Tfh cells were downregulated. Figure 2c indicates the correlation between these differentially expressed types of immune cells. The five types of immune cells were weakly to moderately correlated. Monocytes were negatively correlated with memory B cells and Tfh cells (r = − 0.42 and r = − 0.42, respectively), which indicated that the function of monocytes, Tfh cells, and memory B cells in LN may be antagonistic. However, the relationship between memory B cells and Tfh cells was synergistic.
GSEA-based GO analysis
On the basis of the GO biological process, the top 10 most significantly enriched GO terms are presented in Fig. 3a. Genes in GO terms were primarily associated with “activation of immune response (GO:0002253),” “chemotaxis (GO:0006935),” and “taxis (GO:0042330).” A total of 478 genes were involved in “activation of immune response.” These results confirmed that immune response is very important in LN. Our GO analysis presented numerous important genes associated with this function. The details of GO analysis are shown in Additional file 1: Table S1.
The connection between the most prominent GO terms is shown in Fig. 3b. The network-presented numerous genes, such as RSAD2, C1QA, C1QB, CX3CR1, ITGB2, FCER1G, and CCR1, that were significantly differentially expressed in LN. Moreover, ITGB2, FCER1G, C5AR1, LYN, CD36, and PTPRC were important bridge genes between different biological processes. We used all of the “activation of immune response” gene sets for GSEA, and the gene set enrichment result is presented in Fig. 3c. The enrichment showed that the gene set was enriched at the front of the sequence (ES = 0.61). Over 100 genes were core genes that increased during this process. We obtained the list of all core genes, such as C1QA, RSAD2, C1QB, ITGB2, HCK, C3AR1, FCN1 and FCER1G, for subsequent analysis.
GSEA-based KEGG analysis
A total of 24 prominent KEGG pathways including activated and suppressed pathways were selected (Fig. 4a). Activated pathways, such as “Epstein–Barr virus infection,” “Herpes simplex virus 1 infection,” “Influenza A,” “Human cytomegalovirus infection,” and “Kaposi sarcoma-associated herpesvirus infection,” were related to cellular immunity against viral infection. The result indicated that the activation of signaling pathways in LN is similar to that of viral infection. However, suppressed pathways were mainly concentrated on metabolic process, such as “Biosynthesis of amino acids,” “Valine, leucine and isoleucine degradation,” “Steroid hormone biosynthesis,” and “Oxidative phosphorylation.”
GSEA enrichment plots of representative gene sets on “Epstein–Barr virus infection” and “Biosynthesis of amino acids” are shown in Fig. 4b and 4c, respectively. In the activated pathway, 182 genes participated in the EB virus infection pathway and were concentrated at the front of the sequence. The core genes such as ISG15, OAS1, OAS2, OAS3, LYN, HLA-DQB1, and TLR2 were upregulated. In the suppressed pathway, only 62 functional genes were involved and were enriched at the back of the sequence.
Discovery of core genes
The correlation between core genes came from the GSEA GO term “activation of immune response” and five types of immune infiltrating cells are shown in Fig. 5a. A total of 44 genes showed close connection with immune infiltrating cells. Genes such as RSAD2, PSMB8, PSMA6, and MARCO were negatively correlated with Tfh cells. PLSCR1, ITGB2, HCK, and GBP1 were negatively related to memory B cells. FCN1, PSMB9, PRKCH, and A2M were positively correlated with monocytes. SYK, PYCARD, LPXN, and BTK were positively related to M0 macrophages. However, our analysis only found four genes correlated with activated NK cells.
Validation of core genes
Figure 5b shows the clinical information of GSE113342. The LN grade was mainly concentrated on 3–5 classes. The core gene list was validated in the clinical dataset. Grade, age, and 12-month response were chosen as clinical indicators (Fig. 5c). Through Spearman correlation analysis between core gene list and clinical information, GBP1, CD36, ITGB2, FCER1G, CLEC7A, LILRB4, HLA-DRA, BTK, PYCARD, CFP, CFD, PSMB9, MARCO, CD3D, and PSMB8 were found to be active in both networks, which indicated that these core genes were associated with immune infiltration and affected clinical manifestation. Among them, CLEC7A was positively correlated with age (r = 0.5) but negatively correlated with grade and 12-month response (r = − 0.56 and r = − 0.66, respectively). MARCO was negatively correlated with treatment response (r = − 0.53). FCER1G, PSMB8, and PSMB9 were positively correlated with treatment response (r = 0.58, r = 0.58, and r = 0.54, respectively).
With the development of bioinformatics, increasing attention has been focused on finding hub genes in various diseases, and the collected information on these genes can provide new means for exploring diseases. Multiple susceptibility genes may determine disease occurrence.
In this study, we uncovered different expressional cell patterns of immune infiltration in LN and association with clinical features. Monocytes were the prominent differentially expressed cells. These are important components of the innate immune system; they have an antigen presentation capacity and produce several inflammatory cytokines in SLE . Monocytes accounted for approximately 4% of blood leukocytes in healthy mice and over 50% in lupus-prone mice . Our result also showed that monocytes constituted 30–50% of immune cells in human LN glomeruli. Activated NK cells were also increased in glomeruli. However, reports from other studies showed lower proportions of NK cells in SLE patient blood, especially in patients with LN [21, 22]. However, in rheumatoid arthritis tissue, NK cells were reported to contradict the function of circulating NK cells, which indicated that tissue NK cells may have different effects as compared with blood NK cells in autoimmune disease . Clinical and experimental evidence indicated that aberrant memory B cells and Tfh cells played an important role in the pathogenesis of human SLE [24,25,26]. Resting M0 macrophages can polarize into M1 and M2 macrophages in the presence of the appropriate cytokines . However, no research has explained the function of increased M0 macrophages in LN. The specific role of these immune cells in functional immune responses still remains to be elucidated.
“Activation of immune response” was the top associated pathway under GSEA-based GO analysis. The activation of innate and adaptive immune system triggering immune complex deposition, complement activation, and self-antigen production displayed a toxic effect on renal glomerular and tubular cells, thereby promoting the development of nephritis in patients with SLE [28, 29]. Through KEGG pathway analysis, several kinds of virus infection pathways were associated with LN. The immunoreaction of LN and response to virus may share several common features.
By combining CIBERSORT results and “activation of immune response” GO term, we found many novel commonly expressed genes, some of which were important in autoimmune diseases. For example, FCN1 was proven to be associated with monocytes in patients with microscopic polyangiitis . Another study involving weighted correlation network analysis showed that RSAD2 related to CD4+ T cells may be the most highly ranked hub gene in SLE . BTK mediates TLR signaling in macrophages and may be a promising treatment approach for LN [31,32,33]. These genes were observed to be highly or mildly associated with immune cells in kidney tissues.
Through a review of documents about lupus and related genes [34,35,36,37,38,39,40,41,42,43,44,45,46,47,48], 15 core genes related to clinical manifestation were found to be associated in autoimmune disease (Table 1). FCER1G, CLEC7A, MARCO, CLEC7A, PSMB9, and PSMB8 showed apparent correlation with clinical manifestation. FCER1G, which is associated with multiple leukocyte receptor complexes and mediates signal transduction, plays a negative regulatory role in the B cell responses . CLEC7A, also known as dectin-1, is a type II membrane receptor expressed in the membrane of some leukocytes and likely contributes to the synthesis of pro-inflammatory cytokines in autoimmune conditions . MARCO, a scavenger receptor family, plays important roles in the clearance of apoptotic cells. The presence of anti-MARCO antibodies in SLE patients might contribute to the breakdown of self-tolerance and the pathogenesis of SLE . PSMB8 is involved in antigen processing and presentation in naïve CD4+ T cells, and PSMB9 is induced by interferon stimulation in SLE [41, 48]. All these core genes require additional studies to elucidate the complex interaction with clinical features.
The current work is the first to use CIBERSORT to analyze immune cell infiltration of glomerular tissue in LN. All data were derived from GEO and were therefore reliable. The correlation results of CIBERSORT and GSEA to obtain core genes were validated in clinical data, leading to many new information for our future research. The analytical methods were scientific and novel. However, our study has some limitations. Only a few datasets of LN were available on the GEO database; therefore, the number of samples included in this study was relatively small. However, despite the small sample sizes, we still found some significant differences among groups. In addition, clinical tests need to be conducted to support our results.
Our study provided a new insight into the immune filtration of LN. Five types of immune cells revealed important associations with LN, and monocytes showed the largest differences in the cellular composition of immune infiltration. Fifteen core genes that were related to clinical manifestation were analyzed. These genes may perform crucial functions, and further analysis of these genes in LN may identify targets for immunotherapy.
Microarray data processing
The data in our study came from a public domain. The normalized expression matrix and sample information were downloaded from the GEO database (www.ncbi.nlm.nih.gov/geo). We used “lupus nephritis” as a keyword for searching. The data selection criteria were as follows: (1) the study type was expression profiling by array; (2) the organisms must be Homo sapiens; (3) the samples of each dataset must include glomerular tissue. In accordance with the above criteria, the GSE32591 microarray dataset based on the Affymetrix Human GeneChip U133A (affy) platform was hit and adopted for CIBERSORT. The GSE113342 microarray dataset based on nCounter Nanostring Human Immunology v2 was used to demonstrate the association between selected genes and clinical feature later. Only 500 immune-related genes were detected in this dataset.
Evaluation of immune cell infiltration
Gene expression datasets of GSE32591 were processed to remove the null values. The missing values were supplemented by KNN method in “impute” package , the format was prepared in accordance with the accepted format of CIBERSORT, and then data were uploaded to the CIBERSORT web portal (http://cibersort.stanford.edu/). We used the original CIBERSORT gene signature file LM22, which defines 22 immune cell subtypes, to analyze datasets from human glomerular tissues and normal tissues. CIBERSORT p-value < 0.05 was included.
Differential analysis of immune cell infiltration types
To analyze the significant differential expression of different cell types of immune cells, we used the difference analysis between the disease group and the control group. Limma package and Bayesian method were used to construct a linear model . P-value < 0.05 was the cut-off standard. To further understand the relationship between these different types of immune cell infiltration, Pearson correlation coefficient was used to find the correlation between these differentially expressed types of immune cells.
GSEA is an analytical method for genome-wide expression profile microarray data. It can identify functional enrichment by comparing genes with predefined gene sets. A gene set is a group of genes that shares localization, pathways, functions, or other features. GSEA was conducted using clusterProfiler package (version 3.5) . The fold change of gene expression between LN group and control group was calculated, and the gene list was generated according to the change of |log2FC|. Then, we utilized GSEA-based enriched Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses.
GSEA-based enriched GO analysis
GO analysis includes three categories: molecular function, biological process, and cellular component. In the present study, we only selected biological process to perform GO analysis. GO analysis was performed through gseGO function in clusterProfiler package. The adjusted p-value < 0.05 was set as the cut-off criteria. The connections between the most significant GO terms and participating genes were visualized by GOenrich package with a network diagram.
GSEA-based KEGG pathway analysis
KEGG pathway enrichment analyses were also conducted by gseKEGG function in clusterProfiler package. The adjusted p-value < 0.05 was set as the cut-off criteria.
Core gene list and correlation analysis
The core gene list obtained in the most significant GO term was analyzed by Spearman correlation with the differentially expressed immune cells from CIBERSORT results. Five groups of correlation analysis data were obtained. P-value < 0.05 was used as the cut-off standard, and genes with the top 10 highest absolute values of correlation coefficients were visualized in each group.
Validation of core genes and association with clinical manifestations
In dataset GSE113342 with clinical information, patient part B was excluded because it was data after treatment, and only first renal biopsy data (patient part A), which had approximately 500 immune gene expression values that coincided with the genes obtained in the most significant GO term associated with immune response, were chosen for analysis. Gene intersection was calculated first, and the Spearman correlation analysis between these intersecting genes and clinical information, such as age, grade, and 12-month treatment response, was further applied.
Availability of data and materials
The datasets in the current study come from CEO database: GSE32591 and GSE113342.
Bruton Tyrosine Kinase
Complement C1q A Chain
Complement C1q B Chain
Complement C3a Receptor 1
Complement C5a Receptor 1
C-C Motif Chemokine Receptor 1
Complement Factor D
Complement Factor Properdin
C-Type Lectin Domain Containing 7A
C-Type Lectin Domain Containing 7A
C-X3-C Motif Chemokine Receptor 1
Fc Fragment of IgE Receptor Ig
Glycoprotein Hormone Alpha 2
Major Histocompatibility Complex, Class II, DQ Beta 1
Major Histocompatibility Complex, Class II, DR Alpha
ISG15 Ubiquitin Like Modifier
Integrin Subunit Beta 2
Leukocyte Immunoglobulin Like Receptor B4
Macrophage Receptor with Collagenous Structure
- OAS1, OAS2, OAS3:
2′-5′-Oligoadenylate Synthetase 1, 2, 3
Phospholipid Scramblase 1
Protein Kinase C Eta
Proteasome Subunit Beta 8
Proteasome Subunit Beta 9
Protein Tyrosine Phosphatase Receptor Type C
PYD And CARD Domain Containing
Radical S-Adenosyl Methionine Domain Containing 2
Spleen Associated Tyrosine Kinase
Toll Like Receptor 2
Steinmetz OM, Turner JE, Paust HJ, Lindner M, Peters A, Heiss K, Velden J, Hopfer H, Fehr S, Krieger T, Meyer-Schwesinger C, Meyer TN, Helmchen U, Mittrucker HW, Stahl RA, Panzer U. CXCR3 mediates renal Th1 and Th17 immune response in murine lupus nephritis. J Immunol. 2009;183(7):4693–704.
Yu F, Haas M, Glassock R, Zhao MH. Redefining lupus nephritis: clinical implications of pathophysiologic subtypes. Nat Rev Nephrol. 2017;13(8):483–95.
Dall'Era M. Treatment of lupus nephritis: current paradigms and emerging strategies. Curr Opin Rheumatol. 2017;29(3):241–7.
Forero DA, Guio-Vega GP, Gonzalez-Giraldo Y. A comprehensive regional analysis of genome-wide expression profiles for major depressive disorder. J Affect Disord. 2017;218:86–92.
Khan MI, Debski KJ, Dabrowski M, Czarnecka AM, Szczylik C. Gene set enrichment analysis and ingenuity pathway analysis of metastatic clear cell renal cell carcinoma cell line. Am J Physiol Renal Physiol. 2016;311(2):F424–36.
Bing PF, Xia W, Wang L, Zhang YH, Lei SF, Deng FY. Common marker genes identified from various sample types for systemic lupus erythematosus. PLoS One. 2016;11(6):e0156234.
Sezin T, Vorobyev A, Sadik CD, Zillikens D, Gupta Y, Ludwig RJ. Gene expression analysis reveals novel shared gene signatures and candidate molecular mechanisms between pemphigus and systemic lupus erythematosus in CD4(+) T cells. Front Immunol. 2017;8:1992.
Coit P, Jeffries M, Altorok N, Dozmorov MG, Koelsch KA, Wren JD, Merrill JT, McCune WJ, Sawalha AH. Genome-wide DNA methylation study suggests epigenetic accessibility and transcriptional poising of interferon-regulated genes in naive CD4+ T cells from lupus patients. J Autoimmun. 2013;43:78–84.
Lim EL, Okkenhaug K. PI3Kdelta is a Treg target in cancer immunotherapy. Immunology. 2019.
Ali HR, Provenzano E, Dawson SJ, Blows FM, Liu B, Shah M, Earl HM, Poole CJ, Hiller L, Dunn JA, Bowden SJ, Twelves C, Bartlett JM, Mahmoud SM, Rakha E, Ellis IO, Liu S, Gao D, Nielsen TO, Pharoah PD, Caldas C. Association between CD8+ T-cell infiltration and breast cancer survival in 12,439 patients. Ann Oncol. 2014;25(8):1536–43.
Chen B, Khodadoust MS, Liu CL, Newman AM, Alizadeh AA. Profiling tumor infiltrating immune cells with CIBERSORT. Methods Mol Biol. 2018;1711:243–59.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, Hoang CD, Diehn M, Alizadeh AA. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.
Ali HR, Chlon L, Pharoah PD, Markowetz F, Caldas C. Patterns of immune infiltration in breast Cancer and their clinical implications: a gene-expression-based retrospective study. PLoS Med. 2016;13(12):e1002194.
Rohr-Udilova N, Klinglmuller F, Schulte-Hermann R, Stift J, Herac M, Salzmann M, Finotello F, Timelthaler G, Oberhuber G, Pinter M, Reiberger T, Jensen-Jarolim E, Eferl R, Trauner M. Deviations of the immune cell landscape between healthy liver and hepatocellular carcinoma. Sci Rep. 2018;8(1):6220.
Xiong Y, Wang K, Zhou H, Peng L, You W, Fu Z. Profiles of immune infiltration in colorectal cancer and their clinical significant: a gene expression-based study. Cancer Med. 2018;7(9):4496–508.
Moulton VR, Tsokos GC. T cell signaling abnormalities contribute to aberrant immune cell function and autoimmunity. J Clin Invest. 2015;125(6):2220–7.
Khan SQ, Khan I, Gupta V. CD11b activity modulates pathogenesis of lupus nephritis. Front Med (Lausanne). 2018;5:52.
Bomback AS, Appel GB. Updates on the treatment of lupus nephritis. J Am Soc Nephrol. 2010;21(12):2028–35.
Hirose S, Lin Q, Ohtsuji M, Nishimura H, Verbeek JS. Monocyte subsets involved in the development of systemic lupus erythematosus and rheumatoid arthritis. Int Immunol. 2019.
Santiago-Raber ML, Amano H, Amano E, Baudino L, Otani M, Lin Q, Nimmerjahn F, Verbeek JS, Ravetch JV, Takasaki Y, Hirose S, Izui S. Fcgamma receptor-dependent expansion of a hyperactive monocyte subset in lupus-prone mice. Arthritis Rheum. 2009;60(8):2408–17.
Park YW, Kee SJ, Cho YN, Lee EH, Lee HY, Kim EM, Shin MH, Park JJ, Kim TJ, Lee SS, Yoo DH, Kang HS. Impaired differentiation and cytotoxicity of natural killer cells in systemic lupus erythematosus. Arthritis Rheum. 2009;60(6):1753–63.
Spada R, Rojas JM, Barber DF. Recent findings on the role of natural killer cells in the pathogenesis of systemic lupus erythematosus. J Leukoc Biol. 2015;98(4):479–87.
Schleinitz N, Vely F, Harle JR, Vivier E. Natural killer cells in human autoimmune diseases. Immunology. 2010;131(4):451–8.
Blanco P, Ueno H, Schmitt N. T follicular helper (Tfh) cells in lupus: activation and involvement in SLE pathogenesis. Eur J Immunol. 2016;46(2):281–90.
Iwata YTS. B-cell subsets, signaling and their roles in secretion of autoantibodies. Lupus. 2016;25:850–6.
Kim SJ, Lee K, Diamond B. Follicular helper T cells in systemic lupus erythematosus. Front Immunol. 2018;9:1793.
Zhao YL, Tian PX, Han F, Zheng J, Xia XX, Xue WJ, Ding XM, Ding CG. Comparison of the characteristics of macrophages derived from murine spleen, peritoneal cavity, and bone marrow. J Zhejiang Univ Sci B. 2017;18(12):1055–63.
Z.M. Cristina pamfil, Aurélie de Groof, Gaëlle tilman, Sepideh Babaei, Christine Galant, pauline Montigny, nathalie demoulin, Michel Jadoul,Selda Aydin, ralf Lesche, Fiona Mcdonald, Frédéric A Houssiau, Bernard r Lauwerys. Intrarenal activation of adaptive immune effectors is associated with tubular damage and impaired renal function in lupus nephritis. Ann Rheum Dis. 2018;77:1782–9.
Herrada AA, Escobedo N, Iruretagoyena M, Valenzuela RA, Burgos PI, Cuitino L, Llanos C. Innate immune Cells' contribution to systemic lupus erythematosus. Front Immunol. 2019;10:772.
Muso E, Okuzaki D, Kobayashi S, Iwasaki Y, Sakurai MA, Ito A, Nojima H. Ficolin-1 is up-regulated in leukocytes and glomeruli from microscopic polyangiitis patients. Autoimmunity. 2013;46(8):513–24.
Feng M, Chen JY, Weissman-Tsukamoto R, Volkmer JP, Ho PY, McKenna KM, Cheshier S, Zhang M, Guo N, Gip P, Mitra SS, Weissman IL. Macrophages eat cancer cells using their own calreticulin as a guide: roles of TLR and Btk. Proc Natl Acad Sci U S A. 2015;112(7):2145–50.
Bender AT, Pereira A, Fu K, Samy E, Wu Y, Liu-Bujalski L, Caldwell R, Chen YY, Tian H, Morandi F, Head J, Koehler U, Genest M, Okitsu SL, Xu D, Grenningloh R. Btk inhibition treats TLR7/IFN driven murine lupus. Clin Immunol. 2016;164:65–77.
Chalmers SA, Glynn E, Garcia SJ, Panzenbeck M, Pelletier J, Dimock J, Seccareccia E, Bosanac T, Khalil S, Harcken C, Webb D, Nabozny G, Fine JS, Souza D, Klein E, Herlitz L, Ramanujam M, Putterman C. BTK inhibition ameliorates kidney disease in spontaneous lupus nephritis. Clin Immunol. 2018;197:205–18.
Liu M, Liu J, Hao S, Wu P, Zhang X, Xiao Y, Jiang G, Huang X. Higher activation of the interferon-gamma signaling pathway in systemic lupus erythematosus patients with a high type I IFN score: relation to disease activity. Clin Rheumatol. 2018;37(10):2675–84.
Reiss AB, Wan DW, Anwar K, Merrill JT, Wirkowski PA, Shah N, Cronstein BN, Chan ES, Carsons SE. Enhanced CD36 scavenger receptor expression in THP-1 human monocytes in the presence of lupus plasma: linking autoimmunity and atherosclerosis. Exp Biol Med (Maywood). 2009;234(3):354–60.
Sweet RA, Nickerson KM, Cullen JL, Wang Y, Shlomchik MJ. B cell-extrinsic Myd88 and Fcer1g negatively regulate autoreactive and Normal B cell immune responses. J Immunol. 2017;199(3):885–93.
Salazar-Aldrete C, Galan-Diez M, Fernandez-Ruiz E, Nino-Moreno P, Estrada-Capetillo L, Abud-Mendoza C, Layseca-Espinosa E, Baranda L, Gonzalez-Amaro R. Expression and function of dectin-1 is defective in monocytes from patients with systemic lupus erythematosus and rheumatoid arthritis. J Clin Immunol. 2013;33(2):368–77.
Zimmer J, Nakou M, Bertsias G, Stagakis I, Centola M, Tassiulas I, Hatziapostolou M, Kritikos I, Goulielmos G, Boumpas DT, Iliopoulos D. Gene network analysis of bone marrow mononuclear cells reveals activation of multiple kinase pathways in human systemic lupus erythematosus. PLoS One. 2010;5(10):e13351.
Jensen MA, Patterson KC, Kumar AA, Kumabe M, Franek BS, Niewold TB. Functional genetic polymorphisms in ILT3 are associated with decreased surface expression on dendritic cells and increased serum cytokines in lupus patients. Ann Rheum Dis. 2013;72(4):596–601.
Liu Z, Yu Y, Yue Y, Hearth-Holmes M, Lopez PD, Tineo C, Paulino G, Fu WN, Loyo E, Su K. Genetic alleles associated with SLE susceptibility and clinical manifestations in Hispanic patients from the Dominican Republic. Curr Mol Med. 2019;19(3):164–71.
Nakamura K, Jinnin M, Kudo H, Inoue K, Nakayama W, Honda N, Kajihara I, Masuguchi S, Fukushima S, Ihn H. The role of PSMB9 upregulated by interferon signature in the pathophysiology of cutaneous lesions of dermatomyositis and systemic lupus erythematosus. Br J Dermatol. 2016;174(5):1030–41.
Kong W, Deng W, Sun Y, Huang S, Zhang Z, Shi B, Chen W, Tang X, Yao G, Feng X, Sun L. Increased expression of Bruton's tyrosine kinase in peripheral blood is associated with lupus nephritis. Clin Rheumatol. 2018;37(1):43–9.
Min Sun Shin SK, Kang Y, Wahl ER, Bucala R, Park H-J, Kang I. Macrophage migration inhibitory factor regulates U1 small nuclear RNP immune complex–mediated activation of the NLRP3 Inflammasome. Arthritis Rheumatol. 2019;71:109–20.
Cohen D, Buurma A, Goemaere NN, Girardi G, le Cessie S, Scherjon S, Bloemenkamp KW, de Heer E, Bruijn JA, Bajema IM. Classical complement activation as a footprint for murine and human antiphospholipid antibody-induced fetal loss. J Pathol. 2011;225(4):502–11.
Chougule D, Nadkar M, Venkataraman K, Rajadhyaksha A, Hase N, Jamale T, Kini S, Khadilkar P, Anand V, Madkaikar M, Pradhan V. Adipokine interactions promote the pathogenesis of systemic lupus erythematosus. Cytokine. 2018;111:20–7.
Chen XW, Shen Y, Sun CY, Wu FX, Chen Y, Yang CD. Anti-class a scavenger receptor autoantibodies from systemic lupus erythematosus patients impair phagocytic clearance of apoptotic cells by macrophages in vitro. Arthritis Res Ther. 2011;13(1):R9.
Linden M, Ramirez Sepulveda JI, James T, Thorlacius GE, Brauner S, Gomez-Cabrero D, Olsson T, Kockum I, Wahren-Herlenius M. Sex influences eQTL effects of SLE and Sjogren's syndrome-associated genetic polymorphisms. Biol Sex Differ. 2017;8(1):34.
Renauer P, Coit P, Jeffries MA, Merrill JT, McCune WJ, Maksimowicz-McKinnon K, Sawalha AH. DNA methylation patterns in naive CD4+ T cells identify epigenetic susceptibility loci for malar rash and discoid rash in systemic lupus erythematosus. Lupus Sci Med. 2015;2(1):e000101.
Moorthy K, Jaber AN, Ismail MA, Ernawan F, Mohamad MS, Deris S. Missing-values imputation algorithms for microarray gene expression data. Methods Mol Biol. 2019;1986:255–66.
Diboun I, Wernisch L, Orengo CA, Koltzenburg M. Microarray analysis after RNA amplification can detect pronounced differences in gene expression using limma. BMC Genomics. 2006;7:252.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
This work was supported by the NSFC (grant no. 81270805) and Science and Technology Department of Sichuan province grant (No. 2018SZ0378). The funding bodies had no role in the design of the study and collection, analysis, interpretation of data and writing of the manuscript.
Ethics approval and consent to participate
Consent for publication
On behalf of all authors, the corresponding author states that there is no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Cao, Y., Tang, W. & Tang, W. Immune cell infiltration characteristics and related core genes in lupus nephritis: results from bioinformatic analysis. BMC Immunol 20, 37 (2019). https://doi.org/10.1186/s12865-019-0316-x
- Systemic lupus erythematosus
- Lupus nephritis
- Immune infiltration