Next Article in Journal
DFT Studies of Selected Epoxies with Mesogenic Units–Impact of Molecular Structure on Electro-Optical Response
Next Article in Special Issue
Molecular and Immunological Identification of Low Allergenic Fruits among Old and New Apple Varieties
Previous Article in Journal
Diabetes and Cardiovascular Risk in Renal Transplant Patients
Previous Article in Special Issue
T and B Lymphocyte Transcriptional States Differentiate between Sensitized and Unsensitized Individuals in Alpha-Gal Syndrome
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 7
10.3390/ijms22073417
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Raw Milk-Induced Protection against Food Allergic Symptoms in Mice Is Accompanied by Shifts in Microbial Community Structure

by
Suzanne Abbring
1,
Phillip A. Engen
2,
Ankur Naqib
2,
Stefan J. Green
3,
Johan Garssen
1,4,
Ali Keshavarzian
1,2,5 and
Betty C. A. M. van Esch
1,4,*
1
Division of Pharmacology, Faculty of Science, Utrecht Institute for Pharmaceutical Sciences, Utrecht University, 3584 CG Utrecht, The Netherlands
2
Rush Center for Integrated Microbiome and Chronobiology Research, Rush Medical College, Rush University Medical Center, Chicago, IL 60612, USA
3
Genomics and Microbiome Core Facility, Rush University Medical Center, Chicago, IL 60612, USA
4
Danone Nutricia Research, 3584 CT Utrecht, The Netherlands
5
Departments of Medicine and Physiology, Rush University Medical Center, Chicago, IL 60612, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(7), 3417; https://doi.org/10.3390/ijms22073417
Submission received: 27 February 2021 / Revised: 15 March 2021 / Accepted: 23 March 2021 / Published: 26 March 2021
(This article belongs to the Special Issue Molecular Mechanisms of Allergy and Asthma)
Browse Figures
Review Reports Versions Notes

Abstract

:
The mechanism underlying the allergy-protective effects of raw cow’s milk is still unknown, but the modulation of the gut microbiome may play a role. The effects of consuming raw cow’s milk or processed milk on fecal microbial communities were therefore characterized in an experimental murine model. C3H/HeOuJ mice were treated with raw milk, pasteurized milk, skimmed raw milk, pasteurized milk supplemented with alkaline phosphatase (ALP), or phosphate-buffered saline (PBS) for eight days prior to sensitization and challenge with ovalbumin (OVA). Fecal samples were collected after milk exposure and after OVA sensitization, and microbiomes were characterized using 16S ribosomal RNA gene amplicon sequencing. Treatment with raw milk prior to OVA sensitization increased the relative abundance of putative butyrate-producing bacteria from the taxa Lachnospiraceae UCG-001, Lachnospiraceae UCG-008, and Ruminiclostridium 5 (Clostridial clusters XIVa and IV), while it decreased the relative abundance of Proteobacterial genera such as Parasutterella, a putative pro-inflammatory bacterial genus. This effect was observed after eight days of raw milk exposure and became more pronounced five weeks later, after allergic sensitization in the absence of milk. Similar trends were observed after treatment with skimmed raw milk. Conversely, the feeding of pasteurized milk led to a loss of allergy protection and a putative dysbiotic microbiome. The addition of ALP to pasteurized milk restored the protective effect observed with raw milk and mitigated some of the microbial community alterations associated with milk pasteurization. Raw milk-induced protection against food allergic symptoms in mice is accompanied by an increased relative abundance of putative butyrate-producing Clostridiales and a decreased relative abundance of putative pro-inflammatory Proteobacteria. Given the safety concerns regarding raw milk consumption, this knowledge is key for the development of new, microbiologically safe, preventive strategies to reduce the incidence of allergic diseases.

1. Introduction

The prevalence of allergic diseases such as allergic rhinitis, atopic dermatitis, asthma, and food allergy has increased dramatically in recent decades and has become a major health problem, particularly in affluent and fast-developing societies [1]. Due to the severe impact on the patient’s quality of life and extensive healthcare costs, this vast prevalence has major socio-economic consequences [2]. Even though therapies for allergic diseases have improved over the years with the introduction of allergen-specific immunotherapy and the development of new pharmacological agents aimed at combating inflammatory responses, as well as providing symptomatic relief, these therapies have remained largely noncurative [3,4,5].
Changes in environmental exposures and lifestyle practices such as diet, microbiome, hygiene, pollutants, exercise, and infections are believed to be associated with the increase in and severity of allergic diseases [1]. Numerous studies have emphasized the interaction of such factors with the immune system leading to modified immune responses [6,7,8]. A well-known example is the “farm effect” wherein children growing up on a farm have a lower risk of developing allergic diseases compared to children who are raised in the same rural area but not on a farm [9,10,11,12,13]. This beneficial “farm effect” has been mainly attributed to the consumption of raw, unprocessed cow’s milk [14,15,16,17,18,19]. Since the allergy-protective raw milk effect also has been observed in children not exposed to a farming environment [14,15,16], raw cow’s milk consumption could provide a potential natural solution for allergic diseases.
In previous work, we strengthened the epidemiological findings for an asthma- and allergy-protective effect of raw cow’s milk consumption by demonstrating a causal relationship using preclinical murine models [20,21]. In addition, we showed that the heat-sensitive components, not the fat content, are responsible for the raw milk-induced protection against food allergic symptoms [22]. These studies have increased our knowledge about the protective raw milk constituents, although the bioactive component(s) remain to be elucidated.
Current knowledge on the mechanism(s) underlying the allergy-protective effects of raw cow’s milk is still limited. It is hypothesized that the many immunomodulatory components present in raw milk create a regulatory environment favoring unresponsiveness upon allergen exposure [23]. Potential mechanisms of action are improved intestinal barrier function, the promotion of regulatory T-cell development, and modulation of the gut microbiome [23,24,25].
The gut microbial composition and functionality is influenced by dietary elements and can in turn influence the host immune response with effects on asthma and allergy development [26]. Interestingly, several components present in raw cow’s milk may modulate the gut microbiome. Proteins such as lactoferrin, lactoperoxidase, and lysozyme have antimicrobial activity, and milk oligosaccharides such as sialyllactose can promote the growth of probiotic Bifidobacteria that in turn produce short-chain fatty acids (SCFA), known to prevent the development of allergic diseases [23,26,27]. In addition, raw cow’s milk has been shown to contain significantly higher abundances of viable microbial cells compared to heat-treated milk [18]. Nevertheless, compelling evidence for a raw milk-induced change in the gut microbiome is still lacking.
In the current study, we therefore analyzed the gut microbiota of mice treated with raw milk and then sensitized and challenged with ovalbumin (OVA) to induce food allergy. These mice were previously shown to be protected against food allergic symptoms by raw milk [22]. We also previously demonstrated that the suppression of food allergic symptoms by raw milk was retained after skimming but eliminated after pasteurization of the milk. Supplementing pasteurized milk with alkaline phosphatase (ALP), one of the first bioactive raw milk components that loses its activity upon heat treatment, restored the protective effect [22]. In the current study, we compared fecal microbial communities of mice undergoing different milk treatments to assess the effects on microbial community structure and to assess whether changes in the gut microbiota contribute to observed allergy-protective effects.

2. Results

2.1. Microbiome Sequencing

Fecal samples (n = 100: day −9 (n = 8); day −1 (phosphate-buffered saline (PBS) n = 14, raw, pasteurized, skimmed, pasteurized + ALP n = 8 per group); and day 31 (PBS n = 6, OVA, raw, pasteurized, skimmed, pasteurized + ALP n = 8 per group); Figure 1) were processed for microbiome sequencing. A total of 4,099,786 sequences were generated, with an average depth of sequencing of 40,998 per sample (median = 41,629; min = 7061; max = 63,657). Sequence data were rarefied to a depth of 7000 sequences per sample for diversity analyses (Supplementary Materials Figure S1). At the start of the study (day −9), before tolerance induction, the mouse fecal samples (n = 8) microbial communities were similar, with a Bray–Curtis within-group similarity index average of 73.48 ± 4.318 SD, and baseline samples had Bray–Curtis similarity values greater than 70%.

2.2. Microbial Communities Are Highly Similar between Treatments Directly after Milk Exposure

