Influence of Milk-Feeding Type and Genetic Risk of Developing Coeliac Disease on Intestinal Microbiota of Infants: The PROFICEL Study

Abstract

Interactions between environmental factors and predisposing genes could be involved in the development of coeliac disease (CD). This study has assessed whether milk-feeding type and HLA-genotype influence the intestinal microbiota composition of infants with a family history of CD. The study included 164 healthy newborns, with at least one first-degree relative with CD, classified according to their HLA-DQ genotype by PCR-SSP DQB1 and DQA1 typing. Faecal microbiota was analysed by quantitative PCR at 7 days, and at 1 and 4 months of age. Significant interactions between milk-feeding type and HLA-DQ genotype on bacterial numbers were not detected by applying a linear mixed-model analysis for repeated measures. In the whole population, breast-feeding promoted colonization of C. leptum group, B. longum and B. breve, while formula-feeding promoted that of Bacteroides fragilis group, C. coccoides-E. rectale group, E. coli and B. lactis. Moreover, increased numbers of B. fragilis group and Staphylococcus spp., and reduced numbers of Bifidobacterium spp. and B. longum were detected in infants with increased genetic risk of developing CD. Analyses within subgroups of either breast-fed or formula-fed infants indicated that in both cases increased risk of CD was associated with lower numbers of B. longum and/or Bifidobacterium spp. In addition, in breast-fed infants the increased genetic risk of developing CD was associated with increased C. leptum group numbers, while in formula-fed infants it was associated with increased Staphylococcus and B. fragilis group numbers. Overall, milk-feeding type in conjunction with HLA-DQ genotype play a role in establishing infants' gut microbiota; moreover, breast-feeding reduced the genotype-related differences in microbiota composition, which could partly explain the protective role attributed to breast milk in this disorder.

About the Authors:

Giada De Palma

Affiliation: Instituto de Agroquímica y Tecnología de Alimentos, Consejo Superior de Investigaciones Científicas (IATA-CSIC), Valencia, Spain

Amalia Capilla

Affiliation: Instituto de Biomedicina de Valencia (CSIC), CIBER de Enfermedades Raras (CIBERER), Valencia, Spain

Esther Nova

Affiliation: Department Metabolismo y Nutrición, ICTAN-CSIC, Madrid, Spain

Gemma Castillejo

Affiliation: Unidad de Gastroenterología Pediátrica, Hospital Universitario Sant Joan de Reus, Tarragona, Spain

Vicente Varea

Affiliation: Gastroenterología, Nutrición y Hepatología Pediátrica, Hospital Universitario Sant Joan de Deu and Unidad de Gastroenterología Pediátrica del Institut Dexeus, Barcelona, Spain

Tamara Pozo

Affiliation: Department Metabolismo y Nutrición, ICTAN-CSIC, Madrid, Spain

José Antonio Garrote

Affiliation: Unidad de Gastroenterología Pediátrica, Hospital Clínico Universitario de Valladolid, Valladolid, Spain

Isabel Polanco

Affiliation: Servicio de Gastroenterología y Nutrición Pediátrica, Hospital Universitario La Paz, Madrid, Spain

Ana López

Affiliation: Unidad de Gastroenterología Pediátrica, Hospital Universitario La Fe, Valencia, Spain

Carmen Ribes-Koninckx

Affiliation: Unidad de Gastroenterología Pediátrica, Hospital Universitario La Fe, Valencia, Spain

Ascensión Marcos

Affiliation: Department Metabolismo y Nutrición, ICTAN-CSIC, Madrid, Spain

María Dolores García-Novo

Affiliation: Unidad de Gastroenterología, Hospital Universitario Infantil Niño Jesús, Madrid, Spain

Carmen Calvo

Affiliation: Unidad de Gastroenterología Pediátrica, Hospital Clínico Universitario de Valladolid, Valladolid, Spain

Luis Ortigosa

Affiliation: Unidad de Gastroenterología, Hepatología y Nutrición Pediátrica, Hospital Universitario Nuestra Señora de Candelaria, Santa Cruz de Tenerife, Canarias, Spain

Luis Peña-Quintana

Affiliation: Unidad de Gastroenterología, Hepatología y Nutrición Pediátrica, Hospital Universitario Materno-Infantil de Canarias, Universidad de Las Palmas de Gran Canaria, Las Palmas de Gran Canaria, Spain

Francesc Palau

Affiliation: Instituto de Biomedicina de Valencia (CSIC), CIBER de Enfermedades Raras (CIBERER), Valencia, Spain

Yolanda Sanz

* E-mail: [email protected]

Affiliation: Instituto de Agroquímica y Tecnología de Alimentos, Consejo Superior de Investigaciones Científicas (IATA-CSIC), Valencia, Spain

Introduction

Celiac disease (CD) is a chronic inflammatory disorder of the small intestine that presents in genetically predisposed individuals following gluten consumption [1]. This disease often manifests in early childhood with small intestinal villous atrophy and signs of malabsorption [1]. Currently, a strict gluten-free diet is the only available treatment for patients but compliance with this dietary practice is extremely complex and, therefore, preventive strategies are being investigated [2].

