Vitamin d3-induced tolerogenic dendritic cells modulate the transcriptomic profile of t cd4 cells towards a functional hyporesponsiveness

Introduction

In the last years, tolerogenic dendritic cells (tolDC) have become one of the most promising alternatives for the treatment of autoimmune diseases, such as multiple sclerosis (MS), rheumatoid arthritis, or type 1 diabetes. In fact, several Phase I clinical trials have already finished or are currently ongoing, with positive results regarding the safety and the tolerability of this therapeutic cell-based approach ( 1 ). In general, tolDC are commonly defined as a stable and semi-mature subset of dendritic cells (DC), between antigen-capturing immature DC (iDC) and immunogenic mature DC (mDC)—characterized by their increased expression of MHC class II and co-stimulatory molecules. But most importantly, tolDC are presumably capable to induce immune tolerance towards the peptides these cells are presenting, in an antigen-specific manner ( 2 – 5 ).

TolDC can be generated in vitro from peripheral blood monocytes. In the last years, a wide variety of protocols for their production have been reported, ranging from the use of different drugs and chemical agents to genetic engineering techniques ( 6 , 7 ). In this regard, the use of 1, 25-dyhydroxyvitamin D3, the active form of vitamin D3, constitutes one of the most widely studied approaches for the differentiation of tolDC. Briefly, vitamin D3-induced tolDC (vitD3-tolDC) are thought to develop their regulatory properties through a semi-mature profile, their ability to inhibit or reduce T cell responses, and a switch of the immune response towards a T _H 2 profile ( 8 – 18 ). Furthermore, vitD3-tolDC are characterized by a reduced NF- κ B-mediated activity and an increase of mTOR-mediated glucose metabolism ( 10 , 19 ).

Even though tolDC—and vitD3-tolDC in particular—have been characterized with a developing knowledge over their metabolism, molecular mechanisms, and functional pathways, the specific effect of these cells over the rest of the immune-related components still remains elusive. It is known that tolDC can usually induce either anergy, hyporesponsiveness or depletion over activated T cells, as well as regulatory T cell (Treg) differentiation ( 20 ). However, to our knowledge, so far only one study has focused its attention on the actual processes that autologous T cells might be undergoing upon tolDC interaction—reporting an induction of hyporesponsiveness of CD4 ⁺ memory and naïve T cells towards antigen-specific stimulation mediated by dexamethasone-induced tolDC ( 21 )—but neither at the transcriptomic level nor with vitD3-tolDC in particular.

In previous studies, our group has already extensively characterized vitD3-tolDC phenotypically, functionally, and transcriptomically, evidencing the regulatory potential of these cells both in vitro and in vivo in the animal model of MS, experimental autoimmune encephalomyelitis (EAE) ( 13 , 16 , 22 – 24 ). Consequently, we wanted to take one step further for the elucidation of the mechanisms of immune tolerance induction of vitD3-tolDC. With that aim, here we present a full phenotypical, functional, and transcriptomic characterization of T CD4 ⁺ cells after their interaction with autologous vitD3-tolDC loaded with tetanus toxin (TT), in order to study the antigen-specific effect mediated by these cells compared to TT-loaded immunogenic mDC. The purpose of this study is to identify one or several potential biomarkers of the immune modulation developed by vitD3-tolDC over T cells, which could constitute an interesting tool for the monitoring of patients treated with these cells in clinical trials, and the understanding of the mechanisms of tolerance induction.

Material and Methods

Sample Collection

Buffy coat samples from 16 randomized healthy controls were obtained from the Banc de Sang i Teixits (Barcelona, Spain), according to the institutional Standard Operating Procedures for blood donation, including a signed informed consent. In parallel, whole blood samples from 12 different healthy donors were collected by standard venipuncture in lithium heparin tubes for the allogeneic functional assays (see below).

Monocyte Isolation

Healthy donor buffy coat samples were processed first depleting CD3 ⁺ cells using the RoseetteSep® Human Monocyte Enrichment Kit (StemCell Technologies, Vancouver, Canada) prior to a density gradient separation using ficoll-hypaque (Rafer, Zaragoza, Spain). Afterwards, CD14 ⁺ cells were isolated using the EasySep® Human CD14 Positive Selection Kit (StemCell), according to the manufacturer’s instructions. Cell viability was determined using 7-amino-actinomycin D (7-AAD) (BD Biosciences, Franklin Lakes, NK, USA) and phycoerythrin (PE)-conjugated annexin V (Immunotools, Friesoythe, Germany) staining for 20 min at 4°C, protected from light, and cell counts were quantified simultaneously using PerfectCount beads (Cytognos, Salamanca, Spain). Samples were acquired on a FACSCanto II flow cytometer (BD Biosciences), and monocyte purity was determined using forward and side scatter gating strategies on FACSDiva software (BD Biosciences).

TT-Loaded DC Cultures

The protocol for the generation of antigen-loaded tolDC was adapted from a previous study ( 25 ). Briefly, isolated monocytes were cultured for 6 days in 24-well plates at 37°C at a density of 1 × 10 ⁶ cells/ml in X-VIVO 15 medium, in the presence of 400 U/ml granulocyte-macrophage colony-stimulating factor (GM-CSF) and 500 U/ml IL-4 (both from Peprotech, London, UK). The whole volume of medium and cytokines was replenished on day 4. If no further treatment was performed, monocytes were differentiated into iDC. For the generation of mDC, we further added a maturation cocktail, containing 1, 000 U/ml IL-1β, 1, 000 U/ml TNF-α (both from Peprotech) and 1 µM prostaglandin E2 (PGE2) (Pfizer, New York, NY, USA) on day 4. Finally, in addition to the maturation cocktail, we added 1 nM vitamin D3 (Calcijex, Abbott, Chicago, IL, USA) on days 0 and 4 for the differentiation of vitD3-tolDC. For the generation of TT-loaded mDC (mDC-TT) and TT-loaded vitD3-tolDC (vitD3-tolDC-TT) as antigen-specific experimental conditions, 0. 1 µg/ml of the whole TT protein (Sigma-Aldrich, St. Louis, MO, USA) were added to the mDC and vitD3-tolDC cultures on day 3, 18 h before the addition of the maturation stimulus, while still in an immature status. On day 6, cells were harvested after an accutase (Invitrogen, Carlsbad, CA, USA) detaching treatment for 30 min, and washed twice. As shown above, cell counts and viability were determined by flow cytometry.

Autologous PBMC Isolation, Co-Culture, and Sorting

For the isolation of autologous PBMC, 3 ml of the buffy coat samples from each healthy donor was processed using a ficoll-hypaque density gradient separation and washed twice. Afterwards, cells were counted by flow cytometry, as described above, and plated in round-bottom 96-well plates at a density of 1 × 10 ⁶ cells/ml in RPMI medium supplemented with 10% heat-inactivated fetal bovine serum, 2 mM L-glutamine (Lonza, Basel, Switzerland), 100 U/ml penicillin (Reig Jofre, Sant Joan Despí, Spain) and 100 µg/ml streptomycin (Normon, Tres Cantos, Spain). The plates were then incubated for 6 days at 37°C in a 5% CO ₂ atmosphere. Afterwards, cells were harvested, and cell counts and viability were determined by flow cytometry.