Fecal samples were collected after the tolerance induction period, directly after PBS/milk treatment (day −1; Figure 1) and assessed for differences in microbial community structure. Microbial communities were broadly similar at the day −1 sampling and were dominated by bacteria from the phyla Firmicutes and Bacteroidetes, and bacteria from the genera Muribaculaceae, Lactobacillus, Lachnospiraceae, Bacteroides, and Roseburia (>75% of all sequences at day −1; Figure 2A,B). With the exception of genus-level richness, which was increased in pasteurized milk-treated mice relative to raw milk-treated mice (one-way ANOVA: F(4,41) = 2.702, p = 0.0436; Tukey post hoc: p = 0.028), no significant differences in alpha diversity metrics (i.e., Shannon index, Simpson’s index, richness, and evenness) between groups were observed (Supplementary Materials Figure S2).
In addition, beta diversity analyses (i.e., analysis of similarity (ANOSIM)) were performed to determine if bulk microbial communities were significantly different. At the day −1 sampling, the observed microbial community structures of most treatment groups were not significantly different (ANOSIM R-values from −0.098 to 0.16; p-values from 0.057 to 0.897). However, the fecal microbial community structure of the raw milk treatment group was significantly different from the pasteurized milk group (R = 0.286, p = 0.003) and from the skimmed raw milk group (R = 0.288, p = 0.010; Table 1).
Features analyses were conducted at the taxonomic levels of phylum, family, and genus to identify microbial taxa that were significantly differently abundant between treatment groups at the day −1 sampling. Kruskal–Wallis tests were performed on centered log ratio transformed (CLR) relative abundance data, and false discovery rate (FDR)-corrected p-values and Dunn’s multiple group comparison’s post hoc test adjusted p-values are reported (Supplementary Materials Table S1). No significant differences (FDR p-values < 0.05) in the relative abundance of individual taxa were observed between the treatment groups at the day −1 sampling (Supplementary Materials Table S1).
We further assessed the relative abundance of putative butyrate-producing taxa together, including multiple Lachnospiraceae, Roseburia, and Blautia genera. No significant differences in the combined relative abundance of these taxa between treatment groups were observed at the day −1 sampling (Supplementary Materials Figure S3).

2.3. Significant Differences in Microbial Community Structure Primarily Manifested Themselves after Sensitization

At day 31, after the last sensitization, mouse fecal microbiomes were characterized again. Microbial communities were broadly similar again at the day 31 sampling and dominated by bacteria from the phyla Firmicutes and Bacteroidetes, and bacteria from the taxa Muribaculaceae, Lactobacillus, Lachnospiraceae, Bacteroides, and Roseburia (>75% of all sequences at day 31; Figure 2C,D). However, across all groups, the relative abundance of Bacteroidetes increased relative to day −1 (Figure 2A,C), as did the relative abundances of genera Muribaculaceae and Lactobacillus (Figure 2B,D). At the day 31 sampling, no significant differences in fecal microbiome Shannon index, Simpson’s index, and evenness between groups were observed (Supplementary Materials Figure S4A–C). However, genus-level richness values varied significantly between groups, with OVA-sensitized allergic control mice (OVA group) and pasteurized milk-treated mice having significantly higher richness relative to PBS-sensitized control mice (PBS group) and raw milk-, skimmed raw milk-, and pasteurized milk + ALP-treated mice (Figure 3; one-way ANOVA: F(4,41) = 15.19, p ˂ 0.0001).
Beta-diversity analyses (i.e., ANOSIM) demonstrated that at the day 31 sampling, the fecal microbial communities of raw milk- and skimmed raw milk-treated mice were no longer significantly different (R = 0.122, p = 0.106), while the differences between raw milk and pasteurized milk persisted (R = 0.226, p = 0.013; Table 1). Microbial community structure in all other treatment groups were significantly different from each other (Table 1; R-values from 0.159 to 0.610; p-values from 0.002 to 0.049).
At day 31, nine microbial features were trending as different (FDR-p ≤ 0.07) or significantly different (FDR-p < 0.05) between groups, including the phylum Proteobacteria, family Burkholderiaceae (Proteobacteria), genus Oscillibacter (Firmicutes), genus Ruminiclostridium (Firmicutes), genus Parasutterella (Proteobacteria), and multiple Lachnospiraceae genera (Blautia, NK4A136 group, UCG-001, and UCG-008; Firmicutes; Table 2). The relative abundances of Gram-negative, putative pro-inflammatory phylum Proteobacteria (i.e., Burkholderiaceae and Parasutterella) were significantly higher in pasteurized milk-treated mice and pasteurized milk + ALP-treated mice compared to PBS-treated allergic control mice (OVA group) and raw milk-treated mice (Table 2). The relative abundance of bacteria from the Lachnospiraceae NK4A136 group (putative butyrate producers) was significantly increased in skimmed raw milk-treated mice compared to PBS-treated allergic control mice (OVA group; Table 2).
The relative abundances of Blautia, Oscillibacter, and Lachnospiraceae UCG-008 were increased in OVA-sensitized allergic control mice compared to PBS-sensitized control mice (Table 2). Treating mice with raw milk or skimmed raw milk prior to OVA sensitization significantly increased the relative abundances of Lachnospiraceae UCG-001, Lachnospiraceae UCG-008, and Ruminiclostridium 5 compared to PBS-treated allergic control mice (OVA group; Table 2). Bacteria from the Lachnospiraceae UCG-001 group were elevated in most milk treatments, but not in the pasteurized milk treatment (Table 2). Of the nine microbial trending and significant features, Lachnospiraceae UCG-001 is the only feature significantly different by Dunn’s post hoc test between mice in the pasteurized milk treatment and mice in the pasteurized milk + ALP treatment (Table 2).
The percent relative abundance of putative butyrate-producing genera (Kruskal–Wallis: p = 0.0055), and the ratio of total butyrate-to-total SCFA (Kruskal–Wallis: p = 0.0012) were significantly different across treatment groups at day 31 (Figure 4A,B). Specifically, skimmed raw milk-treated mice had higher relative abundances of putative butyrate producers when compared to PBS-treated allergic control mice (Butyrate (Dunns: p = 0.020); Ratio (Dunns: p = 0.013); Figure 4A,B).
A machine-learning algorithm (Boruta) was used to further identify taxa differentiating microbial communities by treatment at the day 31 sampling. This analysis identified the genera Parasutterella and Bifidobacterium were particularly driving the differences between groups observed at day 31 (Supplementary Materials Figure S5). Parasutterella were identified by using a taxon-by-taxon Kruskal–Wallis test (above), while Bifidobacterium were not significant (FDR-p = 0.162; Supplementary Materials Table S2) by this analysis. We note, however, that Bifidobacterium levels were extremely low in all treatment groups except for OVA (PBS-treated allergic mice) and pasteurized milk treatments (Supplementary Materials Table S2).

3. Discussion