The major genetic risk factor in CD is represented by Human Leukocyte Antigen (HLA)-DQ genes. Several studies have documented that the HLA-DQA1*05 and DQB1*02 alleles, encoding for particular DQ2 molecules, confer high susceptibility to CD [3]. This heterodimer can be encoded both in cis or trans forms. The susceptibility to CD is increased in homozygous subjects with a cis haplotype or possessing a second HLA-DQB1*02 allele [4]. In Europe, approximately 90% of patients have these genetic markers, whereas most of the remaining cases carry the HLA-DQA1*03 and DQB1*0302 alleles coding for DQ8 molecules [5]. Gluten is the main environmental factor responsible for the signs and symptoms of the disease but other environmental elements are also thought to play a role in the disease risk, including the type of milk-feeding, incidence of infections and intestinal dysbiosis [6]–[8].

Gut colonization starts immediately after birth and depends on multiple factors such as the type of delivery, contamination from the environment, the type of milk-feeding and, possibly, the genotype [9]–[11]. An adequate gut microbial colonization process contributes to the physiological development of the gut and the maturation of the immune system, thereby determining the risk of developing disease later in life. The early stage of colonization is characterized by the presence of higher levels of facultative anaerobes (enterobacteria, enterococci and streptococci) than of strictly anaerobic bacteria (e.g. bifidobacteria, bacteroides, clostridia, etc.); however, these proportions are reversed within a week following birth. In breast-fed infants Bifidobacterium spp. predominate, representing up to 90% of the total faecal microbiota, whereas in formula-fed infants the microbiota is more heterogeneous [9], [12], [13]. Breast-milk has been shown to be a continuous source of commensal bacteria to the infant gut, including species of the genera Lactobacillus and Bifidobacterium [14], [15]. It also contains prebiotic substances which are considered the main factors that stimulate the growth of Bifidobacterium spp. [16]. Epidemiological studies suggest that breast-feeding confers a protective effect against the risk of CD development, particularly when gluten is introduced in the diet while the infant is still breast-fed [1], [17]. However, the mechanisms underlying the beneficial effects of breast-milk on CD risk and their relationship with the gut microbiota are unknown. A preliminary study was previously conducted to test whether the HLA-DQ genotype could influence the composition of the gut microbiota, although the population under study was composed by a small number of exclusively breast-fed infants [18] and the aim was to establish the basis for this long-term study, including a representative cohort of infants.

The objective of this study was to assess the gut microbial colonization process during the first 4 months of life in breast-fed and formula-fed babies at risk of developing CD, by using quantitative PCR (qPCR). The ultimate purpose of our research is to gain a better understanding of the effects of early events leading to the acquisition of intestinal microbiota and their interactions with predisposing genes on CD risk.

Methods

Subjects and study design

A prospective observational study was carried out with a cohort of 164 healthy full-term newborns recruited between June 2006 and November 2010, who had a first-degree relative affected by CD. Data of mode of delivery, size, weight, weeks of gestation and type of feeding were recorded at birth and over the study period (Table 1). Infants were grouped for factors that influence the gut microbial colonization process, including age, genetic risk of CD development (HLA-DQ status), and type of milk-feeding (admitted randomly). The infants classified into the breast-fed group were those that received exclusively breast-feeding during the first 7 days of life, during the first month of life or during the first 4 months of life. Infants classified into the formula-fed group were those that received either exclusively formula or both formula and breast-milk at each sampling time. Infants were also grouped according the duration of breast-feeding (never breast-fed, breast-fed less than 1 month, breast-fed more than 1 month but less than 4 months and breast-fed for the 4 months). The study was approved by the ethics committees of Consejo Superior de Investigaciones Científicas (CSIC) and the Hospitals involved in the study, including Hospital Universitario Sant Joan de Reus, Hospital Universitario Sant Joan de Deu, Institut Dexeus, Hospital Clínico Universitario de Valladolid, Hospital Universitario La Paz, Hospital Infantil Universitario La Fe, Hospital Universitario Infantil Niño Jesús, Hospital Universitario Ntra Señora de Candelaria and Hospital Universitario Materno-Infantil de Canarias and conducted in accordance with the Helsinki Declaration of 1975 as revised in 1983. Written informed consent was obtained from the parents of infants included in the study.

[Figure omitted. See PDF.]

Table 1. Demographic data of infants under study.

https://doi.org/10.1371/journal.pone.0030791.t001

HLA-DQ Genotyping

DNA was extracted from yugal mucosa cells by scraping the inner side of the infants' cheek with sterile swabs (Copan innovation, Sarstedt, Germany) and purified according to the DNA IQ™ Casework Sample Kit for Maxwell® 16 protocol (Promega Biotech Iberica, Spain). Low-resolution HLA-DQB1 typing was performed by PCR-SSP (Polymerase Chain Reaction-Sequence Specific Primers) analysis [19]. Each PCR reaction was performed on about 20 ng of extracted DNA, 0.5 U of DNA polymerase (BIOTOOLS B&M S.A, Spain), 1× PCR Master Mix (Dynal AllSet⁺™ SSP or Olerup SSP™) containing nucleotides (200 µmol L⁻¹ each), PCR buffer, 5% glycerol and 100 µg mL⁻¹ cresol red, 0·25 µmol L⁻¹ of each allele- or group-specific primer pair and 0·1 µmol L⁻¹ of internal positive control primer pair matching a segment of the human growth hormone gene in a final volume of 10 µL. An initial denaturation step at 94°C for 2 min was followed by 10 two-temperature cycles (94°C for 10 s and 65°C for 60 s) and 20 three-temperature cycles (94°C for 10 s, 61°C for 50 s and 72°C for 30 s). Detection of amplified alleles was carried out on 2% agarose gels after ethidium bromide staining. HLA-DQA1 alleles were genotyped in a stepwise fashion for high resolution typing to hone the risk classification of each individual.