Subsequently, an antigen-specific proliferation experimental setup was performed in 96-well round-bottom plates with co-cultures of 10 ⁵ autologous PBMC and 5, 000 either mDC-TT or vitD3-tolDC-TT (1: 20 ratio) in a final volume of 200 µl of supplemented RPMI medium. For each condition, 48 replicates were performed. Cells were then incubated for 5 days at 37°C in a 5% CO ₂ atmosphere. Afterwards, cells were harvested and the whole volume of each cell suspension was incubated for 20 min, protected from light, with the adequate amounts of monoclonal antibodies anti-CD3 Violet 450 (V450) and anti-CD4 PerCP-Cyanine dye (Cy)5. 5. Finally, cells were washed and the whole CD3 ⁺ CD4 ⁺ cell subpopulation was isolated and counted using a FACSAria II cell sorter (BD Biosciences). Finally, dry pellets of mDC-TT-conditioned CD3 ⁺ CD4 ⁺ cells (mDC-Tcell) and vitD3-tolDC-TT-conditioned CD3 ⁺ CD4 ⁺ cells (vitD3-Tcell) were obtained by centrifugation and stored at −80°C.

Phenotype Characterization of DC and Autologous PBMC

Surface expression of CD11c, CD14, CD83, CD86 and HLA-DR in iDC, mDC, mDC-TT, vitD3-tolDC and vitD3-tolDC-TT was determined by flow cytometry. In each case, DC suspensions were incubated for 20 min, protected from light, with the appropriate amounts of monoclonal antibodies anti-: CD11c PE-Cy7, CD14 V450, CD83 allophycocyanin (APC), CD86 fluorescein isothiocyanate (FITC) and HLA-DR Violet 500 (V500) (all of them from BD Biosciences). Subsequently, at least 10, 000 CD11c ⁺ events of each condition were acquired using a FACSCanto II flow cytometer and analyzed using FACSDiva software.

For the phenotypical characterization of mDC-Tcell and vitD3-Tcell, cell suspensions of these conditions were incubated for 20 min with the adequate amounts of monoclonal antibodies indicated below. Afterwards, samples were washed twice and acquired on a LSRFortessa flow cytometer, setting the stopping gate at 300, 000 peripheral blood mononuclear cells. The definition of each peripheral blood mononuclear cell subpopulation was determined as specified in Supplementary Table 1 , using several combinations of the following monoclonal antibodies anti-: CXCR3 AlexaFluor (AF)488, CD4 PerCP-Cy5. 5, CCR7 PE, CD45RA PE-Cy7, CD38 APC, CD45 AF700, CD8 APC-H7, CD3 V450, HLA-DR V500, CCR6 Brilliant Violet (BV) 605, CD25 PE, CCR4 PE-Cy7, CD127 AF647, CD45RO APC-H7, CD49b FITC and LAG-3 PE (BD Biosciences). Results were analyzed with FACSDiva software (BD Biosciences). Forward and side scatter gating strategy was used in order to select the desired lymphocyte subpopulations, and their relative percentages were analyzed for each cell subset.

Allogeneic and Autologous Cell Proliferation Assays

For the determination of the reactivity of PBMC from each donor against TT, 2 × 10 ⁵ PBMC were plated in 96-well round bottom plates at day 0 of each culture in supplemented RPMI medium containing 0. 1 µg/ml TT. As control conditions, the same number of cells was cultured with either supplemented RPMI medium only (negative control) or 50 ng/ml phorbol 12-myristate-13-acetate (PMA) and 500 ng/ml ionomycin (positive control). Ten replicates were performed for the negative control and the condition of analysis, and six replicates for the positive control. Cells were then cultured for 5 days at 37°C in a 5% CO ₂ atmosphere. Afterwards, 1 µCi [ ³ H]-thymidine (PerkinElmer, Waltham, MA, USA) was added to each well, and the plate was incubated for further 18 h under the same conditions. Cells were then collected using a HARVESTER96 2M cell harvester (Tomtec Inc, Hamdem, CT, USA) and read on a 1450 MicroBeta TriLux liquid scintillation counter (Wallac, Turku, Finland). Donors were considered positive for TT reactivity when the counts per minute (cpm) of at least five replicates from the condition of analysis were over the mean plus two times the standard deviation (SD) of the negative control.

For the isolation of allogeneic PBMC, whole blood samples of different healthy donors were processed by ficoll-hypaque density gradient separation. Cells were washed twice, and afterwards, their absolute number and viability was determined as shown above. Subsequently, 10 ⁵ either allogeneic or autologous viable PBMC were co-cultured with 5, 000 either iDC, mDC, mDC-TT, vitD3-tolDC or vitD3-tolDC-TT (1: 20 ratio) in 96-well round bottom plates, in a total volume of 200 µl of supplemented RPMI medium. Again, as negative and positive controls, either supplemented RPMI medium or a mix of 50 ng/ml PMA and 500 ng/ml ionomycin was used, respectively. Six replicates of each condition were performed. Cells were then plated for 4 days at 37°C in a 5% CO ₂ atmosphere, and afterwards, 1 µCi [ ³ H]-thymidine was added to each well, and the plates were incubated, harvested, and read as described above.

Cytokine and Soluble Protein Production

The production of granzyme B (GZMB), as well as of IL-1β, IL-6, IL-10, IFN- γ , IL-12p70 and TNF-α cytokines, was quantified in the supernatants of mDC-TT and vitD3-tolDC-TT with autologous PBMC co-cultures, using the Human Soluble Protein CBA Flex Set (BD biosciences) according to manufacturer’s instructions. Samples were acquired on an LSR Fortessa flow cytometer (BD Biosciences), and the results were analyzed using FACSDiva software. The production of TGF- β was determined using the Human/Mouse TGF beta 1 Uncoated ELISA kit (Invitrogen) in 100 µl of the co-culture supernatants after sample activation with HCl 1N, following the manufacturer’s instructions. The optical density of each well was measured at λ = 450 nm, and the optical density at λ = 570 nm was then subtracted as background signal, using a Varioskan Flash Multimode Reader (Thermo Fisher Scientific, Waltham, MA, USA).

RNA Extraction and RNA-seq Analysis

Total RNA of autologous mDC-Tcell and vitD3-Tcell samples was isolated using the automated Maxwell 16 LEV simplyRNA Purification Kit (Promega Biotech, Madison, WI, USA), including a DNAse I digestion step, according to manufacturer’s instructions. Samples were quantified using a Nanodrop ND-1000 spectrophotometer (Thermo Fisher Scientific), and subsequently stored at −80°C in RNAse-free tubes. RNA integrity number (RIN) was determined in an Agilent BioAnalyzer with the RNA6000 Pico assay (Agilent Technologies, Santa Clara, CA, USA). Afterwards, the sequencing libraries were prepared using the TruSeq Stranded Total RNA Sample Preparation kit (Illumina, San Diego, CA, USA) with 200 ng of total RNA per sample as input. Paired-end sequencing (2 × 75 bp) was then performed on a HiSeq-2500 instrument (Illumina). Reads were quality trimmed and adapters removed using Trimmomatic V0. 30. TopHat software v2. 1. 0 was used to map RNA-seq reads to the human reference genome (Ensembl release 78) ( 26 ). FeatureCounts function was used to assign reads to genomic features focusing on RNA biotypes. A matrix with summarized raw counts of reads assigned through mapping to high confidence protein coding genes only was generated (“ golden” annotation label), corresponding to stable and unlikely to change transcripts from the Consensus CDS (CCDS) Project. Data exploration results from hierarchical clustering and principal component analysis (PCA) in R software were used to exclude any outliers and assess sample similarities based on global gene expression patterns, and to guide the modeling design to be used for subsequent analyses.

