Abstract
1,25-Dihydroxyvitamin D3 (1,25(OH)2D3) can modulate immune responses, but whether it directly affects B cell function is unknown. Patients with systemic lupus erythematosus, especially those with antinuclear Abs and increased disease activity, had decreased 1,25(OH)2D3 levels, suggesting that vitamin D might play a role in regulating autoantibody production. To address this, we examined the effects of 1,25(OH)2D3 on B cell responses and found that it inhibited the ongoing proliferation of activated B cells and induced their apoptosis, whereas initial cell division was unimpeded. The generation of plasma cells and postswitch memory B cells was significantly inhibited by 1,25(OH)2D3, although the up-regulation of genetic programs involved in B cell differentiation was only modestly affected. B cells expressed mRNAs for proteins involved in vitamin D activity, including 1α-hydroxylase, 24-hydroxylase, and the vitamin D receptor, each of which was regulated by 1,25(OH)2D3 and/or activation. Importantly, 1,25(OH)2D3 up-regulated the expression of p27, but not of p18 and p21, which may be important in regulating the proliferation of activated B cells and their subsequent differentiation. These results indicate that 1,25(OH)2D3 may play an important role in the maintenance of B cell homeostasis and that the correction of vitamin D deficiency may be useful in the treatment of B cell-mediated autoimmune disorders.
The biologically active form of vitamin D3, 1,25-dihydroxyvitamin D3 (1,25(OH)2D3),3 which was originally described as an essential hormone for bone and mineral homeostasis, has also been shown to have immunomodulatory effects (1). Vitamin D3 can be derived from dietary sources, but the main source of vitamin D3 is obtained through photosynthesis in the skin. Exposure of skin to UV light (270–300 nm) catalyzes the first step in vitamin D3 biosynthesis, converting 7-dehydroxycholesterol into previtamin D3, which is followed by a spontaneous and temperature-dependent isomerization into vitamin D3. To obtain the biologically active metabolite 1,25(OH)2D3, vitamin D3 must first be hydroxylated by D3-25-hydroxylase (CYP2D25) in the liver into 25-hydroxyvitamin D3 (25(OH)D3). Further hydroxylation of 25(OH)D3 by 25(OH)D3-1-α-hydroxylase (CYP27B1) occurs mainly in the proximal convoluted tubule cells of the kidney, resulting in the biologically active 1,25(OH)2D3 (2). In addition, CYP27B1 is expressed by activated macrophages and dendritic cells (3, 4), but their contribution to systemic vitamin D metabolism is unknown. 1,25(OH)2D3 also induces the expression of the major 1,25(OH)2D3-degrading enzyme 24-hydroxylase (CYP24A1) (2). Most of the known biological effects of 1,25(OH)2D3 are mediated through the vitamin D3 receptor (VDR), a member of the superfamily of nuclear hormone receptors (5) that is expressed by different cells of the immune system and elsewhere (5, 6).
The immunoregulatory effects of 1,25(OH)2D3 are mainly thought to be mediated through its action on APCs, with the most potent reported effects on dendritic cells (DCs) (7, 8). The in vitro differentiation of DCs from monocytes or murine bone marrow-derived precursors is inhibited by 1,25(OH)2D3. Moreover, the Ag-presenting function of monocytes and DCs is profoundly inhibited, as is the surface expression of costimulatory molecules and also IL-12 production (9). In addition, expression of the immunosuppressive cytokine IL-10 by DCs, which opposes the Th1-inducing effects of IL-12, is increased by 1,25(OH)2D3 (7, 8).
Direct effects of 1,25(OH)2D3 on T lymphocytes have also been reported. T cell proliferation and cell cycle progression from G1a to G1b are inhibited in vitro by 1,25(OH)2D3. Cytokine production by T cells is also modulated, with inhibition of the Th1 cytokine IFN-γ and an increase of the Th2 cytokines IL-4, IL-5, and IL-10 (10). Thus, 1,25(OH)2D3 is thought to polarize activated T cells toward a Th2 phenotype. In addition, it has been shown that expression of the Fas ligand by activated T cells could be repressed by 1,25(OH)2D3 (11).
The effects of l,25(OH)2D3 on B cell function has not been examined in detail. Published reports have yielded conflicting observations. It has been claimed that 1,25(OH)2D3 may indirectly inhibit B cell function as a consequence of the impairment of CD4 T cell responses or the inhibition of cytokine production by monocytes/macrophages (12). Alternatively, a direct inhibitory effect of 1,25(OH)2D3 on IgE production by human B cells has also been reported (13).
Interestingly, decreased 1,25(OH)2D3 serum concentrations have been reported in many autoimmune diseases, including systemic lupus erythematosus (SLE) (14, 15). SLE is an autoimmune disease characterized by immune dysregulation resulting in overproduction of autoantibodies (16). Although the exact cause of SLE remains unknown, recent studies, including a demonstration of the effectiveness of therapeutic B cell depletion, strongly implicate B cells as central players in the pathogenesis of this disease (17, 18, 19, 20). We therefore, hypothesized that 1,25(OH)2D3 may be important in maintaining B cell homeostasis and that deficiency of 1,25(OH)2D3 might contribute to B cell hyperactivity in SLE.
In this study, we measured 1,25(OH)2D3 and 25(OH)D3 levels in patients with SLE and correlated them with disease activity and the presence of antinuclear Abs (ANA). Subsequently, we investigated the direct effects of 1,25(OH)2D3 on primary human B cells to gain more insight into the potential role of vitamin D on autoimmune disease pathogenesis. Our results provide evidence that vitamin D might be a useful alternative therapy for the B cell hyperactivity characteristic of SLE.
Materials and Methods
Clinical samples
For the analysis of serological levels of 25(OH)D and 1,25(OH)2D3, peripheral blood samples from patients with SLE, patients with rheumatoid arthritis (RA), and demographically comparable healthy controls were obtained from the Shanghai Clinical Center for Rheumatology, Renji Hospital, Shanghai, China after informed consent. The diagnosis of SLE relied on American College of Rheumatology revised criteria (21). Clinical charts of all patients were reviewed and disease activity was scored according to the SLE disease activity index (SLEDAI) (22). For the analysis of B cell function in vitro, blood samples from normal healthy adult donors were obtained from the Warren G. Magnuson Clinical Center Blood Bank (Bethesda, MD). The collection of samples was approved by the National Institute of Arthritis and Musculoskeletal and Skin Diseases/National Institute of Diabetes & Digestive & Kidney Diseases (NIAMS/NIDDK) Institutional Review Board (Bethesda, MD), and informed consent was obtained according to the Declaration of Helsinki.
Determination of serum 25(OH)D and 1,25(OH)2D3 levels
Serum total 25(OH)D and 1,25(OH)2D3 levels were measured in SLE patients, patients with RA, and healthy donors by RIA or ELISA (both from ALPCO Diagnostics) according to manufacturer’s instructions. Both assays yielded similar results.
B cell enrichment, flow cytometry, and cell sorting
B cells were enriched by negative selection from buffy coats or leukapheresis samples using RosetteSep or StemSep B cell purification Ab mixtures (StemCell Technologies). Enriched B cells were >90 and 95% pure, respectively. Purified B cells were stained with various mAb combinations for 20 min on ice in staining buffer (1% BSA and 5% FCS in PBS). The directly conjugated mAb used were anti-IgD-FITC, anti-CD27-PE, anti-CD40-PE, anti-CD86-FITC, anti-IgG-PE, anti-CD19-allophycocyanin, anti-HLA-DR-PE, (BD Pharmingen), anti-CD19-PerCpCy5.5, and anti-CD38-allophycocyanin (clone HB7) (BD Immunocytometry). Stained cells were washed and data were collected immediately using a four-color FACScalibur (BD Biosciences) or fixed in 1% paraformaldehyde and analyzed within 24 h. Data were analyzed using FlowJo software (Tree Star). B cell populations were sorted using the DakoCytomation MoFlo. Sorted B cell populations were >98% pure.
Cell apoptosis and viability assay
To quantify apoptosis, cultured B cells were double-stained with annexin V-FITC conjugate and propidium iodide (PI) using TACS annexin V kits from Oncogene according to the manufacturer’s directions. The stained B cells were immediately analyzed using a FACSCalibur flow cytometer (BD Biosciences).
B cell stimulation
B cells were cultured in 96-well plates at 1 × 106 cells/ml in 100-μl volumes of culture medium (10% FCS in RPMI 1640 medium with l-glutamine and penicillin-streptomycin) alone or supplemented with various cytokines and stimuli: 1 μg/ml anti-CD40 (R&D Systems), 3 μg/ml goat F(ab′)2 anti-IgM (Jackson ImmunoReseach Laboratories), 250 ng/ml IL-21 (Cell Sciences), 50 ng/ml IL-4 (R&D Systems), and 1,25(OH)2D3 (usually 10 nM) or 25(OH) D3 (both from Sigma-Aldrich).
Proliferation assay
Cells were cultured for 3, 4, or 5 days as described above and pulsed for an additional 16 h with [3H]thymidine (1 μCi, 37 kBq) and then harvested. The [3H]thymidine incorporation was measured with a Top Count microplate scintillation counter (Packard Instruments). Alternatively, proliferation was assessed by CFSE (Invitrogen Life Technologies) dilution. Briefly, 1 × 107 B cells were incubated with 1 μM of CFSE for 15 min at room temperature and then quenched with an equal volume of 100% FCS for 2 min and washed. CFSE-labeled cells were cultured for 3 to 6 days as mentioned above and data were collected by FACS analysis.
Ig ELISA
Culture supernatants were incubated on goat antihuman IgG-and IgM-coated (Bethyl Laboratories) Maxisorp plates (Nalge Nunc International) and bound Ig was detected using alkaline phosphatase-conjugated goat anti-human IgG or IgM (Bethyl Laboratories) with Sigma Fast pNPP alkaline phosphate substrate (Sigma-Aldrich).
ELISPOT assay
The number of IgG- and IgM-secreting cells was determined with the ELISPOT assay. MultiScreenHTS plates (polyvinylidene difluoride membrane; Millipore) were coated overnight at 4°C with 50 μl/well goat anti-human IgG or IgM (Bethyl Laboratories), respectively, at a concentration of 5 μg/ml diluted in 0.05 M carbonate PBS (Sigma-Aldrich). After washing, unbound surfaces of wells were blocked with culture medium for 1 h at 25°C. Stimulated cells from 7-day cultures were washed three times and resuspended with medium containing 10% Ultra Low IgG FCS (Invitrogen Life Technologies). Serial dilutions of cells were added to wells in triplicate and incubated at 37°C overnight. Medium containing 10% Ultra Low IgG FCS was included in all experiments as a negative control. Thereafter, plates were washed with TBS plus 0.05% Tween 20 four times and incubated with either biotinylated goat-anti-human IgM or IgG (Invitrogen) diluted 1/2000 in TBS plus 0.05% Tween 20 at 25°C for 3 h. Plates were then washed in TBS plus 0.05% Tween 20 and incubated with 0.5 μg/ml alkaline phosphatase-conjugated streptavidin (Sigma-Aldrich) in TBS plus 0.05% Tween 20 for 1 h at 25°C. All plates were developed with the alkaline phosphatase substrate kit III (Vector Laboratories) after washing in distilled water. Ig-secreting cells appeared as blue spots, which were counted using a CTL Series 3B analyzer (Cellular Technology).
RNA preparation and quantitative RT-PCR
Total RNA from B cells was purified by using the RNeasy mini kit (Qiagen) according to the manufacturer’s manual. Reverse transcription reactions were prepared using the SuperScript one-step PCR system with platinum Taq polymerase and ROX reference dye (Invitrogen Life Technologies). Fifty nanograms of isolated RNA was added per reaction with 1.2 mM MgSO4. TaqMan Assays-on-Demand gene expression primer/probe sets for β2-microglobulin (β2M; assay identifier Hs99999907_ m1), VDR (Hs01045840_m1), CYP27B1 (Hs00168017_m1), CYP24A1 (Hs00167999_m1), Pax5 (paired box gene 5; Hs00277134_m1), Bcl6 (B cell CLL/lymphoma 6, zinc finger protein 51; Hs00153368_m1), BLIMP1(B-lymphocyte-induced maturation protein 1; Hs00153357_m1), IFN regulatory factor 4 (IRF4; Hs00180031_m1), X box-binding protein 1 (XBP1; Hs00231936_m1), endoplasmic reticulum to nucleus signaling 1 (ERN1; Hs00176385_m1), activation-induced cytidine deaminase (AID; Hs00221068_m1), p27 (cyclin-dependent kinase (CDK) inhibitor 1B; Hs00153277_m1), p21 (CDK inhibitor 1A; Hs00355782_m1), and p18 (CDK inhibitor 2C; Hs00176227_m1) were purchased from Applied Biosystems, and quantitative RT-PCR was performed according to the manufacturer’s instructions. mRNA expression for each gene was calculated in triplicate using the comparative cycle threshold method with efficiency calculations, and all mRNA levels were normalized to β2M.
The effects of 1,25(OH)2D3 on the expression of 84 cell cycle regulated genes and 84 genes associated with NF-κB-mediated signal transduction were examined using human cell cycle RT2 Profiler PCR array and human NFκB signaling pathway RT2 Profiler PCR array, respectively (SuperArray Bioscience) according to the manufacturer’s instructions. In brief, cDNA was prepared from 1 μg total RNA by using a RT2 PCR array first strand kit. A total volume of 25 μl of PCR mixture, which included 12.5 μl of RT2 Real-Time SYBR Green/ROX PCR master mix from SuperArray Bioscience (containing HotStart DNA polymerase, SYBR Green dye, and the ROX reference dye), 11.5 μl of double-distilled H2O, and 1 μl of template cDNA, was loaded in each well of the PCR array. PCR amplification was conducted with an initial 10-min step at 95°C followed by 40 cycles of 95°C for 15 s and 60°C for 1 min). The fluorescent signal from SYBR Green was detected immediately after the extension step of each cycle, and the cycle at which the product was first detectable was recorded as the cycle threshold. Data were imported into an Excel database and analyzed using the comparative cycle threshold method with normalization of the raw data to housekeeping genes including β2M, hypoxanthine phosphoribosyltransferase 1, ribosomal protein L13a, GAPDH, and ACTB (β-actin).
Analysis and statistics
Statistical significance of the data was analyzed with a two-tailed, nonpaired or paired Student’s t test, using Microsoft Excel Windows software. Data are plotted as mean ± SEM. An asterisk (∗) denotes p < 0.05, ∗∗ denotes p < 0.01, and ∗∗∗ denotes p < 0.001.
Results
SLE patients have significantly diminished levels of 25(OH) D and 1,25(OH)2D3
A total of 112 SLE patients who fulfilled American College of Rheumatology revised criteria (21) were assessed. Their characteristics are shown in Table I. The serum 25(OH)D level was significantly lower in SLE patients (11.5 ± 1.5 ng/ml, n = 57) than in demographically comparable healthy controls (59.2 ± 6.5 ng/ml, n = 28) and RA patients (54.6 ± 5.2 ng/ml, n = 29) (p < 0.001) (Fig. 1,A). Notably, the level of 25(OH)D in newly diagnosed SLE patients without any treatment was 11.6 ± 2.1 ng/ml (n = 12), which was not statistically different from that in SLE patients with established disease (11.8 ± 1.8 ng/ml, n = 45) (p > 0.05) (Fig. 1,B). No correlation was found between the level of 25(OH)D and various clinical parameters including renal involvement, disease activity assessed by SLEDAI (22), and glucocorticoid use. The serum level of 1,25(OH)2D3 in SLE patients (14.5 ± 1.2 pg/ml, n = 87) was also significantly lower than in normal controls (29.8 ± 1.5 pg/ml, n = 30, p < 0.001) (Fig. 1,C). 1,25(OH)2D3 levels were significantly lower in patients with active SLE as indicated by a SLEDAI score of >4 (12.2 ± 1.6 pg/ml, n = 52) compared with those with a SLEDAI score of ≤4 (19.4 ± 1.9 pg/ml, n = 36, p < 0.001) (Fig. 1,D). Additional analysis demonstrated that the level of 1,25(OH)2D3 was significantly lower in patients with ANA compared with those without ANA (14.4 ± 1.3 pg/ml, n = 71 vs 22.0 ± 4.2 pg/ml, n = 11, p < 0.05), but no significant difference was observed in patients with or without anti-dsDNA (15.0 ± 1.8 pg/ml, n = 39 vs 13.5 ± 1.7 pg/ml, n = 43) (p > 0.05) (Fig. 1, E and F). The data indicated that patients with SLE had decreased levels of 25(OH)2D and 1,25(OH)2D3, and the latter was related to disease activity and ANA production. These results prompted us to examine the impact of vitamin D on B cell function.
Characteristic . | SLE Patients (n = 112) . |
---|---|
Age in years, mean ± SEM (range in years) | 33.8 ± 1.57 (15–66) |
Female (%) | 90.2% (101/112) |
Disease duration in years, mean ± SEM (range in years) | 4.5 ± 0.5 (0.01–20) |
SLEDAI score, mean ± SEM (range) | 6.5 ± 0.9 (0–22) |
Arthritis percentage (no./total) | 31.3% (35/112) |
Rash percentage (no./total) | 23.2% (26/112) |
Oral ulceration percentage (no./total) | 17.9% (20/112) |
Renal involvement percentage (no./total) | 37.5% (42/112) |
Vasculitis percentage (no./total) | 7.1% (8/112) |
Hematological abnormalities percentage (no./total) | 25% (28/112) |
CNS involvement percentage (no./total) | 3.6% (4/112) |
Glucocorticoid use, mean ± SEM mg/day (range) | 23.2 ± 2.9 (0–160) |
Characteristic . | SLE Patients (n = 112) . |
---|---|
Age in years, mean ± SEM (range in years) | 33.8 ± 1.57 (15–66) |
Female (%) | 90.2% (101/112) |
Disease duration in years, mean ± SEM (range in years) | 4.5 ± 0.5 (0.01–20) |
SLEDAI score, mean ± SEM (range) | 6.5 ± 0.9 (0–22) |
Arthritis percentage (no./total) | 31.3% (35/112) |
Rash percentage (no./total) | 23.2% (26/112) |
Oral ulceration percentage (no./total) | 17.9% (20/112) |
Renal involvement percentage (no./total) | 37.5% (42/112) |
Vasculitis percentage (no./total) | 7.1% (8/112) |
Hematological abnormalities percentage (no./total) | 25% (28/112) |
CNS involvement percentage (no./total) | 3.6% (4/112) |
Glucocorticoid use, mean ± SEM mg/day (range) | 23.2 ± 2.9 (0–160) |
1,25(OH)2D3 inhibits proliferation of activated B cells
Initially, the effect of 1,25(OH)2D3 on B cell survival was examined. 1,25(OH)2D3 did not affect the survival of unstimulated B cells in culture (Fig. 2,A). We next assessed the effect of 1,25(OH)2D3 on B cell proliferation. As shown in Fig. 2 B, 1,25(OH)2D3 exerted a significant inhibitory effect on B cell proliferation, especially in cultures activated with a combination of stimuli, e.g., anti-IgM/anti-CD40 (p < 0.05), anti-CD40/IL-21, or anti-IgM/anti-CD40/IL-21 (p < 0.01). B cells stimulated with either anti-CD40 or anti-IgM alone had a limited proliferative capacity and, as a consequence, the effect of 1,25(OH)2D3 was less consistent.
CFSE labeling was also used to examine the effect of 1,25(OH)2D3 on the proliferation of B cells. After 3 days in culture some CFSE dilution was observed. Notably, cultures with and without 1,25(OH)2D3 exhibited a similar frequency of cells with modestly diluted CFSE (Fig. 2 C). As the culture period was prolonged, cells with increasing dilution of CFSE were observed. However, the extent of CFSE dilution was strikingly diminished in 1,25(OH)2D3-containing cultures compared with control cultures. These results were consistent with the conclusion that 1,25(OH)2D3 had little effect on initial cell division, whereas ongoing proliferation of stimulated B cells was suppressed.
1,25(OH)2D3 induces apoptosis of activated B cells
To determine whether the inhibitory effect of 1,25(OH)2D3 on proliferation was associated with the induction of apoptosis, we examined the percentage of apoptotic cells by annexin V/PI staining. During 4 days in culture without stimulation 1,25(OH)2D3 treatment had no significant effect on survival or apoptosis (Fig. 3,A). Stimulation with anti-IgM/anti-CD40 decreased the frequency of apoptotic cells, but in cultures of activated B cells,1,25(OH)2D3 significantly increased the percentage of early apoptotic cells at day 4 (27.5% vs 37.5%, p < 0.01). This was time dependent, as this difference was not apparent at earlier time points (Fig. 3 B).
VDR and 24-hydroxylase can be induced by 1,25(OH)2D3 on B cells
We next examined the expression of the VDR and the known vitamin D responsive gene 24-hydroxylase (CYP24A1) by quantitative PCR. Freshly isolated B cells constitutively expressed very low levels of VDR mRNA. Compared with day 0 expression, the mean fold increase of VDR expression by B cells after 3 days of culture with anti-IgM/anti-CD40, anti-CD40/IL-21, anti-IgM/anti-CD40/IL-21, or without stimulation was 1.45, 2.88, 4.44, and 1.46-fold respectively. 1,25(OH)2D3 increased VDR expression in both activated and unstimulated B cells (Fig. 4,A), but not earlier than day 3 (data not shown). In contrast, CYP24A1 expression was not detectable in either resting or activated B cells. However, it was significantly induced by 1,25(OH)2D3(Fig. 4 B). The results demonstrate that the VDR is expressed and inducible in primary B cells. Activation signals or 1,25(OH)2D3 up-regulated VDR expression, whereas CYP24A1 was only induced by 1,25(OH)2D3 but not by stimulation.
1,25(OH)2D3 specifically increases CD38 expression on B cells
It was previously reported that 1,25(OH)2D3 could increase CD38 expression by lymphocytes (23, 24). To confirm this and to examine the effect of 1,25(OH)2D3 on other B cell markers, we activated purified B cells with various stimuli with and without 1,25(OH)2D3, and then stained the cells for various surface markers. The expression of CD38 on activated B cells was significantly increased by 1,25(OH)2D3. This difference was not observed with unstimulated B cells or cultures activated with anti-CD40 or anti-IgM alone, but it was apparent in B cells stimulated with the combination of anti-IgM/anti-CD40, anti-CD40/IL-21, or anti-IgM/anti-CD40/IL-21. In contrast, the expressions of IgD, CD21, CD23, CD27, CD86, and HLA-DR were altered by activating stimuli but not further affected by 1,25(OH)2D3 (Fig. 5 A). Similarly, 1,25(OH)2D3 had no effect on the expression of CD19 or CD20 (data not shown). To examine the expression of CD40, we used anti-IgM and IL-4 to activate B cells and up-regulate CD40 expression. Using these conditions, 1,25(OH)2D3 again increased CD38 expression (data not shown), but CD40 expression was not altered by 1,25(OH)2D3. These results show that 1,25(OH)2D3 specifically increases CD38 expression by activated B cells but has no effect on the expression of a variety of other lineage, differentiation, or activation markers.
Increased expression of CD38 is commonly used as a phenotypic marker of plasma cells (25). To determine whether the 1,25(OH)2D3-mediated increase in CD38 expression was associated with plasma cell differentiation, we used a method of in vitro stimulation known to induce plasma cells from peripheral B cells (anti-CD40 and IL-21 with or without anti-IgM) and examined the effect of 1,25(OH)2D3 (26). The induction of CD38 resulting from 1,25(OH)2D3 was not associated with an increase in CD27 expression and, indeed, fewer CD38highCD27high plasma cells were generated in the presence of 1,25(OH)2D3 (Fig. 5,B). Both the increase of CD38 expression and the decrease of plasma cell generation occurred in a concentration-dependent manner, with the maximum effect noted reproducibly at 10 nM or less (Fig. 5 C).
1,25(OH)2D3 inhibits plasma cell differentiation and Ig secretion
A more detailed analysis showed that 1,25(OH)2D3 inhibited plasma cell differentiation and the secretion of Ig (Fig. 6). The generation of plasma cells was inhibited throughout the culture period after day 5 and 1,25(OH)2D3 caused significant decreases in both the percentage (Fig. 6,A, left panel) and absolute number of plasma cells (data not shown). In multiple experiments, the generation of phenotypically defined CD19+/−CD38highCD27high plasma cells was inhibited by 23–34% (Fig. 6 A, right panel). Similar results were noted when B cells purified by cell sorting and >98% pure were examined (data not shown).
Analysis of culture supernatants revealed that 1,25(OH)2D3 also significantly decreased the secretion of Ig. Stimulation with anti-CD40 and IL-21 with or without anti-IgM induced considerable Ig secretion by day 5 that increased through day 9 (Fig. 6,B). 1,25(OH)2D3 significantly inhibited IgG secretion. The percent inhibition progressively increased as the cultures were prolonged with the maximal inhibition noted at day 9 (p < 0.001). A significant albeit smaller reduction in IgM production in 1,25(OH)2D3 cultures was only apparent at day 9. The effect of 1,25(OH)2D3 on the number of Ig-secreting cells generated in these cultures was determined by ELISPOT assay. In agreement with the phenotype and Ig production data, 1,25(OH)2D3 significantly reduced the number of IgG-secreting cells in cultures of activated B cells. The number of IgM secreting cells was also significantly decreased by 1,25(OH)2D3 in anti-CD40/IL-21 cultures (Fig. 6 C).
Notably, when 1,25(OH)2D3 was added at day 5 of culture, no effect on the maintenance of plasma cells or on Ig secretion was noted over the next 5 days (Fig. 6 D). This result indicates that 1,25(OH)2D3 inhibits the generation of plasma cells but not their subsequent persistence or Ig secretion.
1,25(OH)2D3 inhibits the generation of postswitched memory B cells from naive B cells
To determine whether the inhibitory effect of 1,25(OH)2D3 was specific for plasma cell generation, we also examined its role in regulating the generation of postswitch memory B cells. To address this, naive B cells (CD19+IgG−CD27−) were purified by sorting (>99% purity) as shown in Fig. 7,A and stimulated in a manner known to induce class switch recombination (26). As shown in Fig. 7,B, the combination of anti-CD40 and IL-21 with or without anti-IgM was able to induce postswitched cells expressing surface IgG. In the absence of IL-21, few cells expressed surface IgG. Notably, the frequency of IgG positive cells was significantly lower in the 1,25(OH)2D3-containing cultures (Fig. 7, B and C), whereas there were more B cells expressing surface IgM (Fig. 7,B). This reduction in class switching correlated with reduced plasma cell numbers and IgG production in parallel 7-day cultures (Fig. 7 D). These data demonstrate that 1,25(OH)2D3 inhibits generation of both postswitch memory cells and plasma cells.
1,25(OH)2D3 down-regulates the expression of XBP1 and ERN1 mRNA but not that of PAX5, B cell lymphoma 6 (BCL6), AID, B lymphocyte-induced maturation protein 1 (BLIMP1), MTA3, and IRF4
Thus far, the results clearly demonstrated that 1,25(OH)2D3 reduced B cell proliferation, the generation of postswitched memory B cells, plasma cell differentiation, and Ig secretion. However, the stage at which 1,25(OH)2D3 exerts its effects remained unclear. To gain further insight into the nature of B cell inhibition by 1,25(OH)2D3, we examined the expression of key transcription factors involved in B cell maturation and plasma cell differentiation. We compared the expression of genes involved in germinal center reactions, including PAX5, BCL6, and AID as well as those involved in plasma cell differentiation, including BLIMP1, IRF4, MTA3, XBP1, and ERN1 (27). All of the analyses were conducted at day 3 before significant effects on the frequency of memory B cells and plasma cells or Ig secretion were evident. Interestingly, 1,25(OH)2D3 had no significant effect on the expression of PAX5, BCL6, AID, BLIMP1, MAT3, or IRF4. However, even at this relatively early time point the expressions of XBP1 and ERN1 were significantly reduced but not eliminated by 1,25(OH)2D3 (Fig. 8). The lack of a consistent effect of 1,25(OH)2D3 on the genetic programs regulating the differentiation of memory cells or plasma cells suggests that it may inhibit their maturation by other means.
1,25(OH)2D3 increases the expression of p27 mRNA but not that of p21 and p18
To determine the mechanism involved in the inhibitory effect of 1,25(OH)2D3 on B cells, we initially screened an array of NF-κB-mediated signal transduction pathway genes. As shown in Table II, a few mRNAs were up- or down-regulated by 1,25(OH)2D3, but no general pattern of gene regulation was observed consistent with a major effect on the NF-κB signaling pathway. We next examined the effect of 1,25(OH)2D3 on cell cycle-related gene expression by B cells (shown in Table III). We found a number of genes, including cyclin D1, cyclin D2, cyclin T1, cyclin T2, CDK4, and CDK6, were significantly down-regulated by 1,25(OH)2D3. In contrast, the CDK inhibitor p27 was significantly up-regulated by 1,25(OH)2D3, whereas most genes, including CDK inhibitor p21, remained unchanged. We verified these findings using quantitative PCR and found that p27, but not p21 or p18 mRNA, was significantly up-regulated by 1,25(OH)2D3 in activated B cells (Fig. 9).
mRNA . | Mean of Fold Change . | p Value . | ||
---|---|---|---|---|
Angiotensinogen (serpin peptidase inhibitor, clade A, member 8) | 0.55 | 0.09 | ||
v-akt murine thymoma viral oncogene homolog 1 | 1.10 | 0.34 | ||
Activating transcription factor 1 | 1.19 | 0.11 | ||
B cell CLL/lymphoma 10 | 1.16 | 0.11 | ||
B cell CLL/lymphoma 3 | 2.12 | 0.16 | ||
B factor, properdin | 0.68 | 0.05 | ||
Baculoviral IAP repeat-containing 2 | 0.73 | 0.12 | ||
Caspase recruitment domain family, member 4 | 0.79 | 0.07 | ||
Caspase 1, apoptosis-related cysteine peptidase (interleukin 1β, convertase) | 0.85 | 0.29 | ||
Caspase 8, apoptosis-related cysteine peptidase | 0.87 | 0.35 | ||
Chemokine (C-C motif) ligand 2 | 0.09 | 0.12 | ||
CD40 antigen (TNF receptor superfamily member 5) | 0.85 | 0.20 | ||
CASP8 and FADD-like apoptosis regulator | 1.29 | 0.12 | ||
Conserved helix-loop-helix ubiquitous kinase | 1.08 | 0.42 | ||
Colony-stimulating factor 2 (granulocyte-macrophage) | 0.51 | 0.27 | ||
Colony stimulating factor 3 (granulocyte) | 0.63 | 0.16 | ||
Solute carrier family 44, member 2 | 1.45 | 0.00 | ||
EDAR-associated death domain | 0.48 | 0.04 | ||
Endothelial differentiation, lysophosphatidic acid G-protein-coupled receptor, 2 | 0.89 | 0.39 | ||
Early growth response 1 | 1.37 | 0.22 | ||
ELK1, member of ETS oncogene family | 0.88 | 0.30 | ||
Coagulation factor II (thrombin) receptor | 0.59 | 0.10 | ||
Fas (TNFRSF6)-associated via death domain | 0.88 | 0.25 | ||
Fas ligand (TNF superfamily, member 6) | 0.84 | 0.37 | ||
v-fos FBJ murine osteosarcoma viral oncogene homolog | 0.83 | 0.32 | ||
Gap junction protein, α1, 43 kDa (connexin 43) | 0.16 | 0.20 | ||
Heme oxygenase (decycling) 1 | 1.36 | 0.21 | ||
5-Hydroxytryptamine (serotonin) receptor 2B | 0.64 | 0.04 | ||
Intercellular adhesion molecule 1 (CD54), human rhinovirus receptor | 0.91 | 0.40 | ||
Interferon, α1 | 0.61 | 0.16 | ||
Interferon, β1, fibroblast | 0.55 | 0.14 | ||
Interferon, γ | 0.80 | 0.39 | ||
Inhibitor of κ light polypeptide gene enhancer in B cells, kinase β | 1.11 | 0.39 | ||
Inhibitor of κ light polypeptide gene enhancer in B cells, kinase ε | 0.86 | 0.33 | ||
Inhibitor of κ light polypeptide gene enhancer in B cells, kinase γ | 0.81 | 0.04 | ||
Interleukin 10 | 0.81 | 0.27 | ||
Interleukin 1, α | 0.70 | 0.01 | ||
Interleukin 1, β | 0.38 | 0.05 | ||
Interleukin 1 receptor, type I | 0.71 | 0.14 | ||
Interleukin 6 (interferon, β2) | 3.27 | 0.01 | ||
Interleukin 8 | 0.76 | 0.22 | ||
Interleukin-1 receptor-associated kinase 1 | 0.66 | 0.19 | ||
Interleukin-1 receptor-associated kinase 2 | 3.09 | 0.09 | ||
v-jun sarcoma virus 17 oncogene homolog (avian) | 0.93 | 0.36 | ||
Lymphotoxin α (TNF superfamily, member 1) | 0.60 | 0.15 | ||
Lymphotoxin β receptor (TNFR superfamily, member 3) | 0.96 | 0.47 | ||
Mucosa-associated lymphoid tissue lymphoma translocation gene 1 | 0.88 | 0.24 | ||
Mitogen-activated protein kinase kinase kinase 1 | 1.25 | 0.23 | ||
Myeloid differentiation primary response gene (88) | 1.01 | 0.48 | ||
NACHT, leucine rich repeat and PYD containing 12 | 0.26 | 0.19 | ||
Nuclear factor of κ light polypeptide gene enhancer in B-cells 1 (p105) | 0.95 | 0.40 | ||
Nuclear factor of κ light polypeptide gene enhancer in B-cells 2 (p49/p100) | 1.00 | 0.50 | ||
(Table continues) |
mRNA . | Mean of Fold Change . | p Value . | ||
---|---|---|---|---|
Angiotensinogen (serpin peptidase inhibitor, clade A, member 8) | 0.55 | 0.09 | ||
v-akt murine thymoma viral oncogene homolog 1 | 1.10 | 0.34 | ||
Activating transcription factor 1 | 1.19 | 0.11 | ||
B cell CLL/lymphoma 10 | 1.16 | 0.11 | ||
B cell CLL/lymphoma 3 | 2.12 | 0.16 | ||
B factor, properdin | 0.68 | 0.05 | ||
Baculoviral IAP repeat-containing 2 | 0.73 | 0.12 | ||
Caspase recruitment domain family, member 4 | 0.79 | 0.07 | ||
Caspase 1, apoptosis-related cysteine peptidase (interleukin 1β, convertase) | 0.85 | 0.29 | ||
Caspase 8, apoptosis-related cysteine peptidase | 0.87 | 0.35 | ||
Chemokine (C-C motif) ligand 2 | 0.09 | 0.12 | ||
CD40 antigen (TNF receptor superfamily member 5) | 0.85 | 0.20 | ||
CASP8 and FADD-like apoptosis regulator | 1.29 | 0.12 | ||
Conserved helix-loop-helix ubiquitous kinase | 1.08 | 0.42 | ||
Colony-stimulating factor 2 (granulocyte-macrophage) | 0.51 | 0.27 | ||
Colony stimulating factor 3 (granulocyte) | 0.63 | 0.16 | ||
Solute carrier family 44, member 2 | 1.45 | 0.00 | ||
EDAR-associated death domain | 0.48 | 0.04 | ||
Endothelial differentiation, lysophosphatidic acid G-protein-coupled receptor, 2 | 0.89 | 0.39 | ||
Early growth response 1 | 1.37 | 0.22 | ||
ELK1, member of ETS oncogene family | 0.88 | 0.30 | ||
Coagulation factor II (thrombin) receptor | 0.59 | 0.10 | ||
Fas (TNFRSF6)-associated via death domain | 0.88 | 0.25 | ||
Fas ligand (TNF superfamily, member 6) | 0.84 | 0.37 | ||
v-fos FBJ murine osteosarcoma viral oncogene homolog | 0.83 | 0.32 | ||
Gap junction protein, α1, 43 kDa (connexin 43) | 0.16 | 0.20 | ||
Heme oxygenase (decycling) 1 | 1.36 | 0.21 | ||
5-Hydroxytryptamine (serotonin) receptor 2B | 0.64 | 0.04 | ||
Intercellular adhesion molecule 1 (CD54), human rhinovirus receptor | 0.91 | 0.40 | ||
Interferon, α1 | 0.61 | 0.16 | ||
Interferon, β1, fibroblast | 0.55 | 0.14 | ||
Interferon, γ | 0.80 | 0.39 | ||
Inhibitor of κ light polypeptide gene enhancer in B cells, kinase β | 1.11 | 0.39 | ||
Inhibitor of κ light polypeptide gene enhancer in B cells, kinase ε | 0.86 | 0.33 | ||
Inhibitor of κ light polypeptide gene enhancer in B cells, kinase γ | 0.81 | 0.04 | ||
Interleukin 10 | 0.81 | 0.27 | ||
Interleukin 1, α | 0.70 | 0.01 | ||
Interleukin 1, β | 0.38 | 0.05 | ||
Interleukin 1 receptor, type I | 0.71 | 0.14 | ||
Interleukin 6 (interferon, β2) | 3.27 | 0.01 | ||
Interleukin 8 | 0.76 | 0.22 | ||
Interleukin-1 receptor-associated kinase 1 | 0.66 | 0.19 | ||
Interleukin-1 receptor-associated kinase 2 | 3.09 | 0.09 | ||
v-jun sarcoma virus 17 oncogene homolog (avian) | 0.93 | 0.36 | ||
Lymphotoxin α (TNF superfamily, member 1) | 0.60 | 0.15 | ||
Lymphotoxin β receptor (TNFR superfamily, member 3) | 0.96 | 0.47 | ||
Mucosa-associated lymphoid tissue lymphoma translocation gene 1 | 0.88 | 0.24 | ||
Mitogen-activated protein kinase kinase kinase 1 | 1.25 | 0.23 | ||
Myeloid differentiation primary response gene (88) | 1.01 | 0.48 | ||
NACHT, leucine rich repeat and PYD containing 12 | 0.26 | 0.19 | ||
Nuclear factor of κ light polypeptide gene enhancer in B-cells 1 (p105) | 0.95 | 0.40 | ||
Nuclear factor of κ light polypeptide gene enhancer in B-cells 2 (p49/p100) | 1.00 | 0.50 | ||
(Table continues) |
mRNA . | Mean of Fold Change . | p Value . |
---|---|---|
Nuclear factor of κ light polypeptide gene enhancer in B-cells inhibitor, α | 1.56 | 0.03 |
Protein phosphatase 1A (formerly 2C), magnesium-dependent, α isoform | 0.91 | 0.28 |
v-raf-1 murine leukemia viral oncogene homolog 1 | 0.92 | 0.38 |
v-rel reticuloendotheliosis viral oncogene homolog (avian) | 0.70 | 0.21 |
v-rel reticuloendotheliosis viral oncogene homolog A, nuclear factor of κ light polypeptide ger | 1.10 | 0.36 |
v-rel reticuloendotheliosis viral oncogene homolog B, nuclear factor of κ light polypeptide ger | 0.84 | 0.23 |
Ret finger protein 2 | 0.92 | 0.24 |
Ras homolog gene family, member A | 1.08 | 0.36 |
Receptor (TNFRSF)-interacting serine-threonine kinase 1 | 1.19 | 0.25 |
Solute carrier family 20 (phosphate transporter), member 1 | 0.80 | 0.20 |
Signal transducer and activator of transcription 1, 91 kDa | 1.01 | 0.49 |
TANK-binding kinase 1 | 0.98 | 0.45 |
Toll-like receptor adaptor molecule 2 | 0.55 | 0.02 |
Toll-like receptor 1 | 0.75 | 0.16 |
Toll-like receptor 2 | 0.30 | 0.10 |
Toll-like receptor 3 | 1.47 | 0.26 |
Toll-like receptor 4 | 1.04 | 0.46 |
Toll-like receptor 6 | 0.83 | 0.07 |
Toll-like receptor 7 | 0.65 | 0.01 |
Toll-like receptor 8 | 0.72 | 0.34 |
Toll-like receptor 9 | 0.69 | 0.22 |
Transmembrane emp24 protein transport domain containing 4 | 0.99 | 0.48 |
Tumor necrosis factor (TNF superfamily, member 2) | 0.48 | 0.07 |
Tumor necrosis factor, α-induced protein 3 | 0.74 | 0.24 |
Tumor necrosis factor receptor superfamily, member 10a | 0.56 | 0.05 |
Tumor necrosis factor receptor superfamily, member 10b | 0.82 | 0.24 |
Tumor necrosis factor receptor superfamily, member 1A | 0.53 | 0.01 |
Tumor necrosis factor receptor superfamily, member 7 | 0.95 | 0.44 |
Tumor necrosis factor (ligand) superfamily, member 10 | 1.16 | 0.28 |
Tumor necrosis factor (ligand) superfamily, member 14 | 1.82 | 0.23 |
TNFRSF1A-associated via death domain | 1.28 | 0.10 |
Toll-like receptor adaptor molecule 1 | 1.39 | 0.10 |
mRNA . | Mean of Fold Change . | p Value . |
---|---|---|
Nuclear factor of κ light polypeptide gene enhancer in B-cells inhibitor, α | 1.56 | 0.03 |
Protein phosphatase 1A (formerly 2C), magnesium-dependent, α isoform | 0.91 | 0.28 |
v-raf-1 murine leukemia viral oncogene homolog 1 | 0.92 | 0.38 |
v-rel reticuloendotheliosis viral oncogene homolog (avian) | 0.70 | 0.21 |
v-rel reticuloendotheliosis viral oncogene homolog A, nuclear factor of κ light polypeptide ger | 1.10 | 0.36 |
v-rel reticuloendotheliosis viral oncogene homolog B, nuclear factor of κ light polypeptide ger | 0.84 | 0.23 |
Ret finger protein 2 | 0.92 | 0.24 |
Ras homolog gene family, member A | 1.08 | 0.36 |
Receptor (TNFRSF)-interacting serine-threonine kinase 1 | 1.19 | 0.25 |
Solute carrier family 20 (phosphate transporter), member 1 | 0.80 | 0.20 |
Signal transducer and activator of transcription 1, 91 kDa | 1.01 | 0.49 |
TANK-binding kinase 1 | 0.98 | 0.45 |
Toll-like receptor adaptor molecule 2 | 0.55 | 0.02 |
Toll-like receptor 1 | 0.75 | 0.16 |
Toll-like receptor 2 | 0.30 | 0.10 |
Toll-like receptor 3 | 1.47 | 0.26 |
Toll-like receptor 4 | 1.04 | 0.46 |
Toll-like receptor 6 | 0.83 | 0.07 |
Toll-like receptor 7 | 0.65 | 0.01 |
Toll-like receptor 8 | 0.72 | 0.34 |
Toll-like receptor 9 | 0.69 | 0.22 |
Transmembrane emp24 protein transport domain containing 4 | 0.99 | 0.48 |
Tumor necrosis factor (TNF superfamily, member 2) | 0.48 | 0.07 |
Tumor necrosis factor, α-induced protein 3 | 0.74 | 0.24 |
Tumor necrosis factor receptor superfamily, member 10a | 0.56 | 0.05 |
Tumor necrosis factor receptor superfamily, member 10b | 0.82 | 0.24 |
Tumor necrosis factor receptor superfamily, member 1A | 0.53 | 0.01 |
Tumor necrosis factor receptor superfamily, member 7 | 0.95 | 0.44 |
Tumor necrosis factor (ligand) superfamily, member 10 | 1.16 | 0.28 |
Tumor necrosis factor (ligand) superfamily, member 14 | 1.82 | 0.23 |
TNFRSF1A-associated via death domain | 1.28 | 0.10 |
Toll-like receptor adaptor molecule 1 | 1.39 | 0.10 |
mRNA . | Mean of Fold Change . | p Value . |
---|---|---|
Cyclin-dependent kinase 6 | 0.83 | 0.12 |
Cyclin-dependent kinase 7 (MO15 homolog, Xenopus laevis, cdk-activating kinase) | 0.99 | 0.43 |
Cyclin-dependent kinase 8 | 0.79 | 0.13 |
Cyclin-dependent kinase inhibitor 1A (p21, Cip1) | 0.74 | 0.08 |
Cyclin-dependent kinase inhibitor 1B (p27, Kip1) | 2.44 | 0.03 |
Cyclin-dependent kinase inhibitor 2A (melanoma, p16, inhibits CDK4) | 1.65 | 0.38 |
Cyclin-dependent kinase inhibitor 2B (p15, inhibits CDK4) | 0.60 | 0.11 |
Cyclin-dependent kinase inhibitor 3 (CDK2-associated dual specificity phosphatase) | 0.96 | 0.35 |
CHK1 checkpoint homolog (Schizosaccharomyces pombe) | 0.72 | 0.19 |
CHK2 checkpoint homolog (S. pombe) | 0.83 | 0.17 |
CDC28 protein kinase regulatory subunit 1B | 0.76 | 0.08 |
CDC28 protein kinase regulatory subunit 2 | 0.86 | 0.04 |
Cullin 1 | 0.98 | 0.36 |
Cullin 2 | 0.74 | 0.14 |
Cullin 3 | 0.71 | 0.10 |
DEAD/H (Asp-Glu-Ala-Asp/His) box polypeptide 11 (CHL1-like helicase homolog, Saccharomyces cerevisiae) | 1.13 | 0.23 |
Dynamin 2 | 1.13 | 0.31 |
E2F transcription factor 4, p107/p130 binding | 0.82 | 0.23 |
Growth arrest and DNA-damage inducible, α | 1.85 | 0.22 |
General transcription factor IIH, polypeptide 1, 62 kDa | 0.96 | 0.26 |
G2 and S phase expressed 1 | 0.93 | 0.15 |
Hect domain and RLD 5 | 0.84 | 0.13 |
HUS1 checkpoint homolog (S. pombe) | 0.85 | 0.18 |
Kinetochore associated 1 | 0.79 | 0.01 |
Karyopherin α2 (RAG cohort 1, importin α1) | 0.74 | 0.05 |
MAD2 mitotic arrest deficient-like 1 (yeast) | 0.92 | 0.27 |
MAD2 mitotic arrest deficient-like 2 (yeast) | 0.81 | 0.05 |
MCM2 minichromosome maintenance deficient 2, mitotin (S. cerevisiae) | 0.90 | 0.14 |
MCM3 minichromosome maintenance deficient 3 (S. cerevisiae) | 1.19 | 0.42 |
MCM4 minichromosome maintenance deficient 4 (S. cerevisiae) | 0.79 | 0.13 |
MCM5 minichromosome maintenance deficient 5, cell division cycle 46 (S. cerevisiae) | 1.21 | 0.32 |
Antigen identified by monoclonal antibody Ki-67 | 0.65 | 0.04 |
Ménage-à-trois 1 (CAK assembly factor) | 0.72 | 0.05 |
MRE11 meiotic recombination 11 homolog A (S. cerevisiae) | 0.73 | 0.02 |
Nibrin | 0.87 | 0.14 |
Proliferating cell nuclear antigen | 0.83 | 0.09 |
RAD1 homolog (S. pombe) | 1.01 | 0.48 |
RAD17 homolog (S. pombe) | 1.07 | 0.12 |
RAD51 homolog (RecA homolog, Escherichia coli) (S. cerevisiae) | 1.56 | 0.22 |
RAD9 homolog A (S. pombe) | 0.96 | 0.44 |
Retinoblastoma 1 (including osteosarcoma) | 0.94 | 0.27 |
Retinoblastoma binding protein 8 | 0.85 | 0.14 |
Retinoblastoma-like 1 (p107) | 1.31 | 0.27 |
Retinoblastoma-like 2 (p130) | 0.97 | 0.37 |
Retinoblastoma protein A3, 14 kDa | 1.36 | 0.22 |
SERTA domain containing 1 | 1.35 | 0.18 |
S-phase kinase-associated protein 2 (p45) | 0.70 | 0.18 |
SMT3 suppressor of mif two 3 homolog 1 (yeast) | 0.72 | 0.02 |
Transcription factor Dp-1 | 0.85 | 0.11 |
Transcription factor Dp-2 (E2F dimerization partner 2) | 0.91 | 0.26 |
Tumor protein p53 (Li-Fraumeni syndrome) | 0.97 | 0.37 |
Ubiquitin-activating enzyme E1 (A1S9T and BN75 temperature sensitivity complementing) | 2.86 | 0.21 |
mRNA . | Mean of Fold Change . | p Value . |
---|---|---|
Cyclin-dependent kinase 6 | 0.83 | 0.12 |
Cyclin-dependent kinase 7 (MO15 homolog, Xenopus laevis, cdk-activating kinase) | 0.99 | 0.43 |
Cyclin-dependent kinase 8 | 0.79 | 0.13 |
Cyclin-dependent kinase inhibitor 1A (p21, Cip1) | 0.74 | 0.08 |
Cyclin-dependent kinase inhibitor 1B (p27, Kip1) | 2.44 | 0.03 |
Cyclin-dependent kinase inhibitor 2A (melanoma, p16, inhibits CDK4) | 1.65 | 0.38 |
Cyclin-dependent kinase inhibitor 2B (p15, inhibits CDK4) | 0.60 | 0.11 |
Cyclin-dependent kinase inhibitor 3 (CDK2-associated dual specificity phosphatase) | 0.96 | 0.35 |
CHK1 checkpoint homolog (Schizosaccharomyces pombe) | 0.72 | 0.19 |
CHK2 checkpoint homolog (S. pombe) | 0.83 | 0.17 |
CDC28 protein kinase regulatory subunit 1B | 0.76 | 0.08 |
CDC28 protein kinase regulatory subunit 2 | 0.86 | 0.04 |
Cullin 1 | 0.98 | 0.36 |
Cullin 2 | 0.74 | 0.14 |
Cullin 3 | 0.71 | 0.10 |
DEAD/H (Asp-Glu-Ala-Asp/His) box polypeptide 11 (CHL1-like helicase homolog, Saccharomyces cerevisiae) | 1.13 | 0.23 |
Dynamin 2 | 1.13 | 0.31 |
E2F transcription factor 4, p107/p130 binding | 0.82 | 0.23 |
Growth arrest and DNA-damage inducible, α | 1.85 | 0.22 |
General transcription factor IIH, polypeptide 1, 62 kDa | 0.96 | 0.26 |
G2 and S phase expressed 1 | 0.93 | 0.15 |
Hect domain and RLD 5 | 0.84 | 0.13 |
HUS1 checkpoint homolog (S. pombe) | 0.85 | 0.18 |
Kinetochore associated 1 | 0.79 | 0.01 |
Karyopherin α2 (RAG cohort 1, importin α1) | 0.74 | 0.05 |
MAD2 mitotic arrest deficient-like 1 (yeast) | 0.92 | 0.27 |
MAD2 mitotic arrest deficient-like 2 (yeast) | 0.81 | 0.05 |
MCM2 minichromosome maintenance deficient 2, mitotin (S. cerevisiae) | 0.90 | 0.14 |
MCM3 minichromosome maintenance deficient 3 (S. cerevisiae) | 1.19 | 0.42 |
MCM4 minichromosome maintenance deficient 4 (S. cerevisiae) | 0.79 | 0.13 |
MCM5 minichromosome maintenance deficient 5, cell division cycle 46 (S. cerevisiae) | 1.21 | 0.32 |
Antigen identified by monoclonal antibody Ki-67 | 0.65 | 0.04 |
Ménage-à-trois 1 (CAK assembly factor) | 0.72 | 0.05 |
MRE11 meiotic recombination 11 homolog A (S. cerevisiae) | 0.73 | 0.02 |
Nibrin | 0.87 | 0.14 |
Proliferating cell nuclear antigen | 0.83 | 0.09 |
RAD1 homolog (S. pombe) | 1.01 | 0.48 |
RAD17 homolog (S. pombe) | 1.07 | 0.12 |
RAD51 homolog (RecA homolog, Escherichia coli) (S. cerevisiae) | 1.56 | 0.22 |
RAD9 homolog A (S. pombe) | 0.96 | 0.44 |
Retinoblastoma 1 (including osteosarcoma) | 0.94 | 0.27 |
Retinoblastoma binding protein 8 | 0.85 | 0.14 |
Retinoblastoma-like 1 (p107) | 1.31 | 0.27 |
Retinoblastoma-like 2 (p130) | 0.97 | 0.37 |
Retinoblastoma protein A3, 14 kDa | 1.36 | 0.22 |
SERTA domain containing 1 | 1.35 | 0.18 |
S-phase kinase-associated protein 2 (p45) | 0.70 | 0.18 |
SMT3 suppressor of mif two 3 homolog 1 (yeast) | 0.72 | 0.02 |
Transcription factor Dp-1 | 0.85 | 0.11 |
Transcription factor Dp-2 (E2F dimerization partner 2) | 0.91 | 0.26 |
Tumor protein p53 (Li-Fraumeni syndrome) | 0.97 | 0.37 |
Ubiquitin-activating enzyme E1 (A1S9T and BN75 temperature sensitivity complementing) | 2.86 | 0.21 |
B cells can directly metabolize 25(OH)D3 into a biologically active molecule
It has been shown that the precursor of 1,25(OH)2D3, 25(OH)D3, is present in serum and is converted to the active form by 1α-hydroxylase (CYP27B1), which is mainly expressed in kidney and also can be found in monocytes/macrophages (3, 4). The expression of 1α-hydroxylase (CYP27B1) mRNA by human peripheral B cells was therefore investigated by RT-PCR. As shown in Fig. 10,A, B cells constitutively expressed 1α-hydroxylase mRNA, and it was up-regulated following stimulation but not further induced by 1,25(OH)2D3. To determine whether this up-regulation had a functional consequence, we added the precursor, 25(OH)D3, to purified B cells stimulated with anti-CD40 and IL-21 and examined the expression of CD38, plasma cell differentiation, and IgG production. 25(OH)D3 increased CD38 expression and inhibited plasma cell generation and IgG production in a concentration-dependent manner (Fig. 10 B). The precursor had effects similar to those of the active form of vitamin D3, but at a 25-fold higher concentration. Because monocytes contaminating the peripheral B cell populations might contribute to these results, B cells that were sorted from tonsils and lacked detectable monocytes/macrophage contamination were assessed and similar results were obtained (data not shown). There results are consistent with the conclusion that human B cells can take up and metabolize 25(OH)D3 to active 1,25(OH)2D3 that can inhibit B cell function.
Discussion
In this study we show that 1,25(OH)2D3 has potent direct effects on B cell responses, inhibiting proliferation, generation of class-switched memory B cells, plasma cell differentiation, and Ig production. Reduced levels of 1,25(OH)2D3 in SLE patients, particularly in those patients with high disease activity scores and ANAs, suggest that vitamin D-dependent B cell regulation may play an important role in maintaining normal B cell homeostasis and that decreased levels of 1,25(OH)2D3 may contribute to B cell hyperactivity in SLE. Finally, because activated B cells express many of the molecules involved in vitamin D metabolism, their characteristic hyperactivity in SLE could contribute to the decreased levels of 1,25(OH)2D3 found in these patients.
The varied functions of 1,25(OH)2D3 are mainly mediated through binding to the VDR. We clearly demonstrate that VDR mRNA is constitutively expressed in human primary B cells at low levels and is up-regulated following stimulation and enhanced in the presence of 1,25(OH)2D3 in a time-dependent manner. Previously, the expression of VDR on B cells has been controversial. It has been claimed that resting B cells do not contain detectable amounts of VDR (13). However, it has also been reported that VDR is constitutively expressed on human tonsil B cells and can be further up-regulated by activation (13, 24) and that VDR is expressed on the B cell lymphoma cell lines SUDHL4 and SUDHL5 (28). Our data analyzing the expression of VDR mRNA and also the functional activity of 1,25(OH)2D3 clearly show that the VDR is expressed by human peripheral B cells. Moreover, our data clearly demonstrate that VDR expression by B cells is regulated both by 1,25(OH)2D3, as it is on other cells (5), but also by activation signals. These results indicate that vitamin D may exert differential effects on activated vs resting B cells and also may have different effects in individuals with different levels of serum 1,25(OH)2D3.
Active 1,25(OH)2D3 is mainly degraded by 24-hydroxylase (CYP24A1) (2). Microarray analysis from human colon carcinoma and ovarian cancer cell lines revealed that CYP24A1 was the most inducible gene responsive to 1,25(OH)2D3 (29, 30). In this report we show that 24-hydroxylase was significantly up-regulated following the incubation of human B cells with 1,25(OH)2D3. As opposed to the VDR, CYP24A1 was not altered by B cell activation. These results further demonstrate that human B cells can respond to 1,25(OH)2D3 directly. In addition, the results suggest that the activity of vitamin D on B cells might be influenced not only by VDR expression but also by the capacity to degrade the active molecule. The increased susceptibility of activated B cells to many of the effects of vitamin D might therefore reflect the up-regulation of VDR but not CYP24A1 by these cells.
25(OH)D3-1-α-hydroxylase (CYP27B1), the enzyme responsible for the final hydroxylation and activation of 25(OH)D3 into 1,25(OH)2D3, is mainly found in the proximal convoluted tubule cells of the kidney (2). Interestingly, we found that CYP27B1 mRNA was also expressed by resting B cells and could be further induced by stimulation, but not by 1,25(OH)2D3. Moreover, we found that the precursor, 25(OH)D3, had similar effects on purified B cells compared with the active form, albeit at higher concentrations. Therefore, 25(OH)D3 might be metabolized to 1,25(OH)2D3 by B cells themselves and may represent a source for the extrarenal synthesis of 1,25(OH)2D3, as is thought to be the case for macrophages and dendritic cells (3, 4). In conditions such as SLE, in which there is diffuse B cell activation (20), systemic vitamin D metabolism might therefore be significantly influenced.
The capacity of 1,25(OH)2D3 to inhibit proliferation has been reported in lymphocytes and in a variety of human cancer cell lines (31). We showed that 1,25(OH)2D3 also exerted an inhibitory effect on B cell proliferation. By using CFSE analysis, initial B cell division was not affected by 1,25(OH)2D3 whereas subsequent ongoing proliferation was progressively inhibited. This effect coincides with activation and 1,25(OH)2D3-mediated up-regulation of the VDR and suggests that a threshold level of VDR engagement might be required for the antiproliferative effect of the occupied VDR to become apparent. Notably, 1,25(OH)2D3-mediated inhibition of proliferation was associated with apoptosis of the activated and dividing B cells. These results suggest that differentiative events that require the initial expansion of B cells might be eliminated as a result of the 1,25(OH)2D3-mediated death of proliferating B cells.
The inhibition of proliferation and the induction of apoptosis are likely to result in the significant reduction in plasma cell differentiation and Ig production that we observed. Previous studies have shown that plasma cell differentiation requires initial B cell proliferation (32) and increases with the number of cell divisions (33). We have previously shown that the combination of IL-21 and anti-CD40 stimulation with or without BCR crosslinking is a potent inducer of proliferation and plasma cell differentiation and that all of the plasma cells generated had undergone extensive proliferation (26). Using this system, we demonstrated that 1,25(OH)2D3 had inhibitory effects on plasma differentiation and Ig production. Interestingly, this effect is not evident when the B cells are treated with 1,25(OH)2D3 after 5 days in culture, indicating that 1,25(OH)2D3 inhibits the generation of plasma cells but not their subsequent persistence and lends support to the conclusion that the inhibition of B cell proliferation by 1,25(OH)2D3 is responsible for the reduction of Ig-secreting cells and Ig production.
Several transcription factors regulate different stages of B cell maturation and plasma cell differentiation (34). The most notable finding of the current study was that expression of most of these transcriptional regulators was unaffected by 1,25(OH)2D3, including PAX5, BCL6, AID, BLIMP1, MTA3, and IRF4. These results provide no support for the possibility that 1,25(OH)2D3 inhibits memory cell or plasma cell differentiation by directly affecting the expression of any of these regulators of B cell maturation. However, the expression of XBP1 was modestly but significantly down-regulated by 1,25(OH)2D3. In addition, mRNA for ERN1, which is required for processing XBP1 mRNA to a spliced form that encodes a more stable and active protein (35), was also down-regulated by 1,25(OH)2D3. The inhibitory effect of 1,25(OH)2D3 on the expressions of XBP1 and ERN1 mRNA levels may explain the greater inhibitory effect of vitamin D on Ig secretion and detection of plasma cells by ELISPOT compared with its more modest inhibition of the differentiation of phenotypically defined plasma cells. Although inhibition of the expressions of XBP1 and ERN1 mRNA by 1,25(OH)2D3 may contribute to a reduced amount of Ig secretion per phenotypically defined plasma cell, it is unlikely to explain the inhibition of the generation of both plasma cells and postswitched memory cells as well as the ongoing proliferation of activated B cells by vitamin D.
Several important intracellular pathways have been reported to be inhibited by 1,25(OH)2D3, such as NF-κB activation and cell cycle progression. The suppressive effect of 1,25(OH)2D3 on the NF-κB signaling pathway was previously observed in T cells, monocytes or macrophages (9, 36, 37), and could have affected the expression of various essential cell surface and secreted molecules. The NF-κB pathway, however, did not appear to be generally suppressed by 1,25(OH)2D3 in activated B cells as assessed by the analysis of expression of various genes influenced by this signaling pathway. Because 1,25(OH)2D3 directly inhibited the proliferation of activated B cells, we therefore examined whether vitamin D might suppress the expression of a variety of cell cycle regulators and found that 1,25(OH)2D3 increased the expression of p27 and decreased the mRNA levels of CDK4 and 6 and cyclin D. Cell cycle progression is under the control of CDK, the activity of which is dependent on the expression of specific CDK inhibitors. Previous studies have shown that 1,25(OH)2D3 inhibits the proliferation of cancer cells by inducing gene transcription of the CDK inhibitors p21 and/or p27, which can inhibit cell cycle progression at the G1/S transition, depending on the cell type (38, 39, 40). In vivo, 1,25(OH)2D3 administration can effectively restore p27 accumulation in cancer cells as detected by immunohistochemistry and reduce the tumor burden (41). Several mechanisms have been reported to be involved, including increasing p27 transcription mediated by VDR-Sp1 interactions with the p27 promoter (42), stabilization of the protein through VDR-induced reduction of CDK2 activity and Skp2 abundance (the main factors responsible for p27 degradation), or induction of PTEN, a phosphatase that dephosphorylates p27 (40). Other CDK inhibitors have been shown to play an essential role in B cell responses. In this regard, p18 has been shown to be required for B cells to terminate proliferation and differentiate into functional plasma cells (38, 43, 44). Previous studies have confirmed that p18 and p27 were involved in regulating different steps of B proliferation. Whereas p27 was thought to regulate cell cycle entry by forming a ternary complex with CDKs and cyclins, p18 was directly involved in induction of the G1 cell cycle arrest by blocking the association of CDK4 and CDK6 with cyclin D (44). Interestingly, our data demonstrated that the mRNA level of p27, but not of p21 or p18, was up-regulated by 1,25(OH)2D3 in activated human B cells. These results suggest the possibility that 1,25(OH)2D3 could inhibit B cell proliferation by up-regulating p27 and thereby inhibit the cell cycle entry of previously cycling B cells. Because 1,25(OH)2D3 also decreased mRNA levels of CDK4 and CDK6 as well as cyclin D, the effect of p27 on ongoing B cell proliferation might be more profound. The 1,25(OH)2D3-mediated induction of p27 may limit ongoing B cell proliferation and thereby play an important role in the control of B cell responses. These findings suggest the possibility that the major effects of vitamin D on plasma cell and memory cell differentiation may result from the suppression of ongoing B cell proliferation, which is required before these differentiation steps can occur (32). Additional experiments will be necessary to test this hypothesis.
Decreased vitamin D concentrations have been reported in many autoimmune diseases such as multiple sclerosis, type 1 diabetes mellitus, RA, fibromyalgia, and SLE (14, 15, 45, 46). We therefore examined the serum levels of 25(OH)D and 1,25(OH)2D3 in Chinese lupus patients. Consistent with some previous findings, both 25(OH)D and 1,25(OH)2D3 levels were significantly lower in patients with SLE compared with healthy donors (14, 46). However, another study reported significantly lower 25(OH)D3 but not 1,25(OH)2D3 levels in SLE patients (47). The explanation for the differences in these studies is not clear but could relate to ethnicity or disease duration (46). Although previous studies suggested a possible relationship between depressed levels of 25(OH)D3 and the presence of lupus renal disease (46), we were able to demonstrate that lower levels of 1,25(OH)2D3 significantly correlated with increased disease activity. Moreover, significantly decreased vitamin D levels were detected in 12 newly diagnosed patients before any treatment was administered, indicating that the deficiency may be present at the onset of lupus and possibly before. This could explain, in part, why African Americans who are often chronically vitamin D deficient (46, 48, 49) also have a higher risk and severity of SLE (50). It is known that SLE may be preceded by autoantibody production that may exist for many years before diagnosis (51), and it is possible that vitamin D deficiency may contribute to this. The current results are consistent with the conclusion that vitamin D is involved in maintaining normal B cell homeostasis and that vitamin D deficiency in SLE patients may contribute to B cell hyperactivity, the breakdown of B cell tolerance, and the generation of some autoantibodies.
Consistent with this, we also found that those SLE patients with positive ANA titers had lower levels of 1,25(OH)2D3 than those without ANA. Interestingly, such a difference was not seen between the patients who were positive or negative for dsDNA Abs. The persistent and therapy-resistant presence of ANA in SLE even during the clinically quiescent phases has been explained by the production of these autoantibodies from long-lived plasma cells (52). In contrast, the presence of the dsDNA Ab is usually correlated with SLE activity and decreases after immunosuppressive or B cell-depletion therapy, suggesting that these autoantibodies originate from short-lived dividing plasmablasts (20, 53). Thus, 1,25(OH)2D3 may play a role in regulating the generation of long-lived plasma cell and be less potent in suppressing the rapid production of autoantibodies from dividing plasmablasts. The kinetics of the in vitro inhibition of B cell responses by 1,25(OH)2D3 is consistent with this possibility.
Besides the inhibitory effect on plasma cell differentiation and Ig production, our data also demonstrated that 1,25(OH)2D3 exerted its inhibitory effects on memory B cell differentiation. Relative increases in memory cells and plasma cells are both characteristic B cell abnormalities in lupus (54). In contrast to naive B cells and plasma cells, memory B cells are not targeted by conventional immunosuppressive therapies and may contribute to disease flares (55). Thus, vitamin D therapy may offer a new opportunity to induce remission in SLE by preventing or decreasing the differentiation of memory cells. This possibility is supported by murine studies in which administration of 1,25(OH)2D3 or its analog has been shown to decrease proteinuria or increase survival in a spontaneous mouse lupus model, MRL/lpr (56).
In conclusion, 1,25-dihydroxyvitamin D3 may play an important role in the maintenance of B cell homeostasis. Correction of vitamin D deficiency may be useful not only to prevent osteopenia but also in the suppression of B cell overactivity in autoimmune disorders such as SLE.
Acknowledgments
We gratefully acknowledge the assistance of the Flow Cytometry Section of the Office of Science and Technology, National Institute of Arthritis and Musculoskeletal and Skin Diseases (Bethesda, MD) for sorting B cell subsets. We are also grateful to Dr. John A. Hardin for many stimulating discussions and to Dr. Nancy Longo for reviewing the manuscript.
Disclosures
The authors have no financial conflict of interest.
Footnotes
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Abbreviations used in this paper: 1,25(OH)2D3, 1,25-dihydroxyvitamin D3; AID, activation-induced cytidine deaminase; ANA, antinuclear Ab; β2M, β2-microglobulin; BCL6, B cell lymphoma 6; BLIMP1, B lymphocyte-induced maturation protein 1; DC, dendritic cell; ERN1, endoplasmic reticulum to nucleus signaling; IRF4, IFN regulatory factor 4; PI, propidium iodide; RA, rheumatoid arthritis; SLE, systemic lupus erythematosus; SLEDAI, SLE disease activity index; VDR, vitamin D3 receptor; XBP1, X box-binding protein.