We previously showed that raw, unprocessed cow’s milk suppressed allergic symptoms in a murine model for food allergy [22]. This protective effect was retained after skimming, but lost upon pasteurization of the milk, indicating that heat-sensitive components could underlie the allergy-protective effects of raw cow’s milk. Additionally, we showed that ALP, a heat-sensitive bioactive raw milk component, was able to partially restore the allergy-protective effect lost with heat treatment. The observed protection against food allergic symptoms was accompanied by a reduction in allergen-specific Th2 responsiveness, and an induction of tolerance-associated cell types, such as CD103+ dendritic cells (DCs) and regulatory T-cells (Tregs). Moreover, SCFA concentrations in the cecum hinted towards a potential immune modulation via the gut microbiome [22]. In the present study, we therefore interrogated mouse fecal samples to determine if changes in the microbial community structure are induced by raw milk treatment and to determine the consequences of milk processing steps on the fecal microbiome.
In this study, analysis of the fecal microbiome of mice indicated that raw milk exposure increased the relative abundance of several putative butyrate-producing bacterial taxa from the order Clostridiales, including Lachnospiraceae UCG-001, Lachnospiraceae UCG-008, and Ruminiclostridium 5, while it decreased the relative abundance of some Proteobacterial taxa (e.g., Parasutterella). This effect was visible after eight days of raw milk exposure (i.e., after the tolerance induction period) but became more pronounced five weeks later, after allergic sensitization in the absence of milk. Similar to animals treated with raw cow’s milk, animals treated with skimmed raw cow’s milk prior to OVA sensitization had higher relative abundances of putative anti-inflammatory bacteria and lower relative abundances of pro-inflammatory bacteria. Moreover, the relative abundance of all identified putative butyrate-producing genera was higher in raw milk and skimmed raw milk treatments relative to all other samples, and significantly higher in the skimmed raw milk treatment relative to PBS and OVA treatments. Mice treated with pasteurized milk had relatively high relative abundances of putative pathobiont inflammatory bacteria and low relative abundances of putative anti-inflammatory bacteria, and this may represent a dysbiotic microbiome. The addition of ALP to pasteurized milk increased the relative abundances of anti-inflammatory bacteria, but did not reduce pro-inflammatory bacteria levels. Remarkably, allergic sensitization to OVA increased bacterial richness compared to sham sensitization. This finding was not observed in any of the other groups, except in the group treated with pasteurized milk prior to OVA sensitization.
Raw cow’s milk consumption has been reported to protect against childhood asthma and allergies [16,17,18,19,20]. If the risk of infections could be overcome, consumption of raw milk would be an attractive preventive strategy for allergic diseases. Further research into the underlying mechanisms is therefore highly desirable. One of the hypothesized working mechanisms of raw cow’s milk is modulation of the gut microbiome and its metabolic activity (i.e., SCFA production) [23,24,25]. Recently, an observational study was conducted in which the impact of unpasteurized dairy products on the human gut microbiome was investigated in participants undertaking a residential 12-week cookery course on an organic farm [28]. This study was one of the first to show that the intake of unpasteurized milk and dairy products increased the relative abundance of putative probiotic bacteria [28]. However, due to the observational nature of this study, authors were unable to disentangle the impact of diet and geographical environment on the gut microbiome.
In our study, fecal samples were collected after milk exposure and after sensitization to OVA. We investigated whether the different milk types induced changes in the microbiota composition, and whether these changes could be associated with the allergy-protective effects previously observed after raw milk, skimmed raw milk, and pasteurized milk + ALP exposure [22]. One of the most striking findings was that raw and skimmed raw milk exposure increased the relative abundance of genera belonging to the butyrate-producing Clostridiales clusters XIVa and IV (informal nomenclature: Clostridium clusters), including the taxa Lachnospiraceae UCG-001, Lachnospiraceae UCG-008, and Ruminiclostridium 5 [29,30,31]. Both raw milk types showed the strongest increase in relative abundance of these bacteria at day 31 (after sensitization), five weeks after the last milk exposure. This could indicate that bacteria present in raw milk needed time to colonize in the gut or that raw milk components promoting the growth of probiotic bacteria, such as sialyllactose [23,32] and lactoferrin [33,34,35], needed time to exert their effect.
Bacteria within the Clostridium clusters XIVa and IV are known for their high level of metabolite production. They can produce acetic, propionic, and butyric acids through the degradation of non-digestible oligosaccharides. The organisms typically reside in close proximity to the host epithelium, and they can exert a strong influence on the host immune system [36,37]. For example, the colonization of mice with a defined mixture of Clostridial strains belonging to clusters XIVa and IV induced the accumulation and differentiation of colonic FoxP3+ Treg cells [38]. The presence of Clostridium induced the release of active TGF-β and other Treg-inducing factors from colonic epithelial cells and thereby presumably conditioned CD103+ DCs to drive Treg differentiation. Furthermore, Clostridium species have also been shown to induce IL-10-expressing Tregs in the colon [38]. Importantly, the oral inoculation of a mixture of these species resulted in resistance to allergy, as demonstrated by reduced Th2 responsiveness in the spleen and lower systemic IgE responses [38].
Although we did not assess Treg numbers in the colon, we did observe an increase in TGF-β- and IL-10-producing Treg cells in lymphoid organs after raw milk treatment [22,39]. In addition, raw milk exposure induced tolerogenic CD103+ DCs [22]. Since it is well known that Treg cells travel between compartments [40] and that bacterial metabolites such as SCFAs can enter the systemic circulation [41], we speculate that the induction of SCFA-producing Clostridiales by raw milk underlies the observed increase in tolerance-associated cell types in lymphoid organs which may have contributed to the observed allergy-protective effect. This hypothesis is supported by the estimation of the total relative abundance of all identified putative butyrate-producing taxa in mouse feces. These measurements indicated a significantly higher relative abundance of total butyrate-producing genera in skimmed raw milk-treated mice compared to PBS-treated allergic control mice at the day 31 sampling. For raw milk-treated mice, a similar pattern was observed although the difference was not significant at the p < 0.05 level (p = 0.233). Interestingly, bacterial communities with the potential of producing butyrate were also found to be associated with the protective farm effect on childhood asthma in the PASTURE birth cohort study [42].
The current literature particularly speculates about a raw milk-induced outgrowth of Bifidobacteria and Lactobacilli [23,24,25]. In line with this hypothesis, the previously mentioned observational study demonstrated an increase in the relative abundance of Lactobacilli following the consumption of unpasteurized milk and dairy products [28]. However, in our study, the relative abundance of both probiotic bacterial genera was low in raw milk-treated mice. Our findings are consistent with another study which showed that the relative abundance of Bifidobacterium members was lower in the feces of infants raised in farming versus nonfarming environments [43]. Solid evidence is therefore still lacking on the role of these putative probiotic bacterial taxa in allergic responses. Strikingly, in our study, the relative abundances of Bifidobacteria and Lactobacilli were highest in PBS-treated allergic control mice and pasteurized milk-treated mice, while these mice showed the strongest allergic symptoms.
Another remarkable observation in this study is of the higher microbial richness in PBS-treated allergic control mice (i.e., OVA group) and pasteurized milk-treated mice when compared to all other groups at the day 31 sampling. A high microbial diversity is usually linked to improved gut health [44], and several studies have found an association between low microbial diversity and allergic diseases such as eczema and allergic rhinitis [45]. Elsewhere, bacterial diversity was found to be increased in sputum samples of asthmatic patients, compared to non-asthmatic subjects [46,47]. These results do not suggest a straightforward relationship between fecal microbial diversity and allergic responses.
Besides increasing the relative abundance of anti-inflammatory bacteria, raw milk treatment also decreased the relative abundance of pro-inflammatory bacteria. Mice in the raw milk treatments had both higher relative abundance of putative anti-inflammatory bacteria and lower relative abundance of putative pro-inflammatory bacteria, especially relative to animals in the pasteurized milk treatment. In particular, Proteobacteria levels were greatly reduced in raw milk treatments, relative to levels in pasteurized milk-treated mice. Proteobacteria are generally recognized as the microbial signature of gut dysbiosis, a condition that has been characterized by a microbial imbalance leading to metabolic disorders and intestinal inflammation [48]. In this study, levels of bacteria from the putative pathobiont inflammatory bacterial genus Parasutterella were most strongly associated with pasteurization. Parasutterella levels have been previously associated with food allergy in infancy [49]. Although we did not observe an increase in the relative abundance of Parasutterella in OVA-sensitized allergic control mice compared to PBS-sensitized control mice, the relative abundance of Parasutterella did increase significantly in pasteurized milk-treated mice compared to raw milk-treated mice.
The addition of ALP to pasteurized milk was previously shown to partially restore the allergy-protective effect lost upon heat treatment [22]. Interestingly, the fecal microbiota community structure of pasteurized milk + ALP-treated mice also differed significantly from mice treated with pasteurized milk alone. By adding ALP to pasteurized milk, the relative abundance of anti-inflammatory bacteria (e.g., Lachnospiraceae UCG-001) increased, and even reached similar bacterial levels as observed in raw milk-treated mice. Spiking pasteurized milk with ALP did not, however, lower the relative abundance of Proteobacteria. Since mice treated with pasteurized milk + ALP were still protected against food allergic symptoms, these results suggest that inducing butyrate-producing, anti-inflammatory bacteria is more important than reducing pro-inflammatory bacteria levels in preventing allergic diseases.
In summary, we demonstrated that the raw milk-induced protection against food allergic symptoms coincides with an increased relative abundance of putative butyrate-producing genera and a decreased relative abundance of Proteobacterial taxa in the feces of mice long after raw milk exposure. Similar microbial responses to raw milk and skimmed raw milk were observed. The previously observed loss of allergy protection after pasteurization was accompanied by a putative dysbiosis-associated microbial phenotype that could be partly restored by the addition of ALP. The effect of the origin of the milk (e.g., organic vs. conventional) on the allergy-protective capacity and associated microbiota changes needs to be assessed in future studies.
By showing a raw milk-induced modulation of the gut microbiome, this study contributes to the current knowledge on mechanisms underlying the allergy-protective effects of raw cow’s milk. Since the consumption of raw cow’s milk is not recommended due to the potential presence of pathogens, this knowledge is crucial to promote the development of microbiologically safe alternatives.

4. Materials and Methods

4.1. Mice

Three-week-old, specific-pathogen-free, female C3H/HeOuJ mice (Charles River Laboratories, Sulzfeld, Germany) were housed at the animal facility of Utrecht University (Utrecht, The Netherlands) in filter-topped Makrolon cages (one cage/group, n = 6–8/cage) with standard chip bedding, Kleenex tissues, and a plastic shelter on a 12 h light/dark cycle with access to food (“Rat and Mouse Breeder and Grower Expanded”; Special Diet Services, Witham, UK) and water ad libitum. Upon arrival, mice were randomly assigned to the control or experimental groups and were acclimatized to the laboratory conditions for one week prior to experimentation. All animal procedures were approved by the Ethical Committee for Animal Research of the Utrecht University and complied with the European Directive 2010/63/EU on the protection of animals used for scientific purposes (AVD1008002015346).