Differential Gene Expression Analysis

Transcriptional changes at the gene level between mDC-Tcell and vitD3-Tcell were assessed using the Bioconductor DeSeq2 package in R ( 27 ). A paired sample comparison design, factoring in inter-individual differences, was applied. The results were considered statistically significant with an adjusted p-value (padj) < 0. 05. We set a 20% fold change (FC) cutoff as the threshold for relevant biological effects (| FC| > 1. 2).

Gene Ontology Enrichment Analysis

Unranked lists of the significant differentially upregulated (FC > 1. 2 and padj < 0. 05) and downregulated genes (FC < −1. 2 and padj < 0. 05) were tested for enrichment in Gene Ontology (GO) functional categories using the GOrilla web tool, applying the default settings for comparison to the background list of genes found in the dataset ( 28 ). We tested for enrichment in three types of GO categories: “ biological process” (GOPROCESS), “ molecular function” (GOFUNCTION), and “ cellular component” (GOCOMPONENT). Enrichment score (ES) was defined as ES = (b/n)/(B/N), where “ N” is the total number of genes in the background list, “ B” is the total number of genes in N associated with a specific GO term, “ n” is the number of differentially expressed genes being tested for enrichment and “ b” is the number of n intersecting with B. Enrichment p-value is computed according to the hypergeometric (HG) model. False discovery rate (FDR) q-value is the Benjamini and Hochberg multiple testing correction adjusted p-value. For the i ^th term (ranked according to p-value), the FDR q-value is the p-value multiplied by the number of GO terms assessed and divided by i.

Statistical Analysis

All the statistical analyses were performed with either parametric or non-parametric tests depending on the normality of each compared data set, as determined by the D’Agostino & Pearson test using Prism 6. 0 software (GraphPad, La Jolla, CA, USA). For multiple comparisons, either the non-parametric Friedman test with Dunn’s correction or the one-way ANOVA test with Geisser–Greenhouse correction were used, and analogously, either paired t tests or Wilcoxon tests for the comparisons between two groups if they were normally distributed or not, respectively. Results were expressed as mean ± SD, unless noted otherwise, and they were considered statistically significant when p < 0. 05.

Results

Functional and Phenotypical Characteristics of TT-Loaded vitD3-tolDC

Monocytes from 16 healthy donor samples were isolated (94. 4 ± 2. 8% purity) with viabilities of CD14 ⁺ cells above 95%. After their differentiation into DC, with or without exposition to TT, cells were harvested and their purity, viability, and phenotype were determined by flow cytometry, as shown previously ( 25 ). In all cases, purity was > 90%, as determined by the percentage of CD11c ⁺ cells, with a mean viability of 94. 2 ± 3. 3%, which was not affected by the addition of TT ( Supplementary Table 2 ). The study of the phenotype of vitD3-tolDC-TT showed significant reductions in the surface expression of CD86 (77. 2 ± 8. 7%) and HLA-DR (79. 5 ± 7. 7%) compared to mDC, but more importantly, evidenced that the exposure of DC to TT on day 3 of the culture did not have an effect per se over the expression of these molecules, neither in vitD3-tolDC nor in mDC, since there were no relevant differences on the percentages of reduction ( Supplementary Figure 1 ). The same could be observed regarding the functionality of these cells. On the one hand, as also shown in Supplementary Figure 1 , both vitD3-tolDC and vitD3-tolDC-TT exhibited a similar and strongly reduced induction of allogeneic proliferation compared to mDC (vitD3-tolDC: 50. 6 ± 30. 7, p < 0. 001; vitD3-tolDC-TT: 49. 2 ± 36. 7, p = 0. 001). On the other hand, there were no statistically significant differences in the mean induction of allogeneic proliferation induced by mDC-TT compared to mDC (p = 0. 916). Altogether, our results evidence that vitD3-tolDC-TT show the same tolerogenic properties as vitD3-tolDC, thus demonstrating that loading these cells with TT does not affect their phenotype nor their functionality.

VitD3-tolDC-TT Induce an Antigen-Specific Response Over Autologous PBMC

In order to test the antigen-specific functionality of vitD3-tolDC-TT in an autologous setup, we assessed the baseline reactivity of each donor against the TT itself to measure their potential to respond under these conditions. As shown in Figure 1A , we were able to assess the TT reactivity in all of our healthy donors, but only nine of them resulted positive, according to the criteria described in the Material and Methods section—when the mean proliferation of at least 5 out of 10 replicates was over the mean plus two times the SD of the control condition—and reaching statistical significance (p < 0. 05).

FIGURE 1

PBMC reactivity and antigen-specific induction of autologous proliferation mediated by DC against tetanus toxin.(A)Induction of proliferation of PBMC without stimuli (C−) and against tetanus toxin (TT) after 5 days of culture (n = 16). Data presented as counts per minute (cpm), measured as tritiated thymidine incorporation after 18 h. Ten replicated measurements of each condition were performed.(B)Induction of antigen-specific autologous proliferation against TT mediated by immature DC (iDC), mature DC (mDC), TT-loaded mDC (mDC-TT), vitamin D3-induced tolerogenic DC (vitD3-tolDC) and TT-loaded vitD3-tolDC (vitD3-tolDC-TT), as well as a negative control (C−), without any stimuli (n = 16) and(C)comparison of autologous antigen-specific proliferation against TT mediated by mDC-TT and vitD3-tolDC-TT on each donor. Data presented as relative percentage of induced proliferation compared to mDC-TT, measured as tritiated thymidine incorporation after 18 h. Six replicated measurements of each condition were performed. Error bars corresponding to SEM. ns, not significant; *p < 0. 05. One-way ANOVA test with Geisser–Greenhouse correction or paired t test.

Subsequently, we analyzed the capability of our cells to induce proliferation over autologous PBMC. As shown in Figure 1B , a significant proliferation was only induced by mDC-TT, as evidenced by the statistically significant differences observed with the remaining conditions. Specifically, reductions of a 38. 4 ± 44. 3% (p = 0. 020), a 40. 0 ± 21. 0% (p < 0. 001), a 56. 9 ± 19. 2% (p < 0. 001) and a 37. 3 ± 17. 4% (p < 0. 001) were observed in iDC, mDC, vitD3-tolDC and vitD3-tolDC-TT, respectively, compared to mDC-TT. Our results therefore evidence that autologous proliferation is only primed if an antigenic peptide is presented by an immunogenic DC condition, such as mDC-TT, confirming the antigen-specific modulation developed by our cells. Furthermore, reduced autologous proliferation mediated by vitD3-tolDC-TT was observed in all donors ( Figure 1C ).

VitD3-tolDC-TT Drive a Reduction of T _H 1 CD4 ⁺ Cell Subpopulations

