VDR binding sites on Chr17q12 and q21.1 overlap with open chromatin signatures.

A, Significant disease and trait associations in 17q12-21.1 (y-axis), retrieved from the NHGRI EBI GWAS Catalog, with at least one significant SNP (x-axis) in Chr17q12-21.1. Colored bars denote autoimmune diseases mapped to Th1 (red bars) or Th2 (blue bars) driven immunity. B, University of California Santa Cruz (UCSC) Genome Browser tracks showing chromosomal location, common dbSNPs, VDR-, RXR-, CTCF– and IKZF3-binding sites and H3K27Ac marks present in the Chr17q12-21.1 locus. VDR binding sites overlapping with active regulatory elements are highlighted by colored boxes. C, Linkage disequilibrium (LD) between the significant SNPs in the Chr17q12-21.1 (magenta) and 21.2 (green) region, highlighting Th1/Th2 associated SNPs (red and blue, respectively), eQTLs in the VDR binding regions in Chr17q12 21.1 and 21.2 (yellow), and functional asthma SNPs from the literature (cyan).

VDR expression is elevated in Th2 cells and vitamin D regulated expression of genes on Chr1q12 and Chr17q21.

A, Scheme of the house dust mite (HDM) induced airway inflammation protocol. B, Flow cytometric analysis and cell frequencies of VDR expression in CD45+ CD4+ T cells in the respective groups (n=3-5). C, Flow cytometric analyses and cell frequencies of VDR expression in Th1 (Tbet+), Th2 (Gata3+) and Treg (Foxp3+) cells in the respective groups (n=4-8). Quantitative RT-PCR analysis of relative mRNA expression levels of genes encoded on Chr17q12-21.1 in (D) Th2 cultures of WT and vitamin D deficient mice and (E) control and calcitriol stimulated cultures. Quantitative RT-PCR analysis of relative mRNA expression levels of genes encoded on Chr17q21.2 in (F) Th2 cultures of WT and vitamin D deficient mice and (G) control and calcitriol stimulated cultures. Gene expression levels were normalized to the house keeping gene L32 and are expressed relative to the expression level in WT (D, F) or WT-vehicle (E, G) Th2 cells. (n ≥ 4 per group). Each symbol represents one mouse. Numbers in flow plots indicate percentages. Error bars indicate the s.e.m. Statistical tests: two-tailed Mann-Whitney U test (B, C). *P < 0.05, **P < 0.01, ***P < 0.001. Data summarize results from two independent experiments.

Vitamin D deficiency augments allergic airway inflammation.

A, Schematic presentation of the different mouse strains used. B, Representative H&E-stained lung sections from saline and HDM-exposed wild type (C57/BL6; WT), vitamin D deficient and VDR-/-mice (x 10 magnification). Black arrows point to the peri bronchial infiltrate. C, Representative flow plots assessing VDR expression in lung CD4+ T cells from the indicated groups. D, Total IgE levels detected in the sera of respective groups. E-G, Absolute numbers of (E) leukocytes, (F) eosinophils, and (G) lymphocytes in the airways of the respective groups. Data summarize results from two or three independent experiments with n≥8 per group. H, Frequencies of IL-13 producing Th2 in the lung. Data summarize results from two or three independent experiments with n>7 per group. I, Heat map showing differentially expressed genes between WT and vitamin D deficient Th2 cultures identified by RNA-Sequencing analysis. Analysis and visualization were performed with DESeq2. J, Overrepresentation analysis of Hallmark gene sets in the differentially expression gene set. The x-axis represents the gene ratio in the respective gene set. Dot sizes denote the number of genes in the respective pathway. The dot color indicates the FDR-adjusted p-value. Analysis and visualization were performed with the clusterProfiler package. K, Quantitative RT-PCR analysis of Il-2 mRNA expression in Th2 cultures from the respective groups (n=4 per group). Gene expression levels were normalized to the house keeping gene L32 and are expressed as fold-induction compared to WT Th2 cells. L, Flow cytometric analysis and cell frequencies of IL-2 production in WT and VDR-/-Th2 cultures (n=5 per group). Each symbol represents one mouse. Error bars indicate the s.e.m. Statistical significance was determined with: Mixed-effect analysis with Holm-Šidák’s post-hoc analysis (D –H), two-tailed Mann-Whitney U test (K, L). *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001. Data summarize results from two or three independent experiments with n>8 per group.

Vitamin D stimulation of polarizing Th2 cells suppresses the IL-2/Stat5 pathway.

A, Heat map of differentially expressed genes between vehicle (EtOH [% v/v]) and calcitriol treated WT Th2 cells. Analysis and visualization were performed with DESeq2. B, Gene set enrichment analysis (GSEA) of differentially expressed genes. Analysis and visualization were performed with the clusterProfiler package. C, Quantitative RT-PCR for Il-2, Stat5a and Stat5b mRNA expression in Th2 cells exposed to either vehicle or calcitriol. Gene expression levels were normalized to the house keeping gene L32 and are expressed as fold-induction compared to control cells. D, Flow cytometric analysis and mean fluorescence intensities of Itk and Il2rβ expression in indicated groups. E, Flow cytometric analysis and quantification of IL-13 expression in respective groups. F, Network analysis to visualize vitamin D induced molecular changes in the molecular interactome of calcitriol stimulated CD4+ Th2. Predefined set of cell membrane receptors were set as starting points (orange nodes) and annotated transcription factors were set as ending points (blue). All intermediate proteins on the PPI are visualized as grey nodes. Molecular interactions only overrepresented in calcitriol stimulated cells are shown as orange arrows. Each symbol in the bar graphs represents one individual sample and data summarize the results from two or three independent experiments with n≥4. Error bars indicate the s.e.m. Statistical test: two-tailed Student’s t-test (C, E) and two-tailed Mann-Whitney U-test (D). *P < 0.05, **P < 0.01, ***P < 0.001.

Prenatal vitamin D supplementation protects from asthma development.

A, H&E staining of representative lung sections of indicated experimental groups. Arrowheads point to peri bronchial infiltrates. B-E, Absolute numbers of total BAL (B) leukocytes, (C) eosinophils, (D) lymphocytes, and (E) lung IL-13+ Th2 cells (n≥8 per group). F, Flow cytometric analysis of Aiolos expression in differentiating Th2 cells at indicated time points. G, Flow cytometric analysis of Aiolos expression in control and calcitriol stimulated Th2 cells. H, Flow cytometric analysis and frequencies of IL-2+ cells in WT and Ikzf3-/-cultures post indicated treatment. Each symbol represents and independent sample (n≥9 per group). Flow cytometric analysis and quantification of VDR expression in lung CD4+ T cells of the respective groups. Each symbol represents one individual mouse and data summarize the results from two to three independent experiments with n>8 per group. Error bars indicate the s.e.m. Statistical significance was determined with Mixed-effect analysis with Holm-Šidák’s post-hoc analysis (B-E), two-tailed Mann-Whitney U-test (G) and One-way ANOVA with Holm-Šidák’s post-hoc analysis (H). *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001.

A, Significant disease and trait associations in 17q21.2 (y-axis), retrieved from the NHGRI-EBI GWAS Catalog, with at least one significant SNPs (x-axis) in Chr17q21.2. Colored bars denote autoimmune diseases mapped to Th1 (red bars) or Th2 (blue bars) driven immunity. B, University of California Santa Cruz (UCSC) Genome Browser tracks showing chromosomal location, common dbSNPs, VDR-, RXR-, CTCF– and IKZF3-binding sites and H3K27Ac marks present in the Ch17q21.2 locus. VDR binding sites overlapping with active regulatory elements are highlighted by colored boxes. C, Overview of the analyses conducted on the VDAART transcriptomic and clinical data.

Baseline VDR expression is confined to the CD4+ Foxp3-Teff cell population.