4.2. Milk Types

Raw, unprocessed cow’s milk was collected from a dairy farm (Macroom, Ireland) and subsequently divided into three aliquots. Aliquot 1 was stored directly at −20 °C with no further treatment (raw milk). Aliquot 2 was heated at 78 °C for 15 s, cooled on ice to 4 °C and then stored at −20 °C (pasteurized milk) and aliquot 3 was skimmed using a centrifugal separator at 55 °C to remove milk fat and then cooled on ice to 4 °C and stored at −20 °C (skimmed milk; 0.1% fat). All milk types were solely produced for experimental purposes (Danone Nutricia Research, Macroom, Ireland). Before use, milk samples were thawed at room temperature and a portion of the pasteurized milk was supplemented with bovine intestinal ALP (pasteurized milk + ALP; 3 units/0.5 mL pasteurized milk; kindly provided by Prof. Dr. W. Seinen (Alloksys Life Sciences BV/AMRIF BV, Wageningen, The Netherlands)).

4.3. Experimental Design—Tolerance Induction, Sensitization and Challenges

The experimental design of the study is illustrated in Figure 1. Mice were orally treated by gavage once daily with 0.5 mL raw milk, pasteurized milk, skimmed raw milk, pasteurized milk supplemented with ALP, or PBS (as a control) for eight consecutive days (experimental days −9 to −2) after acclimatization. After this tolerance induction period, mice were orally sensitized with 20 mg of the chicken egg protein OVA (grade V; Sigma–Aldrich, Zwijndrecht, The Netherlands) dissolved in 0.5 mL of PBS containing 10 µg of cholera toxin (CT; List Biological Laboratories, Campbell, CA, USA) as an adjuvant (days 0, 7, 14, 21, and 28; n = 8/group). PBS-sensitized control mice (n = 6) received CT alone (10 µg/0.5 mL of PBS). On days −9, −1, and 31 (before tolerance induction, after tolerance induction, and after sensitization), fecal samples were collected for microbiota analysis. On day 33, five days after the last sensitization, all mice were intradermally and orally challenged with OVA (10 µg of OVA/20 µL of PBS and 50 mg of OVA/0.5 mL of PBS, respectively) to assess the allergic response (results previously published [22]). Sixteen hours after the oral challenge (day 34), mice were killed by cervical dislocation.

4.4. Microbiota Profiling and Bioinformatics Analysis