Once determined that an antigen-specific modulation was established by TT-loaded DC, we studied which changes were being induced over the autologous T lymphocytes. Therefore, we characterized the phenotype of CD3 ⁺ CD4 ⁺ and CD3 ⁺ CD8 ⁺ cells using an exhaustive multiparametric flow cytometry panel, described in previous studies ( 29 ). First, our results evidenced a reduction in the prevalence of activated T CD4 ⁺ cells, determined by HLA-DR and/or CD38 staining, in vitD3-Tcell compared to mDC-Tcell (Activated CD4 ⁺ _mDC-Tcell : 23. 57 ± 15. 81 vs Activated CD4 ⁺ _vitD3-Tcell : 18. 52 ± 14. 16; p = 0. 002). The same effect was observed over T CD8 ⁺ cells (Activated CD8 ⁺ _mDC-Tcell : 15. 94 ± 12. 48 vs Activated CD8 ⁺ _vitD3-Tcell : 11. 33 ± 9. 81; p = 0. 002). Furthermore, we found a reduction in the relative percentages of CD4 ⁺ T _H 1 Central Memory (CM) and Effector Memory (EM) subpopulations in vitD3-T cell (T _H 1 CM _mDC-Tcell : 33. 98 ± 6. 44 vs T _H 1 CM _vitD3-Tcell : 30. 23 ± 7. 48; p = 0. 013; T _H 1 EM _mDC-Tcell : 44. 46 ± 8. 72 vs T _H 1 EM _vitD3-Tcell : 40. 95 ± 8. 08; p = 0. 001). All these results are shown in Figures 2A, B . Thus, our data suggest that vitD3-tolDC-TT are driving an antigen-specific switch towards a more anti-inflammatory—or less T _H 1-like—profile over T CD4 ⁺ lymphocytes. We could not detect any significant changes over any other T cell subpopulation, nor Treg nor Tr1 subpopulations (data not shown).

FIGURE 2

Phenotype and cytokine production of DC-co-cultured autologous PBMC.(A)Gating strategy.(B)Relative percentages of CD4 ⁺ CD38 ⁺ HLA-DR ⁺, CD8 ⁺ CD38 ⁺ HLA-DR ⁺, T _H 1 Central Memory and T _H 1 Effector Memory subpopulations on PBMC co-cultured with either tetanus toxin (TT)-loaded mature DC (mDC-TT) or TT-loaded vitamin D3-induced tolerogenic DC (vitD3-tolDC-TT) after 5 days of culture (n = 10). Data presented as the relative percentage of each lymphocyte subpopulation within its respective parent subpopulation, measured by multiparametric flow cytometry, in a box and whiskers representation. Error bars corresponding to the maximum and minimum values of each condition.(C)Analysis of the secretion of granzyme B (GZMB) (n = 12), IFN- γ (n = 12), IL-1β (n = 8), IL-6 (n = 8), IL-10 (n = 10), TGF- β (n = 7), and TNF-α (n = 8) in the supernatants of PBMC co-cultured with either mDC-TT or vitD3-tolDC-TT after 5 days of culture using either ELISA (TGF- β ) or cytometric bead array (GZMB, IFN- γ , IL-1β, IL-6, IL-10, and TNF-α) techniques. One ELISA experiment was performed for all samples, with duplicated measurements for each sample. One single cytometric bead array experiment was performed for the analysis of all the samples, with one single measurement for each sample. Sample sizes vary due to some measurements being negative or under the detection limit of the technique. Error bars corresponding to SEM. ns, not significant; * p < 0. 05. Paired t test or Wilcoxon test.

Next, we analyzed the cytokine secretion profile present in the autologous co-cultures of mDC-TT and vitD3-tolDC-TT. Our results, as shown in Figure 2C , evidenced a statistically significant increased secretion of the cytokine IL-6 and lower levels of IFN- γ in the co-culture of autologous PBMC with vitD3-tolDC-TT compared to mDC-TT (IL-6 _mDC-TT : 61. 4 ± 84. 3 pg/ml vs IL-6 _{vitD3-tolDC-TT} : 77. 7 ± 94. 5 pg/ml; p = 0. 039; and IFN- γ _mDC-TT : 3. 3 ± 2. 9 pg/ml vs IFN-γ _{vitD3-tolDC-TT} : 2. 2 ± 2. 8 pg/ml; p = 0. 002). Therefore, the reduction in the production of IFN- γ , combined with the increase of IL-6, again suggest a reduction of the T _H 1-like cytokine profile, in line with the phenotype results. No statistically significant changes could be found in the production of GZMB, IL-1β, IL-10, TGF-β nor TNF-α.

VitD3-tolDC Induce a General Transcriptomic Repression Over T CD4 ⁺ Cells

For all the 16 donors, at least 700, 000 CD3 ⁺ CD4 ⁺ cells in both conditions (mDC-Tcell and vitD3-Tcell) were successfully isolated by flow cytometry cell sorting. The gating strategy is shown in Supplementary Figure 1 . Afterwards, we extracted their RNA and selected 10 donors that showed sufficient nucleic acid concentration and integrity for the RNA-seq analysis (RIN > 7) in both mDC-Tcell and vitD3-Tcell conditions. Consequently, donors HD4, HD5, HD9, HD10, HD11, and HD12 were discarded from downstream studies.

After processing the samples through the RNA-seq analysis, 39% of total reads could be assigned to different known RNA classes ( Supplementary Figure 2A ), and out of them, around 47% of these assigned reads could be related to protein coding genes ( Supplementary Figure 2B ). Interestingly, the hierarchical clustering analysis revealed that our samples tended to cluster by individual rather than by treatment ( Figure 3A ), but also that there is a consistent pattern by which vitD3-Tcell samples ranked higher on both axes from the PCA ( Figure 3B ). These results led to choose a paired comparative analysis approach for the differential expression analysis.

FIGURE 3

Exploratory analysis of the RNA-seq study of T CD4+ cells co-cultured with autologous antigen-specific DC.(A)Hierarchical clustering analysis by gene expression of the 20 samples of CD4 ⁺ T cells of the RNA-seq study.(B)Representation of the first two principal component analysis (PCA) components at the gene level. Each color corresponds to a different sample of CD4+ T cells, as depicted in the legend, and the co-culture condition of each sample is indicated with a prefix, either “ m” for tetanus toxin (TT)-loaded mature DC or “ v” for TT-loaded vitamin D3-induced tolerogenic DC.(C)Volcano plot showing the significant differentially expressed genes. Axis is the log2 fold change. Color code: green, significantly regulated genes (padj < 0. 05; | FC| > 1. 2) considered in the Gene Ontology enrichment analysis; orange, genes with | FC| > 1. 2 below the significance threshold; red, genes with padj < 0. 05 below the relevant fold change cutoff.

After the subsequent filtering process described in the Material and Methods section, a total of 16, 333 protein coding genes with detectable reads were tested for differential expression. Among all of them, 546 genes showed a statistically significant change in their expression (adjusted p < 0. 05) in vitD3-Tcell compared to mDC-Tcell, and only 373 also presented an absolute value of FC superior to 1. 20 (| FC _{vitD3-Tcell vs mDC-Tcell} | > 1. 20). While only 29 of these genes were up-modulated in vitD3-Tcell compared to mDC-Tcell, the majority of them, 344 genes, were down-modulated, indicating a strong transcriptomic repression induced by vitD3-tolDC-TT over these cells ( Figure 3C ).

T CD4 ⁺ Cells Selectively Undergo a Strong Functional and Immune-Related Transcriptomic Down-Modulation Upon Interaction With vitD3-tolDC

When we studied those differentially expressed genes that appeared up-modulated (FC _{vitD3-Tcell vs mDC-Tcell} > 1. 20) in our analysis ( Table 1 ), we did not find many relevant or immune-related genes. Specifically, 18 of these 29 genes did not have any GO annotation, and among the rest, we could only find the genes encoding the JUNB and SCML1 transcription factors and several other genes encoding different molecule transporters ( ABCC2 and SLC10A1 ), G-protein modulators ( GRTP1 and RASA4 ) and kinases ( AK5 and CKMT2 ).

TABLE 1

Up-modulated genes in vitD3-Tcell compared to mDC–Tcell.

However, as mentioned above, the study of the down-modulated genes (FC _{vitD3-Tcell vs mDC-Tcell} < 1. 20) yielded many more results. Within the 50 most down-regulated results (FC _{vitD3-Tcell vs mDC-Tcell} < 1. 78) we found several genes encoding proteins involved in the immune response ( CCL17 , CCL22 , EBI3 , IL13 and LIF ), antigen presentation ( HLA-DQA1 , HLA-DQA2 and HLA-DRB5 ) and microtubule binding ( KIF4A , KIF15 , KIF18B and KIFC1 ), among others ( Table 2 ). Furthermore, when we analyzed the whole list, we could find not only several more genes included in these categories, but also many other genes encoding proteins related to cytoskeleton and cell adhesion ( ARPC1B , CAPG , CTNNA1 , LGALS1 , MYL6B , LGALS9 or SDC4 ), actin related functionalities ( ACTB , ACTG1 , PARVB or TPM4 ), G-proteins and modulators ( GBP2 , GBP4 , GNA15 , GNG4 , IQGAP3 , MYO1G , MYO1E or SRGAP3 ), nucleic acid binding ( ASF1B , DEPDC1 , EXO1 or FEN1 ), histones ( HIST1H2BL , HIST2H2BF or, HIST1H4H ), the pro-inflammatory transcription factor STAT1 and other pro-inflammatory mediators ( TNFSF4 ), different kinase activators and modulators, proteases and protease inhibitors, oxydases, oxygenases, transferases and many other metabolic mediators. The whole list is shown in Supplementary Table 3 . Altogether, these results indicate that vitD3-tolDC-TT mediate a strong down-modulation of metabolic and immune-related functions over vitD3-Tcell.

TABLE 2

Top 50 down-modulated genes in vitD3-Tcell compared to mDC-Tcell.

VitD3-Tcell Present Decreased Cell Cycle and Mitotic Activity

In this regard, the GO enrichment analysis further supported the results observed in the differential gene expression (DGE) study. Thus, first, the enrichment analysis produced a total of 482 protein sets and pathways with p < 0. 001, four of them up-modulated—although none of them showed an FDR value below 0. 25—and the remaining 478 down-modulated ( Supplementary Table 4 ). We further filtered the results to analyze only those GO terms that presented a much more significant enrichment (p < 10 ^-9 ). This process left us with a total of 66 down-modulated GO terms, but none up-modulated. These 66 elements, ordered by decreasing ES, are shown in Table 3 . Interestingly, among the most significantly enriched down-modulated pathways, we found several GO annotations referring to immune-related functionality (for instance Interferon-Gamma-Mediated Signaling Pathway , ES: 10. 65; Cytokine-Mediated Signaling Pathway , ES: 5. 09; or Immune Response , ES: 3. 03), class II-related antigen presentation (like MHC Class II Protein Complex , ES: 30. 23; Antigen Processing And Presentation Of Exogenous Peptide Antigen Via MHC Class II , ES: 8. 08; or Antigen Processing And Presentation Of Exogenous Peptide Antigen , ES: 6. 74), cell response to different stimuli ( Cell Surface Receptor Signaling Pathway , ES: 2. 21; Response To Stress , ES: 1. 93; or Cellular Response To Stimulus , ES: 1. 80) and, specially, to cell cycle and mitotic division (for instance Condensed Chromosome Outer Kinetochore , ES: 21. 98; Mitotic Spindle Organization , ES: 8. 79; Microtubule Cytoskeleton Organization Involved In Mitosis , ES: 7. 48; or Cell Cycle Checkpoint , ES: 6. 05). Our results, in line with the DGE analysis, would suggest that vitD3-Tcell are undergoing a process of transcriptomic down-modulation leading to reduced immune-related, metabolic and proliferative functionalities.

TABLE 3

Most significantly down-regulated Gene Ontology terms in vitD3-Tcell compared to mDC-Tcell.

Discussion

In this study we analyzed the specific effect of vitD3-tolDC over autologous CD4 ⁺ T cells. Thus, we switched the attention from the study of tolDC themselves—widely studied so far—to focus on the study of the functional effect that these cells develop over T cells upon their interaction. In homeostatic conditions, either depletion, inactivation and/or induction of anergy is often induced on T cells due to a lack of one or more of the three immunogenic activation signals. This causes T cells to become hyporesponsive or to die ( 20 , 30 ). However, in the case of autoimmunity, where T cells are already activated and developing an immunogenic response, an antigen-specific process of tolerance induction is required. In this regard, previous in vivo studies with vitD3-tolDC in the EAE model showed that an antigen-specific setup—and therefore an active process—is required, provided that a beneficial effect of this therapy was only observed when vitD3-tolDC were pulsed with the adequate immunogenic peptide ( 23 , 24 ). Consequently, we developed an experimental model for the generation of autologous antigen-specific vitD3-tolDC and T cells from healthy donors, using an immunogenic peptide presented via class II MHC with the aim to reproduce antigen presentation in the context of CD4 ⁺-mediated autoimmune diseases. In this regard, we selected TT for its compliance with this feature—since the vaccination against TT is included in European health systems—which also allowed us to use healthy donors instead of patients of a determined autoimmune disease, thus eliminating disease-conditioned immune variations.

All in all, our approach aimed to be as versatile as possible, and to serve as a preliminary study for future research, with the idea that the immunogenic peptide/s might be replaced depending on the disease of interest, as well as, of course, using patient cells instead of healthy donors. On the one hand, for autoimmune diseases with identified autoantigens, this decision would be trivial. On the other hand, for those conditions in which autoimmune antigens are yet to be identified, the experimental design should be considered case by case. For instance, using autologous synovial fluid as a source for autoantigens to load tolDC has provided promising results in rheumatoid arthritis ( 31 ), so analog workaround solutions could be taken into account.

After validating our experimental setup—meaning that vitD3-tolDC-TT were able to induce an antigen-specific response—we focused on the study of the actual phenotypic, functional and transcriptomic modulations induced by vitD3-tolDC. First, the analysis of the phenotype of T CD4 ⁺ cells evidenced that their interaction with vitD3-tolDC-TT caused a relative reduction in the activation of these cells. More importantly, a switch in the immune response of these cells towards a more immunoregulatory profile was induced, with a reduction in the prevalence of T _H 1 memory subpopulations. These results were further supported by the decrease of IFN- γ production in the autologous co-culture supernatants, consequently supporting that vitD3-tolDC were inducing a switch towards a more anti-inflammatory immune profile, in line with previous findings in the literature ( 16 – 18 ).

