- 1School of Life Sciences and Technology, Tongji University, Shanghai, China
- 2Biomedical Research Center, Tongji University Suzhou Institute, Suzhou, Jiangsu, China
- 3Department of Neurology, Tongji Hospital, Tongji University, Shanghai, China
Introduction: Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease. Anti-B-cell-activating factor (BAFF) therapy effectively depletes B cells and reduces SLE disease activity. This research aimed to evaluate the effect of BAFF blockade on B cell receptor (BCR) repertoire and gene expression.
Methods: Through next-generation sequencing, we analyzed gene expression and BCR repertoire in MRL/lpr mice that received long-term anti-BAFF therapy. Based on gene expression profiles, we predicted the relative proportion of immune cells using ImmuCellAI-mouse, validating our predictions via flow cytometry and FluoroSpot.
Results: The loss of BCR repertoire diversity and richness, along with increased clonality and differential frequency distribution of the immunoglobulin heavy chain variable (IGHV) segment gene usage, were observed in BAFF-blockade mice. Meanwhile, the distribution of complementarity-determining region 3 (CDR3) length and CDR3 amino acid usage remained unaffected. BAFF blockade resulted in extensive changes in gene expression, particularly that of genes related to B cells and immunoglobulins. Besides, the tumor necrosis factor (TNF)-α responses and interferon (IFN)-α/γ were downregulated, consistent with the decrease in IFN-γ and TNF-α serum levels following anti-BAFF therapy. In addition, BAFF blockade significantly reduced B cell subpopulations and plasmacytoid dendritic cells, and caused the depletion of antibody-secreting cells.
Discussion: Our comparative BCR repertoire and transcriptome analyses of MRL/lpr mice subjected to BAFF blockade provide innovative insights into the molecular pathophysiology of SLE.
1 Introduction
Systemic lupus erythematosus (SLE) is a multisystem autoimmune disease characterized by abnormal immune activity that causes tissue damage and extensive inflammation in multiple organs (1). Disruption of immune tolerance mechanisms leads to the massive proliferation and activation of B and T cells, which in turn produce autoantibodies and proinflammatory cytokines (2). B-cell-activating factor (BAFF, BLyS), which belongs to the tumor necrosis factor (TNF) ligand family, plays a key role in B-cell proliferation, differentiation, and activation (3). Elevated serum levels of BAFF correlate positively with anti-double-stranded DNA antibody (anti-dsDNA) titers and an increased present-day risk of autoimmunity (4, 5). To date, belimumab and telitacicept are the only two biologics against BAFF that have been approved for SLE therapy worldwide (6, 7).
Considering the central role of B cells in SLE pathogenesis, several researchers have explored changes in the B cell receptor (BCR) repertoire. Tipton and Odendah reported that the plasma cells of SLE patients exhibited increased usage of the immunoglobulin heavy chain variable (IGHV) 4 family, especially IGHV4-34, which encodes 9G4+ antibodies that target double-stranded DNA (8, 9). More diverse complementarity-determining region 3 (CDR3) sequences have been observed in the B cell receptor (BCR) repertoire of patients with SLE than in healthy controls (10). Additionally, a significant shortening of CDR3 length in peripheral blood BCR has been reported in the former (11).
Transcriptome analyses of whole blood samples have been used to identify the biomarkers and pathogenic drivers of SLE. The expression of interferon (IFN)-stimulated genes (ISGs), more specifically type I IFN-stimulated genes, has been reported to correlate with SLE disease activity (12, 13). High ISG expression leads to autoantibody production together with the abnormal proliferation of B cells (14). A whole-blood transcriptome study of 158 pediatric patients with SLE found an enrichment of plasmablast and neutrophil transcripts (15). Based on gene expression profiles, Akhtar et al. estimated immune cell proportions using CIBERSORTx and found that patients with SLE had more naïve CD4 T cells and fewer memory-activated CD4 T cells than patients exposed to mycophenolate mofetil (16).
The complex interplay between BAFF, BAFF receptor, and BCR, can reshape the BCR repertoire and result in a skewed usage of the IGHV gene (17–19).. However, changes in the BCR repertoire and transcriptome of SLE murine models treated with biologics against BAFF remain largely unknown. MRL/lpr mice develop an autoimmune disease resembling SLE caused by a mutation in the Fas gene, which promotes the survival of self-reactive lymphocytes (20). In this study, we performed a thorough analysis of the BCR repertoire and transcriptome of MRL/lpr mice treated with the BAFF receptor fusion protein (BAFF-R-Fc). Our results demonstrate that the frequency distribution of the IGHV gene usage was distinct from those of controls. We also observed a decrease in BCR diversity and richness, and an increase in clonality in the treatment group concurrently. Furthermore, blocking BAFF downregulates genes related to B cell proliferation and antibody production, and alters the composition of splenic immune cells, especially by reducing the number of B cell subsets. Overall, our findings advance a better knowledge of the mechanisms of BAFF blockade against SLE.
2 Methods
2.1 Mice and mouse BAFF-R-Fc fusion protein
Female 6-week-old MRL/lpr mice were purchased from Shanghai Laboratory Animal Co., Ltd. (Shanghai, China). All mice were raised under specific pathogen-free conditions. Mouse BAFF-R-Fc fusion protein was provided by RemeGen Co., Ltd. (Yantai, China). All experimental procedures were performed following the regulations of the Animal Ethics Committee of Tongji University (No. TJLAC-018-032).
2.2 Animal experiment
The treatment group received BAFF-R-Fc at a dose of 10 mg/kg, whereas the control group received normal saline at the age of 8 weeks. All drugs were injected subcutaneously once per day for a total of 28 injections.
2.3 BCR repertoire construction, sequencing, and analysis
Total RNA was extracted from the MRL/lpr mice spleens after treatment using Trizol according to the user manual. 5’RACE was performed for reverse-transcription, and 5’RACE template switch oligos contain unique molecular identifiers (HiScript-TS 5’/3’ RACE Kit, Vazyme, China). Reverse transcription products were amplified via PCR using specific primers designed for IgG and IgM sequences (primers were provided by Azenta Life Sciences, Chelmsford, MA, USA). PCR products were purified by DNA magnetic beads. PCR purified product is terminal repair (including 5 ‘terminal phosphorylation and 3’ terminal plus’ A ‘) by the End Prep Enzyme Mix, with sequencing adapters at both ends. Then PCR products were amplified with P5 and P7 primers, cleaned up using beads, validated using a Qsep100 (Bioptic, Taiwan, China), and quantified by Qubit3.0 Fluorometer (Invitrogen, Carlsbad, CA, USA). Sequencing was performed with an Illumina MiSeq instrument and carried out using a 2×300 paired-end configuration, according to the manufacturer’s instructions (Illumina, San Diego, CA, USA).
Raw fastq files were first subject to quality assessment. Adapters and bases with poor quality scores (Q value lower than 20) were removed using Cutadapt (V1.9.1) to generate clean data. Paired-end reads were merged using FLASH (version 2.2.00). The merged sequences were BLASTed against the IMGT reference database to identify the best match between germline V(D)J genes and sequences in the CDR1, CDR2, and CDR3 regions using MIXCR (V3.0.13). To assess BCR repertoire diversity, IGHV, and IGHJ gene usage, distributions of CDR3 region length, and clone information, clonotype sequences were analyzed using the immunarch R package (21).
2.4 RNA-seq and data analysis
RNA-seq was performed using 1 μg total RNA (Azenta Life Sciences, Chelmsford, MA, USA). The poly(A) mRNA isolation was performed using Oligo(dT) beads. The mRNA fragmentation was performed using divalent cations and high temperature. Priming was performed using Random Primers. First strand cDNA and the second-strand cDNA were synthesized. The purified double-stranded cDNA was then treated to repair both ends and add a dA-tailing in one reaction, followed by a T-A ligation to add adaptors to both ends. Size selection of Adaptor-ligated DNA was then performed using DNA Clean Beads. Each sample was then amplified by PCR using P5 and P7 primers and the PCR products were validated. Sequencing was performed on an Illumina HiSeq instrument and carried out using a 2×150 paired-end configuration, according to the manufacturer’s instructions (Illumina, San Diego, CA, USA).
The DESeq2 Bioconductor package was performed to identify differentially expressed genes (DEGs) (22, 23), which were defined as those with an adjusted p-value < 0.05 and |log2 fold change| ≥ 1. ClusterProfiler and pathview R packages were used to identify Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway terms and Gene Ontology (GO) biological process, with a significant adjusted p-value < 0.05 as the threshold (24). Gene set enrichment analysis (GSEA) was used to identify differentially enriched hallmark gene sets downloaded from the MSigDB database using the ClusterProfiler R package. p-value < 0.05 and |Normalized enrichment score| > 1 were set as the cutoffs.
Immune-related DEGs were submitted to the STRING database to obtain interaction information, with the cutoff interaction score set at 0.7. The visualization of the protein-protein interaction (PPI) network was performed using Cytoscape software. The betweenness centrality for each gene was estimated using the cytoNCA plugin of Cytoscape. DEGs were arranged according to the value of betweenness centrality. The ImmuCellAI-mouse online tool was used to predict the relative proportion of immune cells based on RNA-seq data (25). The TRRUST v2 database, which contains 6552 transcription factors (TFs)-target interactions between mouse TFs, was used to identify key TFs among DEGs (26).
2.5 Measurement of immunoglobulin isotype, anti-dsDNA antibody titer, and immunoregulatory cytokines levels in serum
Anti-dsDNA antibody titers were measured using an ELISA kit (Alpha Diagnostic Intl. Inc. Texas, USA). Immunoglobulin (IgA, IgM, IgG1, IgG2a, IgG2b, IgG3) and cytokine (IL-10, IFN-γ, TNF-α) levels in serum were measured with a LEGENDplex™ CBA kit (BioLegend, San Diego, CA, USA), using a CytoFLEX LX (Beckman Coulter, Brea, CA, USA), according to the manufacturer’s instructions.
2.6 Analysis of B cell subsets in the spleen via flow cytometry
MRL/lpr mice were sacrificed on the day after the 28th treatment to obtain the spleens. Forty-µm mesh filters (BioFIL, Shanghai, China) were used to obtain single-cell suspensions in phosphate-buffered saline (BioFIL, Shanghai, China). The sorting strategy for splenic lymphocytes was based on OMIP 076 (27). Subsets of B cell lymphocytes were analyzed using CytoFLEX LX (Beckman Coulter, Brea, CA, USA) after incubation with the following antibody combinations: PerCP/Cyanine5.5 anti-mouse CD45R/B220, TruStain FcX™ (anti-mouse CD16/32) Antibody, APC/Cyanine7 anti-mouse CD19, FITC anti-mouse TCRβ chain, PE/Cyanine7 anti-mouse CD138 (Syndecan-1), APC anti-mouse IgM, Brilliant Violet 785™ anti-mouse I-A/I-E, APC anti-mouse CD1d (CD1.1, Ly-38), Zombie Aqua™, PE anti-mouse CD5, PE anti-mouse PD-L2. All antibodies were purchased from BioLegend (San Diego, CA, USA).
2.7 Detection of antibody-secreting cells in the spleen via FluoroSpot assay
The number of IgM, IgA, and IgG ASCs was determined using FluoroSpot flex (MABTECH AB, Sweden) following the manufacturer’s instructions. In brief, 5×104 cells were added per well in 200 µl cell culture medium containing 10% fetal bovine serum. The cells were then incubated overnight in a 37°C humidified incubator with 5% CO2. Cells were then removed and detection antibodies were added to per well and incubated for 2 hours at 20°C. The spots were counted using a FluoroSpot reader (MABTECH AB, Sweden).
2.8 Statistical analysis
The statistical significance of difference between the two groups was analyzed using the unpaired two-tailed t-test or Mann-Whitney U-test. P-value adjusting is done using the Holm-Bonferroni correction in multiple testing. Pearson correlation coefficient was performed to indicate the extent of linear correlation between two arbitrary variables. Statistical analysis was performed using GraphPad Prism version 9.5 and R software version 4.3.0. Data were presented as the mean ± standard deviation. A p-value < 0.05 was regarded as statistically significant. *p < 0.05; **p < 0.01; ***p < 0.001.
3 Results
3.1 Changes in BCR repertoire after BAFF blockade
To analyze BCR repertoire clonotypes in MRL/lpr mice, we performed next-generation sequencing of the BCR transcripts expressed by splenic B cells. First, we determined the ratio of unique/total BCR sequences. Notably, the ratio of unique/total BCR clonotypes in the treatment group was significantly lower than in the control group (Figure 1A). Next, we used ecological parameters to quantitatively measure repertoire diversity and complexity (28). The Chao1 index is a non-parametric asymptotic estimator of clonotype richness. Compared with the control, the treated MRL/lpr mice had a lower Chao1 index for the BCR repertoire (Figure 1B). To measure repertoire diversity, we used the Shannon index. A higher Shannon index indicates greater diversity of the BCR repertoire. A significant reduction in the Shannon index was noted compared to the control (Figure 1C). We also evaluated the inverse Simpson’s index for a more precise assessment. Consistent with the Shannon index results, the inverse Simpson’s index was reduced in the treatment group, although no statistical difference was observed (Figure 1D). Collectively, after BAFF-R-Fc treatment, the BCR repertoire of MRL/lpr mice exhibited a lower ratio of unique/total clonotypes, together with an obvious reduction in richness and diversity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1307392/fimmu-14-1307392-HTML-r1/image_m/fimmu-14-1307392-g001.jpg)
Figure 1 Analysis of BCR repertoire in the spleen of MRL/lpr mice. (A) The number of unique clonotypes. (B) Clonotypes richness estimated by Chao1 index. (C) Quantification of the BCR repertoire diversity using the Shannon index and (D) the inverse-Simpson index. (E) The cumulative relative abundance of small-expanded (0-0.01%), medium-expanded (0.01%-0.1%), large-expanded (0.1%-0.5%), and hyper-expanded (0.5%-100%) clonotypes in two groups (mean ± S.D, Mann-Whitney U-test and Holm-Bonferroni correction). (F) The bar plot of the D50 diversity index, D50 represents the number of clonotypes occupying 50% of repertoires. (G) Visualization of overlap coefficient among samples to measure repertoire similarity. (H) Multi-dimensional scaling (MDS) is applied to overlap coefficient, then clusters the MDS resulting components using K-means. Statistical analysis was performed with unpaired t-test comparing untreated to treated. Bars represent mean ± S.D. *p < 0.05; **p < 0.01; ***p < 0.001: ns, not significant.
According to the clonotype proportions in total sequences (Supplementary Figure 1A), we classified these into four groups to assess the clonal space homeostasis: small-expanded (0-0.01%), medium-expanded (0.01%-0.1%), and large-expanded (0.1%-0.5%) and hyper-expanded (0.5%-100%). The cumulative relative abundances of these clonotypes were estimated (Figure 1E; Supplementary Figure 1B), as shown, the proportion of hyper-expanded clonotypes in the treatment group was significantly higher than that in the control, whereas the proportions of small- and medium-expanded clonotypes were significantly reduced. We also considered the most abundant cell clonotypes. Top 11-100 (p = 0.011) clonotypes occupied more repertoire space in the treatment group than those in the control, while 1001-3000 (p = 0.011), and 3001-10000 (p = 0.011) clonotypes occupied less repertoire space (Supplementary Figures 1C, D). Additionally, we estimated the relative proportion of the least prolific clonotypes. The proportions of clonotypes count 1, 2-3, 4-10, 11-30, and 31-100 were significantly lower in the treatment group than those in control, while clonotypes count more than 100 occupied more proportion (Supplementary Figures 1E, F). The diversity 50 index represents the number of clonotypes that occupy 50% of all clonotypes observed in the BCR repertoire. We found that the diversity 50 index in the treatment group was markedly less than that in the control group (Figure 1F). In conclusion, these outcomes indicate that fewer clonotypes tended to expand massively in treated mice.
To measure repertoire similarity, we took advantage of repertoire overlap by computing the overlap coefficient, which was defined as the ratio of the intersection’s size to the smaller size of the two sets. (Figure 1G). With this result, multidimensional scaling (MDS) clustered using K-means was applied to visualize repertoire similarity. MDS partially segregated the treated mice from the control group, as four of six samples from the treated group were highly clustered (Figure 1H). This distinction existed in the repertoire of drug-treated and control MRL/lpr mice.
3.2 Comparison of IGHV and immunoglobulin heavy chain joining gene usage
To determine whether BAFF blockade alters the usage of IGHV and IGHJ genes, we calculated the proportion of V and J alleles. The IGHV genes with frequency > 0.1% are shown in Figure 2A. IGHV1-14, IGHV1S61, and IGHV1-82 showed extremely high abundance (> 5%) in both groups (Supplementary Table 1). We compared the proportion of IGHV allele usage between the two groups, and found no significant difference in VH gene usage between them (Supplementary Table 1). We then analyzed the frequency of IGHV gene families in two groups. IGHV1 family was the most commonly used gene family in both groups (Figure 2B). And the frequency of IGHV1 family in BAFF-R-fc group did not reach statistical significance (64% ± 8%) compared to the control (57% ± 5%). The frequencies of IGHJ alleles are shown in Figure 2C. No statistical difference was observed in the usage of IGHJ alleles, and IGHJ4 was the most commonly used gene in both groups (Figure 2D). To fully understand the differences in gene usage, we estimated the Jensen-Shannon divergence between the two samples (Figure 2E), a method of measuring similarity between two probability distributions. A higher Jensen-Shannon divergence indicated a larger difference in the distribution of IGHV gene usage. Based on this result, PCA clustered using K-means was applied to map IGHV gene usage. As shown in Figure 2F, the PCA significantly segregated four samples in the treatment group from the control group, whereas two samples could not be distinguished. Together, these results suggest that the usage of IGHV alleles in treated mice differed from that in control mice.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1307392/fimmu-14-1307392-HTML-r1/image_m/fimmu-14-1307392-g002.jpg)
Figure 2 Differential usage of IGHV/J genes in BCR repertoires of MRL/lpr mice. (A) Heat map of IGHV gene usage showing those abundances > 0.1% among IGH sequences from each sample. The color indicates clone abundance per sample, as indicated in the legend. (B) The frequency of IGHV gene families in the BAFF-R-Fc treated and control group. (C) Heat map of IGHJ gene usage of each sample. (D) IGHJ gene usage is shown in the boxplot. (E) Visualization of Jensen-Shannon divergence between samples for IGHV genes usage analysis. (F) Principal components analyses (PCA) of Jensen-Shannon divergence, then clusters the PCA resulting components using K-means. Statistical analysis was performed with Mann-Whitney U-test and Holm-Bonferroni correction comparing untreated to treated. Bars represent mean ± S.D.
3.3 Characteristics of CDR3 sequences
CDR3 region of the immunoglobulin heavy chain is a key region for recognizing and binding antigens (29). Therefore, we analyzed the distribution of CDR3 length. As shown in Figure 3A, the most common CDR3 length in the control mice was 15 amino acids. In BAFF-R-Fc-treated mice, the most common CDR3 length decreased to 14 amino acids. However, no statistical differences in the usage of amino acids with different lengths were observed. We investigated the mean length of CDR3 sequences. The CDR3 mean length observed in BAFF-R-Fc-treated mice (13.71± 0.27) did not reach statistical significance compared to the control (13.93 ± 0.52) (Figure 3B). We also estimated the 20-amino-acid composition of CDR3 (Figure 3C). We then determined the cumulative relative abundances of hydrophilic, neutral, and hydrophobic amino acids. No statistically significant differences were observed in the usage of these three amino acid types (Figure 3D). We identified the public CDR3 sequences of BCR heavy chains shared between two samples and exhibited in Figure 3E, these shared CDR3 sequences had different unique molecular identifiers. For all samples, only three CDR3 sequences, CAKTTRATSAMDYW, CARRVPVYFDYW and CASHRELWFAYW, were used in all samples (Supplementary Table 2). And we tracked seven CDR3 sequences shared between most samples (Figure 3F). Noticeably, a few CDR3 sequences exhibited extremely high proportions in a specific sample.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1307392/fimmu-14-1307392-HTML-r1/image_m/fimmu-14-1307392-g003.jpg)
Figure 3 Characteristics of the CDR3 sequences in BCR repertoires of MRL/lpr mice. (A) Distribution of the length of CDR3 region in treated and control groups. (B) CDR3 region mean length in two groups (mean ± S.D., unpaired t-test). (C) Stacked plot displaying CDR3 region amino acid composition in each sample. (D) The usage of hydrophilic, neutral, and hydrophobic amino acids in two groups. (E) Heatmap of public CDR3 sequences shared between samples. (F) Tracking the seven CDR3 amino acid sequences shared between most samples. Statistical analysis was performed with Mann-Whitney U-test and Holm-Bonferroni correction comparing untreated to treated. Bars represent mean ± S.D.
3.4 Distinct transcriptional signatures linked to BAFF-R-Fc treatment responses
We first performed PCA on the gene expression data, as shown in Figure 4A. PCA segregated treated from control mice, suggesting significant differences in gene expression profiles between the two groups. We identified 1059 DEGs in total (Supplementary Table 3), among which 757 were significantly downregulated, including IGHD, EBF1, and CR2, while 302 genes were upregulated (Figure 4B) in the treated group. The expression levels of DEGs are shown in Figure 4C.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1307392/fimmu-14-1307392-HTML-r1/image_m/fimmu-14-1307392-g004.jpg)
Figure 4 Distinct transcriptional signatures in the spleen of MRL/lpr mice received BAFF-R-Fc. (A) Principal components analyses (PCA) of spleen gene expression profiles from BAFF-R-Fc treated group and control. The first two principal components PCA1 and PCA2 are evaluated for each sample vector and plotted. (B) Volcano plot and (C) heatmap of differential expression genes between treatment and control group. Red points denote upregulated genes, blue points denote downregulated genes, and gray points denote genes with no significant differences. (D) Venn diagram of the overlap among immune-related genes (IRGs) and upregulated/downregulated DEGs. IRGs are obtained from the Immport database. The protein-protein interaction (PPI) network of 87 immune-related DEGs is constructed with STRING database (E) and Cytoscape software (F). Color from red to yellow and circular sizes from large to small indicate the betweenness centrality of genes from high to low. Betweenness centrality is estimated by the cytoNCA plugin.
To fully understand the biological mechanism of BAFF blockade in SLE, we used multiple external databases to interpret DEGs. We identified key genes encoding cytokines that play important roles in autoimmune disease (30). CCL12, CXCL11, CXCL3, IL-10, IL-1A, IL-21, TNFSF15, and TNFSF18 were downregulated after treatment (Supplementary Figures 2A, B). Among these, IL-21 and CXCL13 drive B cell activation, differentiation, and antibody production (31, 32). The expression of CXCL5 (Supplementary Figures 2A, B), which suppresses neutrophil and myeloid activation, degranulation, and proliferation, was significantly higher in the treatment group (33). Next, we analyzed ISGs that correlate with SLE disease activity (34). Changes in the expression levels of ISGs were observed in the two groups (Supplementary Figure 2C). However, no ISGs meet the criteria for the DEGs, even though the treatment group showed an improvement in disease activity.
Next, we focused on immune-related genes (IRGs). Among these, 70 DEGs were downregulated, and 17 were upregulated after BAFF-R-Fc treatment (Figure 4D; Supplementary Table 4). To determine the interactions between these immune-related DEGs, a PPI network was constructed using the STRING database (Figure 4E). We then evaluated the betweenness centrality of each gene using Cytoscape and arranged the genes based on betweenness centrality (Figure 4F). IL-10, CD19, and CD40 were the top three highest-ranking IRGs. The downregulation of IL-10 expression level is consistent with the serum level we detected after treatment (Supplementary Figure 3B). Immune complexes can induce IL-10 production from peripheral blood mononuclear cells through the Fcγ receptor, further enhancing autoantibody production and maintaining the hyperactivity of B cells in SLE (35). Further, interactions between the CD40 ligand on T cells and CD40 on B cells lead to the over-activation of T, and B cells in SLE (36).
Finally, we used the TRRUST v2 database to identify key TFs. Among the downregulated DEGs, we observed 58 TFs, 33 of which with a p-value < 0.05, suggesting that these might play important roles during treatment (Supplementary Table 5). The functions of these TFs, such as PAX5, XBP1, CEBPB, NF-κb1, EBF1, MEF2C, and POU2AF1 were strongly associated with B cell activation and proliferation (Supplementary Figure 2D) (37–41), further providing more evidence that BAFF blockade suppresses B cells proliferation in SLE.
3.5 Functional and gene set enrichment analysis
We performed KEGG pathway and GO biological process enrichment analyses to reveal the underlying mechanism of BAFF blockade in MRL/lpr mice based on DEGs. KEGG pathway enrichment analysis demonstrated that DEGs were strongly enriched in the BCR signaling pathway, antibody production, cytokine-cytokine receptor interaction, and several autoimmune diseases, as shown in Figure 5A. These pathways are closely linked (Figure 5B). Enriched GO biological process terms included BCR signaling pathway, membrane invagination, B cell activation, and immune recognition (Figure 5C). Considering that some genes with important biological functions were not significantly differentially expressed, we performed GSEA to confirm the enrichment of hallmark pathways using whole gene expression profiles. Noticeably, the IFN-α/γ and TNFα responses were significantly downregulated, consistent with serum levels of IFN-γ and TNFαwe detected in the two groups (Supplementary Figures 3C, D). In addition, E2F targets, and G2M checkpoint gene sets, two hallmark gene sets related to cell cycle progression and cell proliferation, were activated in the treatment compared to the control group (Figure 5D; Supplementary Table 6). These results identify the inhibited signal pathways in MRL/lpr mice that received anti-BAFF therapy.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1307392/fimmu-14-1307392-HTML-r1/image_m/fimmu-14-1307392-g005.jpg)
Figure 5 Functional enrichment analysis of DEGs and GSEA. (A) Enrichment analysis of DEGs using the KEGG pathway database. The y-axis represents pathway names, the x-axis represents enrichment score in -log10(p value). The number of genes enriched on pathways is represented by the plot size, larger sizes indicate more significant DEGs enrichment on pathways. The color gradient ranges from blue to red, reflecting decreasing p values from high to low. (B) Pathway-pathway interaction network of the ten KEGG pathways with the smallest p values is shown. (C) GO biological process enrichment analysis of DEGs. (D) Activated or inhibited pathways were identified using GESA based on hallmark gene sets downloaded from the MSigDB database. The enrichment score line graph highest peak at the bottom an enrichment score < 0 represents the pathway is inhibited, the highest peak at the top and an enrichment score > 0 represents the pathway is activated. Five pathways with the smallest p values were represented by different colors.
3.6 Immune cell composition of spleen analysis
Using gene expression profiles, we computed the relative proportion of immune cells in MRL/lpr mice spleens (Figure 6A). Distinct changes were observed after BAFF-R-Fc treatment, as shown in Figure 6B. Significant decreases in the proportion of B cell subsets, including memory B, B1, plasma, and marginal zone B cells, as well as plasmacytoid dendritic cells (pDCs) were observed. Meanwhile, the proportion of monocytes, basophils, eosinophils, mast cells, neutrophils, M2 macrophages, naïve CD4 T cells, T helper cells, and exhausted CD8 T cells significantly increased. Principal component analysis (PCA) revealed significant group-bias clustering in the proportion of immune cells in samples from treated mice and controls (Figure 6C). The decrease in B cell subpopulations was consistent with the downregulation of TFs related to B cell proliferation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1307392/fimmu-14-1307392-HTML-r1/image_m/fimmu-14-1307392-g006.jpg)
Figure 6 Using gene expression profiles to predict the relative proportion of immune cells. (A) The relative abundance of 29 immune cells in the spleen of each sample was estimated by ImmuCellAI-mouse. (B) Comparison of each immune cell type in control and BAFF blockade groups. (C) Principal components analyses (PCA) of immune cells’ proportions. The first two principal components PCA1 and PCA2 are evaluated for each sample vector and plotted. cDC, conventional dendritic cell. CD8 Tc, CD8 cytotoxic T cell. pDC, plasmacytoid dendritic cell. NKT, natural killer T cell. CD8 Tc, CD8 cytotoxic T cell. T gamma delta, γδ T cell. MoDC, monocyte-derived dendritic cell. NK, natural killer cell. Treg, regulatory T cell. CD8 Tex, exhausted CD8 T cell. CD4 Tm, memory CD4 T cell. CD8 Tcm, CD8 central memory T cell. CD8 Tem, CD8 effector memory T cell. Statistical analysis was performed with Mann-Whitney U-test. Bars represent mean ± S.D. *p < 0.05; **p < 0.01.
3.7 The correlations between immunoglobulins, B cell proportion, BCR repertoire diversity, and richness
To investigate the correlation between BCR diversity, richness, B cell proportion in the spleen, and serum immunoglobulins levels, we conducted a correlation analysis of these variables in the two groups. First, we measured the serum levels of immunoglobulins. All immunoglobulin isotypes, including IgA, IgM, IgG1, IgG2a, IgG2b, and IgG3, exhibited lower serum levels following treatment (Figure 7A; Supplementary Figure 3A). We also measured anti-dsDNA antibody titers in MRL/lpr mice, which are highly specific to SLE and positively correlated with disease activity. As expected, BAFF-R-Fc decreased anti-dsDNA antibody titers compared to that in the control (Figure 7B). Correlation analysis revealed significant positive correlations between immunoglobulin levels and BCR diversity and richness (Figures 7C, D). Besides, B cell proportion was also positively correlated with BCR diversity and richness (Figures 7E, F). These data suggest that anti-BAFF therapy in MRL/lpr mice could impact the BCR diversity and richness by depleting B cells and preventing antibody production, which is relevant to prognostic outcomes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1307392/fimmu-14-1307392-HTML-r1/image_m/fimmu-14-1307392-g007.jpg)
Figure 7 Correlation analysis of the BCR diversity, richness, immunoglobulins levels, and B cells proportion. (A) Serum levels of total immunoglobulins in the two groups. (B) Anti-dsDNA antibodies titers of BAFF blockade and control mice. Titers were shown in activity units. Correlation between total immunoglobulins levels and the Shannon index (C) and CHAO-1 (D). Correlation between B cells proportion and the Shannon index (E) and CHAO-1 (F). Pearson correlation coefficient was performed to indicate the extent of linear correlation between two arbitrary variables. Confidence intervals and line of regression were shown. Statistical analysis was performed with unpaired t-test comparing untreated to treated. Bars represent mean ± S.D. **p < 0.01; ***p < 0.001.
3.8 BAFF-R-Fc regulated splenic B cell subset development and reduced ASC numbers
We validated B cell depletion caused by BAFF blockade by flow cytometry and FluoroSpot. We first determined the proportion of B cells (TCRβ-CD45R+), which was obviously reduced in the spleen of mice treated with BAFF-R-Fc (Figures 8A, B). Consistently, memory B cells (PD-L2+IgM+MHCII+) and B10 cells (CD1d+CD5+), a regulatory B cell subset producing IL10 to control T cell-dependent inflammatory responses (42), were also reduced (Figures 8A, B). Meanwhile, we found that the proportion of T cells (TCRβ+) was higher than that in controls (Figures 8A, B). We confirmed that BAFF-R-Fc suppressed autoantibody production by decreasing the number of ASCs, including plasma cell precursors (plasmablasts, IgM+MHCII+CD138+) and plasma cells (CD19-CD138+) (Figures 8A, B). We then investigated whether BAFF blockade only leads to the depletion of certain ASCs, subjecting them to FluoroSpot analysis. The numbers of IgM-, IgA-, and IgG-secreting cells were consistently reduced (Figures 8C, D). These immune cell proportions obtained by flow cytometry were consistent with the predictions using ImmuCellAI-mouse, with Pearson’s correlation coefficients ranged from 0.66 to 0.89 (Supplementary Figures 3E–H). Our results suggest that BAFF blockade contributes to the depletion of these B cell subsets.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1307392/fimmu-14-1307392-HTML-r1/image_m/fimmu-14-1307392-g008.jpg)
Figure 8 BAFF blockade regulated splenic T cells, B cell subsets, and ASCs. B cell subsets and T cells were obtained from MRL/lpr mice spleens, and detected by flow cytometry. IgM-, IgA-, and IgG-secreting cells were detected using the FluoroSpot assay. (A) Representative FACS analysis of T cells (TCRβ+), B cells (TCRβ-CD45R+), and B cell subsets including plasma cells (CD19-CD138+) plasmablast (IgM+MHCII+CD138+), B10 cells (CD1d+CD5+) and memory B cells (PD-L2+IgM+MHCII+). (B) The box plot shows the percent of B cells, T cells, and plasma cells in splenic lymphocytes and the percent of memory B cells plasmablasts, and B10 cells in B cells. (C) Images of BAFF blockade and control’s IgA/IgM/IgG FluoroSpot with three-color spots: IgA-green spot, IgM-red spots, IgG-yellow spot, and composite picture of three-color spots. (D) Number of IgA, IgM, IgG, and total ASCs spots in 5×104 cells detected by FluoroSpot assay. Statistical analysis was performed with unpaired t-test and Holm-Bonferroni correction comparing untreated to treated. Bars represent mean ± S.D. **p < 0.01; ***p < 0.001.
4 Discussion
In this study, we performed transcriptome and BCR repertoire analyses on the spleens of MRL/lpr mice treated with the BAFF-F-R-Fc fusion protein. We noted lower BCR diversity, richness, and a decrease in the number of unique clonotypes in treated mice. IGHV gene usage was also different from that in controls. Additionally, our study successfully detected changes in the proportions of splenic lymphocyte subpopulations and gene expression. We observed downregulation of IFN-α/γ and TNF-α responses, with most DEGs related to B cell proliferation and antibody production. These results imply that BAFF blockade dramatically remodeled the immune microenvironment in the spleen of MRL/lpr mice and regulated the B cell receptor repertoire.
The comprehensive characterization of the BCR repertoire has improved our knowledge of autoimmune disease pathogenesis (43, 44). The presence of abnormal specific clonotype expansion and diversity as well as the abnormal usage of IGHV genes in SLE indicate a broad dysregulation of the BCR repertoire (45, 46). Thus, we studied the effect of BAFF blockade on the BCR repertoire. Tipton et al. previously noted that the plasma and plasmablast populations of patients with SLE expand considerably as a result of prolonged immune activation, which resulting in an increase in the proportion of hyperexpanded clonotypes and a decrease in BCR repertoire diversity (8). Interestingly, instead of returning to normal levels, a further decrease in BCR repertoire diversity and richness were observed in the BAFF-R-Fc treatment group compare to the control. Besides, hyperexpanded clonotypes occupied more BCR repertoire space in the treatment group. So, we hypothesized that BAFF blockade reduces the number of plasma cells and plasmablasts, while also sharply reducing the number of other B cell subpopulations, which causes dramatic declines in BCR repertoire diversity and richness. As we found the BCR diversity and richness were positively correlated with B cells proportion. Similar outcomes have been reported in patients receiving stable immunosuppressive therapy such as mycophenolate mofetil and leflunomide (47).
Huang et al. reported a notable loss of IGHV4-34 among plasmablasts and a mutated IgM BCR repertoire in chronic belimumab-treated individuals when compared to controls (48). Here, we found no significant difference in the usage of IGHV genes after BAFF-R-Fc therapy in MRL/lpr mice. Notably, some VH genes such as VH1-14, VH1S61, and VH1-82 showed extremely high frequencies in both treated and control mice. Further studies are required to confirm whether these genes are involved in the production of autoreactive antibodies. Nonetheless, we observed a difference in the frequency distribution of VH genes between the treated mice and the controls, resulting in larger Jensen-Shannon divergence between them. We noticed that our analysis of variations in IGHV gene usage successfully distinguished between the treatment and control groups.
Wardemann and Mefre previously showed that a longer immunoglobulin heavy chain CDR3 correlates with antibody autoreactivity in SLE (49, 50). We found that no significant difference in CDR3 mean length of IGHV in BAFF-R-Fc-treated mice compared to those observed in controls, consistent with a previous study (47). However, the distribution of CDR3 lengths was altered, with the most commonly used length decreasing from 15 amino acids to 14 amino acids.
Here, we systematically analyzed the transcriptome profiles of spleens from MRL/lpr mice and identified genes regulated by anti-BAFF therapy. Similar to a previous study on gene expression changes in patients treated with tabalumab (51), most of the downregulated genes were associated with B cell proliferation and the B cell receptor signaling pathway, confirming the cell type-specific effects of BAFF. Consistent with the decrease in B cell subsets we observed in the spleens of treated mice using flow cytometry, a similar decline in B cell subsets was observed in patients with SLE treated with belimumab (52).
Elevated ISGs expression in the peripheral blood cells of patients with SLE is associated with disease activity (53, 54). However, in two phase III trials of BAFF blockade with tabalumab, patients with improved disease activity showed stable IFN response gene expression (51). We analyzed downregulated DEGs and confirmed that no significant changes in the expression levels of ISGs in MRL/lpr mice treated with BAFF-R-Fc compared to the control. We then performed GSEA to identify the integrated expression levels of the IFN signature. The GSEA results successfully demonstrated a significant downregulation of IFN-α/γ and TNF-α responses after anti-BAFF therapy, consistent with the observed decrease in IFN-γ and TNF-α serum levels.
In our study, we detected a decrease in splenic IL-10 mRNA levels and serum IL-10 protein levels following BAFF blockade. The PPI network of IRGs suggested that IL-10 is involved in the regulation of B cell proliferation. While the role of IL-10 in SLE is controversial, significantly increased serum levels of IL-10 have been reported in patients, being positively correlated with SLE disease activity and dsDNA antibody titers (55, 56). In the early stages of disease, IL-10 inhibits IFN-γ-mediated autoantibody production and renal inflammation. However, in the later stages, excess IL-10 enhances major histocompatibility complex II expression on B cells and their differentiation into plasmablasts that secrete IgM and IgG (57–59). Thus, the downregulation of IL-10 expression level may contribute to mitigating disease activity.
The disruption of immune tolerance mechanisms leads to the expansion of B cell subpopulations and is a characteristic of SLE. Herein, BAFF blockade reduced the number of B cell subpopulations, including plasma cells, plasmablasts, and marginal zone B cells, in the spleen, which is consistent with previous studies (60, 61). We examined whether blocking BAFF resulted in the preferential reduction of certain plasma cell subpopulations. Our findings revealed significant reductions in the numbers of IgM, IgA, and IgG ASCs, consistent with the decreased serum immunoglobulin isotype levels. Our results noted that anti-BAFF therapy depleted memory B cells in MRL/lpr mice, in contrast to a previous study where patients treated with belimumab retained memory B cells (61). One possible reason is that the depletion of B cells occurs in the early phase of the disease in mice. Germinal center B cells were preserved following anti-BAFF therapy, in accordance with observations in patients receiving long-term belimumab treatment (60). In the absence of BAFF, germinal center responses were reportedly attenuated, yet still allowed antibody production (62). The depletion of CD1D+CD5+ B10 Bregs was also observed after BAFF blockade. Iwata et al. reported a significant increase in blood B10 cell frequency among patients and mice with autoimmune diseases compared to healthy controls (42, 63). The survival of CD19+-/- NZB/W mice greatly improved following the adoptive transfer of B10 cells (64), with these cells suppressing the generation of IFN-γ, TNF-α, and pathogenic autoantibodies in MRL/lpr mice, independently from IL-10 (65). The above-described findings raise the question of whether preserving the B10 subpopulation while blocking BAFF to deplete B cells will achieve better therapeutic outcomes.
Early studies reported that pDCs were chronically activated and continuously released IFN-α in patients with SLE (66, 67). DNA-containing immune complexes from SLE patient sera promote the synthesis of chemokines and cytokines by pDCs through cooperative interactions between CD32 and Toll-like receptor 9 (68, 69). Here, we observed a decrease in the proportion of pDCs and the downregulation of IFN-α response in MRL/lpr mice after BAFF blockade. However, more experimental research is needed to fully understand the mechanism.
Our study had several limitations. First, we only obtained the bulk transcriptome and BCR repertoire. Thus, we were unable to examine gene expression levels in different immune cell types and determine interactions between immune cells. Second, the lack of databases with information on immune receptors in mice prevented us from annotating our data and identifying autoimmunity-associated receptors. We could not confirm whether the highly expanded clonotypes were related to autoantibody. Third, the prediction of pDC and other myeloid cells by cell deconvolution algorithms was not accurate enough (70, 71), and required experimental verification. Finally, a large-scale clinical trial is required to confirm the potential of BCR sequencing in predicting disease prognosis.
Data availability statement
The RNA-seq and BCR repertoire data during the current study are available in the SRA repository, accession: PRJNA1017616.
Ethics statement
The animal study was approved by Animal Ethics Committee of Tongji University (No. TJLAC-018-032). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
TH: Conceptualization, Data curation, Formal Analysis, Writing – original draft. CP: Methodology, Writing – original draft. XX: Formal Analysis, Writing – original draft. YF: Methodology, Writing – original draft. JZ: Validation, Writing – original draft. HG: Funding acquisition, Supervision, Writing – review & editing. JF: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work is supported by the National Basic Research Program of China (973 Program) (2015CB553706), the National Natural Science Foundation of China (NSFC82172604), and the Shanghai Science & Technology Basic Research Program (18JC1414400). The authors declare that this study received funding from RemeGen Co., Ltd. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.
Acknowledgments
We would like to thank Editage for editing and reviewing this manuscript for the English language.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1307392/full#supplementary-material
References
1. Kiriakidou M, Ching CL. Systemic lupus erythematosus. Ann Internal Med (2020) 172(11):Itc81–itc96. doi: 10.7326/AITC202006020
2. Arbuckle MR, McClain MT, Rubertone MV, Scofield RH, Dennis GJ, James JA, et al. Development of autoantibodies before the clinical onset of systemic lupus erythematosus. New Engl J Med (2003) 349(16):1526–33. doi: 10.1056/NEJMoa021933
3. Moore PA, Belvedere O, Orr A, Pieri K, LaFleur DW, Feng P, et al. BLyS: member of the tumor necrosis factor family and B lymphocyte stimulator. Sci (New York NY). (1999) 285(5425):260–3. doi: 10.1126/science.285.5425.260
4. Cheema GS, Roschke V, Hilbert DM, Stohl W. Elevated serum B lymphocyte stimulator levels in patients with systemic immune-based rheumatic diseases. Arthritis Rheumatol (2001) 44(6):1313–9. doi: 10.1002/1529-0131(200106)44:6<1313::AID-ART223>3.0.CO;2-S
5. Steri M, Orrù V, Idda ML, Pitzalis M, Pala M, Zara I, et al. Overexpression of the cytokine BAFF and autoimmunity risk. New Engl J Med (2017) 376(17):1615–26. doi: 10.1056/NEJMoa1610528
6. Navarra SV, Guzmán RM, Gallacher AE, Hall S, Levy RA, Jimenez RE, et al. Efficacy and safety of belimumab in patients with active systemic lupus erythematosus: a randomised, placebo-controlled, phase 3 trial. Lancet (London England) (2011) 377(9767):721–31. doi: 10.1016/S0140-6736(10)61354-2
7. Dhillon S. Telitacicept: first approval. Drugs (2021) 81(14):1671–5. doi: 10.1007/s40265-021-01591-1
8. Tipton CM, Fucile CF, Darce J, Chida A, Ichikawa T, Gregoretti I, et al. Diversity, cellular origin and autoreactivity of antibody-secreting cell population expansions in acute systemic lupus erythematosus. Nat Immunol (2015) 16(7):755–65. doi: 10.1038/ni.3175
9. Odendahl M, Jacobi A, Hansen A, Feist E, Hiepe F, Burmester GR, et al. Disturbed peripheral B lymphocyte homeostasis in systemic lupus erythematosus. J Immunol (2000) 165(10):5970–9. doi: 10.4049/jimmunol.165.10.5970
10. Zhang Y, Lee TY. Revealing the immune heterogeneity between systemic lupus erythematosus and rheumatoid arthritis based on multi-omics data analysis. Int J Mol Sci (2022) 23(9):5166. doi: 10.3390/ijms23095166
11. Liu S, Hou XL, Sui WG, Lu QJ, Hu YL, Dai Y. Direct measurement of B-cell receptor repertoire’s composition and variation in systemic lupus erythematosus. Genes immunity. (2017) 18(1):22–7. doi: 10.1038/gene.2016.45
12. Chiche L, Jourde-Chiche N, Pascual V, Chaussabel D. Current perspectives on systems immunology approaches to rheumatic diseases. Arthritis Rheumatol (2013) 65(6):1407–17. doi: 10.1002/art.37909
13. Chiche L, Jourde-Chiche N, Whalen E, Presnell S, Gersuk V, Dang K, et al. Modular transcriptional repertoire analyses of adults with systemic lupus erythematosus reveal distinct type I and type II interferon signatures. Arthritis Rheumatol (Hoboken NJ). (2014) 66(6):1583–95. doi: 10.1002/art.38628
14. Rai R, Chauhan SK, Singh VV, Rai M, Rai G. RNA-seq analysis reveals unique transcriptome signatures in systemic lupus erythematosus patients with distinct autoantibody specificities. PloS One (2016) 11(11):e0166312. doi: 10.1371/journal.pone.0166312
15. Banchereau R, Hong S, Cantarel B, Baldwin N, Baisch J, Edens M, et al. Personalized immunomonitoring uncovers molecular networks that stratify lupus patients. Cell (2016) 165(3):551–65. doi: 10.1016/j.cell.2016.03.008
16. Akthar M, Nair N, Carter LM, Vital EM, Sutton E, McHugh N, et al. Correction: Deconvolution of whole blood transcriptomics identifies changes in immune cell composition in patients with systemic lupus erythematosus (SLE) treated with mycophenolate mofetil. Arthritis Res Ther (2023) 25(1):160. doi: 10.1186/s13075-023-03160-1
17. Thapa M, Tedesco D, Gumber S, Elrod EJ, Han JH, Kitchens WH, et al. Blockade of BAFF reshapes the hepatic B cell receptor repertoire and attenuates autoantibody production in cholestatic liver disease. J Immunol (2020) 204(12):3117–28. doi: 10.4049/jimmunol.1900391
18. Smulski CR, Eibel H. BAFF and BAFF-receptor in B cell selection and survival. Front Immunol (2018) 9:2285. doi: 10.3389/fimmu.2018.02285
19. Metzler G, Kolhatkar NS, Rawlings DJ. BCR and co-receptor crosstalk facilitate the positive selection of self-reactive transitional B cells. Curr Opin Immunol (2015) 37:46–53. doi: 10.1016/j.coi.2015.10.001
20. Balomenos D, Rumold R, Theofilopoulos AN. Interferon-gamma is required for lupus-like disease and lymphoaccumulation in MRL-lpr mice. J Clin Invest. (1998) 101(2):364–71. doi: 10.1172/JCI750
21. Team I. immunarch: an R package for painless bioinformatics analysis of T-cell and B-cell immune repertoires. Zenodo (2019). 0.4. 1. doi: 10.5281/zenodo.3383240
22. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol (2010) 11(10):R106.
23. Ge SX, Son EW, Yao R. iDEP: an integrated web application for differential expression and pathway analysis of RNA-Seq data. BMC Bioinf (2018) 19(1):534. doi: 10.1186/s12859-018-2486-6
24. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics: J Integr Biol (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118
25. Miao YR, Xia M, Luo M, Luo T, Yang M, Guo AY. ImmuCellAI-mouse: a tool for comprehensive prediction of mouse immune cell abundance and immune microenvironment depiction. Bioinf (Oxford England). (2022) 38(3):785–91. doi: 10.1093/bioinformatics/btab711
26. Han H, Cho J-W, Lee S, Yun A, Kim H, Bae D, et al. TRRUST v2: an expanded reference database of human and mouse transcrip tional regulatory interactions. Nucleic Acids Res (2018) 46(D1):D380–D6. doi: 10.1093/nar/gkx1013
27. Mincham KT, Young JD, Strickland DH. OMIP 076: High-dimensional immunophenotyping of murine T-cell, B-cell, and antibody secreting cell subsets. Cytometry Part A: J Int Soc Analytical Cytology. (2021) 99(9):888–92. doi: 10.1002/cyto.a.24474
28. Chaudhary N, Wesemann DR. Analyzing immunoglobulin repertoires. Front Immunol (2018) 9:462. doi: 10.3389/fimmu.2018.00462
29. D’Angelo S, Ferrara F, Naranjo L, Erasmus MF, Hraber P, Bradbury ARM. Many routes to an antibody heavy-chain CDR3: necessary, yet insufficient, for specific binding. Front Immunol (2018) 9:395. doi: 10.3389/fimmu.2018.00395
30. Nakano M, Ota M, Takeshima Y, Iwasaki Y, Hatano H, Nagafuchi Y, et al. Distinct transcriptome architectures underlying lupus establishment and exacerbation. Cell (2022) 185(18):3375–89.e21. doi: 10.1016/j.cell.2022.07.021
31. Sasaki T, Bracero S, Keegan J, Chen L, Cao Y, Stevens E, et al. Longitudinal immune cell profiling in patients with early systemic lupus erythematosus. Arthritis Rheumatol (Hoboken NJ). (2022) 74(11):1808–21. doi: 10.1002/art.42248
32. Wang S, Wang J, Kumar V, Karnell JL, Naiman B, Gross PS, et al. IL-21 drives expansion and plasma cell differentiation of autoreactive CD11c hi T-bet + B cells in SLE. Nat Commun (2018) 9(1):1758. doi: 10.1038/s41467-018-03750-7
33. Fan X, Ng CT, Guo D, Lim F, Tan JC, Law A, et al. Dampened inflammation and improved survival after CXCL5 administration in murine lupus via myeloid and neutrophil pathways. Arthritis Rheumatol (Hoboken NJ). (2023) 75(4):553–66. doi: 10.1002/art.42383
34. Chaussabel D, Quinn C, Shen J, Patel P, Glaser C, Baldwin N, et al. A modular analysis framework for blood genomics studies: application to systemic lupus erythematosus. Immunity (2008) 29(1):150–64. doi: 10.1016/j.immuni.2008.05.012
35. Rönnelid J, Tejde A, Mathsson L, Nilsson-Ekdahl K, Nilsson B. Immune complexes from SLE sera induce IL10 production from normal peripheral blood mononuclear cells by an FcgammaRII dependent mechanism: implications for a possible vicious cycle maintaining B cell hyperactivity in SLE. Ann Rheum Dis (2003) 62(1):37–42. doi: 10.1136/ard.62.1.37
36. Ramanujam M, Steffgen J, Visvanathan S, Mohan C, Fine JS, Putterman C. Phoenix from the flames: Rediscovering the role of the CD40-CD40L pathway in systemic lupus erythematosus and lupus nephritis. Autoimmun Rev (2020) 19(11):102668. doi: 10.1016/j.autrev.2020.102668
37. Li Y, Xie X, Jie Z, Zhu L, Yang JY, Ko CJ, et al. DYRK1a mediates BAFF-induced noncanonical NF-κB activation to promote autoimmunity and B-cell leukemogenesis. Blood (2021) 138(23):2360–71. doi: 10.1182/blood.2021011247
38. Somasundaram R, Jensen CT, Tingvall-Gustafsson J, Åhsberg J, Okuyama K, Prasad M, et al. EBF1 and PAX5 control pro-B cell expansion via opposing regulation of the Myc gene. Blood (2021) 137(22):3037–49. doi: 10.1182/blood.2020009564
39. Xu C, Xu J, Lu L, Tian W, Ma J, Wu M. Identification of key genes and novel immune infiltration-associated biomarkers of sepsis. Innate immunity. (2020) 26(8):666–82. doi: 10.1177/1753425920966380
40. Canté-Barrett K, Meijer MT, Cordo V, Hagelaar R, Yang W, Yu J, et al. MEF2C opposes Notch in lymphoid lineage decision and drives leukemia in the thymus. JCI Insight (2022) 7(13):e150363. doi: 10.1172/jci.insight.150363
41. Zhao J, Li L, Feng X, Gao C, Gao L, Zhan Y, et al. TIGIT-Fc fusion protein alleviates murine lupus nephritis through the regulation of SPI-B-PAX5-XBP1 axis-mediated B-cell differentiation. J autoimmunity. (2023) 139:103087. doi: 10.1016/j.jaut.2023.103087
42. Yanaba K, Bouaziz JD, Haas KM, Poe JC, Fujimoto M, Tedder TF. A regulatory B cell subset with a unique CD1dhiCD5+ phenotype controls T cell-dependent inflammatory responses. Immunity (2008) 28(5):639–50. doi: 10.1016/j.immuni.2008.03.017
43. Sáez Moya M, Gutiérrez-Cózar R, Puñet-Ortiz J, Rodríguez de la Concepción ML, Blanco J, Carrillo J, et al. Autoimmune B cell repertoire in a mouse model of sjögren’s syndrome. Front Immunol (2021) 12:666545. doi: 10.3389/fimmu.2021.666545
44. Lee YN, Frugoni F, Dobbs K, Tirosh I, Du L, Ververs FA, et al. Characterization of T and B cell repertoire diversity in patients with RAG deficiency. Sci Immunol (2016) 1(6):eaah6109. doi: 10.1126/sciimmunol.aah6109
45. Clarke T, Du P, Kumar S, Okitsu SL, Schuette M, An Q, et al. Autoantibody repertoire characterization provides insight into the pathogenesis of monogenic and polygenic autoimmune diseases. Front Immunol (2023) 14:1106537. doi: 10.3389/fimmu.2023.1106537
46. Bashford-Rogers RJM, Bergamaschi L, McKinney EF, Pombal DC, Mescia F, Lee JC, et al. Analysis of the B cell receptor repertoire in six immune-mediated diseases. Nature (2019) 574(7776):122–6. doi: 10.1038/s41586-019-1595-3
47. Wu M, Pan W, Jia C, He Z, Zhao M, Tang C, et al. Systemic lupus erythematosus patients contain B-cell receptor repertoires sensitive to immunosuppressive drugs. Eur J Immunol (2022) 52(4):669–80. doi: 10.1002/eji.202149596
48. Huang W, Quach TD, Dascalu C, Liu Z, Leung T, Byrne-Steele M, et al. Belimumab promotes negative selection of activated autoreactive B cells in systemic lupus erythematosus patients. JCI Insight (2018) 3(17):e122525. doi: 10.1172/jci.insight.122525
49. Wardemann H, Yurasov S, Schaefer A, Young JW, Meffre E, Nussenzweig MC. Predominant autoantibody production by early human B cell precursors. Sci (New York NY). (2003) 301(5638):1374–7. doi: 10.1126/science.1086907
50. Meffre E, Milili M, Blanco-Betancourt C, Antunes H, Nussenzweig MC, Schiff C. Immunoglobulin heavy chain expression shapes the B cell receptor repertoire in human B cell development. J Clin Invest. (2001) 108(6):879–86. doi: 10.1172/JCI13051
51. Hoffman RW, Merrill JT, Alarcón-Riquelme MM, Petri M, Dow ER, Nantz E, et al. Gene expression and pharmacodynamic changes in 1,760 systemic lupus erythematosus patients from two phase III trials of BAFF blockade with tabalumab. Arthritis Rheumatol (Hoboken NJ). (2017) 69(3):643–54. doi: 10.1002/art.39950
52. Wallace DJ, Stohl W, Furie RA, Lisse JR, McKay JD, Merrill JT, et al. randomized, double-blind, placebo-controlled, dose-ranging study of belimumab in patients with active systemic lupus erythematosus. Arthritis Rheumatol (2009) 61(9):1168–78. doi: 10.1002/art.24699
53. Deng Y, Zheng Y, Li D, Hong Q, Zhang M, Li Q, et al. Expression characteristics of interferon-stimulated genes and possible regulatory mechanisms in lupus patients using transcriptomics analyses. EBioMedicine (2021) 70:103477. doi: 10.1016/j.ebiom.2021.103477
54. Nehar-Belaid D, Hong S, Marches R, Chen G, Bolisetty M, Baisch J, et al. Mapping systemic lupus erythematosus heterogeneity at the single-cell level. Nat Immunol (2020) 21(9):1094–106. doi: 10.1038/s41590-020-0743-0
55. Chun HY, Chung JW, Kim HA, Yun JM, Jeon JY, Ye YM, et al. Cytokine IL-6 and IL-10 as biomarkers in systemic lupus erythematosus. J Clin Immunol (2007) 27(5):461–6. doi: 10.1007/s10875-007-9104-0
56. Gateva V, Sandling JK, Hom G, Taylor KE, Chung SA, Sun X, et al. A large-scale replication study identifies TNIP1, PRDM1, JAZF1, UHRF1BP1 and IL10 as risk loci for systemic lupus erythematosus. Nat Genet (2009) 41(11):1228–33. doi: 10.1038/ng.468
57. Heine G, Drozdenko G, Grün JR, Chang HD, Radbruch A, Worm M. Autocrine IL-10 promotes human B-cell differentiation into IgM- or IgG-secreting plasmablasts. Eur J Immunol (2014) 44(6):1615–21. doi: 10.1002/eji.201343822
58. Yin Z, Bahtiyar G, Zhang N, Liu L, Zhu P, Robert ME, et al. IL-10 regulates murine lupus. J Immunol (2002) 169(4):2148–55. doi: 10.4049/jimmunol.169.4.2148
59. Amend A, Wickli N, Schäfer AL, Sprenger DTL, Manz RA, Voll RE, et al. Dual role of interleukin-10 in murine NZB/W F1 lupus. Int J Mol Sci (2021) 22(3):1347. doi: 10.3390/ijms22031347
60. Jacobi AM, Huang W, Wang T, Freimuth W, Sanz I, Furie R, et al. Effect of long-term belimumab treatment on B cells in systemic lupus erythematosus: extension of a phase II, double-blind, placebo-controlled, dose-ranging study. Arthritis Rheumatol (2010) 62(1):201–10. doi: 10.1002/art.27189
61. Stohl W, Hiepe F, Latinis KM, Thomas M, Scheinberg MA, Clarke A, et al. Belimumab reduces autoantibodies, normalizes low complement levels, and reduces select B cell populations in patients with systemic lupus erythematosus. Arthritis Rheumatol (2012) 64(7):2328–37. doi: 10.1002/art.34400
62. Rahman ZS, Rao SP, Kalled SL, Manser T. Normal induction but attenuated progression of germinal center responses in BAFF and BAFF-R signaling-deficient mice. J Exp Med (2003) 198(8):1157–69. doi: 10.1084/jem.20030495
63. Iwata Y, Matsushita T, Horikawa M, Dilillo DJ, Yanaba K, Venturi GM, et al. Characterization of a rare IL-10-competent B-cell subset in humans that parallels mouse regulatory B10 cells. Blood (2011) 117(2):530–41. doi: 10.1182/blood-2010-07-294249
64. Watanabe R, Ishiura N, Nakashima H, Kuwano Y, Okochi H, Tamaki K, et al. Regulatory B cells (B10 cells) have a suppressive role in murine lupus: CD19 and B10 cell deficiency exacerbates systemic autoimmunity. J Immunol (2010) 184(9):4801–9. doi: 10.4049/jimmunol.0902385
65. Yu M, Song Y, Zhu MX, Liang W, Long Q, Ding PW, et al. B10 cells ameliorate the progression of lupus nephritis by attenuating glomerular endothelial cell injury. Cell Physiol biochemistry: Int J Exp Cell physiology biochemistry Pharmacol (2015) 36(6):2161–9. doi: 10.1159/000430182
66. Rönnblom L, Alm GV. Systemic lupus erythematosus and the type I interferon system. Arthritis Res Ther (2003) 5(2):68–75. doi: 10.1186/ar625
67. Chaudhary V, Ah Kioon MD, Hwang SM, Mishra B, Lakin K, Kirou KA, et al. Chronic activation of pDCs in autoimmunity is linked to dysregulated ER stress and metabolic responses. . J Exp Med (2022) 219(11):e20221085. doi: 10.1084/jem.20221085
68. Means TK, Latz E, Hayashi F, Murali MR, Golenbock DT, Luster AD. Human lupus autoantibody-DNA complexes activate DCs through cooperation of CD32 and TLR9. J Clin Invest. (2005) 115(2):407–17. doi: 10.1172/JCI23025
69. Båve U, Magnusson M, Eloranta ML, Perers A, Alm GV, Rönnblom L. Fc gamma RIIa is expressed on natural IFN-alpha-producing cells (plasmacytoid dendritic cells) and is required for the IFN-alpha production induced by apoptotic cells combined with lupus IgG. J Immunol (2003) 171(6):3296–302. doi: 10.4049/jimmunol.171.6.3296
70. Jiménez-Sánchez A, Cast O, Miller ML. Comprehensive benchmarking and integration of tumor microenvironment cell estimation methods. Cancer Res (2019) 79(24):6238–46. doi: 10.1158/0008-5472.CAN-18-3560
Keywords: systemic lupus erythematosus, B cell receptor (BCR) repertoire, transcriptome, B cell-activating factor (BAFF), MRL/lpr
Citation: Huang T, Pi C, Xu X, Feng Y, Zhang J, Gu H and Fang J (2024) Effect of BAFF blockade on the B cell receptor repertoire and transcriptome in a mouse model of systemic lupus erythematosus. Front. Immunol. 14:1307392. doi: 10.3389/fimmu.2023.1307392
Received: 04 October 2023; Accepted: 18 December 2023;
Published: 09 January 2024.
Edited by:
Jing Jing Wang, Flinders University, AustraliaReviewed by:
Lisa K. Blum, Consultant, South San Francisco, CA, United StatesKirsten Fairfax, University of Tasmania, Australia
Copyright © 2024 Huang, Pi, Xu, Feng, Zhang, Gu and Fang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hua Gu, gu_hua@tongji.edu.cn; Jianmin Fang, jfang@tongji.edu.cn