Total genomic DNA was extracted from mouse fecal samples utilizing the QIAamp DNA Microbiome Kit (QIAGEN, Germantown, MD, USA), and extract DNA concentrations were measured using fluorometric quantitation (Qubit, Life Technologies, Grand Island, NY, USA). Primers 515F/806R (515F: GTGYCAGCMGCCGCGGTAA; 806R: GGACTACNVGGGTWTCTAAT), targeting the V4 variable region of microbial small subunit (SSU or 16S) ribosomal RNA (rRNA) genes, were used for PCR (e.g., Caporaso et al. [50]), and prepared for high-throughput amplicon sequencing using a two-stage PCR method, as described previously [51]. Negative controls were used with each set of amplifications, which indicated no background contamination. Sequencing was performed using an Illumina MiSeq (Illumina, San Diego, CA, USA), with a V2 kit and paired-end 250 base reads at the University of Illinois at Chicago (UIC) Genome Research Core (GRC). Raw sequence data (FASTQ files) were deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under the BioProject identifier PRJNA576858.
Basic processing of the raw data was performed by the Research Informatics Core (RIC) at UIC. Raw FASTQ forward and reverse read files for each sample were merged using the software package PEAR (Paired-End-reAd mergeR; v0.9.8; (http://www.exelixis-lab.org/web/software/pear (accessed on 1 February 2021)) [52,53]. Sequences were then processed using the software package QIIME, v1.9 [54]. Briefly, merged reads were quality trimmed and primer sequences were removed. Sequences shorter than 225 base pairs were discarded. Chimeric sequences were identified and removed using the USEARCH algorithm [55], along with the annotation being performed against the SILVA (silva_132_16S.97 [56]) reference database at a similarity. A biological observation matrix (BIOM) [57] was generated at each taxonomic level from phylum to species. Downstream statistical analyses were performed using the software packages Primer7, as well as vegan (https://cran.r-project.org/package=vegan (accessed on 1 February 2021)) and boruta (http://www.jstatsoft.org/v36/i11/ (accessed on 1 February 2021)), both within the R programming environment [58,59,60,61].

4.5. Statistics

Analyses of alpha diversity and beta diversity were used to examine changes in fecal microbial community structure between animals and between experimental groups. Alpha diversity indices (i.e., Shannon, Simpson, richness, and Pielou’s evenness) were generated using the package “vegan” implemented in the R programming environment (https://cran.r-project.org (accessed on 1 February 2021), https://github.com/vegandevs/vegan accessed date: 30 December 2019). Alpha diversity indices were calculated on datasets rarefied to 7000 sequences per sample [62]. Differences in alpha diversity indices between groups were assessed for significance using one-way analysis of variance (ANOVA) tests with Tukey’s post hoc test adjusted p-values reported.
To examine differences in microbial community structure between samples, a pairwise Bray–Curtis dissimilarity (non-phylogenetic) metric was generated using the Primer7 software package (PRIMER-E Ltd., Lutton, United Kingdom), and analysis of similarity (ANOSIM) calculations were performed on pair-wise distance matrices to determine if differences in microbial community structure between treatments were significant. ANOSIM was performed at the taxonomic level of genus on square-root transformed data. Before tolerance induction (day −9), the fecal microbial community structure was examined using within-group Bray–Curtis similarity scores. Random forest (RF) models were used to predict allergy status based on the genus-level bacterial profiles using the default parameters of the Boruta algorithm (“randomForest” package within R; (http://www.jstatsoft.org/v36/i11/ (accessed on 1 February 2021)) [61]. Differences in the relative abundance of individual taxa across treatment groups were assessed for significance using the Kruskal–Wallis test on centered log ratio transformed (CLR) data, with false discovery rate (FDR)-corrected p-values, plus Dunn’s multiple group comparison’s post hoc test. A predictive model depicting the relative abundances of butyrate-producing taxa and total butyrate-to-total SCFA ratio was examined in the mice groups, as previously described [63,64]. Figures were created using GraphPad Prism (v8.00; GraphPad Software Inc, La Jolla, San Jose, CA, USA).

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms22073417/s1, Figure S1: Rarefaction curves of alpha diversity metrics, Figure S2: Pasteurized milk-treated mice had increased bacterial richness compared to raw milk-treated mice at day −1, Figure S3: No differences in butyrate-producing genera between treatment groups were observed at day −1, Figure S4: Alpha diversity indices in control and experimental groups at day 31, Figure S5: Taxonomic feature selection using a machine learning random forest model (Boruta) at day 31, Table S1: Average total count and percent relative abundance of sequences of individual fecal taxa between the different groups at day −1, Table S2: Average total count and percent relative abundance of sequences of individual fecal taxa between the different groups at day 31.

Author Contributions

Conceptualization, S.A., J.G. and B.C.A.M.v.E.; methodology; S.A., P.A.E., A.N., S.J.G., A.K. and B.C.A.M.v.E.; formal analysis, P.A.E., A.N. and S.J.G.; investigation, S.A., P.A.E., A.N. and S.J.G.; writing—original draft preparation, S.A.; writing—review and editing, P.A.E., A.N., S.J.G., J.G., A.K. and B.C.A.M.v.E.; visualization, S.A., P.A.E., A.N. and S.J.G.; project administration, S.A., P.A.E. and B.C.A.M.v.E.; funding acquisition, J.G. and B.C.A.M.v.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by Danone Nutricia Research.

Institutional Review Board Statement

All animal procedures are approved by the Ethical Committee for Animal Research of the Utrecht University and complied with the European Directive 2010/63/EU on the protection of animals used for scientific purposes (AVD1008002015346).

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this published article (and its supplementary information files). Additional datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. Raw sequence data are deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive under the BioProject identifier PRJNA576858.

Acknowledgments

The authors would like to thank W. Seinen (Alloksys Life Sciences BV/AMRIF BV, Wageningen, The Netherlands) for kindly providing bovine intestinal alkaline phosphatase.

Conflicts of Interest

J.G. and B.C.A.M.v.E. are (partly) employed at Danone Nutricia Research. All other authors declare no conflict of interest.

References

  1. Platts-Mills, T.A. The allergy epidemics: 1870–2010. J. Allergy Clin. Immunol. 2015, 136, 3–13. [Google Scholar] [CrossRef] [Green Version]
  2. Zuberbier, T.; Lotvall, J.; Simoens, S.; Subramanian, S.V.; Church, M.K. Economic burden of inadequate management of allergic diseases in the European Union: A GA(2)LEN review. Allergy 2014, 69, 1275–1279. [Google Scholar] [CrossRef]
  3. Galli, S.J.; Tsai, M.; Piliponsky, A.M. The development of allergic inflammation. Nature 2008, 454, 445–454. [Google Scholar] [CrossRef] [Green Version]
  4. Barnes, P.J. New therapies for asthma: Is there any progress? Trends Pharmacol. Sci. 2010, 31, 335–343. [Google Scholar] [CrossRef]
  5. Burks, A.W.; Sampson, H.A.; Plaut, M.; Lack, G.; Akdis, C.A. Treatment for food allergy. J. Allergy Clin. Immunol. 2018, 141, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Hong, X.; Wang, X. Early life precursors, epigenetics, and the development of food allergy. Semin. Immunopathol. 2012, 34, 655–669. [Google Scholar] [CrossRef] [Green Version]
  7. Van Esch, B.; Porbahaie, M.; Abbring, S.; Garssen, J.; Potaczek, D.P.; Savelkoul, H.F.J.; van Neerven, R.J.J. The impact of milk and its components on epigenetic programming of immune function in early life and beyond: Implications for allergy and asthma. Front. Immunol. 2020, 11, 2141. [Google Scholar] [CrossRef]
  8. Acevedo, N.; Alashkar Alhamwe, B.; Caraballo, L.; Ding, M.; Ferrante, A.; Garn, H.; Garssen, J.; Hii, C.S.; Irvine, J.; Llinas-Caballero, K.; et al. Perinatal and early-life nutrition, epigenetics, and allergy. Nutrients 2021, 13, 724. [Google Scholar] [CrossRef]
  9. Braun-Fahrlander, C.; Gassner, M.; Grize, L.; Neu, U.; Sennhauser, F.H.; Varonier, H.S.; Vuille, J.C.; Wuthrich, B. Prevalence of hay fever and allergic sensitization in farmer’s children and their peers living in the same rural community. SCARPOL team. Swiss Study on Childhood Allergy and Respiratory Symptoms with Respect to Air Pollution. Clin. Exp. Allergy 1999, 29, 28–34. [Google Scholar] [CrossRef]
  10. Kilpelainen, M.; Terho, E.O.; Helenius, H.; Koskenvuo, M. Farm environment in childhood prevents the development of allergies. Clin. Exp. Allergy 2000, 30, 201–208. [Google Scholar] [CrossRef] [PubMed]
  11. Von Ehrenstein, O.S.; von Mutius, E.; Illi, S.; Baumann, L.; Bohm, O.; von Kries, R. Reduced risk of hay fever and asthma among children of farmers. Clin. Exp. Allergy 2000, 30, 187–193. [Google Scholar] [CrossRef]
  12. Ernst, P.; Cormier, Y. Relative scarcity of asthma and atopy among rural adolescents raised on a farm. Am. J. Respir. Crit. Care Med. 2000, 161, 1563–1566. [Google Scholar] [CrossRef]
  13. Riedler, J.; Eder, W.; Oberfeld, G.; Schreuer, M. Austrian children living on a farm have less hay fever, asthma and allergic sensitization. Clin. Exp. Allergy 2000, 30, 194–200. [Google Scholar] [CrossRef] [PubMed]
  14. Riedler, J.; Braun-Fahrlander, C.; Eder, W.; Schreuer, M.; Waser, M.; Maisch, S.; Carr, D.; Schierl, R.; Nowak, D.; von Mutius, E.; et al. Exposure to farming in early life and development of asthma and allergy: A cross-sectional survey. Lancet 2001, 358, 1129–1133. [Google Scholar] [CrossRef]
  15. Waser, M.; Michels, K.B.; Bieli, C.; Floistrup, H.; Pershagen, G.; von Mutius, E.; Ege, M.; Riedler, J.; Schram-Bijkerk, D.; Brunekreef, B.; et al. Inverse association of farm milk consumption with asthma and allergy in rural and suburban populations across Europe. Clin. Exp. Allergy 2007, 37, 661–670. [Google Scholar] [CrossRef]
  16. Perkin, M.R.; Strachan, D.P. Which aspects of the farming lifestyle explain the inverse association with childhood allergy? J. Allergy Clin. Immunol. 2006, 117, 1374–1381. [Google Scholar] [CrossRef]
  17. Ege, M.J.; Frei, R.; Bieli, C.; Schram-Bijkerk, D.; Waser, M.; Benz, M.R.; Weiss, G.; Nyberg, F.; van Hage, M.; Pershagen, G.; et al. Not all farming environments protect against the development of asthma and wheeze in children. J. Allergy Clin. Immunol. 2007, 119, 1140–1147. [Google Scholar] [CrossRef] [PubMed]
  18. Loss, G.; Apprich, S.; Waser, M.; Kneifel, W.; Genuneit, J.; Buchele, G.; Weber, J.; Sozanska, B.; Danielewicz, H.; Horak, E.; et al. The protective effect of farm milk consumption on childhood asthma and atopy: The GABRIELA study. J. Allergy Clin. Immunol. 2011, 128, 766–773. [Google Scholar] [CrossRef] [PubMed]
  19. Brick, T.; Schober, Y.; Bocking, C.; Pekkanen, J.; Genuneit, J.; Loss, G.; Dalphin, J.C.; Riedler, J.; Lauener, R.; Nockher, W.A.; et al. Omega-3 fatty acids contribute to the asthma-protective effect of unprocessed cow’s milk. J. Allergy Clin. Immunol. 2016, 137, 1699–1706. [Google Scholar] [CrossRef] [Green Version]
  20. Abbring, S.; Verheijden, K.A.T.; Diks, M.A.P.; Leusink-Muis, A.; Hols, G.; Baars, T.; Garssen, J.; van Esch, B.C.A.M. Raw cow’s milk prevents the development of airway inflammation in a murine house dust mite-induced asthma model. Front. Immunol. 2017, 8, 1045. [Google Scholar] [CrossRef] [Green Version]
  21. Abbring, S.; Wolf, J.; Ayechu Muruzabal, V.; Diks, M.A.P.; Alashkar Alhamwe, B.; Alhamdan, F.; Harb, H.; Renz, H.; Garn, H.; Garssen, J.; et al. Raw cow’s milk reduces allergic symptoms in a murine model for food allergy—A potential role for epigenetic modifications. Nutrients 2019, 11, 1721. [Google Scholar] [CrossRef] [Green Version]
  22. Abbring, S.; Ryan, J.T.; Diks, M.A.P.; Hols, G.; Garssen, J.; van Esch, B.C.A.M. Suppression of food allergic symptoms by raw cow’s milk in mice is retained after skimming but abolished after heating the milk—A promising contribution of alkaline phosphatase. Nutrients 2019, 11, 1499. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Van Neerven, R.J.; Knol, E.F.; Heck, J.M.; Savelkoul, H.F. Which factors in raw cow’s milk contribute to protection against allergies? J. Allergy Clin. Immunol. 2012, 130, 853–858. [Google Scholar] [CrossRef]
  24. Abbring, S.; Hols, G.; Garssen, J.; van Esch, B.C.A.M. Raw cow’s milk consumption and allergic diseases—The potential role of bioactive whey proteins. Eur. J. Pharmacol. 2019, 843, 55–65. [Google Scholar] [CrossRef]
  25. Sozanska, B. Raw cow’s milk and its protective effect on allergies and asthma. Nutrients 2019, 11, 469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. McKenzie, C.; Tan, J.; Macia, L.; Mackay, C.R. The nutrition-gut microbiome-physiology axis and allergic diseases. Immunol. Rev. 2017, 278, 277–295. [Google Scholar] [CrossRef]
  27. Tan, J.; McKenzie, C.; Potamitis, M.; Thorburn, A.N.; Mackay, C.R.; Macia, L. Chapter three—The role of short-chain fatty acids in health and disease. In Advances in Immunology; Alt, F.W., Ed.; Academic Press: Cambridge, MA, USA, 2014; Volume 121, pp. 91–119. [Google Scholar]
  28. Butler, M.I.; Bastiaanssen, T.F.S.; Long-Smith, C.; Berding, K.; Morkl, S.; Cusack, A.M.; Strain, C.; Busca, K.; Porteous-Allen, P.; Claesson, M.J.; et al. Recipe for a healthy gut: Intake of unpasteurised milk is associated with increased Lactobacillus abundance in the human gut microbiome. Nutrients 2020, 12, 1468. [Google Scholar] [CrossRef] [PubMed]
  29. Duncan, S.H.; Barcenilla, A.; Stewart, C.S.; Pryde, S.E.; Flint, H.J. Acetate utilization and butyryl coenzyme A (CoA):acetate-CoA transferase in butyrate-producing bacteria from the human large intestine. Appl. Environ. Microbiol. 2002, 68, 5186–5190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Barcenilla, A.; Pryde, S.E.; Martin, J.C.; Duncan, S.H.; Stewart, C.S.; Henderson, C.; Flint, H.J. Phylogenetic relationships of butyrate-producing bacteria from the human gut. Appl. Environ. Microbiol. 2000, 66, 1654–1661. [Google Scholar] [CrossRef] [Green Version]
  31. Biddle, A.; Stewart, L.; Blanchard, J.; Leschine, S. Untangling the Genetic Basis of Fibrolytic Specialization by Lachnospiraceae and Ruminococcaceae in Diverse Gut Communities. Diversity 2013, 5, 627–640. [Google Scholar] [CrossRef]
  32. Perdijk, O.; van Baarlen, P.; Fernandez-Gutierrez, M.M.; van den Brink, E.; Schuren, F.H.J.; Brugman, S.; Savelkoul, H.F.J.; Kleerebezem, M.; van Neerven, R.J.J. Sialyllactose and galactooligosaccharides promote epithelial barrier functioning and distinctly modulate microbiota composition and short chain fatty acid production in vitro. Front. Immunol. 2019, 10, 94. [Google Scholar] [CrossRef] [Green Version]
  33. Giansanti, F.; Panella, G.; Leboffe, L.; Antonini, G. Lactoferrin from milk: Nutraceutical and pharmacological properties. Pharmaceuticals 2016, 9, 61. [Google Scholar] [CrossRef] [Green Version]
  34. Oda, H.; Wakabayashi, H.; Yamauchi, K.; Abe, F. Lactoferrin and bifidobacteria. Biometals 2014, 27, 915–922. [Google Scholar] [CrossRef] [PubMed]
  35. Tian, H.; Maddox, I.S.; Ferguson, L.R.; Shu, Q. Influence of bovine lactoferrin on selected probiotic bacteria and intestinal pathogens. Biometals 2010, 23, 593–596. [Google Scholar] [CrossRef]
  36. Nagano, Y.; Itoh, K.; Honda, K. The induction of Treg cells by gut-indigenous Clostridium. Curr. Opin. Immunol. 2012, 24, 392–397. [Google Scholar] [CrossRef]
  37. Lopetuso, L.R.; Scaldaferri, F.; Petito, V.; Gasbarrini, A. Commensal Clostridia: Leading players in the maintenance of gut homeostasis. Gut Pathog. 2013, 5, 23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Atarashi, K.; Tanoue, T.; Shima, T.; Imaoka, A.; Kuwahara, T.; Momose, Y.; Cheng, G.; Yamasaki, S.; Saito, T.; Ohba, Y.; et al. Induction of colonic regulatory T cells by indigenous Clostridium species. Science 2011, 331, 337–341. [Google Scholar] [CrossRef] [Green Version]
  39. Abbring, S.; Xiong, L.; Diks, M.A.P.; Baars, T.; Garssen, J.; Hettinga, K.; van Esch, B. Loss of allergy-protective capacity of raw cow’s milk after heat treatment coincides with loss of immunologically active whey proteins. Food Funct. 2020, 11, 4982–4993. [Google Scholar] [CrossRef]
  40. Siegmund, K.; Feuerer, M.; Siewert, C.; Ghani, S.; Haubold, U.; Dankof, A.; Krenn, V.; Schon, M.P.; Scheffold, A.; Lowe, J.B.; et al. Migration matters: Regulatory T-cell compartmentalization determines suppressive activity in vivo. Blood 2005, 106, 3097–3104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Koh, A.; De Vadder, F.; Kovatcheva-Datchary, P.; Backhed, F. From dietary fiber to host physiology: Short-chain fatty acids as key bacterial metabolites. Cell 2016, 165, 1332–1345. [Google Scholar] [CrossRef] [Green Version]
  42. Depner, M.; Taft, D.H.; Kirjavainen, P.V.; Kalanetra, K.M.; Karvonen, A.M.; Peschel, S.; Schmausser-Hechfellner, E.; Roduit, C.; Frei, R.; Lauener, R.; et al. Maturation of the gut microbiome during the first year of life contributes to the protective farm effect on childhood asthma. Nat. Med. 2020, 26, 1766–1775. [Google Scholar] [CrossRef] [PubMed]
  43. Thorsen, J.; McCauley, K.; Fadrosh, D.; Lynch, K.; Barnes, K.L.; Bendixsen, C.G.; Seroogy, C.M.; Lynch, S.V.; Gern, J.E. Evaluating the effects of farm exposure on infant gut microbiome. J. Allergy Clin. Immunol. 2019, 143, 299. [Google Scholar] [CrossRef] [Green Version]
  44. Fujimura, K.E.; Slusher, N.A.; Cabana, M.D.; Lynch, S.V. Role of the gut microbiota in defining human health. Expert Rev. Anti Infect. Ther. 2010, 8, 435–454. [Google Scholar] [CrossRef] [Green Version]
  45. Melli, L.C.; do Carmo-Rodrigues, M.S.; Araujo-Filho, H.B.; Sole, D.; de Morais, M.B. Intestinal microbiota and allergic diseases: A systematic review. Allergol. Immunopathol. 2016, 44, 177–188. [Google Scholar] [CrossRef]
  46. Marri, P.R.; Stern, D.A.; Wright, A.L.; Billheimer, D.; Martinez, F.D. Asthma-associated differences in microbial composition of induced sputum. J. Allergy Clin. Immunol. 2013, 131, 346–352. [Google Scholar] [CrossRef] [Green Version]
  47. Huang, Y.J.; Nelson, C.E.; Brodie, E.L.; Desantis, T.Z.; Baek, M.S.; Liu, J.; Woyke, T.; Allgaier, M.; Bristow, J.; Wiener-Kronish, J.P.; et al. Airway microbiota and bronchial hyperresponsiveness in patients with suboptimally controlled asthma. J. Allergy Clin. Immunol. 2011, 127, 372–381. [Google Scholar] [CrossRef] [Green Version]
  48. Shin, N.R.; Whon, T.W.; Bae, J.W. Proteobacteria: Microbial signature of dysbiosis in gut microbiota. Trends Biotechnol. 2015, 33, 496–503. [Google Scholar] [CrossRef]
  49. Rachid, R.A.; Gerber, G.; Li, N.; Umetsu, D.T.; Bry, L.; Chatila, T.A. Food allergy in infancy is associated with dysbiosis of the intestinal microbiota. J. Allergy Clin. Immunol. 2016, 137, 235. [Google Scholar] [CrossRef]
  50. Caporaso, J.G.; Lauber, C.L.; Walters, W.A.; Berg-Lyons, D.; Huntley, J.; Fierer, N.; Owens, S.M.; Betley, J.; Fraser, L.; Bauer, M.; et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012, 6, 1621–1624. [Google Scholar] [CrossRef] [Green Version]
  51. Naqib, A.; Poggi, S.; Wang, W.; Hyde, M.; Kunstman, K.; Green, S.J. Making and sequencing heavily multiplexed, high-throughput 16S ribosomal RNA gene amplicon libraries using a flexible, two-stage PCR protocol. Methods Mol. Biol. 2018, 1783, 149–169. [Google Scholar] [CrossRef] [PubMed]
  52. Zhang, J.; Kobert, K.; Flouri, T.; Stamatakis, A. PEAR: A fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics 2014, 30, 614–620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Schmieder, R.; Edwards, R. Quality control and preprocessing of metagenomic datasets. Bioinformatics 2011, 27, 863–864. [Google Scholar] [CrossRef] [Green Version]
  54. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Peña, A.G.; Goodrich, J.K.; Gordon, J.I.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [Green Version]
  56. Glockner, F.O.; Yilmaz, P.; Quast, C.; Gerken, J.; Beccati, A.; Ciuprina, A.; Bruns, G.; Yarza, P.; Peplies, J.; Westram, R.; et al. 25 years of serving the community with ribosomal RNA gene reference databases and tools. J. Biotechnol. 2017, 261, 169–176. [Google Scholar] [CrossRef]
  57. McDonald, D.; Clemente, J.C.; Kuczynski, J.; Rideout, J.R.; Stombaugh, J.; Wendel, D.; Wilke, A.; Huse, S.; Hufnagle, J.; Meyer, F.; et al. The Biological Observation Matrix (BIOM) format or: How I learned to stop worrying and love the ome-ome. Gigascience 2012, 1, 7. [Google Scholar] [CrossRef] [Green Version]
  58. R Development Core Team. R: A Language and Environment for Statistical Computing. Available online: https://www.R-project.org/ (accessed on 30 December 2019).
  59. Clarke, K.R. Non-parametric multivariate analyses of changes in community structure. Aust. J. Ecol. 1993, 18, 117–143. [Google Scholar] [CrossRef]
  60. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package. Available online: https://CRAN.R-project.org/package=vegan (accessed on 30 December 2019).
  61. Kursa, M.B.; Rudnicki, W.R. Feature selection with the Boruta package. J. Stat. Softw. 2010, 36, 1–13. [Google Scholar] [CrossRef] [Green Version]
  62. Gihring, T.M.; Green, S.J.; Schadt, C.W. Massively parallel rRNA gene sequencing exacerbates the potential for biased community diversity comparisons due to variable library sizes. Environ. Microbiol. 2012, 14, 285–290. [Google Scholar] [CrossRef] [PubMed]
  63. Bishehsari, F.; Engen, P.A.; Preite, N.Z.; Tuncil, Y.E.; Naqib, A.; Shaikh, M.; Rossi, M.; Wilber, S.; Green, S.J.; Hamaker, B.R.; et al. Dietary fiber treatment corrects the composition of gut microbiota, promotes SCFA production, and suppresses colon carcinogenesis. Genes 2018, 9, 102. [Google Scholar] [CrossRef] [Green Version]
  64. Lordan, C.; Thapa, D.; Ross, R.P.; Cotter, P.D. Potential for enriching next-generation health-promoting gut bacteria through prebiotics and other dietary components. Gut Microbes 2020, 11, 1–20. [Google Scholar] [CrossRef] [Green Version]
Ijms 22 03417 g001 550
Figure 1. Schematic overview of the study design. Female C3H/HeOuJ mice were grouped as depicted: PBS (PBS-sensitized control mice; n = 6), OVA (OVA-sensitized allergic mice; n = 8), raw milk (raw milk-treated mice; n = 8), pasteurized milk (pasteurized milk-treated mice; n = 8), skimmed milk (skimmed raw milk-treated mice; n = 8), and pasteurized milk + ALP (pasteurized milk + ALP-treated mice; n = 8). On experimental days 0, 7, 14, 21, and 28, mice were orally sensitized to the hen’s egg protein OVA using cholera toxin as an adjuvant (20 mg of OVA + 10 µg of CT/0.5 mL of PBS). PBS-sensitized control mice received cholera toxin alone. Prior to sensitization, mice were orally treated with 0.5 mL of PBS (as a control), raw milk, pasteurized milk, skimmed raw milk, or pasteurized milk supplemented with ALP for eight consecutive days (day −9 to −2). On days −9, −1, and 31 (before tolerance induction, after tolerance induction, and after sensitization), fecal samples were collected for microbiota analysis. On day 33, all mice were intradermally (10 µg of OVA/20 µL of PBS) and orally (50 mg of OVA/0.5 mL of PBS) challenged with OVA. Mice were killed on day 34 (as indicated by †). The allergic response in each group is displayed with a “+” for an allergic response or a “–” for no allergic response (results previously published [22]). PBS, phosphate-buffered saline; ALP, alkaline phosphatase; OVA, ovalbumin; CT, cholera toxin.
Figure 1. Schematic overview of the study design. Female C3H/HeOuJ mice were grouped as depicted: PBS (PBS-sensitized control mice; n = 6), OVA (OVA-sensitized allergic mice; n = 8), raw milk (raw milk-treated mice; n = 8), pasteurized milk (pasteurized milk-treated mice; n = 8), skimmed milk (skimmed raw milk-treated mice; n = 8), and pasteurized milk + ALP (pasteurized milk + ALP-treated mice; n = 8). On experimental days 0, 7, 14, 21, and 28, mice were orally sensitized to the hen’s egg protein OVA using cholera toxin as an adjuvant (20 mg of OVA + 10 µg of CT/0.5 mL of PBS). PBS-sensitized control mice received cholera toxin alone. Prior to sensitization, mice were orally treated with 0.5 mL of PBS (as a control), raw milk, pasteurized milk, skimmed raw milk, or pasteurized milk supplemented with ALP for eight consecutive days (day −9 to −2). On days −9, −1, and 31 (before tolerance induction, after tolerance induction, and after sensitization), fecal samples were collected for microbiota analysis. On day 33, all mice were intradermally (10 µg of OVA/20 µL of PBS) and orally (50 mg of OVA/0.5 mL of PBS) challenged with OVA. Mice were killed on day 34 (as indicated by †). The allergic response in each group is displayed with a “+” for an allergic response or a “–” for no allergic response (results previously published [22]). PBS, phosphate-buffered saline; ALP, alkaline phosphatase; OVA, ovalbumin; CT, cholera toxin.
Ijms 22 03417 g001
Ijms 22 03417 g002 550
Figure 2. Microbial profiles of the control and experimental groups at day −1 and day 31. (A) Day −1 stacked column plots depicting the percent mean relative abundance (>1%) of bacterial phyla. (B) Day −1 stacked column plots depicting the percent mean relative abundance (>0.1%) of bacterial genera. (C) Day 31 stacked column plots depicting the percent mean relative abundance (>1%) of bacterial phyla. (D) Day 31 stacked column plots depicting the percent mean relative abundance (>0.1%) of bacterial genera. PBS, phosphate-buffered saline; OVA, ovalbumin; raw, raw cow’s milk; pasteurized, pasteurized cow’s milk; skimmed, skimmed raw cow’s milk; pasteurized + ALP, pasteurized milk spiked with alkaline phosphatase.
Figure 2. Microbial profiles of the control and experimental groups at day −1 and day 31. (A) Day −1 stacked column plots depicting the percent mean relative abundance (>1%) of bacterial phyla. (B) Day −1 stacked column plots depicting the percent mean relative abundance (>0.1%) of bacterial genera. (C) Day 31 stacked column plots depicting the percent mean relative abundance (>1%) of bacterial phyla. (D) Day 31 stacked column plots depicting the percent mean relative abundance (>0.1%) of bacterial genera. PBS, phosphate-buffered saline; OVA, ovalbumin; raw, raw cow’s milk; pasteurized, pasteurized cow’s milk; skimmed, skimmed raw cow’s milk; pasteurized + ALP, pasteurized milk spiked with alkaline phosphatase.
Ijms 22 03417 g002
Ijms 22 03417 g003 550
Figure 3. Ovalbumin sensitization and pasteurized milk treatment significantly increased bacterial richness at day 31. Microbial richness was examined at the taxonomic level of genus. Measurements were based on 16S rRNA gene amplicon sequence datasets rarefied to 7000 sequences per sample. One-way ANOVA results and Tukey post hoc test adjusted p-values shown: * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. PBS, phosphate-buffered saline; OVA, ovalbumin; raw, raw cow’s milk; pasteurized, pasteurized cow’s milk; skimmed, skimmed raw cow’s milk; pasteurized + ALP, pasteurized milk spiked with alkaline phosphatase.
Figure 3. Ovalbumin sensitization and pasteurized milk treatment significantly increased bacterial richness at day 31. Microbial richness was examined at the taxonomic level of genus. Measurements were based on 16S rRNA gene amplicon sequence datasets rarefied to 7000 sequences per sample. One-way ANOVA results and Tukey post hoc test adjusted p-values shown: * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. PBS, phosphate-buffered saline; OVA, ovalbumin; raw, raw cow’s milk; pasteurized, pasteurized cow’s milk; skimmed, skimmed raw cow’s milk; pasteurized + ALP, pasteurized milk spiked with alkaline phosphatase.
Ijms 22 03417 g003
Ijms 22 03417 g004 550
Figure 4. Skimmed raw milk is associated with increased butyrate-producing genera at day 31. The percent mean relative abundance of (A) putative butyrate-producing taxa and (B) total butyrate-to-total SCFA ratio examined is shown for the different control and experimental groups at day 31. Kruskal–Wallis results and Dunn’s post hoc test adjusted p-values shown: * p < 0.05, ** p < 0.01. PBS, phosphate-buffered saline; OVA, ovalbumin; raw, raw cow’s milk; pasteurized, pasteurized cow’s milk; skimmed, skimmed raw cow’s milk; pasteurized + ALP, pasteurized milk spiked with alkaline phosphatase; KW, Kruskal–Wallis; SCFA, short-chain fatty acids.
Figure 4. Skimmed raw milk is associated with increased butyrate-producing genera at day 31. The percent mean relative abundance of (A) putative butyrate-producing taxa and (B) total butyrate-to-total SCFA ratio examined is shown for the different control and experimental groups at day 31. Kruskal–Wallis results and Dunn’s post hoc test adjusted p-values shown: * p < 0.05, ** p < 0.01. PBS, phosphate-buffered saline; OVA, ovalbumin; raw, raw cow’s milk; pasteurized, pasteurized cow’s milk; skimmed, skimmed raw cow’s milk; pasteurized + ALP, pasteurized milk spiked with alkaline phosphatase; KW, Kruskal–Wallis; SCFA, short-chain fatty acids.
Ijms 22 03417 g004
Table
Table 1. Group analysis of similarity (ANOSIM) results for overall fecal microbiota compositions between control and experimental groups at the taxonomic level of genus.
Table 1. Group analysis of similarity (ANOSIM) results for overall fecal microbiota compositions between control and experimental groups at the taxonomic level of genus.
Day −1Rp-Value
PBS vs. Raw0.0250.311
PBS vs. Pasteurized0.1600.057
PBS vs. Skimmed0.0730.157
PBS vs. Pasteurized + ALP−0.0980.897
Raw vs. Pasteurized0.2860.003
Raw vs. Skimmed0.2880.010
Raw vs. Pasteurized + ALP0.0980.114
Pasteurized vs. Skimmed−0.0750.854
Pasteurized vs. Pasteurized + ALP0.0830.151
Skimmed vs. Pasteurized + ALP−0.0180.541
Day 31
PBS vs. OVA0.3890.002
PBS vs. Raw0.3650.005
PBS vs. Pasteurized0.3880.004
PBS vs. Skimmed0.6100.002
PBS vs. Pasteurized + ALP0.4530.004
OVA vs. Raw0.3730.003
OVA vs. Pasteurized0.1740.026
OVA vs. Skimmed0.4640.003
OVA vs. Pasteurized + ALP0.3130.005
Raw vs. Pasteurized0.2660.013
Raw vs. Skimmed0.1220.106
Raw vs. Pasteurized + ALP0.1590.049
Pasteurized vs. Skimmed0.3260.008
Pasteurized vs. Pasteurized + ALP0.2140.038
Skimmed vs. Pasteurized + ALP0.2270.034
Global R comparison was based on ANOSIM performed within the software package Primer7 (PRIMER-E Ltd., Lutton, UK); p-values were calculated based on a permutational analysis vs. employing 999 permutations; square-root transformed. p-values < 0.05 are considered significant. PBS, phosphate-buffered saline; OVA, ovalbumin; raw, raw cow’s milk; pasteurized, pasteurized cow’s milk; skimmed, skimmed raw cow’s milk; pasteurized + ALP, pasteurized milk spiked with alkaline phosphatase.
Table
Table 2. Significant average total count and percent relative abundance of sequences of individual fecal taxa between the different groups at day 31.
Table 2. Significant average total count and percent relative abundance of sequences of individual fecal taxa between the different groups at day 31.
Taxonomic LevelsPBS 1
Count
(% RA)		OVA 2
Count
(% RA)		Raw 3
Count
(% RA)		Pasteurized 4 Count
(% RA)		Skimmed 5 Count
(% RA)		Pasteurized + ALP 6 Count
(% RA)		p-Value *FDR
p-Value *		Dunn’s Post Hoc Tests #
Phylum
Proteobacteria312.17
(0.76)		575.13
(1.33)		179.88
(0.44)		910.38
(2.06)		244.00
(0.60)		811.50
(1.93)		0.0000.0081 v 4, 1 v 6, 2 v 4, 2 v 6, 3 v 4, 3 v 6, 4 v 5, 5 v 6
Family
Burkholderiaceae280.33
(0.68)		170.63
(0.37)		55.13
(0.14)		786.25
(1.78)		132.13
(0.32)		657.63
(1.56)		0.0000.0051 v 4, 2 v 4, 2 v 6, 3 v 4, 3 v 6, 4 v 5, 5 v 6
Genus
Blautia22.00
(0.08)		160.38
(0.37)		216.00
(0.54)		388.88
(0.87)		99.25
(0.23)		121.13
(0.29)		0.0030.0661 v 2, 1 v 3, 1 v 4, 1 v 6, 3 v 5
Lachnospiraceae NK4A136 group1608.33
(5.87)		4399.75
(9.77)		6321.13
(15.59)		6045.00
(13.78)		9674.00
(24.60)		5868.88
(14.00)		0.0010.0381 v 3, 1 v 4, 1 v 5, 1 v 6, 2 v 5
Lachnospiraceae UCG-0011175.67
(4.40)		796.13
(1.80)		2478.63
(6.10)		787.38
(1.79)		2000.75
(5.06)		2200.00
(5.40)		0.0030.0662 v 3, 2 v 5, 2 v 6, 3 v 4, 4 v 5, 4 v 6
Lachnospiraceae UCG-00831.50
(0.13)		95.38
(0.21)		162.38
(0.40)		110.25
(0.25)		172.75
(0.42)		157.75
(0.38)		0.0020.0541 v 2, 1 v 3, 1 v 4, 1 v 5, 1 v 6, 2 v 3
Oscillibacter42.00
(0.15)		144.88
(0.32)		179.63
(0.43)		220.63
(0.50)		308.50
(0.78)		130.63
(0.32)		0.0030.0701 v 2, 1 v 3, 1 v 4, 1 v 5, 5 v 6
Parasutterella273.67
(0.66)		165.63
(0.36)		54.00
(0.14)		766.00
(1.73)		127.75
(0.31)		637.88
(1.51)		0.0000.0041 v 3, 2 v 4, 2 v 6, 3 v 4, 3 v 6, 4 v 5, 5 v 6
Ruminiclostridium 534.00
(0.13)		51.13
(0.12)		116.88
(0.28)		96.88
(0.22)		164.25
(0.41)		83.38
(0.20)		0.0020.0541 v 3, 1 v 4, 1 v 5, 2 v 3, 2 v 4, 2 v 5
Individual taxa differences between groups were assessed for significance using Kruskal–Wallis test on centered log ratio transformed (CLR) data, with false discovery rate (FDR)-corrected p-values reported. Adjusted p-values were considered significant at: * FDR-p ˂ 0.05; trends at * p ˂ 0.05 are reported. Dunn’s multiple group comparison’s post hoc tests significance indicated # p ˂ 0.05. 1 PBS, phosphate-buffered saline; 2 OVA, ovalbumin; 3 raw, raw cow’s milk; 4 pasteurized, pasteurized cow’s milk; 5 skimmed, skimmed raw cow’s milk; 6 pasteurized + ALP, pasteurized milk spiked with alkaline phosphatase; count = mean number of sequences in defined group; % RA, percent relative abundance in defined group; NS, non-significant.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abbring, S.; Engen, P.A.; Naqib, A.; Green, S.J.; Garssen, J.; Keshavarzian, A.; van Esch, B.C.A.M. Raw Milk-Induced Protection against Food Allergic Symptoms in Mice Is Accompanied by Shifts in Microbial Community Structure. Int. J. Mol. Sci. 2021, 22, 3417. https://doi.org/10.3390/ijms22073417

AMA Style

Abbring S, Engen PA, Naqib A, Green SJ, Garssen J, Keshavarzian A, van Esch BCAM. Raw Milk-Induced Protection against Food Allergic Symptoms in Mice Is Accompanied by Shifts in Microbial Community Structure. International Journal of Molecular Sciences. 2021; 22(7):3417. https://doi.org/10.3390/ijms22073417

Chicago/Turabian Style

Abbring, Suzanne, Phillip A. Engen, Ankur Naqib, Stefan J. Green, Johan Garssen, Ali Keshavarzian, and Betty C. A. M. van Esch. 2021. "Raw Milk-Induced Protection against Food Allergic Symptoms in Mice Is Accompanied by Shifts in Microbial Community Structure" International Journal of Molecular Sciences 22, no. 7: 3417. https://doi.org/10.3390/ijms22073417

APA Style

Abbring, S., Engen, P. A., Naqib, A., Green, S. J., Garssen, J., Keshavarzian, A., & van Esch, B. C. A. M. (2021). Raw Milk-Induced Protection against Food Allergic Symptoms in Mice Is Accompanied by Shifts in Microbial Community Structure. International Journal of Molecular Sciences, 22(7), 3417. https://doi.org/10.3390/ijms22073417

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop