Skip to main content

Vitamin D and allergic airway disease shape the murine lung microbiome in a sex-specific manner

Abstract

Background

Vitamin D is under scrutiny as a potential regulator of the development of respiratory diseases characterised by chronic lung inflammation, including asthma and chronic obstructive pulmonary disease. It has anti-inflammatory effects; however, knowledge around the relationship between dietary vitamin D, inflammation and the microbiome in the lungs is limited. In our previous studies, we observed more inflammatory cells in the bronchoalveolar lavage fluid and increased bacterial load in the lungs of vitamin D-deficient male mice with allergic airway disease, suggesting that vitamin D might modulate the lung microbiome. In the current study, we examined in more depth the effects of vitamin D deficiency initiated early in life, and subsequent supplementation with dietary vitamin D on the composition of the lung microbiome and the extent of respiratory inflammation.

Methods

BALB/c dams were fed a vitamin D-supplemented or -deficient diet throughout gestation and lactation, with offspring continued on this diet post-natally. Some initially deficient offspring were fed a supplemented diet from 8 weeks of age. The lungs of naïve adult male and female offspring were compared prior to the induction of allergic airway disease. In further experiments, offspring were sensitised and boosted with the experimental allergen, ovalbumin (OVA), and T helper type 2-skewing adjuvant, aluminium hydroxide, followed by a single respiratory challenge with OVA.

Results

In mice fed a vitamin D-containing diet throughout life, a sex difference in the lung microbial community was observed, with increased levels of an Acinetobacter operational taxonomic unit (OTU) in female lungs compared to male lungs. This effect was not observed in vitamin D-deficient mice or initially deficient mice supplemented with vitamin D from early adulthood. In addition, serum 25-hydroxyvitamin D levels inversely correlated with total bacterial OTUs, and Pseudomonas OTUs in the lungs. Increased levels of the antimicrobial murine ß-defensin-2 were detected in the bronchoalveolar lavage fluid of male and female mice fed a vitamin D-containing diet. The induction of OVA-induced allergic airway disease itself had a profound affect on the OTUs identified in the lung microbiome, which was accompanied by substantially more respiratory inflammation than that induced by vitamin D deficiency alone.

Conclusion

These data support the notion that maintaining sufficient vitamin D is necessary for optimal lung health, and that vitamin D may modulate the lung microbiome in a sex-specific fashion. Furthermore, our data suggest that the magnitude of the pro-inflammatory and microbiome-modifying effects of vitamin D deficiency were substantially less than that of allergic airway disease, and that there is an important interplay between respiratory inflammation and the lung microbiome.

Background

Vitamin D has been proposed as an important regulator of both inflammation [1] and the microbiome [2, 3]. Indeed, it is hypothesised that there is a protective role for vitamin D in regulating the gut microbiome, which may shape predisposition for the development of inflammatory autoimmune and allergic diseases [2–4]. Vitamin D is commonly acquired through skin exposure to ultraviolet B radiation found in sunlight and through the diet. Circulating 25-hydroxyvitamin D (25(OH)D) is typically used as a measure of vitamin D status, although 1,25-dihydroxyvitamin D (1,25(OH)2D) is the most active vitamin D metabolite.

Our understanding of the effects of vitamin D on the microbiome and its associations with reduced tissue inflammation is currently limited to the gastrointestinal tract. Colitis severity and bacterial numbers in the colons of mice were increased during vitamin D deficiency [5]. Dietary-induced vitamin D deficiency altered the composition of the fecal microbiome of C57Bl/6 mice, which had increased relative quantities of Bacteroidetes, Firmicutes, Actinobacteria, and Gammaproteobacteria [6]. These mice also had increased colonic injury in comparison to mice fed a vitamin D-sufficient diet, but exhibited signs of hypocalcemia [6]. In other studies, the colons of 21-day-old vitamin D-deficient mice were enriched for Bacteroides/Prevotella colony forming units, although this difference disappeared with age [7]. CYP27B1−/− mice, which lack the expression of the 1α-hydroxylase enzyme, responsible for converting circulating 25(OH)D to active 1,25(OH)2D, had an increased fecal burden of Proteobacterium phylum (including Helicobacteraceae species) in a colitis model [8]. Treatment of CYP27B1−/− mice with 1,25(OH)2D (1.25 μg/100 g diet) suppressed colitis severity and Helicobacteraceae numbers [8]. Vitamin D supplementation (980 IU vitamin D3/kg per week for 4 weeks) of 16 healthy young adults reduced the relative abundance of Gammaproteobacteria such as Pseudomonas spp. and Escherichia/Shigella spp., increasing bacterial richness in the upper gastrointestinal tract, but not in ileum, colon nor in stool samples [9]. Changes in Bacteroides spp. were identified in the stool samples of African American men (n = 115) with prediabetes and hypovitaminosis D after weekly supplementation with a vitamin D analogue (ergocalciferol; 50,000 IU) [10]. Together, these data suggest that dietary vitamin D can alter the gut microbiome, reducing the abundance of potentially pathogenic species with these changes linked to reduced gastrointestinal inflammation and injury.

Vitamin D may shape the microbiome through a number of interdependent mechanisms. Firstly, vitamin D regulates innate immune responses. Vitamin D can induce the expression of antimicrobial peptides and proteins, such as cathelicidins and ß-defensins, which are produced by monocytes, macrophages and epithelial cells in the skin and lung [1, 3]. Vitamin D-deficient mice had reduced colonic expression of the antimicrobial protein, angiogenin-4, which was associated with increased bacterial load, and tissue inflammation [5]. Such antimicrobials may directly kill microbiota or be involved in activating innate immune processes such as autophagy in macrophages, promoting the ingestion of microbes into phagolysosomes for neutralization [11]. Vitamin D may promote tolerant adaptive immune responses by modulating the gut microbiome. As highlighted above, there is evidence linking the modified gut microbiome and increased gut inflammation of CYP27B1−/− mice, with fewer tolerogenic CD103+ dendritic cell in lamina propria [8]. These are cells responsible for shaping the regulatory T cell repertoire of the gut [12]. Ongoing presentation of bacterial antigens by dendritic cells may be required for tolerance towards gut microbes, with an essential role for tolerogenic dendritic cells and T cells to limit inflammation [13, 14]. Finally, vitamin D may up-keep epithelial integrity. Colonic epithelial cells from CYP27B1−/− mice expressed reduced levels of the cell-to-cell adhesion protein, E-cadherin [8]. Assa et al. [6] demonstrated that vitamin D-deficient mice had reduced colonic epithelial barrier function, with increased permeability that associated with increased proinflammatory cytokine expression, signs of colitis and relative quantities of Bacteroidetes, Firmicutes, Actinobacteria and Gammaproteobacteria in the intestinal microbiome.

In healthy human lungs, the bronchial tree holds approximately 2000 distinct bacterial genomes per cm2 [15]. There are diverse microbial communities, with Bacteroidetes, Proteobacteria and Firmicutes the mostly commonly identified at the phylum level [16]. Similarly, the healthy mouse lung microbiome is dominated by species of these phyla and also Actinobacteria and Cyanobacteria [17]. The lung microbiome changes from infancy through to adulthood, with the ratio of Bacteroidetes:Firmicutes/Gammproteobacteria increasing with age in specific-pathogen-free mice [18]. However, dysregulation of the microbiome of the lungs and associated tissues could contribute towards respiratory inflammation.

Indeed, microbial communities shift in the respiratory tract during respiratory disease [15, 19–23]. Hilty et al. [15] demonstrated that pathogenic proteobacteria were more common in bronchial brushings of the left upper lung lobe of adult asthmatics (n = 11) or patients with chronic obstructive pulmonary disease (COPD) (n = 5) than controls (n = 8), with similar findings in the bronchoalveolar lavage fluid (BALF) of asthmatic children (n = 13). Lung tissue samples from patients with COPD (n = 8) had increased bacteria of the Lactobacillus genus with bacterial communities distinct to healthy controls (n = 8), smokers without COPD (n = 8) or patients with cystic fibrosis (n = 8) [20]. Lactobacillus, Pseudomonas, and Rickettsia species were enriched in endobronchial brush samples from patients with chronic persistent asthma (n = 39) in comparison to healthy controls (n = 19) [21]. The bronchial brushings of patients with severe asthma (n = 30) were enriched with Actinobacteria and a member of the Klebsiella genus in comparison to healthy controls (n = 8) and mild-to-moderate asthmatics (n = 41) [22]. Together, these findings suggest that pathogenic bacteria are more common in the lungs and associated respiratory tissue of patients with respiratory disease, which may vary depending on the site examined.