When we deepened into the analysis of the vitD3-tolDC-mediated transcriptomic profile of T cells, we observed several genes and GO terms regulated in line with the abovementioned phenotypical and functional switch towards a less activated an more immunoregulatory profile; for instance, a down-modulation of STAT1 gene and the interferon-gamma-mediated signaling pathway was observed, which on the other hand supported the robustness of our RNA-seq study. Beyond this, the results pointed towards a generalized down-modulation of the transcriptomic profile of vitD3-Tcell that could either respond to an induction of T CD4 ⁺ cell hyporesponsiveness or even to a process of clonal deletion. On the one hand, the down-modulation of genes and pathways involved in crucial cellular mechanisms—in particular those related to cell proliferation, mitosis, cell cycle and response to immune stimuli, some of them never reported before—could be explained by both of these processes. However, on the other hand, the lack of induction of cell death and apoptotic-related pathways makes clonal deletion very unlikely to be happening. Therefore, our results suggest that the antigen-specific interaction of vitD3-tolDC with autologous T CD4 ⁺ cells is mediating, in fact, an induction of hyporesponsiveness over these cells. Furthermore, previous studies from our group already pointed in this direction ( 16 ), a biological situation that, potentially, might lead to the abrogation of an autoimmune immunogenic response in patients. Moreover, our current results provide evidence that these modulations are taking place at the transcriptomic level in T CD4 ⁺ cells, indicating that the antigen-specific modulation induced by vitD3-tolDC is deeper than expected and, in consequence, probably also long-lasting.

Unfortunately, the lack of strongly up-regulated genes among the protein-coding RNA transcripts did not allow us to point towards many clear candidate biomarkers that might become indicators of the response of T cells upon their interaction with vitD3-tolDC to monitor patients in clinical trials. One of the most relevant exceptions came given by JUNB gene, encoding a member of the AP-1 family of transcription factors. In experimental models, this gene has been reported to be crucial in maintaining Treg suppressive function ( 32 ), although it is also apparently involved in the induction and maintenance of IL-23-related pathogenicity of T _H 17 cells ( 33 , 34 ). However, neither of these functionalities seems to fit in our model based on our results, since neither T _H 17 nor Treg induction was evidenced. Consequently, it would be interesting to elucidate the specific role of JunB in our experimental setting, and to what extent small changes in its expression can actually be sufficient or not for the induction and maintenance of immune tolerance. Furthermore, if this hypothesis proves to be valid, either JUNB and/or other related genes might also constitute potential biomarkers of response to vitD3-tolDC treatment in the clinic. It is also worth noting that we may have overlooked other potential biomarkers that might be found among the non-protein-coding and alternate splicing RNA transcripts. Although this possibility, if true, would have a limited functional value in our experimental model, it could be addressed in future studies.

In addition, our results did not allow us to reach any conclusion regarding a potential induction of anergy, and, as discussed above, they also rule out any kind of Treg or Tr1 response. Even though our previous in vivo experiments with the murine EAE model pointed towards an induction of Treg mediated by vitD3-tolDC ( 23 , 24 ), we have observed that, at least in this experimental setting, this is not the case with human cells. These results are in line with what previous studies from both our group and other authors have already reported ( 13 , 16 , 35 ), although there seems to be some controversy ( 36 , 37 ). However, it is also worth mentioning that in these reports, Treg induction was only observed after two rounds of stimulation of T cells, which might explain why we have not detected it. This is definitely something to be taken into account, since Treg induction is, undoubtedly, one of the main mechanisms for the induction of immune tolerance of tolDC and other antigen presenting cell approaches ( 38 , 39 ). Indeed, Treg themselves, when expanded in vitro , present a huge therapeutic potential as a cell therapy for autoimmune diseases in humans ( 40 ). Consequently, the transcriptomic study of vitD3-tolDC-induced Treg should probably be addressed separately in future studies, since two rounds of T cell stimulation might have masked some of the results that we have reported here.

Our current study presents some limitations. First, since we focused on the study of CD4 ⁺ T cells alone, we were naturally omitting the potential modulation that vitD3-tolDC might be mediating through other subpopulations, such as regulatory B cells or regulatory NK cells. Furthermore, the election of the timepoint for the RNA-seq analysis intrinsically establishes another limitation, which is the status of the transcriptomic profile at different timepoints of the co-culture. However, our selection came based on the phenotypical and functional results shown in the study, which evidence that, by day 5 of the co-culture, there is a significant and differential modulation mediated by vitD3-tolDC over T CD4 ⁺ cells. Consequently, even though it is true that other timepoints might provide valuable additional information, we think that our election provided the best compromise, and a full time-course characterization of the antigen-specific transcriptomic changes induced by vitD3-tolDC will be addressed in future studies. On the other hand, we cannot fully discard the presence of non-antigen-specific CD4 ⁺ cells by the time the cell sorting was performed. However, even with a residual amount of non-antigen-specific T cells, the obtained results were consistent not only within the different techniques, but also with the literature, supporting our findings.

In conclusion, our results evidence that vitD3-tolDC are inducing a strong antigen-specific transcriptomic down-modulation over autologous T CD4 ⁺ cells, with a reduced ability to respond to immune- and non-immune-related stimuli. Consequently, it constitutes one of the first attempts to understand the changes that T cells are undergoing at the transcriptomic level upon an antigen-specific interaction with a tolerogenic cell product, such as vitD3-tolDC. In that regard, we identified several specific genes and pathways selectively down-modulated, as well as the induction of JUNB , which might constitute a putative biomarker of the modulation mediated by vitD3-tolDC over CD4 ⁺ T cells and, consequently, a potential biomarker to monitor the effect of vitD3-tolDC after their administration to patients. Therefore, the results presented in this article allowed us to better understand the process of T cell hyporesponsiveness at the molecular level and, more importantly, to set the path for future studies to fully elucidate the specific processes that are taking place in one of the most important mechanisms that the promising tolDC-based therapies can trigger in order to restore tolerance in autoimmune diseases.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www. ncbi. nlm. nih. gov/geo/ , GSE128816.

Ethics Statement

This study was reviewed and approved by the Germans Trias i Pujol Hospital Ethical Committee. The participants provided their written informed consent to participate in this study.

Author Contributions

EM-C, JN-B and MM conceived the experiments. JN-B and MM performed the cell cultures and the cell sortings, and prepared the samples for the RNA-seq analysis. JN-B and MM analyzed the results. EM-C, JN-B and MM interpreted the results. JN-B wrote the manuscript. AT-S, BQ-S, CR-T, EM-C, JN-B, and MM reviewed the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was partially supported by projects PI14/01175, PI16/01737 and PI17/01521, integrated in the Plan Nacional de I+D+I and co-supported by the ISCIII- Subdirección General de Evaluación and the Fondo Europeo de Desarrollo Regional (FEDER), and by grant 140191 from IWT-TBM (Belgium). JN-B was beneficiary of a FI-DGR PhD contract from AGAUR ( Agència de Gestió d’Ajuts Universitaris i de Recerca ), supported by the Government of Catalonia (02/2015-01/2018). The contracts of JN-B and MM have been partially supported by grant 140191 from IWT-TBM. The authors are members of a consolidated group recognized by AGAUR (2017 SGR 103), supported by the Government of Catalonia. This work has been supported by positive discussion through A FACTT network (COST Action BM1305: www. afactt. eu ) and project 779316–ReSToRe ( www. h2020restore. eu ). COST and ReSToRe are supported by the EU Framework Programme Horizon 2020.

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.

Acknowledgments