A, Gating strategy for the analysis of VDR expression in different leukocyte populations in the lung. B, Flow cytometric analysis and cell frequencies of VDR expression in CD45+, CD4+, CD8+ and CD19+ in the respective groups (n=3-5). C, Representative flow plots of lung CD45+ CD4+ T cells isolated from saline and HDM-immunized and –challenged mice for the co-expression of Foxp3 and VDR. Each symbol represents one mouse. Numbers in flow plots indicate percentages. Error bars indicate the s.e.m. Statistical tests: two-tailed Mann-Whitney U test. *P < 0.05. Data summarize results from two independent experiments.

VDR expression is induced during CD4+ Th2 cell differentiation.

A, Representative immunoblot for VDR expression in freshly isolated (0 hrs) and activated CD4+ CD62L+ CD44-naïve T cells (48 hrs.) and in vitro polarized Th1, Th2 and iTreg cells. Kidney samples of WT and VDR-/-mice were used as controls. B, Representative immunoblot for VDR expression in polarizing Th2 at the indicated days of in vitro culture. C, Representative image of VDR localization in in vitro polarized Th2 cells stimulated with vehicle (EtOH, [% v/v]) or calcitriol.

Vitamin D deficiency does not alter CD4+ IL-10 and IL-17A expression allergic airway inflammation.

Frequencies of BAL neutrophils (A), IL-10 (B) and IL-17A (C) CD4+ T cells. Each symbol represents one mouse. Error bars indicate the s.e.m. Data summarize results from two or three independent experiments with n>4 per group.

Vitamin D deficiency increases IL-13 expression in Th2 cells.

Flow cytometric analysis and frequencies of IL-13 expression in WT, vitamin D deficient and VDR-/-Th2 cultures (n=5-9 per group). Each symbol represents one individual sample and data summarize the results from two or three independent experiments. Error bars indicate the s.e.m. Statistical tests: Kruskal-Wallis Test with Dunn’s post-hoc analysis. *P < 0.05, **P < 0.01, ****P < 0.0001.

KEGG pathways analysis post vitamin D stimulation.

List of KEGG pathways enriched (p-adj. < 0.05) in DE genes after calcitriol stimulation. The x-axis indicates the number of genes annotated to each pathway. The bar colors indicate the FDR-adjusted p-value based on the colormap on the right side of the figure. The R package clusterProfiler was used for the enrichment analysis and visualization.

Active signaling pathways in Th2 cells.

A, Schematic representation of the network analysis employed. B, Network analysis to visualize active signaling paths in control stimulated Th2 cultures. Predefined set of cell membrane receptors were set as starting points (orange nodes) and annotated transcription factors were set as ending points (blue). All intermediate proteins on the PPI are visualized as grey nodes. Molecular interactions only overrepresented in control stimulated cells are shown as purple arrows.

Prenatal vitamin D supplementation increases VDR expression in CD4+ T cells.

A, Schematic presentation of the dietary intervention protocol and breeding scheme with vitamin D. B, Absolute neutrophil numbers in BAL (n≥8 per group). C, Flow cytometric analysis and quantification of VDR expression in lung CD4+ T cells of the respective groups. Each symbol represents one individual mouse and data summarize the results from two to three independent experiments with n≥5 per group. Error bars indicate the s.e.m. Statistical significance was determined with One-way ANOVA with Holm-Šidák’s post-hoc analysis. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001.

Prenatal vitamin D supplementation mitigates Th2 development and cytokine production in vitro.

A-E, Representative flow plots and cell frequencies of (A) Gata3, (B) VDR, (C), IL-2, (D) IL-13, (E) IL-5 expression in indicated groups. Each symbol represents one individual sample and data summarize the results from two to three independent experiments (n= 4-9 per group). Error bars indicate the s.e.m. Statistical test: two-tailed Mann-Whitney U-test (a – e). *P < 0.05, **P < 0.01, ***P < 0.001.

Schematic summary of vitamin D molecular genetics in the 17q12-21 region.