Faecal sampling and DNA extractions

Stool samples were collected from every subject at 7 days, 1 month and 4 months of age and frozen at −20°C immediately. Samples (1 g) were diluted 1∶10 (w/v) in PBS (pH 7.2) and homogenized by thorough agitation in a vortex. Aliquots were used for DNA extraction using the QIAamp DNA stool Mini kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. DNA extractions from different pure cultures of reference strains were done following the same protocol.

Quantitative PCR (qPCR) analysis of faecal bacteria

qPCR was used to quantify the different bacterial groups of the faecal microbiota using genus-, group- and species-specific primers as previously described [20], [21]. Briefly, PCR amplification and detection were performed with an ABI PRISM 7000-PCR sequence detection system (Applied Biosystems, UK) using SYBR® Green PCR Master Mix (SuperArray Bioscience Corporation, USA). The bacterial concentration from each sample was calculated by comparing the Ct values obtained from standard curves of reference strains. Standard curves were created using serial 10-fold dilutions of pure culture DNA corresponding to 10² to 10⁹ cells, as determined by microscopy counts after staining with 4′, 6 diamino-2-phenylindole in an epifluorescence microscope (Olympus BX51, Tokio, Japan).

Statistical analyses

Data were analysed using the SPSS 19.0 software for Windows (SPSS Inc, Chicago, IL, USA). The demographic characteristics of the study subjects are given as mean values (standard deviations [SDs]) for continuous variables and as numbers and proportions for categorical variables. Differences in demographic characteristic measures between the study groups were assessed using the Chi-Square Pearson test for categorical variables and ANOVA and post hoc LSD test for continuous variables. Microbiological data were transformed from exponential numbers into logarithms to adjust to normal and, in the tables, are expressed as mean values of log cells/g faeces and standard deviations (SDs). A mixed model with repeated measures with three fixed factors [genetic risk, type of milk-feeding and age (repeated measure)] was applied to determine the effects of genetic risk and type of feeding on bacterial counts. Interactions (magnitudinal) between these three factors were also analysed and not detected. Interactions between faecal microbial counts and type of delivery were not found and, therefore, data from infants with different type of delivery were grouped for statistical analyses. Differences between bacterial numbers at each sampling time (age) were analysed by ANOVA and post hoc LSD test. Correlations between factors (bacterial counts, age and genetic risk) were determined by Pearson correlation coefficients and correlations between bacterial counts and type of milk-feeding were tested by applying the Chi-Square Pearson test. In all cases, p-values less than 0.050 were considered statistically significant.

Results

Subjects and genetic risk of CD

The demographic characteristics of the infants under study are shown in Table 1. All newborns were full-term and most were delivered naturally (115 of 164). The size and weight of the infants at the moment of delivery were within standard parameters.

Infants were classified into three main groups according to their HLA-DQ genotype, and probabilities of developing CD were estimated according to previous studies [22], [23]. The first group included those individuals carrying the DQ2 haplotype in both cis (DQA1*0501-DQB1*0201 in homozygosis) and trans conformations (DQA1*0201-DQB1*0202 with DQA1*0505-DQB1*0301 in heterozygosis). The second group included those subjects carrying the DQ2 haplotype in cis conformation along with any other haplotype, as well as subjects carrying the DQ8 haplotype (DQA1*0301-DQB1*0302) in homozygosis. The third group included those individuals with other common genotypes not associated with CD. Of the 164 infants under study, 28.40% were classified in the first risk group, with the highest probability (>20%) of developing CD (the high genetic risk group) and 42.59% in the second group, with a >7% probability of developing CD (the intermediate genetic risk group). The remaining 29.01% of infants comprised the third group, with the lowest probability (<1%) of developing CD (the low genetic risk group).

Influence of milk-feeding type on faecal microbiota of infants at risk of developing CD

The effect of milk-feeding type on the composition of the faecal microbiota, irrespectively of genotype, was evaluated over the study period. The type of milk feeding significantly influence different bacterial group counts by analysing data with the linear mixed model with time sampling as the repeated measure. C. leptum group (P = 0.005), B. longum (P = 0.050), and B. breve (P = 0.008) numbers were significantly higher in breast-fed than in formula-fed infants, whereas Bacteroides fragilis group (P = 0.004), C. coccoides-E. rectale group (P<0.001), E. coli (P = 0.026), and B. lactis (P = 0.002) numbers were higher in formula-fed infants than in breast-fed infants.

Correlations between milk-feeding type and several bacterial group counts at each sampling time were also analysed. Increased numbers of C. leptum group and B. breve correlated with breast-feeding at 7 days of age (r = −0.255, P = 0.012; r = −0.197, P = 0.036, respectively). Moreover, formula-feeding correlated with increased numbers of E. coli, C. coccoides-E. rectale group and B. lactis at 1 month of age (r = 0.218, P = 0.006; r = 0.217, P = 0.011; r = 0.574, P = 0.002, respectively), and with increased numbers of C. coccoides-E. rectale and Bacteroides fragilis groups at 4 months of age (r = 0.280, P = 0.001; r = 0.267, P = 0.004, respectively).

When analysing the cumulative effect of breast-feeding on the microbiota composition of 4-month-old infants, statistically significant negative correlations were established between increased Bacteroides fragilis and C. coccoides-E. rectale group counts and longer breast-feeding duration (r = −0.218, P = 0.020; r = −0.245, P = 0.003, respectively).

Influence of HLA-DQ genotype on the fecal microbiota of infants at risk of developing CD

The influence of HLA-DQ genotype, irrespective of milk-feeding type, on fecal microbiota composition was established by applying a linear mixed-model analysis with sampling time as the repeated measure. According to this analysis, the effect of the genetic risk on bacterial numbers was significant for Bifidobacterium spp. (P<0.001) and B. longum (P<0.001), whose numbers increased when genetic risk of CD decreased. The effect of genetic risk was also significant for Staphylococcus spp. (P = 0.010) and Bacteroides fragilis group (P = 0.050), whose counts were higher when infants' genetic risk was greater. The influence of HLA-DQ genotype on fecal microbiota composition at each sampling time (7 days, 1 and 4 months) is also shown in Tables 2, 3 and 4.

[Figure omitted. See PDF.]

Table 2. Faecal microbiota of infants with different HLA-DQ genotype at 7 days of age analysed by qPCR.

https://doi.org/10.1371/journal.pone.0030791.t002

[Figure omitted. See PDF.]

Table 3. Faecal microbiota of infants with different HLA-DQ genotype at 1 month of age analysed by qPCR.

https://doi.org/10.1371/journal.pone.0030791.t003

[Figure omitted. See PDF.]

Table 4. Faecal microbiota of infants with different HLA-DQ genotype at 4 months of age analysed by qPCR.

https://doi.org/10.1371/journal.pone.0030791.t004

Correlations between genetic risk and several bacterial group counts at different infant ages (sampling time) were also established. Statistically significant correlations were found between increased numbers of both Bifidobacterium spp. (r = 0.235, P = 0.007 at 7 days; r = 0.268, P = 0.001 at 1 month; r = 0.257, P = 0.001 at 4 months) and B. longum (r = 0.209, P = 0.013 at 7 days; r = 0.207, P = 0.011 at 1 month; r = 0.200, P = 0.016 at 4 months) and reduced genetic risk during the whole sampling period. In contrast, increased numbers of Staphylococcus spp. at 1 and 4 months (r = −0.465, P = 0.001 at 1 month; r = −0.278, P = 0.038 at 4 months) and B. lactis at 4 months (r = −0.332, P = 0.026) correlated with increased genetic risk.

Influence of the HLA-DQ genotype in the faecal microbiota of either breast-fed or formula fed infants

The microbiota composition of infants grouped according to milk-feeding type was also analysed as a function of HLA-DQ genotype over the study period to eliminate the effects of milk-feeding type. In breast-fed infants, the effect of genetic risk of CD was significant on bacterial counts of Bifidobacterium spp. (P = 0.046), which decreased when the genetic risk was increased, and on bacterial counts of Staphylococcus spp. (P = 0.030), C. leptum group (P = 0.047), B. adolescentis ((P = 0.028) and B. dentium (P = 0.009), which increased when the genetic risk was also increased, according to linear mixed-model analysis with sampling time (age) as the repeated measure. The differences in mean bacterial numbers according to the genetic risk of developing CD in breast-fed infants analysed at each sampling point are also shown in Tables 5, 6 and 7.

[Figure omitted. See PDF.]

Table 5. Faecal microbiota of breast-fed infants with different HLA-DQ genotype at 7 days of age analysed by qPCR.

https://doi.org/10.1371/journal.pone.0030791.t005

[Figure omitted. See PDF.]

Table 6. Faecal microbiota of breast-fed infants with different HLA-DQ genotype at 1 month of age analysed by qPCR.

https://doi.org/10.1371/journal.pone.0030791.t006

[Figure omitted. See PDF.]

Table 7. Faecal microbiota of breast-fed infants with different HLA-DQ genotype at 4 months of age analysed by qPCR.

https://doi.org/10.1371/journal.pone.0030791.t007

Correlations between genetic risk and bacterial counts were also analysed and these were significant between increased genetic risk of CD development and increased C. leptum group counts in breast-fed infants (r = −0.408, P = 0.004).

In formula-fed infants, the effect of genetic risk on Bifidobacterium spp. and B. longum counts was found significant (P<0.001 and P<0.001, respectively) by applying a linear mixed model analysis with sampling time as the repeated measure, and those infants with reduced genetic risk had increased numbers of these bacterial groups. The effect of genetic risk on numbers of B. fragilis group and Staphylococcus spp. was also significant, following the opposite trend (P = 0.008 and P = 0.004, respectively). The differences in mean bacterial numbers according to the genetic risk of developing CD in formula-fed infants analysed at each sampling point are also shown in Tables 8, 9 and 10.

[Figure omitted. See PDF.]

Table 8. Faecal microbiota of formula-fed infants with different HLA-DQ genotype at 7 days of age analysed by qPCR.

https://doi.org/10.1371/journal.pone.0030791.t008

[Figure omitted. See PDF.]

Table 9. Faecal microbiota of formula-fed infants with different HLA-DQ genotype at 1 month of age analysed by qPCR.

https://doi.org/10.1371/journal.pone.0030791.t009

[Figure omitted. See PDF.]

Table 10. Faecal microbiota of formula-fed infants with different HLA-DQ genotype at 4 months of age analysed by qPCR.

https://doi.org/10.1371/journal.pone.0030791.t010

Statistically significant correlations were also established between the genetic risk of CD development and bacterial group numbers in formula-fed infants at different sampling times. At 7 days, reduced numbers of Bifidobacterium spp. correlated with increased genetic risk of developing CD (r = 0.362, P = 0.028). At 1 month, reduced numbers of Bifidobacterium spp. and B. longum correlated with increased genetic risk of developing CD (r = 0.511, P<0.001 and r = 0.454, P = 0.001, respectively), and increased numbers of Staphylococcus spp. correlated with increased genetic risk of developing CD (r = −0.573, P = 0.013). At 4 months, reduced numbers of Bifidobacterium spp. and B. longum correlated with increased genetic risk of developing CD (r = 0.412, P<0.001 and r = 0.336, P = 0.002, respectively), and increased numbers of Staphylococcus spp. and B. lactis correlated with increased genetic risk (r = −0.345, P = 0.040 and r = −0.442, P = 0.010, respectively).

The cumulative effect of breast-feeding on microbiota composition was also evaluated at 4 months of age and some significant correlations were established. When infants had never been breast-fed, significant correlations between increased genetic risk and low numbers of Bifidobacterium spp., and B. longum were detected (r = 0.343, P = 0.026; r = 0.455, P = 0.003, respectively). Similarly, when infants were breast-fed for more than 1 month and less than 4 months, low numbers of Bifidobacterium spp. correlated with high genetic CD risk (r = 0.408, P = 0.031), whereas high B. lactis numbers correlated with high genetic CD risk (r = −0.697, P = 0.006). On the other hand, when breast-feeding was exclusive until the fourth month of age, the high genetic risk correlated with increased numbers of C. leptum group (r = −0.454, P = 0.001).

Discussion

This is the first report on the effects of both milk feeding type and HLA-DQ genotype on the gut microbial colonization process of healthy full-term infants with a family history of CD. In the present study, milk-feeding type influenced the composition of the microbiota, partially in agreement with previous studies [24]. Breast-feeding favoured C. leptum group, B. longum and B. breve gut colonization, while formula-feeding favoured that of B. lactis, E. coli, C. coccoides-E. rectale group and B. fragilis group. Similar trends were detected when considering the cumulative effect of breast-feeding particularly for C. coccoides-E. rectale and B. fragilis groups.

In addition, specific features of fecal microbiota were associated with the genetic risk of developing CD, based on HLA-DQ genotype, when considered irrespectively of milk-feeding type. Increased numbers of Bifidobacterium spp. and B. longum were characteristic of microbiota of infants with the lowest genetic risk, whereas increased numbers of Staphylococcus spp. and B. fragilis group were characteristic of that of infants with the highest genetic CD risk. To date, there is limited evidence of a correlation between genotype and intestinal microbiota composition in humans. In previous human studies, monozygotic twins were demonstrated to have more similar faecal bacterial DNA profiles than unrelated individuals [25] and monozygotic twins more so than dizygotic twins [10]. More recently a strong bond between genotype, phenotype and changes in gut microbiota has been reported in patients with Cohn's disease, demonstrating that specific gene alterations (e.g. NOD2 and ATG16L1) can have an impact on intestinal microbiota composition [26]. Animal and human studies indicate that the host genotype may influence factors such as the repertoire of mucins, which act as bacterial adhesion sites in the intestinal mucosa, as well as the immune responses. Together, these can contribute to modulating the colonization of certain microorganisms [27]. In this context, our results suggest that HLA-DQ genotype influences the microbial colonization pattern early in life and, therefore, could be an additional factor influencing the risk of developing CD later in life. Enterocytes, which are in close proximity with intestinal content and bacteria, can express HLA class II molecules of the MHC to a certain extent and are able to act as antigen presenting cells [28]. Moreover, HLA class II molecules are primarily expressed by dendritic cells present in the lamina propria that can sample the mucosal surface for microbial antigens, which can be presented to naïve B and T cells after processed to peptides that are loaded on MHC class I and class II molecules. This is a critical step in triggering the mucosal innate immune response, which could restrict bacterial colonization and influence disease risk [29]. In particular, the increased Staphylococcus counts detected in the high CD genetic risk group of infants could favour the activation of a robust T-cell response by the preferential interaction of certain superantigens, such as staphylococcal superantigen A, with HLA-DQ molecules, thereby enhancing the risk of T-cell mediated diseases [30], such as CD.

The influence of the HLA-DQ genotype was also analysed in subgroups of either breast-fed or formula fed infants to gather more information about the respective effects of each variable (genotype and milk-feeding type) on the microbiota of the infants under study. In the whole infant population, formula feeding favoured the presence of increased numbers of B. fragilis group, which were also higher in infants with higher genetic risk of developing CD in the whole population and in the subgroup of formula-fed infants, suggesting that the colonization of this bacterial group is greatly defined by type of milk-feeding. Nevertheless, increased counts of Staphylococcus spp. were associated with increased genetic risk of developing CD in the whole population and in both breast- and formula-fed infant subgroups, but their colonization was not favoured by formula feeding, suggesting that the HLA-DQ genotype plays a more prominent role in the colonization of this bacterial group. Notably, reduced numbers of Bifidobacterium spp. were associated with an increased risk of developing CD in the whole population and in both breast- and formula-fed infants, and the colonization of species of this genus was also favoured by breast-feeding. Therefore, the findings suggest that bifidobacterial numbers can be influenced by both the HLA-DQ genotype and the milk-feeding type. Moreover, low counts of the species B. longum were found in infants of higher risk to develop CD in the whole population and in the subgroup of formula-fed infants but not in the subgroup of breast-fed infants indicating that the breast-feeding is providing to the infant's gut microbiota certain Bifidobacterium spp. [31], [16], which could partially explain the protective role attributed to breast feeding in the risk of developing CD in previous epidemiological studies [17].

These and other commensal bacteria are recognized as constituting major stimuli for the adequate development of immune functions and oral tolerance [6], which could also be related to the risk of developing CD. In another prospective study, a reduced ratio of Bifidobacterium to Clostridium counts in the faecal microbiota of infants was shown to precede the development of atopic diseases later in life, indicating that the relative proportions of these bacterial groups may favour or protect against the development of immune-related disorders [32]. Moreover, Bifidobacterium spp. and B. longum levels in both biopsies and faeces have been reported lower in CD patients than in healthy controls [31], [33]. However, breast feeding also favoured the presence of increased numbers of C. leptum group when this factor was considered alone, and this bacterial group was associated with an increased risk of developing CD in the subgroup of breast-fed infants. C. leptum group numbers were also reported to be higher in the fecal and duodenal microbiota of CD patients than in healthy controls [34].

Overall, this study demonstrates that the milk-feeding type and the HLA-DQ genotype differently influence the bacterial colonization pattern of the newborn intestine during the first 4 months of life and, therefore, could also influence the risk of developing CD in later life. Breast-feeding reduced the genotype-related differences in microbiota composition, which could partly explain the protective role attributed to breast milk in this disorder. Further studies are underway to reveal additional evidence of the role played by early intestinal colonization patterns in CD development in this cohort of infants.

Acknowledgments

We thank Dr. Laura Barrios and Fernando Santoveña for statistical advice.

Author Contributions

Conceived and designed the experiments: YS. Performed the experiments: GDP AC. Analyzed the data: GDP AC YS. Contributed reagents/materials/analysis tools: FP EN GC VV TP AJG IP AL CR-K AM ADGN CC LO LPQ. Wrote the paper: GDP AC YS.

Citation: De Palma G, Capilla A, Nova E, Castillejo G, Varea V, Pozo T, et al. (2012) Influence of Milk-Feeding Type and Genetic Risk of Developing Coeliac Disease on Intestinal Microbiota of Infants: The PROFICEL Study. PLoS ONE 7(2): e30791. https://doi.org/10.1371/journal.pone.0030791

References

1. Fasano A, Catassi C (2005) Coeliac disease in children. Best Pract Res Clin Gastroenterol 19: 467–478.A. FasanoC. Catassi2005Coeliac disease in children.Best Pract Res Clin Gastroenterol19467478

2. Sanz Y (2009) Novel Perspectives in Celiac Disease Therapy. Mini-Reviews in Medicinal Chemistry 9: 359–367.Y. Sanz2009Novel Perspectives in Celiac Disease Therapy.Mini-Reviews in Medicinal Chemistry9359367

3. Sollid LM (2002) Coeliac disease: dissecting a complex inflammatory disorder. Nat Rev Immunol 2: 647–655.LM Sollid2002Coeliac disease: dissecting a complex inflammatory disorder.Nat Rev Immunol2647655

4. Dubois PC, van Heel DA (2008) Translational mini-review series on the immunogenetics of gut disease: immunogenetics of coeliac disease. Clin Exp Immunol 153: 162–173.PC DuboisDA van Heel2008Translational mini-review series on the immunogenetics of gut disease: immunogenetics of coeliac disease.Clin Exp Immunol153162173

5. Kagnoff MF (2007) Celiac disease: pathogenesis of a model immunogenetic disease. J Clin Invest 117: 41–49.MF Kagnoff2007Celiac disease: pathogenesis of a model immunogenetic disease.J Clin Invest1174149

6. Sanz Y, Sánchez E, De Palma G, Medina M, Marcos A, et al. (2008) Indigenous gut microbiota, probiotics, and coeliac disease. In: Overton LT, Ewente MR, editors. Child Nutrition & Physiology. Nova Science Publishers Inc, NY. pp. 210–224.Y. SanzE. SánchezG. De PalmaM. MedinaA. Marcos2008Indigenous gut microbiota, probiotics, and coeliac disease.LT OvertonMR EwenteChild Nutrition & PhysiologyNova Science Publishers Inc, NY210224

7. Plot L, Amital H (2009) Infectious associations of Celiac disease. Autoimmunity Reviews 8: 316–319.L. PlotH. Amital2009Infectious associations of Celiac disease.Autoimmunity Reviews8316319

8. Plot L, Amital H, Barzilai O, Ram M, Nicola B, et al. (2009) Infections May Have a Protective Role in the Etiopathogenesis of Celiac Disease.670–674.L. PlotH. AmitalO. BarzilaiM. RamB. Nicola2009Infections May Have a Protective Role in the Etiopathogenesis of Celiac Disease.670674

9. Salminen S, Isolauri E (2006) Intestinal colonization, microbiota, and probiotics. Journal of Pediatrics 149: S115–S120.S. SalminenE. Isolauri2006Intestinal colonization, microbiota, and probiotics.Journal of Pediatrics149S115S120

10. Stewart JA, Chadwick VS, Murray A (2005) Investigations into the influence of host genetics on the predominant eubacteria in the faecal microflora of children. J Med Microbiol 54: 1239–1242.JA StewartVS ChadwickA. Murray2005Investigations into the influence of host genetics on the predominant eubacteria in the faecal microflora of children.J Med Microbiol5412391242

11. Bezirtzoglou E, Romond C (1990) Effect of the feeding practices on the establishment of bacterial interactions in the intestine of the newborn delivered by cesarean section. J Perinat Med 18: 139–143.E. BezirtzoglouC. Romond1990Effect of the feeding practices on the establishment of bacterial interactions in the intestine of the newborn delivered by cesarean section.J Perinat Med18139143

12. Bezirtzoglou E (1997) The intestinal microflora during the first weeks of life. Anaerobe 3: 173–177.E. Bezirtzoglou1997The intestinal microflora during the first weeks of life.Anaerobe3173177S1075-9964(97)90102-5 [pii];10.1006/anae.1997.0102 [doi]. S1075-9964(97)90102-5 [pii];10.1006/anae.1997.0102 [doi].

13. Bezirtzoglou E, Tsiotsias A, Welling GW (2011) Microbiota profile in feces of breast- and formula-fed newborns by using fluorescence in situ hybridization (FISH). Anaerobe. E. BezirtzoglouA. TsiotsiasGW Welling2011Microbiota profile in feces of breast- and formula-fed newborns by using fluorescence in situ hybridization (FISH).AnaerobeS1075-9964(11)00033-3 [pii];10.1016/j.anaerobe.2011.03.009 [doi]. S1075-9964(11)00033-3 [pii];10.1016/j.anaerobe.2011.03.009 [doi].

14. Martin R, Jimenez E, Heilig H, Fernandez L, Marin ML, et al. (2009) Isolation of bifidobacteria from breast milk and assessment of the bifidobacterial population by PCR-denaturing gradient gel electrophoresis and quantitative real-time PCR. Appl Environ Microbiol 75: 965–969.R. MartinE. JimenezH. HeiligL. FernandezML Marin2009Isolation of bifidobacteria from breast milk and assessment of the bifidobacterial population by PCR-denaturing gradient gel electrophoresis and quantitative real-time PCR.Appl Environ Microbiol75965969

15. Reviriego C, Eaton T, Martin R, Jimenez E, Fernandez L, et al. (2005) Screening of virulence determinants in Enterococcus faecium strains isolated from breast milk. J Hum Lact 21: 131–137.C. ReviriegoT. EatonR. MartinE. JimenezL. Fernandez2005Screening of virulence determinants in Enterococcus faecium strains isolated from breast milk.J Hum Lact21131137

16. Martin R, Langa S, Reviriego C, Jiminez E, Marin ML, et al. (2003) Human milk is a source of lactic acid bacteria for the infant gut. J Pediatr 143: 754–758.R. MartinS. LangaC. ReviriegoE. JiminezML Marin2003Human milk is a source of lactic acid bacteria for the infant gut.J Pediatr143754758

17. Ivarsson A, Hernell O, Stenlund H, Persson LA (2002) Breast-feeding protects against celiac disease. Am J Clin Nutr 75: 914–921.A. IvarssonO. HernellH. StenlundLA Persson2002Breast-feeding protects against celiac disease.Am J Clin Nutr75914921

18. De Palma G, Capilla A, Nadal I, Nova E, Pozo T, et al. (2010) Interplay Between Human Leukocyte Antigen Genes and the Microbial Colonization Process of the Newborn Intestine. Current Issues in Molecular Biology 12: 1–10.G. De PalmaA. CapillaI. NadalE. NovaT. Pozo2010Interplay Between Human Leukocyte Antigen Genes and the Microbial Colonization Process of the Newborn Intestine.Current Issues in Molecular Biology12110

19. Olerup O, Aldener A, Fogdell A (1993) HLA-DQB1 and -DQA1 typing by PCR amplification with sequence-specific primers (PCR-SSP) in 2 hours. Tissue Antigens 41: 119–134.O. OlerupA. AldenerA. Fogdell1993HLA-DQB1 and -DQA1 typing by PCR amplification with sequence-specific primers (PCR-SSP) in 2 hours.Tissue Antigens41119134

20. Matsuki T, Watanabe K, Fujimoto J, Miyamoto Y, Takada T, et al. (2002) Development of 16S rRNA-gene-targeted group-specific primers for the detection and identification of predominant bacteria in human feces. Appl Environ Microbiol 68: 5445–5451.T. MatsukiK. WatanabeJ. FujimotoY. MiyamotoT. Takada2002Development of 16S rRNA-gene-targeted group-specific primers for the detection and identification of predominant bacteria in human feces.Appl Environ Microbiol6854455451

21. Malinen E, Kassinen A, Rinttila T, Palva A (2003) Comparison of real-time PCR with SYBR Green I or 5′-nuclease assays and dot-blot hybridization with rDNA-targeted oligonucleotide probes in quantification of selected faecal bacteria. Microbiology 149: 269–277.E. MalinenA. KassinenT. RinttilaA. Palva2003Comparison of real-time PCR with SYBR Green I or 5′-nuclease assays and dot-blot hybridization with rDNA-targeted oligonucleotide probes in quantification of selected faecal bacteria.Microbiology149269277

22. Donat E, Planelles D, Capilla-Villanueva A, Montoro JA, Palau F, et al. (2009) Allelic distribution and the effect of haplotype combination for HLA type II loci in the celiac disease population of the Valencian community (Spain). Tissue Antigens 73: 255–261.E. DonatD. PlanellesA. Capilla-VillanuevaJA MontoroF. Palau2009Allelic distribution and the effect of haplotype combination for HLA type II loci in the celiac disease population of the Valencian community (Spain).Tissue Antigens73255261

23. Bourgey M, Calcagno G, Tinto N, Gennarelli D, Margaritte-Jeannin P, et al. (2007) HLA related genetic risk for coeliac disease. Gut 56: 1054–1059.M. BourgeyG. CalcagnoN. TintoD. GennarelliP. Margaritte-Jeannin2007HLA related genetic risk for coeliac disease.Gut5610541059

24. Penders J, Thijs C, Vink C, Stelma FF, Snijders B, et al. (2006) Factors influencing the composition of the intestinal microbiota in early infancy. Pediatrics 118: 511–521.J. PendersC. ThijsC. VinkFF StelmaB. Snijders2006Factors influencing the composition of the intestinal microbiota in early infancy.Pediatrics118511521

25. Zoetendal EG (2001) The Host Genotype Affects the Bacterial Community in the Human Gastronintestinal Tract. Microb Ecol Health Dis 13: 129–134.EG Zoetendal2001The Host Genotype Affects the Bacterial Community in the Human Gastronintestinal Tract.Microb Ecol Health Dis13129134

26. Frank DN, Robertson CE, Hamm CM, Kpadeh Z, Zhang T, et al. (2011) Disease phenotype and genotype are associated with shifts in intestinal-associated microbiota in inflammatory bowel diseases. Inflamm Bowel Dis 17: 179–184.DN FrankCE RobertsonCM HammZ. KpadehT. Zhang2011Disease phenotype and genotype are associated with shifts in intestinal-associated microbiota in inflammatory bowel diseases.Inflamm Bowel Dis17179184

27. Sanz Y, Nadal I, Sanchez E (2007) Probiotics as drugs against human gastrointestinal infections. Recent Pat Antiinfect Drug Discov 2: 148–156.Y. SanzI. NadalE. Sanchez2007Probiotics as drugs against human gastrointestinal infections.Recent Pat Antiinfect Drug Discov2148156

28. Kaiserlian D, Vidal K (1993) Antigen presentation by intestinal epithelial cells. Immunol Today 14: 144.D. KaiserlianK. Vidal1993Antigen presentation by intestinal epithelial cells.Immunol Today14144

29. Hapfelmeier S, Muller AJ, Stecher B, Kaiser P, Barthel M, et al. (2008) Microbe sampling by mucosal dendritic cells is a discrete, MyD88-independent step in DeltainvG S. Typhimurium colitis. J Exp Med 205: 437–450.S. HapfelmeierAJ MullerB. StecherP. KaiserM. Barthel2008Microbe sampling by mucosal dendritic cells is a discrete, MyD88-independent step in DeltainvG S. Typhimurium colitis.J Exp Med205437450

30. Rajagopalan G, Polich G, Sen MM, Singh M, Epstein BE, et al. (2008) Evaluating the role of HLA-DQ polymorphisms on immune response to bacterial superantigens using transgenic mice. Tissue Antigens 71: 135–145.G. RajagopalanG. PolichMM SenM. SinghBE Epstein2008Evaluating the role of HLA-DQ polymorphisms on immune response to bacterial superantigens using transgenic mice.Tissue Antigens71135145

31. Newburg DS (2005) Innate immunity and human milk. J Nutr 135: 1308–1312.DS Newburg2005Innate immunity and human milk.J Nutr13513081312

32. Kalliomaki M, Kirjavainen P, Eerola E, Kero P, Salminen S, et al. (2001) Distinct patterns of neonatal gut microflora in infants in whom atopy was and was not developing. J Allergy Clin Immunol 107: 129–134.M. KalliomakiP. KirjavainenE. EerolaP. KeroS. Salminen2001Distinct patterns of neonatal gut microflora in infants in whom atopy was and was not developing.J Allergy Clin Immunol107129134

33. Nadal I, Donat E, Ribes-Koninckx C, Calabuig M, Sanz Y (2007) Imbalance in the composition of the duodenal microbiota of children with coeliac disease. J Med Microbiol 56: 1669–1674.I. NadalE. DonatC. Ribes-KoninckxM. CalabuigY. Sanz2007Imbalance in the composition of the duodenal microbiota of children with coeliac disease.J Med Microbiol5616691674

34. Collado MC, Donat E, Ribes-Koninckx C, Calabuig M, Sanz Y (2009) Specific duodenal and faecal bacterial groups associated with paediatric coeliac disease. J Clin Pathol 62: 264–269.MC ColladoE. DonatC. Ribes-KoninckxM. CalabuigY. Sanz2009Specific duodenal and faecal bacterial groups associated with paediatric coeliac disease.J Clin Pathol62264269