We want to thank Marco A. Fernández-Sanmartín and Gerard Requena, from the IGTP Flow Cytometry Platform, Edurne Pedrosa, Verónica Guirao and Ingrid Liarte, from the IGTP Biobank, Lauro Sumoy and Raquel Pluvinet, from the IGTP High Content Genomics and Bioinformatics Unit, and, specifically, Jorge Carrascal and all the members of the IGTP Immunology Division for their support and technical assistance.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www. frontiersin. org/articles/10. 3389/fimmu. 2020. 599623/full#supplementary-material

Supplementary Figure 1 | Phenotype and functionality of DC.(A)Gating strategy.(B)Surface expression of CD86 (FITC) and HLA-DR (V500) on immature DC (iDC), mature DC (mDC), tetanus toxin (TT)-loaded mDC (mDC-TT), vitamin D3-induced tolerogenic DC (vitD3-tolDC) and TT-loaded vitD3-tolDC (vitD3-tolDC-TT) (n = 8). Data presented as relative percentage of median fluorescence intensity (MFI) normalized versus mDC, measured by flow cytometry. In all cases, DC populations were CD11c ⁺(C)Allogeneic proliferation of PBMC from healthy donors co-cultured with either iDC, mDC, mDC-TT, vitD3-tolDC or vitD3-tolDC-TT (n = 14) in a DC : PBMC 1: 20 ratio. Data presented as relative percentage of induced proliferation compared to mDC, measured as tritiated thymidine incorporation after 18 h. Six replicated measurements of each condition were performed. Error bars corresponding to SEM. ns = not significant; * p < 0. 05. Friedman test with Dunn’s correction.(D)Gating strategy for the sorting of autologous T CD4 ⁺ cells.

Supplementary Figure 2 | Assignment and functional classification of the total reads from the RNA-seq study of T CD4 ⁺ cells co-cultured with autologous antigen-specific DC.(A)Assignment of the total reads from the RNA-seq analysis to known RNA classes.(B)Classification of the assigned reads into known RNA functionalities.

Supplementary Table 3 | Up and down-modulated genes in vitD3-Tcell compared to mDC-Tcell.

Supplementary Table 4 | Differentially enriched Gene Ontology terms in vitD3-Tcell compared to mDC-Tcell.

References

1. Ten Brinke A, Hilkens CMU, Cools N, Geissler EK, Hutchinson JA, Lombardi G, et al. Clinical Use of Tolerogenic Dendritic Cells-Harmonization Approach in European Collaborative Effort. Mediators Inflammation (2015) 2015: 471719. doi: 10. 1155/2015/471719

2. Banchereau J, Briere F, Caux C, Davoust J, Lebecque S, Liu Y-J, et al. Immunobiology of Dendritic Cells. Annu Rev Immunol (2000) 18: 767–811. doi: 10. 1146/annurev. immunol. 18. 1. 767

3. Mueller DL. Mechanisms maintaining peripheral tolerance. Nat Immunol (2010) 11: 21–7. doi: 10. 1038/ni. 1817

4. Ganguly D, Haak S, Sisirak V, Reizis B. The role of dendritic cells in autoimmunity. Nat Rev Immunol (2013) 13: 566–77. doi: 10. 1038/nri3477

5. Hopp A-K, Rupp A, Lukacs-Kornek V. Self-Antigen Presentation by Dendritic Cells in Autoimmunity. Front Immunol (2014) 5: 55. doi: 10. 3389/fimmu. 2014. 00055

6. Schinnerling K, García-González P, Aguillón JC. Gene expression profiling of human monocyte-derived dendritic cells – searching for molecular regulators of tolerogenicity. Front Immunol (2015) 6: 528. doi: 10. 3389/fimmu. 2015. 00528

7. Navarro-Barriuso J, Mansilla MJ, Martínez-Cáceres EM. Searching for the Transcriptomic Signature of Immune Tolerance Induction—Biomarkers of Safety and Functionality for Tolerogenic Dendritic Cells and Regulatory Macrophages. Front Immunol (2018) 9: 2062. doi: 10. 3389/fimmu. 2018. 02062

8. Adorini L, Penna G, Giarratana N, Roncari A, Amuchastegui S, Daniel KC, et al. Dendritic cells as key targets for immunomodulation by Vitamin D receptor ligands. J Steroid Biochem Mol Biol (2004) 89–90: 437–41. doi: 10. 1016/j. jsbmb. 2004. 03. 013

9. Penna G, Roncari A, Amuchastegui S, Daniel KC, Berti E, Colonna M, et al. Expression of the inhibitory receptor ILT3 on dendritic cells is dispensable for induction of CD4+Foxp3+ regulatory T cells by 1, 25-dihydroxyvitamin D3. Blood (2005) 106: 3490–7. doi: 10. 1182/blood-2005-05-2044

10. Adorini L, Penna G. Induction of tolerogenic dendritic cells by vitamin D receptor agonists. Handb Exp Pharmacol (2009) 188: 251–73. doi: 10. 1007/978-3-540-71029-5_12

11. Pedersen AW, Holmstrøm K, Jensen SS, Fuchs D, Rasmussen S, Kvistborg P, et al. Phenotypic and functional markers for 1α, 25-dihydroxyvitamin D3-modified regulatory dendritic cells. Clin Exp Immunol (2009) 157: 48–59. doi: 10. 1111/j. 1365-2249. 2009. 03961. x

12. Széles L, Keresztes G, Töröcsik D, Balajthy Z, Krenács L, Póliska S, et al. 1, 25-dihydroxyvitamin D3 is an autonomous regulator of the transcriptional changes leading to a tolerogenic dendritic cell phenotype. J Immunol Baltim Md 1950 (2009) 182: 2074–83. doi: 10. 4049/jimmunol. 0803345

13. Naranjo-Gómez M, Raïch-Regué D, Oñate C, Grau-López L, Ramo-Tello C, Pujol-Borrell R, et al. Comparative study of clinical grade human tolerogenic dendritic cells. J Transl Med (2011) 9: 89. doi: 10. 1186/1479-5876-9-89

14. Ferreira GB, Kleijwegt FS, Waelkens E, Lage K, Nikolic T, Hansen DA, et al. Differential protein pathways in 1, 25-dihydroxyvitamin d(3) and dexamethasone modulated tolerogenic human dendritic cells. J Proteome Res (2012) 11: 941–71. doi: 10. 1021/pr200724e

15. Raïch-Regué D, Naranjo-Gómez M, Grau-López L, Ramo C, Pujol-Borrell R, Martínez-Cáceres E, et al. Differential effects of monophosphoryl lipid A and cytokine cocktail as maturation stimuli of immunogenic and tolerogenic dendritic cells for immunotherapy. Vaccine (2012) 30: 378–87. doi: 10. 1016/j. vaccine. 2011. 10. 081

16. Raïch-Regué D, Grau-López L, Naranjo-Gómez M, Ramo-Tello C, Pujol-Borrell R, Martínez-Cáceres E, et al. Stable antigen-specific T-cell hyporesponsiveness induced by tolerogenic dendritic cells from multiple sclerosis patients. Eur J Immunol (2012) 42: 771–82. doi: 10. 1002/eji. 201141835

17. Ferreira GB, Vanherwegen A-S, Eelen G, Gutiérrez ACF, Van Lommel L, Marchal K, et al. Vitamin D3 Induces Tolerance in Human Dendritic Cells by Activation of Intracellular Metabolic Pathways. Cell Rep (2015) 10: 711–25. doi: 10. 1016/j. celrep. 2015. 01. 013