Animal models offer the opportunity of dissecting various factors that change the lung microbiome to initiate and/or exacerbate inflammatory or allergic lung diseases. The effect of vitamin D deficiency on the lung microbiome has not been specifically investigated in vivo. We previously observed that inflammatory cell numbers in the BALF of vitamin D-deficient male mice with allergic airway disease were associated with increased bacterial load in the lungs [24]. Low levels of circulating 25(OH)D were associated with increased numbers of eosinophils and neutrophils in BALF in male mice with allergic airway disease [24]. Increased levels of BALF inflammation and microbial load in the lungs of initially vitamin D-deficient mice were reversed by dietary vitamin D supplementation [24]. In the current study, we hypothesised that vitamin D deficiency would alter the lung microbiome of naïve mice, potentially contributing towards increased respiratory inflammation, prior to the initiation of allergic airway disease.

Methods

Mice and diets

All experiments were performed according to the ethical guidelines of the National Health and Medical Research Council of Australia and with approval from the Telethon Kids Institute Animal Ethics Committee (AEC#238). Mice were purchased from the Animal Resources Centre, Western Australia. Mice were maintained under specific-pathogen-free conditions. Female 3 week-old BALB/c mice were fed semi-pure diets, which were either supplemented with vitamin D3 (2,280 IU vitamin D3/kg, SF05-34, Specialty Feeds, Perth, Western Australia) or not (0 IU vitamin D3/kg, SF05-033, Specialty Feeds) as previously described [24, 25]. From 8 weeks of age, female mice were mated with male BALB/c mice for up to 2 weeks. These male mice were fed standard mouse chow until breeding (Specialty Feeds, containing 2,000 IU vitamin D3/kg). Litter sizes from dams fed either diet were a mean of 5 pups per litter, with equal proportions of females and males. Offspring born were fed the vitamin D-replete or -deficient diets for the rest of the experiment, except in experiments when initially vitamin D-deficient mice were switched to a vitamin D-replete diet at 8 weeks of age (Fig. 1). Experiments were performed over a 12-month period between November 2011 and November 2012.

Fig. 1
figure 1

The dietary intake of mice in each treatment. In (A), female BALB/c mice (dams) were fed vitamin D-containing (VitD+) or vitamin D-null (VitD−) diets from 3 weeks of age and used to produce offspring. A subgroup of the vitamin D-deficient offspring was fed a vitamin D-supplemented diet from 8 weeks of age (VitD−/+)

Measurement of serum and BALF levels of 25-hydroxyvitamin D (25(OH)D)

Serum and BALF 25(OH)D levels were measured using IDS EIA ELISA kits (Immunodiagnostic Systems Ltd, Fountain Hills, AZ) as described by the manufacturer (limit of detection was 5–7 nmol/L).

Determining bacterial loads in lung samples

Following the BALF procedure (see below), the left lobe of each lung was minced using a sterile scalpel blade. A 2 mm3 sample of the lung was snap-frozen in liquid nitrogen. DNA was extracted from lung samples using the DNeasy Blood and Tissue DNA extraction kit as described by the manufacturer (Qiagen). Universal 16S rRNA primers (F primer = 5′-TCCTACGGGAGGCAGCAG T-3′; R primer = 5′-GGACTACCAGGGTATCTAATCCTGTT-3′ [26]) were used to detect bacteria in DNA samples using the Power SYBR Green PCR Master Mix as described by the manufacturer (Applied Biosystems, Calsbad, USA). These ‘universal primers’ have broad specificity to detect conserved regions of 16S rDNA from 34 bacterial species encompassing most bacterial groups [26]. An 18S rRNA endogenous control including a FAM-MGB probe was used as the internal control for this assay using conditions described by the manufacturer (Applied Biosystems).

Characterising the bacterial lung microbiome

To profile the microbial community of the lungs, we amplified the variable region V3-V4 of the 16S rRNA gene as previously reported [17], with modifications, using the primer pair of 341 F and 806R. The PCR mix was composed of 5 μl 5x Phusion buffer HF (7.5 mM MgCl2, Finnzymes, Finland), 0.5 μl 10 mM dNTPs, 1.25 μl 10 μM of each primer, 0.25 μl DNA polymerase (Hotstart Phusion 540 L, 1 unit/μl Finnzymes) and 3.5 μl template. The PCR program started with 98 °C (2 min), followed by 35 cycles of denaturation at 98 °C (5 s), annealing at 55 °C (30 s) and strand elongation at 60 °C (60 s). The PCR was finalised by a single elongation step at 72 °C for 5 min and than cooled to 4 °C. The size of the PCR product was evaluated using gel electrophoresis. The fragment was then excised and purified using the Montage Gel Extraction Kit (Merck Millipore) with negative controls also excised at the position of the expected fragment size. Adaptors were added to the amplicons in a second PCR run under the same conditions as PCR I with a reduced cycle number of 15. The sequences were generated with GS FLX Titanium (454 Life Sciences, Roche). Two sequencing runs were performed, with the machine cleaned between runs (using standard protocols) and the second run then proceeding, with each run containing samples from the ‘naïve’ and ‘OVA-induced allergic airway disease’ datasets. The 1,475,452 reads were trimmed for low quality (minimal quality score = 25) using the Qiime pipeline version 1.5.0 [27]. Only sequences with a minimal length of 200 bp were considered, the sequences were denoised using the Ampliconnoise algorithm [28]. Chimeras were removed using the Uchime algorithm [29]. Operational taxonomic units (OTUs) were picked de novo from quality-checked reads and clustered at 97 % sequence similarity using Uclust. Taxonomy was assigned using the RDP classifier (version 2.2) method and Greengenes as reference database [30]. After data treatment, 280,699 reads from the ‘naïve dataset’, and 330,733 reads from the ‘OVA-allergic airway disease’ dataset, were used for down stream analysis. To adjust for differences in sequencing depth between samples the OTUs were normalised using the cumulative sum-scaling method available in the metagenomeSeq Bioconductor package [31].

Bronchoalveolar lavage fluid for assessment of ß-defensin-2 (mBD2)

A total of 1 ml BALF was collected for each mouse as described previously [32]. Levels of mBD2 in BALF were detected using an ELISA kit and method supplied by USCN Life Science Inc (Wuhan, China).

Sensitisation and challenge of mice with ovalbumin

Ovalbumin (OVA) (Sigma Chemical Company, St Louis, MO, USA) was mixed with an aluminium hydroxide suspension (Alum, Serva, Heidelberg, Germany). This OVA/Alum solution was diluted in 0.9 % saline to sensitise and boost mice intraperitoneally (200 μl) with 1 μg OVA and 0.2 mg Alum as previously described, with mice sensitised at 12 weeks of age and boosted at 14 days of age [24]. For respiratory challenge, mice were placed in a Perspex box and inhaled a 1 % OVA-in-saline (1 mg/ml) aerosol delivered using an ultrasonic nebuliser (UltraNebs, DeVilbiss, Somerset, PA, USA) for 30 min. This was performed once, 7 days after the OVA/Alum boost. Lungs and BALF was obtained 24 h after the aerosol challenge.

Bronchoalveolar lavage fluid for assessment of inflammatory cells

BALF cells (5 × 105) were spun onto glass slides using a cytocentrifuge and differential counts of inflammatory cells performed after staining cells with the DIFF-Quik Stain Set 64851 (Lab Aids, Narrabeen, NSW, Australia) as per the manufacturer’s instructions. At least 300 cells were counted for each sample from ≥3 independent fields of view (x100).

Measurement of circulating immune cells

Blood was obtained from mice and immediately added to K2EDTA-coated Microtainer tubes (BD, Franklin Lakes, NJ) to prevent clotting. The proportions of circulating neutrophils and monocytes within the white blood cell population was determined using the ADVIA® 120 Haematology System (Siemens Healthcare Diagnostics Inc, Tarrytown, NY).

Statistical analyses

Data were compared using one-way ANOVA or an unpaired two-way student’s t test or correlated (Spearman’s rank correlation) using the Prism 5 for Mac OS X statistical analysis program. Differences were considered significant with a p-value <0.05. We used the non-parametric Kruskal-Wallis test (p < 0.05) to test for significant differences in OTU abundance. Treatment effects on the complex microbial communities in the lung were observed by applying the vegdist function (part of the R vegan package) on the raw metagenomeSeq normalised OTU abundance table to generate the Bray-Curtis dissimilarity between samples.

The lungs of specific-pathogen-free mice contain little microbial DNA, and as recently described for microbiome studies [33] are of high risk of contamination by bacterial DNA present in molecular biology reagents used for high throughput sequencing. These can especially contribute towards incorrect interpretations of cluster analyses. However, simply removing OTUs from the abundance table might introduce bias by potentially removing ‘real’ OTUs present in the sample (e.g. Escherichia spp.). Therefore, the entire cluster analysis that was based on the metagenomeSeq normalised OTU table was performed twice; first with, and then without 126 OTUs detected in the blank reagents and listed in the OTU table that contained afterwards 5,635 OTUs for downstream analysis. Most (96 %) of the contaminating OTUs were taxonomically assigned to the halophilic marine genus Halomonas. The rest were Shewanella, Delftia and Stenotrophomonas. Unless stated differently we used the cleaned microbial OTU table describing the lung microbiome according to [34].

The Bray-Curtis dissimilarity test takes microbial diversity and relative abundance into account. Grouping of samples was visualised by using ordination applying non-metric multidimensional scaling (NMDS) generated in the R vegan package [35]. Microbial clustering was further evaluated with the analysis of similarity (ANOSIM), which takes the ranked dissimilarities (Bray-Curtis) of samples belonging to the same treatment group and compares it to the distance of samples of a different treatment [36]. The closer the ANOSIM generated R-value was to 1, the larger the variation between microbial treatments, whereas 0 indicates no difference. The results were than tested for significance by 999 permutations with a 5 % significance level. Individual variation of selected microbes was displayed with the Euclidean distance in a one-sided dendrogram (Heatmap) based on the log-transformed metagenomeSeq normalised OTU counts. The percentage of each OTU was relative to the number of OTUs for each sample.

Results

Dietary vitamin D increased circulating but not BALF 25(OH)D levels

Circulating serum levels of 25(OH)D were measured in 8-week old offspring born to vitamin D-replete and -deficient dams (Fig. 2a). Serum 25(OH)D levels were <20 or ≥50 nmol.L−1 in offspring fed a vitamin D-deficient or -replete diet, respectively. When vitamin D-deficient offspring were fed a vitamin D-containing diet for four weeks, their serum 25(OH)D increased to levels equivalent to mice fed a vitamin D-containing diet throughout the experiment (Fig. 2a). BALF 25(OH)D levels were low, although greater than the detection limit of the ELISA and not affected by dietary vitamin D (Fig. 2b).

Fig. 2
figure 2

Vitamin D deficiency reduced serum 25(OH)D levels but had no effect on BALF 25(OH)D levels. Female BALB/c mice (dams) were fed vitamin D-containing (VitD+) or vitamin D-null (VitD-) diets from 3 weeks of age and used to produce offspring. A subgroup of the vitamin D-deficient offspring was fed a vitamin D-supplemented diet from 8 weeks of age (VitD−to vitD+). In (a), serum 25(OH)D levels in offspring at 8 and 12 weeks of age (mean ± SEM for ≥5 mice per group). In (b) BALF 25(OH)D levels in offspring at 12 weeks of age (mean ± SEM for ≥3 mice per group, the broken line indicates the level of detection (LOD) for 25(OH)D; 7 pg/ml). (*p < 0.05)

The effects of dietary vitamin D on bacterial loads in the lungs of naive mice

There was a trend similar to our previously published findings (24) for increased bacterial load in the lungs of vitamin D-deficient male mice, although this was not statistically significant (One-way ANOVA, F = 1.481, p = 0.245; Fig. 3). Bacterial load was quantified using a PCR with 16S rRNA primers capable of amplifying DNA of most bacterial species.

Fig. 3
figure 3

The effects of vitamin D deficiency on lung bacterial load in naïve mice. Female BALB/c mice (dams) were fed vitamin D-containing (VitD+) or vitamin D-null (VitD−) diets from 3 weeks of age and used to produce offspring. A subgroup of the vitamin D-deficient offspring was fed a vitamin D-supplemented diet from 8 weeks of age (VitD−to VitD+). Bacterial DNA levels were determined in the lungs of naïve mice using a PCR with universal primers for detection of bacterial 16S rRNA gene (mean + SEM for 10–15 mice/treatment)

Dietary vitamin D had a sex-specific effect on the lung microbiome of naive mice

We used high-throughput sequencing to determine how vitamin D might affect the composition of bacteria in the lung microbiome. After data treatment (see methods) there was an average read distribution of 4,803 sequences per naïve animal. We first evaluated if dietary vitamin D had a significant impact on the microbial composition enumerating the dissimilarity of the OTUs. As shown in NMDS plots, there was no clustering by diet (Fig. 4a), sex (Fig. 4b), or the combination of both (Fig. 4c). We then evaluated the relative abundance of OTUs with 23–112 OTUs identified/mouse (Fig. 5). Neither dietary vitamin D nor sex significantly affected bacterial diversity in the lungs as defined by the relative number (Fig. 5, Kruskal-Wallis Test, p > 0.05) of OTUs. In Fig. 6, a heat map was used to show the most frequent OTUs identified from the lungs of male mice. There was very little effect of dietary vitamin D on the relative proportion of each OTU (Fig. 6), with similar findings in female mice (data not shown), and little difference in bacterial variation between treatments identified using ANOSIM analysis, except in mice fed a vitamin D-containing diet throughout the experiment (Table 1).

Fig. 4
figure 4

The lung microbiomes of vitamin D-replete and -deficient naïve mice were similar. Female BALB/c mice (dams) were fed vitamin D-containing (+) or vitamin D-null (−) diets from 3 weeks of age and used to produce offspring. A subgroup of the vitamin D-deficient offspring was fed a vitamin D-supplemented diet from 8 weeks of age (−/+). Non-metric multidimensional scaling (NMDS) plots depict the dissimilarity of the detected operational taxonomic units (OTUs) from the lungs, with results shown for clustering by (a) diet, (b) sex, or (c) all (M = male, F = female; n = 9-10/treatment)

Fig. 5
figure 5

Dietary vitamin D did not affect the diversity of the lung microbiome. Female BALB/c mice (dams) were fed vitamin D-containing (+) or vitamin D-null (−) diets from 3 weeks of age and used to produce male (M) and female (F) offspring. A subgroup of the vitamin D-deficient offspring was fed a vitamin D-supplemented diet from 8 weeks of age (−/+). Shown are the number of operational taxonomic units (OTUs) observed at a sequencing depth of 1,200 per animal, with no significant effect of sex or vitamin D. Data is shown with the median as a blue line, the first or third quartiles as the upper and lower boxes (respectively) and the upper and lower whiskers show the 1.5 times of the interquartile distance and minimum and maximum data (n = 9–10/treatment, with an outlier denoted as a circle)

Fig. 6
figure 6

The abundance of bacterial operational taxonomic units was minimally affected by vitamin D deficiency. Female BALB/c mice (dams) were fed vitamin D-containing (+) or vitamin D-null (−) diets from 3 weeks of age and used to produce offspring. A subgroup of the vitamin D-deficient offspring was fed a vitamin D-supplemented diet from 8 weeks of age (−/+). Shown is a heat map of the of operational taxonomic units (OTUs) significantly detected in the lungs of male (M) naïve offspring. OTU counts were log-transformed; with red denoting increased frequency of a given OTU, and dark blue no occurrence (n = 9-10/treatment)

Table 1 Comparison of the lung microbiome within groups with analysis of similarity (ANOSIM) examining the effects of sex or dietary vitamin D in naïve mice

We observed a significant separation (NMDS cluster) between female and male mice that were fed the vitamin D-supplemented diet throughout the experiment (Fig. 7). This difference was attributed to increased representation of an Acinetobacter OTU in female mice, which varied significantly between the mice (Wilcoxon Rank Sum test, p = 0.02), together with a single low-abundant OTU that could only be annotated as bacteria. Reduced carriage of the Acinetobacter OTU was observed in male mice (2/10 mice; ≤5 % of sequences in the observed sequences; mean = 0.7 %), compared to female mice (7/10 mice; ≤40 % of sequences in the observed sequences; mean = 7.0 %) fed a vitamin D-containing diet throughout the experiment.

Fig. 7
figure 7

The lung microbiomes of vitamin D-replete and -deficient naïve mice were similar. Female BALB/c mice (dams) were fed vitamin D-containing (+) or vitamin D-null (−) diets from 3 weeks of age and used to produce offspring. A subgroup of the vitamin D-deficient offspring was fed a vitamin D-supplemented diet from 8 weeks of age (−/+). Non-metric multidimensional scaling (NMDS) plots depict the dissimilarity of the detected operational taxonomic units (OTUs) from the lungs, with results shown for clustering for mice fed the vitamin D+ diet only. (M = male, F = female; n = 9–10/treatment)

Antimicrobial levels in the lungs were enhanced by dietary vitamin D in both male and female mice

In both male and female mice, dietary vitamin D increased ß-defensin-2 (mBD2) protein levels in BALF, with the effects of early-life vitamin D-deficiency reversed by dietary vitamin D (Fig. 8a). There was no effect of dietary vitamin D on serum mBD2 levels (Fig. 8b). Levels of mBD2 were significantly greater in serum than BALF. These results suggest that the sex-dependent effects of vitamin D on Acinetobacter in the lungs were not due to differences in the expression of mBD2 as dietary vitamin D increased levels to a similar degree in the BALF of male and female mice.

Fig. 8
figure 8

Vitamin D deficiency reduced ß-defensin-2 (mBD2) protein levels in the BALF of naïve mice. Female BALB/c mice (dams) were fed vitamin D-containing (VitD+) or vitamin D-null (VitD-) diets from 3 weeks of age and used to produce offspring. A subgroup of the vitamin D-deficient offspring was fed a vitamin D-supplemented diet from 8 weeks of age (VitD− to VitD+). ß-defensin-2 protein levels were determined in the BALF (a) (n = 10 mice/treatment) and serum (b) (n ≥ 3 mice/treatment) of naïve mice (mean + SEM) (*p < 0.05, the broken line indicates the level of detection (LOD) for ß-defensin-2; 3 pg/ml)

Serum 25(OH)D levels inversely correlated with the most prominent OTU of Pseudomonas in the lung

In a separate analysis, we identified an inverse correlation between serum 25(OH)D levels and the total number of OTUs (Fig. 9, Spearman’s rho = −0.592, p = 0.001). When each sex was considered separately, some evidence for an inverse correlation was observed between total OTUs in the lungs and serum 25(OH)D levels in females (Spearman’s rho = −0.642, p = 0.052) but not males (Spearman’s rho = 0.463, p = 0.238). When examining individual OTUs in the lungs, there was a strong negative correlation between 25(OH)D levels in serum and an unspecified Pseudomonas OTU (Spearman’s rho = −0.556, p = 0.016), with a marginal evidence observed in females (Spearman’s rho = −0.611, p = 0.081) but not males (Spearman’s rho = −0.343, p = 0.366). Some evidence for a negative correlation was observed between serum 25(OH)D and all Pseudomonas OTUs (Spearman’s rho = −0.429, p = 0.07). However, there was no evidence of a relationship between total OTUs or the unspecified Pseudomonas OTUs in the lung and mBD2 levels in BALF (total OTUs; Spearman’s rho = 0.109, p = 0.429; Pseudomonas OTU Spearman’s rho = 0.202, p = 0.138 for BALF) or serum (total OTUs; Spearman’s rho = 0.114, p = 0.673, Pseudomonas OTU Spearman’s rho = 0.114, p = 0.673).

Fig. 9
figure 9

Number of operational taxonomic units in lungs inversely correlated with serum 25(OH)D. Female BALB/c mice (dams) were fed vitamin D-containing (+) or vitamin D-null (−) diets from 3 weeks of age and used to produce offspring. A subgroup of the vitamin D-deficient male and female offspring was fed a vitamin D-supplemented diet for 4 weeks from 8 weeks of age (−/+). A significant inverse correlation was observed between the number of operational taxonomic units (OTUs) in the lungs and circulating 25(OH)D levels (for n = 18 mice for whom both measures were obtained)

The lung microbiomes of naïve mice, and, ovalbumin (OVA)-sensitised and -challenged mice clustered separately

We hypothesised that allergic sensitisation and challenge with OVA (to induce allergic airway disease) would significantly change the microbiota composition of the lungs. The relative effects of the induction of OVA-induced allergic airway versus that of vitamin D deficiency (alone) on the lung microbiome are unknown. To better understand the magnitude of the observed effect of dietary vitamin D on the lung microbiome, we compared its capacity to regulate the lung microbiome to the effects of the induction of allergic airway disease. Mice were intraperitoneally sensitised and boosted with ovalbumin (OVA) (1 μg) and Alum hydroxide suspension (Alum, adjuvant, 0.2 mg) at 12 and 14 weeks of age (respectively), which was followed by a respiratory challenge at 15 weeks of age with aerosolised OVA-in-saline (1 mg/ml) [24]. Lungs were sampled 24 h after the challenge. We compared the overall microbial community from the lungs of the naïve mice characterised above (n = 58; 4,803 sequences per naïve mouse) with the lungs of mice with OVA-induced allergic airway disease (n = 68; 4,469 sequences per mouse). As shown in Fig. 10a, the lung samples of naïve mice and mice with allergic airway disease clustered significantly apart, which was confirmed using a compositional analysis to compare the difference in numbers of observed bacterial OTUs (Bray-Curtis, R = 0.334, p = 0.001). A NMDS analysis confirmed that observed differences were not induced by performing two sequencing runs (data not shown). Bacterial diversity was not affected with the mean number of OTUs for naïve mice or mice with allergic airway disease (naïve, 138 + 5.1; allergic airway disease, 128 + 5.4; mean + SEM, p = 0.132), respectively.

Fig. 10
figure 10

Ovalbumin-induced allergic airway disease substantially modified the lung microbiome. The lung microbiomes of naïve mice (naïve, n = 56) were compared with those with ovalbumin (OVA)-induced allergic airway disease (OVA, n = 68), 24 h after respiratory challenge with OVA, with results combined for all naïve or all OVA mice. To induce allergic airway disease, mice were injected with OVA and Alum at 12 weeks of age, boosted with OVA and Alum at 14 weeks of age, and then administered a respiratory challenge with OVA at 15 weeks of age. Lungs were obtained 24 h later. In (a), a non-metric multidimensional scaling (NMDS) plot depicts the dissimilarity of the detected lung of operational taxonomic units (OTUs) for all mice, and in (b) for mice only fed the vitamin D-supplemented diet (+). (M = male, F = female; n = 9–10/treatment)

The lung microbiome was significantly altered by allergic airway disease

The strength of the difference between the lung microbiomes of naïve mice, and, mice with allergic airway disease was quantified by ANOSIM. A significant difference was identified (Table 2). As for naïve mice (Fig. 7), a significant separation (NMDS cluster) was observed between female and male mice that were fed the vitamin D-supplemented diet throughout the experiment with OVA-induced allergic airway disease (Fig. 10b; Table 3 for ANOSIM bacterial variance). There was a significant effect of allergic airway disease on the relative proportion of OTUs detected in mice. The most abundant OTUs detected, which each accounted for >1 % of the sequences after normalization, are shown for the naïve mice, or, mice with allergic airway disease in Fig. 11a. Proportions of the environmental microbial group GNO2 (Naïve, mean = 9.5 %; OVA, mean = 0.0 %) and an unassigned bacterial OTU (Naïve, mean = 1.1 %; OVA, mean = 0.0 %) were reduced in mice with OVA-induced allergic airway disease, while Acinetobacter (Naïve, mean = 1.29 %; OVA, mean = 1.8 %), Chloroplast (Naïve, mean = 0.1 %, OVA, mean = 5.0), Comamonadaceae (Naïve, mean = 0.1 %; OVA, mean = 4.0.%), Cyanobacteria (Naive, mean = 0.1 %; OVA, mean = 1.0 %) Staphylococcus (Naïve, mean = 0.5 %; OVA, mean = 1.2 %), Micrococcaceae (Naïve, mean = 0.1 %; OVA, mean = 1.8 %) and Verrucomicrobiaceae (Naïve, mean = 0.1 %; OVA, mean = 7.6 %) OTUs, and an OTU assigned to Peptostreptococcaceae (Naïve, mean = 1.0 %; OVA, mean = 1.5 %) were increased by the induction of allergic airway disease. The Acinetobacter OTU significantly increased in the OVA samples was not the same Acinetobacter OTU identified as responsible for largely driving the difference between male and female mice fed a vitamin D-containing diet throughout the experiment. There was no significant difference in the number of any Pseudomonas OTUs detected in the lungs of naïve mice, or, mice with allergic airway disease. Together, these results suggest that the effects of allergic airway disease on the lung microbiome of mice far exceeded those of vitamin D deficiency.

Table 2 Group comparison of the lung microbiomes of naïve mice, with those from mice with OVA-induced allergic airway disease within treatments for analysis of similarity (ANOSIM)
Table 3 Comparison of the lung microbiomes of mice with OVA-induced allergic airway disease with analysis of similarity (ANOSIM) examining the effects of sex or dietary vitamin D
Fig. 11
figure 11

Ovalbumin-induced allergic airway disease substantially modified the lung microbiome. The lung microbiomes of naïve mice (naïve, n = 56) were compared with those with ovalbumin (OVA)-induced allergic airway disease (OVA, n = 68), 24 h after respiratory challenge with OVA, with results combined for all naïve or all OVA mice. To induce allergic airway disease, mice were injected with OVA and Alum at 12 weeks of age, boosted with OVA and Alum at 14 weeks of age, and then administered a respiratory challenge with OVA at 15 weeks of age. Lungs were obtained 24 h later. In (a), a heat map shows OTUs significantly detected in the lungs of naïve and OVA-sensitised and -challenged mice. OTU counts were log-transformed; with red denoting increased frequency of a given OTU, and dark blue no occurrence. In (b), a heat map compares the sex-specific differences in detected Pseudomonas and Acinetobacter OTUs in naïve mice and mice with OVA-induced allergic airway disease (M = male, F = female), with a significant difference (*p < 0.05) denoting significant more of the most abundant Acinetobacter OTU in male mice with OVA-induced allergic airway disease (for mice fed a vitamin D-containing diet throughout the experiment)

Even though OVA-induced allergic airway disease substantially altered the lung microbiome, the sex-specific effect on the abundance of Acinetobacter OTUs was still apparent in mice fed exclusively a vitamin D-containing diet (Table 3). Acinetobacter OTUs were more frequent in OVA-induced allergic airway disease than naïve mice (Fig. 11b). There were a further 11 OTUs varied in their frequency when the lung microbiomes of female and male mice with OVA-induced allergic airway disease were compared. However, these were present at a very low frequency (<0.01 % of the total normalised sequences). There was no effect of OVA-induced allergic airway disease on the proportions of the dominant Acinetobacter OTU in female mice; however, induction of OVA-induced allergic airway disease increased the proportions of this OTU in male mice. For example, for male mice fed the vitamin D-containing diet throughout the experiment, the percentage of this OTU increased from 0.69 + 0.51 (mean + SEM) in naïve mice to 4.7 + 1.5 in mice with OVA-induced allergic airway disease (*p < 0.05, n = 10–13 mice/treatment). There was no significant correlation between the number of Acinetobacter OTUs and 25(OH)D or mBD2 levels in BALF or serum (Spearman p > 0.05, data not shown).

OVA-induced allergic airway disease had a substantially greater effect on inflammatory cells in BALF than vitamin D deficiency

The effects of vitamin D deficiency and the induction of OVA-induced allergic airway disease on numbers of neutrophils and macrophages/monocytes in the BALF or blood were then compared. Substantially increased BALF neutrophil numbers were observed in both male (Fig. 12a) and female (Fig. 12e) mice with OVA-induced allergic airway disease as expected; however, there was only a marginal effect of vitamin D deficiency. There was also a significant increase in neutrophil numbers in the blood of female mice with OVA-induced allergic airway disease, with no effect of vitamin D deficiency (Fig. 12f). Neither treatment significantly modulated macrophage levels in BALF (Fig. 12b, f), nor monocytes levels in blood (Fig. 12d, h). These data suggest that there is an important interplay between respiratory inflammation and the lung microbiome.

Fig. 12
figure 12

Allergic airway disease had a substantially greater effect on bronchoalveolar lavage neutrophil numbers than vitamin D deficiency. Neutrophil and macrophage numbers in the bronchoalveolar lavage (BALF) and blood of naïve mice fed a vitamin D (VitD)-replete (+) or -deficient (−) or initially vitamin D-deficient then replete diet (−/+) diet were compared with those of vitamin D-replete mice with OVA-induced allergic airway disease, 24 h after respiratory challenge with OVA. For OVA-sensitisation, mice were injected with OVA and Alum at 12 weeks of age, boosted with OVA and Alum at 14 weeks of age, and then administered a respiratory challenge with OVA at 15 weeks of age. BALF or blood was obtained 24 h later. In (a-d) results from males, and (e-h) results from female mice. In (a) and (e) neutrophils, and (b) and (f) macrophages in BALF are shown (n ≥ 14 mice/treatment). In (c) and (g) neutrophils, and (d) and (h) monocytes per ml of blood are depicted (n = 5–10 mice/treatment). Data is shown as mean + SEM (*p < 0.05)

Discussion

We observed modest effects of dietary vitamin D on the bacterial composition of the lung microbiome. Serum 25(OH)D levels inversely correlated with the number of OTUs detected, and more specifically the presence of an unidentified Pseudomonas OTU in the lungs of naïve mice. These results suggest that vitamin D sufficiency limited the number of respiratory pathobionts (like Pseudomonas) rather than increasing the number of protective commensals in the lungs. We also observed a sex difference in mice fed the vitamin D-containing diet throughout the experiment, with differences observed largely limited to a single Acintobacter OTU. The induction of allergic airway disease altered the lung microbiome composition to a degree that far exceeded any of the more subtle effects of vitamin D deficiency. More external environmental microbes, such as Cyanobacteria and Chloroplast, were detected in mice with allergic airways disease as well as OTUs previously identified in the lung microbiome, such as Staphylococcus, Peptostreptococcaceae and Acinetobacter. BALF neutrophil numbers were increased in mice of both sexes by ~10-fold with the induction of OVA-induced allergic airway disease, while the effects of vitamin D deficiency on these cell numbers in BALF were modest.

There were sex-specific differences in the lung microbiomes of mice fed a vitamin D-containing diet throughout the experiment. The Acinetobacter OTU, which accounted these effects, was distinct from the Acinetobacter OTU modulated by OVA-induced allergic airway disease. The nature of the sex-specific effect of vitamin D on the Acinetobacter OTUs differed, depending on whether mice were naïve or had OVA-induced allergic airway disease. We have previously shown sex-specific clustering of the lung microbiome in a different mouse strain (C57Bl/6) using denaturing gradient gel electrophoresis [37]. Acinetobacter may be a common commensal of human [38] and mouse lungs [17]. This genus may colonise the lungs early in life and has been linked to protection from allergy [39–41]. The switching of the direction of the sex-specific effects may reflect differences in the severity of OVA-induced allergic airway disease in male and female mice that are regulated by vitamin D, as previously observed [24]. The sex-specific effects on the lung microbiome were not reproduced by supplementing initially deficient mice with vitamin D, supporting the hypothesis that colonization with Acinetobacter species in the lungs occurs early in life.

As we have previously observed, 25(OH)D levels were reduced in male mice as compared to female mice fed a vitamin D-containing diet (see also [24, 25, 42, 43]). This observation is explained by increased renal expression of 24-hydroxylase (CYP24A1), the enzyme responsible for breaking down active 1,25(OH)2D, in male mice [42]. Hormones influence other effects of vitamin D, with reports of functional synergies between 1,25(OH)2D and 17-β-estradiol in T cells [44, 45], and differential regulation of 1,25(OH)2D-modulated pathways by testosterone in males and females [45, 46]. Testosterone increased Foxp3 (regulatory gene) expression in T cells from women, but not men [46] and 1,25(OH)2D more potently induced CD4 + CD25 + Foxp3+ (regulatory) T cells from PBMCs of female than male subjects [45], although we observed similar reductions in TReg cell percentages in the skin-draining lymph nodes of vitamin D-deficient male and female mice [47]. These effects of estrogen and testosterone on immune cell function suggest that females may be more resilient to the effects of vitamin D deficiency than males.

The inverse relationship between lung Pseudomonas and serum 25(OH)D could be related to increased levels of mBD2 in BALF of vitamin D-sufficient mice. mBD2 is a member of the defensin family of peptides, which have well-established killing effects on Gram-negative bacteria like Pseudomonas. mBD2 exhibits significant protein sequence homology with human ß-defensins 1 and 2 [48] and is expressed by a variety of epithelial cells (including tracheal) as well some immune cells (reviewed in [49]). mBD2 expression is induced during infection with Gram-negative bacteria, their products (e.g. lipopolysaccharide), and various proinflammatory cytokines (e.g. tumour necrosis factor) [48, 50]. Lipopolysaccharide may also drive the production of active 1,25(OH)2D from circulating 25(OH)D, resulting in synthesis of ß-defensins [51]. mBD2 is an immunoregulatory molecule, and its expression improves bacterial resistance by promoting proinflammatory cytokine expression through activation of toll-like receptor-4 [50]. Using immunohistochemistry, other researchers have detected only minimal concentrations of mBD2 in the lungs of BALB/c mice with OVA-induced allergic airway disease [52]. In addition, reduced levels of cathelin-related antimicrobial peptide (the mouse equivalent of cathelicidin) were detected in mice with OVA-induced allergic airway disease in comparison to non-sensitised mice [53]. These findings suggest that OVA-induced allergic airway disease may reduce the expression of antimicrobials like mBD2 levels in the lungs.

Regulation of mBD2 by vitamin D may significantly inhibit specific bacteria like Pseudomonas. Wu et al. demonstrated that mBD2 was required for host resistance against corneal infection with Pseudomonas aeruginosa [49]. These bacteria were more frequently detected in the sputum of vitamin D-deficient patients with bronchiestasis than -replete patients [54]. However, in another study there was no difference in the incidence of P. aeruginosa infection in children with cystic fibrosis that were vitamin D-deficient or -sufficient (25(OH)D > 30 μg/L) [55]. Active 1,25(OH)2D (10 nM) induced antimicrobial activity in bronchial epithelial cells against P. aeruginosa [56]. While we observed a significant inverse relationship between total OTUs or numbers of an unspecified Pseudomonas OTU and serum 25(OH)D, there was no inverse association with mBD2 in BALF. It is likely that another explanation, such as compromised epithelial integrity [6, 8], is responsible for the observed inverse relationship between serum 25(OH)D and lung OTUs.

In addition to Pseudomonas, we identified microbes in the lungs of our naïve mice that were previously detected in the lungs of humans [57] and other mice [17] including OTUs of Micrococcus, Staphylococcus, Cupriavidus and Streptococcus. The microbiome dataset from the lung tissue of naive mice (in particular) was very sparse, with few OTUs found among all lung samples (~130–140 OTUs/mouse). We obtained lung samples directly from the thoracic cavity of euthanised specific-pathogen-free mice, limiting possible contamination from the upper respiratory tract and oral cavity. It is possible that the modest effects of dietary vitamin D were due to the scarcity of the collected DNA. However, we believe this is not likely, due to the rigorous use of negative controls to identify OTUs of contaminating bacteria [33], our observations of association between circulating 25(OH)D levels and lung OTUs, and the significant microbial shift induced by the induction of allergic airway disease.

Methodological considerations are important, as we have previously shown differences in the microbial diversity of BALF and lung tissue [17]. With the small biomass of bacteria expected in the lungs of specific-pathogen-free naïve mice [58], there could also be problems associated with low yield DNA introduced through bacterial contamination of commercially purchased products (eg. in molecular-grade water or PCR reagents, [34]). As described by Salter et al. [34], reagents and extractions kits used for 16S amplicon library preparation can mask real biological effects and or increase false results. This is particularly a problem with lung samples as the bacterial yield is low, Therefore we analyzed data both with and without OTUs considered as contamination by sequences at low levels in our negative controls, allowing us to observe the sex-specific differences in mice fed only a diet containing vitamin D.

A limitation of the current study was that the naïve mice did not receive a placebo treatment (i.e. Alum bolus and nebulisation). We also did not examine the microbiome of other locations including the BALF [15, 17], or the oro- [59] or naso-pharynx [60], which could limit comparisons with human studies. We could not discern cause and effect in our modeling: for example, with allergic airway disease did the induced inflammation change the microbiome (or vice versa)? A detailed time-course analysis and use of further control groups (i.e. Alum only) could help determine the cause and effect relationships in future studies.

Another consideration is the model of asthma used. This well-characterised model [61] involved sensitisation of mice with a ‘low dose’ of the allergen OVA (1 μg) with Alum (0.2 mg), which caused methacholine-induced airway hyperresponsiveness, airway eosinophilia and neutrophilia, increases in BALF levels of IL-5, and circulating allergen-specific IgE and IgG [24]. This model induces a phenotype similar to allergic asthma typified by a T helper type-2 (Th2) immune response, which is induced by a single well-defined allergen [61]. There is no ideal animal model of allergic asthma [62]; commonly known as allergic airway disease in mice. House dust mite (HDM) models offer an advantage of using a human allergen; however, HDM preparations can be contaminated with varying quantities of bacterial-derived lipopolysaccharide. The induction of neutrophilic and eosinophilic inflammation in HDM-induced allergic airway disease models are dependent on toll-like receptor-4 [63], complicating the interpretation of the inflammatory response [62], especially with regards to examining the lung microbiome. We have shown here that the induction of OVA-induced allergic airway disease significantly altered the lung microbiome and induced lung inflammation. We hypothesise that factors that modulate airway inflammation per se, including those induced by other allergens, such as OVA with Alum or HDM extract [61], will also modulate the lung microbiome to perpetuate pulmonary inflammation and disease.

In previous studies, we have shown that the effects of vitamin D deficiency on airway inflammation in male mice were dependent on the dose of OVA and Alum used to sensitise mice [24]. However, vitamin D deficiency did not significantly modify airway resistance, tissue elastance or damping in male mice with OVA-induced allergic airway disease induced by this low-dose sensitisation [24]. In other studies we examined the lung function of naïve vitamin D-deficient and -replete BALB/c mice and observed increased airway resistance and tissue damping in female (but not male) vitamin D-deficient but otherwise naive mice [43]. These observations suggest that any protective effect of vitamin D on lung function is limited to female mice, and combined with a trend for an inverse correlation between total OTUs in the lungs and serum 25(OH)D levels in female mice, may suggest that the lung microbiome could regulate lung function in a sex-dependent fashion. In addition, vitamin D may maintain optimal lung function by preventing airway remodeling through a process dependent on transforming growth factor-β [43].

In addition to the sex-dependent effects of vitamin D deficiency on the severity of lung inflammation in mice with allergic airway disease [24], we have previously shown that vitamin D deficiency increased the capacity of airway-draining lymph node cells from male and female mice to proliferate and produce Th2 cytokines [24, 25]. The effects of deficiency were reversed by subsequent supplementation with dietary vitamin D3 [24]. Vitamin D deficiency increased the influx of lymphocytes into BALF in response to exposure to HDM; however, deficiency was protective and reduced airway smooth muscle mass and airway resistance induced by HDM [65]. Increased OVA-specific IgE and IgG1 were detected in vitamin D-deficient female BALB/c mice following OVA/Alum sensitisation (without further respiratory challenge) [66]. Increased lung eosinophil numbers and CD4 + T1ST2+ cells (Th2 cells), and reduced CD4 + IL-10+ (regulatory) cells were observed in young adult offspring born to vitamin D-deficient BALB/c dams following chronic intranasal instillation of HDM to offspring from 3 days of age [67]. However, as observed in our studies with OVA-induced allergic airway disease, there was no effect of vitamin D deficiency on airway hyperresponsiveness [67]. Collectively, these studies suggest that vitamin D deficiency promotes Th2 responses and the accumulation of eosinophils and neutrophils in the lungs of mice with allergic airway disease, without further compromising lung function.

We detected a significant negative relationship between circulating 25(OH)D and total OTUs detected in the lungs, suggestive of reduced bacterial diversity with increasing serum 25(OH)D. Others have noted that bacterial diversity is increased in healthy lungs when compared to diseased lungs [23]. Patients with poorly controlled asthma (n = 30) had reduced bacterial diversity and species richness, with increased Haemophilus influenzae in sputum samples, particularly in younger males with increased neutrophils (n = 7) [23]. Changes in bacterial diversity were associated with oral corticosteroid intake, airway obstruction, and eosinophilia in lung lavage fluid [21]. Conversely, bacterial diversity in bronchial epithelial brushings was positively correlated with bronchial hyperresponsiveness of adults with sub-optimally controlled asthma (n = 65) [19]. However, those with increased baseline diversity had greater improvements in bronchial hyperresponsiveness in response to clathriomycin treatment [19].

Conclusion

In these studies the capacity of OVA-induced allergic airway disease to modulate the lung microbiome and induce respiratory inflammation far exceeded any more subtle effects of vitamin D. We also show that the lung microbiome is differentially modified by sex in mice fed a vitamin D-containing diet, through specific effects on Acinetobacter OTU. It would be interesting to investigate the influence of other dietary or environmental modifiers on the lung microbiome, such as dietary fibre and dust from houses with dogs, which both change the gut microbiome, and have protective effects on the expression of allergic airway disease in mice [64, 68]. We identified a negative association between circulating 25(OH)D and OTUs of an unidentified Pseudomonas, which could be most clinically relevant for patients where Pseudomonas plays a central role in disease progression, such as for those with cystic fibrosis or COPD [69], and suggest that vitamin D status (or circulating 25(OH)D) may reduce the presence of pathobionts in the lungs. However, further studies in humans are required to reproduce the preclinical findings reported here.

Abbreviations

1,25(OH)2D:

1,25-dihydroxyvitamin D

25(OH)D:

25-hydroxyvitamin D

Alum:

Aluminium hydroxide suspension

ANOSIM:

Analysis of similarity

BALF:

Bronchoalveolar lavage fluid

COPD:

Chronic obstructive pulmonary disease

HDM:

House dust mite

mBD2:

Mouse ß-defensin-2

NMDS:

Non-metric multidimensional scaling

OTU:

Operational taxonomic units

OVA:

Ovalbumin

Th2:

T helper type-2

References

  1. Muehleisen B, Gallo RL. Vitamin D in allergic disease: shedding light on a complex problem. J Allergy Clin Immunol. 2013;131:324–9.

    Article  CAS  PubMed  Google Scholar 

  2. Weiss ST. Bacterial components plus vitamin D: the ultimate solution to the asthma (autoimmune disease) epidemic? J Allergy Clin Immunol. 2011;127:1128–30.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Lucas RM, Gorman S, Geldenhuys S, Hart PH. Vitamin D and immunity. F1000 Prime Rep. 2014;6:118.

    Article  Google Scholar 

  4. Weiss ST, Litonjua AA. Vitamin D, the gut microbiome, and the hygiene hypothesis. How does asthma begin? Am J Respir Crit Care Med. 2015;191:492–3.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Lagishetty V, Misharin AV, Liu NQ, Lisse TS, Chun RF, Ouyang Y, McLachlan SM, Adams JS, Hewison M. Vitamin D deficiency in mice impairs colonic antibacterial activity and predisposes to colitis. Endocrinology. 2010;151:2423–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Assa A, Vong L, Pinnell LJ, Avitzur N, Johnson-Henry KC, Sherman PM. Vitamin D deficiency promotes epithelial barrier dysfunction and intestinal inflammation. J Infect Dis. 2014;210:1296–305.

    Article  CAS  PubMed  Google Scholar 

  7. Jahani R, Fielding KA, Chen J, Villa CR, Castelli LM, Ward WE, Comelli EM. Low vitamin D status throughout life results in an inflammatory prone status but does not alter bone mineral or strength in healthy 3-month-old CD-1 male mice. Mol Nutr Food Res. 2014;58:1491–501.

    Article  CAS  PubMed  Google Scholar 

  8. Ooi JH, Li Y, Rogers CJ, Cantorna MT. Vitamin D regulates the gut microbiome and protects mice from dextran sodium sulfate-induced colitis. J Nutr. 2013;143:1679–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Bashir M, Prietl B, Tauschmann M, Mautner SI, Kump PK, Treiber G, Wurm P, Gorkiewicz G, Hogenauer C, Pieber TR. Effects of high doses of vitamin D on mucosa-associated gut microbiome vary between regions of the human gastrointestinal tract. Eur J Nutr. 2015. doi:10.1007/s00394-015-0966-2.

    PubMed  PubMed Central  Google Scholar 

  10. Ciubotaru I, Green SJ, Kukreja S, Barengolts E. Significant differences in fecal microbiota are associated with various stages of glucose tolerance in African American male veterans. Transl Res. 2015;166:401–11.

    Article  CAS  PubMed  Google Scholar 

  11. Lagishetty V, Liu NQ, Hewison M. Vitamin D metabolism and innate immunity. Mol Cell Endocrinol. 2011;347:97–105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Scott CL, Aumeunier AM, Mowat AM. Intestinal CD103+ dendritic cells: master regulators of tolerance? Trends Immunol. 2011;32:412–9.

    Article  CAS  PubMed  Google Scholar 

  13. Kelly D, Delday MI, Mulder I. Microbes and microbial effector molecules in treatment of inflammatory disorders. Immunol Rev. 2012;245:27–44.

    Article  CAS  PubMed  Google Scholar 

  14. Kim KS, Hong SW, Han D, Yi J, Jung J, Yang BG, Lee JY, Lee M, Surh CD. Dietary antigens limit mucosal immunity by inducing regulatory T cells in the small intestine. Science. 2016;351:858–63.

    Article  CAS  PubMed  Google Scholar 

  15. Hilty M, Burke C, Pedro H, Cardenas P, Bush A, Bossley C, Davies J, Ervine A, Poulter L, Pachter L, et al. Disordered microbial communities in asthmatic airways. PLoS One. 2010;5:e8578.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Beck JM, Young VB, Huffnagle GB. The microbiome of the lung. Transl Res. 2012;160:258–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Barfod KK, Roggenbuck M, Hansen LH, Schjorring S, Larsen ST, Sorensen SJ, Krogfelt KA. The murine lung microbiome in relation to the intestinal and vaginal bacterial communities. BMC Microbiol. 2013;13:303.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Gollwitzer ES, Saglani S, Trompette A, Yadava K, Sherburn R, McCoy KD, Nicod LP, Lloyd CM, Marsland BJ. Lung microbiota promotes tolerance to allergens in neonates via PD-L1. Nat Med. 2014;20:642–7.

    Article  CAS  PubMed  Google Scholar 

  19. Huang YJ, Nelson CE, Brodie EL, Desantis TZ, Baek MS, Liu J, Woyke T, Allgaier M, Bristow J, Wiener-Kronish JP, et al. Airway microbiota and bronchial hyperresponsiveness in patients with suboptimally controlled asthma. J Allergy Clin Immunol. 2011;127:372–81. e371-373.

    Article  PubMed  Google Scholar 

  20. Sze MA, Dimitriu PA, Hayashi S, Elliott WM, McDonough JE, Gosselink JV, Cooper J, Sin DD, Mohn WW, Hogg JC. The lung tissue microbiome in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2012;185:1073–80.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Denner DR, Sangwan N, Becker JB, Hogarth DK, Oldham J, Castillo J, Sperling AI, Solway J, Naureckas ET, Gilbert JA, White SR. Corticosteroid therapy and airflow obstruction influence the bronchial microbiome, which is distinct from that of bronchoalveolar lavage in asthmatic airways. J Allergy Clin Immunol. 2015. doi:10.1016/j.jaci.2015.10.017.

    PubMed  Google Scholar 

  22. Huang YJ, Nariya S, Harris JM, Lynch SV, Choy DF, Arron JR, Boushey H. The airway microbiome in patients with severe asthma: Associations with disease features and severity. J Allergy Clin Immunol. 2015;136:874–84.

    Article  PubMed  Google Scholar 

  23. Simpson JL, Daly J, Baines KJ, Yang IA, Upham JW, Reynolds PN, Hodge S, James AL, Hugenholtz P, Willner D, Gibson PG. Airway dysbiosis: Haemophilus influenzae and Tropheryma in poorly controlled asthma. Eur Respir J. 2016;47:792–800.

    Article  PubMed  Google Scholar 

  24. Gorman S, Weeden CE, Tan DH, Scott NM, Hart J, Foong RE, Mok D, Stephens N, Zosky G, Hart PH. Reversible control by vitamin D of granulocytes and bacteria in the lungs of mice: an ovalbumin-induced model of allergic airway disease. PLoS One. 2013;8:e67823.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Gorman S, Tan DH, Lambert MJ, Scott NM, Judge MA, Hart PH. Vitamin D(3) deficiency enhances allergen-induced lymphocyte responses in a mouse model of allergic airway disease. Pediatr Allergy Immunol. 2012;23:83–7.

    Article  PubMed  Google Scholar 

  26. Nadkarni MA, Martin FE, Jacques NA, Hunter N. Determination of bacterial load by real-time PCR using a broad-range (universal) probe and primers set. Microbiology. 2002;148:257–66.

    Article  CAS  PubMed  Google Scholar 

  27. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Quince C, Lanzen A, Davenport RJ, Turnbaugh PJ. Removing noise from pyrosequenced amplicons. BMC Bioinformatics. 2011;12:38.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Liu Z, DeSantis TZ, Andersen GL, Knight R. Accurate taxonomy assignments from 16S rRNA sequences produced by highly parallel pyrosequencers. Nucleic Acids Res. 2008;36:e120.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Paulson JN, Stine OC, Bravo HC, Pop M. Differential abundance analysis for microbial marker-gene surveys. Nat Methods. 2013;10:1200–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. McGlade JP, Gorman S, Zosky GR, Larcombe AN, Sly PD, Finlay-Jones JJ, Turner DJ, Hart PH. Suppression of the asthmatic phenotype by ultraviolet B-induced, antigen-specific regulatory cells. Clin Exp Allergy. 2007;37:1267–76.

    Article  CAS  PubMed  Google Scholar 

  33. Weiss S, Amir A, Hyde ER, Metcalf JL, Song SJ, Knight R. Tracking down the sources of experimental contamination in microbiome studies. Genome Biol. 2014;15:564.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Salter SJ, Cox MJ, Turek EM, Calus ST, Cookson WO, Moffatt MF, Turner P, Parkhill J, Loman NJ, Walker AW. Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol. 2014;12:87.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Oksanen JF, Blanchet G, Kindt R, Legendre P, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H. Vegan: Community Ecology Package. 2015.

    Google Scholar 

  36. Clarke KR. Non-parametric multivariate analysis of changes in community structure. Aust J Ecol. 1993;18:117–43.

    Article  Google Scholar 

  37. Barfod KK, Vrankx K, Mirsepasi-Lauridsen HC, Hansen JS, Hougaard KS, Larsen ST, Ouwenhand AC, Krogfelt KA. The Murine Lung Microbiome Changes During Lung Inflammation and Intranasal Vancomycin Treatment. Open Microbiol J. 2015;9:167–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Yu G, Gail MH, Consonni D, Carugno M, Humphrys M, Pesatori AC, Caporaso NE, Goedert JJ, Ravel J, Landi MT. Characterizing human lung tissue microbiota and its relationship to epidemiological and clinical features. Genome Biol. 2016;17:163.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Debarry J, Hanuszkiewicz A, Stein K, Holst O, Heine H. The allergy-protective properties of Acinetobacter lwoffii F78 are imparted by its lipopolysaccharide. Allergy. 2010;65:690–7.

    Article  CAS  PubMed  Google Scholar 

  40. Hanski I, von Hertzen L, Fyhrquist N, Koskinen K, Torppa K, Laatikainen T, Karisola P, Auvinen P, Paulin L, Makela MJ, et al. Environmental biodiversity, human microbiota, and allergy are interrelated. Proc Natl Acad Sci U S A. 2012;109:8334–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Lohmann P, Luna RA, Hollister EB, Devaraj S, Mistretta TA, Welty SE, Versalovic J. The airway microbiome of intubated premature infants: characteristics and changes that predict the development of bronchopulmonary dysplasia. Pediatr Res. 2014;76:294–301.

    Article  CAS  PubMed  Google Scholar 

  42. Gorman S, Scott NM, Tan DH, Weeden CE, Tuckey RC, Bisley JL, Grimbaldeston MA, Hart PH. Acute erythemal ultraviolet radiation causes systemic immunosuppression in the absence of increased 25-hydroxyvitamin D3 levels in male mice. PLoS One. 2012;7:e46006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Foong RE, Shaw NC, Berry LJ, Hart PH, Gorman S, Zosky GR. Vitamin D deficiency causes airway hyperresponsiveness, increases airway smooth muscle mass, and reduces TGF-beta expression in the lungs of female BALB/c mice. Physiol Rep. 2014;2:e00276.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Nashold FE, Spach KM, Spanier JA, Hayes CE. Estrogen controls vitamin D3-mediated resistance to experimental autoimmune encephalomyelitis by controlling vitamin D3 metabolism and receptor expression. J Immunol. 2009;183:3672–81.

    Article  CAS  PubMed  Google Scholar 

  45. Correale J, Ysrraelit MC, Gaitan MI. Gender differences in 1,25 dihydroxyvitamin D3 immunomodulatory effects in multiple sclerosis patients and healthy subjects. J Immunol. 2010;185:4948–58.

    Article  CAS  PubMed  Google Scholar 

  46. Walecki M, Eisel F, Klug J, Baal N, Paradowska-Dogan A, Wahle E, Hackstein H, Meinhardt A, Fijak M. Androgen receptor modulates Foxp3 expression in CD4+ CD25+ Foxp3+ regulatory T-cells. Mol Biol Cell. 2015;26:2845–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Gorman S, Geldenhuys S, Judge MA, Weeden CE, Waithman J, Hart PH. Dietary vitamin D increases percentages and function of regulatory T cells in the skin-draining lymph nodes and suppresses dermal inflammation. J Immunol Res. 2016. In press.

  48. Morrison GM, Davidson DJ, Dorin JR. A novel mouse beta defensin, Defb2, which is upregulated in the airways by lipopolysaccharide. FEBS Lett. 1999;442:112–6.

    Article  CAS  PubMed  Google Scholar 

  49. Wu M, McClellan SA, Barrett RP, Hazlett LD. Beta-defensin-2 promotes resistance against infection with P. aeruginosa. J Immunol. 2009;182:1609–16.

    Article  CAS  PubMed  Google Scholar 

  50. Biragyn A, Ruffini PA, Leifer CA, Klyushnenkova E, Shakhov A, Chertov O, Shirakawa AK, Farber JM, Segal DM, Oppenheim JJ, Kwak LW. Toll-like receptor 4-dependent activation of dendritic cells by beta-defensin 2. Science. 2002;298:1025–9.

    Article  CAS  PubMed  Google Scholar 

  51. Lagishetty V, Chun RF, Liu NQ, Lisse TS, Adams JS, Hewison M. 1alpha-hydroxylase and innate immune responses to 25-hydroxyvitamin D in colonic cell lines. J Steroid Biochem Mol Biol. 2010;121:228–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Guo S, Wu LX, Jones CX, Chen L, Hao CL, He L, Zhang JH. Allergic airway inflammation disrupts interleukin-17 mediated host defense against streptococcus pneumoniae infection. Int Immunopharmacol. 2016;31:32–8.

    Article  CAS  PubMed  Google Scholar 

  53. Beisswenger C, Kandler K, Hess C, Garn H, Felgentreff K, Wegmann M, Renz H, Vogelmeier C, Bals R. Allergic airway inflammation inhibits pulmonary antibacterial host defense. J Immunol. 2006;177:1833–7.

    Article  CAS  PubMed  Google Scholar 

  54. Chalmers JD, McHugh BJ, Docherty C, Govan JR, Hill AT. Vitamin-D deficiency is associated with chronic bacterial colonisation and disease severity in bronchiectasis. Thorax. 2013;68:39–47.

    Article  PubMed  Google Scholar 

  55. McCauley LA, Thomas W, Laguna TA, Regelmann WE, Moran A, Polgreen LE. Vitamin D deficiency is associated with pulmonary exacerbations in children with cystic fibrosis. Ann Am Thorac Soc. 2014;11:198–204.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Yim S, Dhawan P, Ragunath C, Christakos S, Diamond G. Induction of cathelicidin in normal and CF bronchial epithelial cells by 1,25-dihydroxyvitamin D(3). J Cyst Fibros. 2007;6:403–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Marsland BJ, Gollwitzer ES. Host-microorganism interactions in lung diseases. Nat Rev Immunol. 2014;14:827–35.

    Article  CAS  PubMed  Google Scholar 

  58. Yun Y, Srinivas G, Kuenzel S, Linnenbrink M, Alnahas S, Bruce KD, Steinhoff U, Baines JF, Schaible UE. Environmentally determined differences in the murine lung microbiota and their relation to alveolar architecture. PLoS One. 2014;9:e113466.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Cardenas PA, Cooper PJ, Cox MJ, Chico M, Arias C, Moffatt MF, Cookson WO. Upper airways microbiota in antibiotic-naive wheezing and healthy infants from the tropics of rural Ecuador. PLoS One. 2012;7:e46803.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Teo SM, Mok D, Pham K, Kusel M, Serralha M, Troy N, Holt BJ, Hales BJ, Walker ML, Hollams E, et al. The infant nasopharyngeal microbiome impacts severity of lower respiratory infection and risk of asthma development. Cell Host Microbe. 2015;17:704–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Martin RA, Hodgkins SR, Dixon AE, Poynter ME. Aligning mouse models of asthma to human endotypes of disease. Respirology. 2014;19:823–33.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Kumar RK, Herbert C, Foster PS. Mouse models of acute exacerbations of allergic asthma. Respirology. 2016;21:842–9.

    Article  PubMed  Google Scholar 

  63. McAlees JW, Whitehead GS, Harley IT, Cappelletti M, Rewerts CL, Holdcroft AM, Divanovic S, Wills-Karp M, Finkelman FD, Karp CL, Cook DN. Distinct Tlr4-expressing cell compartments control neutrophilic and eosinophilic airway inflammation. Mucosal Immunol. 2015;8:863–73.

    Article  CAS  PubMed  Google Scholar 

  64. Vital M, Harkema JR, Rizzo M, Tiedje J, Brandenberger C. Alterations of the Murine Gut Microbiome with Age and Allergic Airway Disease. J Immunol Res. 2015;2015:892568.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Foong RE, Bosco A, Troy NM, Gorman S, Hart PH, Kicic A, Zosky GR. Identification of genes differentially regulated by vitamin D deficiency that alter lung pathophysiology and inflammation in allergic airways disease. Am J Physiol Lung Cell Mol Physiol. 2016. doi:10.1152/ajplung.00026.2016.

    PubMed  Google Scholar 

  66. Heine G, Tabeling C, Hartmann B, Gonzalez Calera CR, Kuhl AA, Lindner J, Radbruch A, Witzenrath M, Worm M. 25-hydroxvitamin D3 promotes the long-term effect of specific immunotherapy in a murine allergy model. J Immunol. 2014;193:1017–23.

    Article  CAS  PubMed  Google Scholar 

  67. Vasiliou JE, Lui S, Walker SA, Chohan V, Xystrakis E, Bush A, Hawrylowicz CM, Saglani S, Lloyd CM. Vitamin D deficiency induces Th2 skewing and eosinophilia in neonatal allergic airways disease. Allergy. 2014;69:1380–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Trompette A, Gollwitzer ES, Yadava K, Sichelstiel AK, Sprenger N, Ngom-Bru C, Blanchard C, Junt T, Nicod LP, Harris NL, Marsland BJ. Gut microbiota metabolism of dietary fiber influences allergic airway disease and hematopoiesis. Nat Med. 2014;20:159–66.

    Article  CAS  PubMed  Google Scholar 

  69. Cullen L, McClean S. Bacterial Adaptation during Chronic Respiratory Infections. Pathogens. 2015;4:66–89.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This research was supported by: the Asthma Foundation of Western Australia, BrightSpark Foundation, Health Department of Western Australia, Raine Medical Research Foundation, Telethon Kids Institute, the University of Western Australia, Innovation Fund Denmark, The National Research Centre for the Working Environment, and the Staten Serum Institut (Dupont). None of these funding bodies had any role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

Please contact author for data requests.

Competing interests

Prof Feelisch is a member of the Scientific Advisory Board of AOBiome LLC, a company commercializing ammonia-oxidizing bacteria for use in inflammatory skin disease. However, this membership is not related to the subject of this manuscript. We have no further disclosures or competing interests to declare.

Authors’ contributions

SG conceived and designed this study with input from all authors. MR acquired and analysed the microbiome data for the study with help from DA, KB and SG. All authors have contributed towards the interpretation of findings from this study, have played a role in drafting the article or revising it critically for its intellectual content and have given their final approval for this version of the paper to be published.

Consent for publication

Not applicable.

Ethics approval and consent to participate

All experiments were performed according to the ethical guidelines of the National Health and Medical Research Council of Australia and with approval from the Telethon Kids Institute Animal Ethics Committee (AEC#238).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Shelley Gorman.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Roggenbuck, M., Anderson, D., Barfod, K.K. et al. Vitamin D and allergic airway disease shape the murine lung microbiome in a sex-specific manner. Respir Res 17, 116 (2016). https://doi.org/10.1186/s12931-016-0435-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12931-016-0435-3

Keywords