18. Malaguarnera L, Marsullo A, Zorena K, Musumeci G, Di Rosa M. Vitamin D3 regulates LAMP3 expression in monocyte derived dendritic cells. Cell Immunol (2017) 311: 13–21. doi: 10. 1016/j. cellimm. 2016. 09. 013

19. Xing N, Maldonado ML L, Bachman LA, McKean DJ, Kumar R, Griffin MD. Distinctive dendritic cell modulation by vitamin D(3) and glucocorticoid pathways. Biochem Biophys Res Commun (2002) 297: 645–52. doi: 10. 1016/S0006-291X(02)02262-3

20. Maldonado RA, von Andrian UH. Chapter 4 – How Tolerogenic Dendritic Cells Induce Regulatory T Cells. In: Alt FW, Austen KF, Honjo T, Melchers F, Uhr JW, Unanue ER, editors. Advances in Immunology . USA: Academic Press (2010). p. 111–65. doi: 10. 1016/B978-0-12-380995-7. 00004-5

21. Maggi J, Schinnerling K, Pesce B, Hilkens CM, Catalán D, Aguillón JC. Dexamethasone and Monophosphoryl Lipid A-Modulated Dendritic Cells Promote Antigen-Specific Tolerogenic Properties on Naive and Memory CD4+ T Cells. Front Immunol (2016) 7: 359. doi: 10. 3389/fimmu. 2016. 00359

22. Navarro-Barriuso J, Mansilla MJ, Naranjo-Gómez M, Sánchez-Pla A, Quirant-Sánchez B, Teniente-Serra A, et al. Comparative transcriptomic profile of tolerogenic dendritic cells differentiated with vitamin D3, dexamethasone and rapamycin. Sci Rep (2018) 8: 14985. doi: 10. 1038/s41598-018-33248-7

23. Mansilla MJ, Sellès-Moreno C, Fàbregas-Puig S, Amoedo J, Navarro-Barriuso J, Teniente-Serra A, et al. Beneficial effect of tolerogenic dendritic cells pulsed with MOG autoantigen in experimental autoimmune encephalomyelitis. CNS Neurosci Ther (2015) 21: 222–30. doi: 10. 1111/cns. 12342

24. Mansilla MJ, Contreras-Cardone R, Navarro-Barriuso J, Cools N, Berneman Z, Ramo-Tello C, et al. Cryopreserved vitamin D3-tolerogenic dendritic cells pulsed with autoantigens as a potential therapy for multiple sclerosis patients. J Neuroinflamm (2016) 13: 113. doi: 10. 1186/s12974-016-0584-9

25. Navarro-Barriuso J, Mansilla MJ, Quirant-Sánchez B, Ardiaca-Martínez A, Teniente-Serra A, Presas-Rodríguez S, et al. MAP7 and MUCL1 are biomarkers of vitamin D3-induced tolerogenic dendritic cells in multiple sclerosis patients. Front Immunol (2019) 10: 1251. doi: 10. 3389/fimmu. 2019. 01251

26. Zerbino DR, Achuthan P, Akanni W, Amode MR, Barrell D, Bhai J, et al. Ensembl 2018. Nucleic Acids Res (2018) 46: D754–61. doi: 10. 1093/nar/gkx1098

27. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol (2014) 15: 550. doi: 10. 1186/s13059-014-0550-8

28. Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z. GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinf (2009) 10: 48. doi: 10. 1186/1471-2105-10-48

29. Maecker HT, McCoy JP, Nussenblatt R. Standardizing immunophenotyping for the Human Immunology Project. Nat Rev Immunol (2012) 12: 191–200. doi: 10. 1038/nri3158

30. Alberts B, Johnson A, Lewis J, Raff M, Roberts K, Walter P. Lymphocytes and the Cellular Basis of Adaptive Immunity. Mol Biol Cell (2002). Available at: https://www. ncbi. nlm. nih. gov/books/NBK26921/ (Accessed February 21, 2019).

31. Bell GM, Anderson AE, Diboll J, Reece R, Eltherington O, Harry RA, et al. Autologous tolerogenic dendritic cells for rheumatoid and inflammatory arthritis. Ann Rheum Dis (2016) 76: 227–34. doi: 10. 1136/annrheumdis-2015-208456

32. Koizumi S, Sasaki D, Hsieh T-H, Taira N, Arakaki N, Yamasaki S, et al. JunB regulates homeostasis and suppressive functions of effector regulatory T cells. Nat Commun (2018) 9: 5344. doi: 10. 1038/s41467-018-07735-4

33. Hasan Z, Koizumi S, Sasaki D, Yamada H, Arakaki N, Fujihara Y, et al. JunB is essential for IL-23-dependent pathogenicity of Th17 cells. Nat Commun (2017) 8: 15628. doi: 10. 1038/ncomms15628

34. Carr TM, Wheaton JD, Houtz GM, Ciofani M. JunB promotes Th17 cell identity and restrains alternative CD4+ T-cell programs during inflammation. Nat Commun (2017) 8: 301. doi: 10. 1038/s41467-017-00380-3

35. Sochorová K, Budinský V, Rožková D, Tobiasová Z, Dusilová-Sulková S, Špíšek R, et al. Paricalcitol (19-nor-1, 25-dihydroxyvitamin D2) and calcitriol (1, 25-dihydroxyvitamin D3) exert potent immunomodulatory effects on dendritic cells and inhibit induction of antigen-specific T cells. Clin Immunol (2009) 133: 69–77. doi: 10. 1016/j. clim. 2009. 06. 011

36. Unger WWJ, Laban S, Kleijwegt FS, van der Slik AR, Roep BO. Induction of Treg by monocyte-derived DC modulated by vitamin D3 or dexamethasone: differential role for PD-L1. Eur J Immunol (2009) 39: 3147–59. doi: 10. 1002/eji. 200839103

37. Vanherwegen A-S, Eelen G, Ferreira GB, Ghesquière B, Cook DP, Nikolic T, et al. Vitamin D controls the capacity of human dendritic cells to induce functional regulatory T cells by regulation of glucose metabolism. J Steroid Biochem Mol Biol (2019) 187: 134–45. doi: 10. 1016/j. jsbmb. 2018. 11. 011

38. Boks MA, Kager-Groenland JR, Haasjes MSP, Zwaginga JJ, van Ham SM, ten Brinke A. IL-10-generated tolerogenic dendritic cells are optimal for functional regulatory T cell induction–a comparative study of human clinical-applicable DC. Clin Immunol Orlando Fla (2012) 142: 332–42. doi: 10. 1016/j. clim. 2011. 11. 011

39. Riquelme P, Haarer J, Kammler A, Walter L, Tomiuk S, Ahrens N, et al. TIGIT + iTregs elicited by human regulatory macrophages control T cell immunity. Nat Commun (2018) 9: 2858. doi: 10. 1038/s41467-018-05167-8

40. Marek-Trzonkowska N, Myśliwiec M, Dobyszuk A, Grabowska M, Derkowska I, Juścińska J, et al. Therapy of type 1 diabetes with CD4(+)CD25(high)CD127-regulatory T cells prolongs survival of pancreatic islets – results of one year follow-up. Clin Immunol Orlando Fla (2014) 153: 23–30. doi: 10. 1016/j. clim. 2014. 03. 016