Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2008, p. 4686-4694, Vol. 74, No. 15
0099-2240/08/$08.00+0 doi:10.1128/AEM.00122-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Eline S. Klaassens,1,
*
Erja Malinen,3
Willem M. de Vos,1,3 and
Elaine E. Vaughan1,4
Laboratory of Microbiology, Wageningen University, Dreijenplein 10, 6703 HB Wageningen, The Netherlands,1 Departamento de Sistemas Biológicos, Universidad Autónoma Metropolitana-Xochimilco, Apartado Postal 23-181, México 16000, México,2 Department of Basic Veterinary Sciences, Helsinki University, FI-00014 Helsinki, Finland,3 Unilever Food and Health Research Institute, Unilever Research and Development, Olivier van Noortlaan 120, 3133 AT Vlaardingen, The Netherlands4
Received 15 January 2008/ Accepted 15 May 2008
|
|
|---|
|
|
|---|
The molecular mechanisms that underlie the influence of breast-feeding on the development of the intestinal bifidobacterial community are not yet completely understood. The publication of the genome sequence from Bifidobacterium longum NCC2705 revealed that the chromosome encodes numerous genes for carbohydrate utilization collectively termed the glycobiome (38). Sequence data indicated a preference for di-, tri-, and oligosaccharides, pointing toward a bias for complex oligosaccharides (33, 34) complemented with transporters for a variety of disaccharides and oligosaccharides (33). This suggests adaptation to a special colonic niche and utilization of specific dietary components (38). Recently, a bifidobacterium-specific microarray platform comprising clone library-containing inserts of seven bifidobacterial species was constructed (9). In the present study the effects of human milk and formula milk on the global gene expression of B. longum were investigated to gain some insight into the impact of the milks on the development and colonization of intestinal bifidobacteria within the infant gut. Transcriptomics was performed on B. longum LMG 13197 during growth in human and formula milk with GOS-lcFOS and in comparison to semisynthetic medium with glucose. This was complemented with a study of the metabolic effect of semisynthetic medium with GOS on the B. longum transcriptome.
|
|
|---|
Growth in milk and semisynthetic medium with prebiotics.
Human milk was aseptically collected from two volunteers with a sterile pump and immediately frozen at –20°C. Written consent was received from all subjects who provided human breast milk. Commercial formula milk (Nutrilon 1 with GOS and polyfructose; Nutricia, Zoetermeer, The Netherlands) was dissolved to 10% in demineralized water. Both milks were centrifuged for 5 min at 5,000 rpm (4°C) and filtered to remove fat. Filtered-sterilized cysteine was added (0.5 g/liter), and the pH was adjusted to 6.5. Fifty milliliters of milk was added to sterile flasks with stoppers; flasks were sealed with aluminum rings and placed in a water bath for 20 min at 80°C. After flasks were cooled to room temperature, they were flushed with CO2/N2 using sterile filters. They were incubated a second time in a water bath for 20 min at 80°C and cooled to room temperature overnight. A semisynthetic medium was prepared with a minimal carbon source (6.7 g/liter yeast nitrogen base [Difco]), 0.5 g/liter cysteine, 1 ml of Tween 80, and a 40% [vol/vol] salt solution consisting of 0.2 g/liter CaCl2, 1 g/liter K2HPO4, 1 g/liter KH2PO4, 10 g/liter NaHCO3, 2 g/liter NaCl, 0.5 g/liter casein enzymatic hydrolysate, 0.05 g/liter sodium thioglycolate) such that it could be supplemented with the different carbon sources to be monitored (19). Transcriptomics of B. longum in human or formula milks were compared to the semisynthetic medium containing 5 g/liter glucose as a sole carbon source (hereafter, glucose medium). The number of CFU was determined by plating on MRS agar medium (Difco, France) with cysteine (0.5 g/liter), followed by incubation at 37°C anaerobically. The semisynthetic medium supplemented with GOS (composed of 33% disaccharides, 39% trisaccharides, 18% tetrasaccharides, 7% pentasaccharides, and 3% hexa-, hepta-, and octasaccharides; Friesland Foods, The Netherlands) (hereafter, GOS medium) as a sole carbon source was used to study the transcriptomics of B. longum in prebiotics. For comparative purposes B. longum was grown in semisynthetic medium with glucose or no carbon source. Carbon solutions were sterilized separately and added prior to fermentation at a final concentration of 5 g/liter. All media and carbon solutions were flushed with CO2/N2; flasks were closed with rubber stoppers, sealed with aluminum rings, and sterilized at 120°C for 15 min. Growth was followed by measuring the optical density at 600 nm (OD600). All tests were performed in duplicate at 37°C with agitation (200 rpm) for 10 h, and samples were taken every 2 h for pH measurements and determination of the number of CFU by plating on MRS medium with cysteine (0.5 g/liter), as described above.
RNA isolation and purification.
Cells for RNA isolation (50 ml) were harvested at mid-exponential phase (pH 5; OD600 of 0.5) and concentrated immediately by centrifugation at 12,000 rpm for 10 min at 4°C. After centrifugation, the supernatant was discarded, and the cell pellet was snap frozen in liquid nitrogen and stored at –80°C. Samples for RNA isolation were taken from milk and glucose medium cultures after 4 h and from GOS cultures after 6 h (Fig. 1 and 2). RNA isolation was performed using the Macaloid clay method as previously described (35). Briefly, cell pellets were resuspended in 0.5 ml of Tris-EDTA buffer (10 mM; pH 7.5) added to tubes containing 0.8 g of 0.1-mm zirconia beads (Biospec Products, Bartlesville, OK), 0.18 g of Macaloid clay in Tris-EDTA buffer (Kronos Titan GmbH, Leverkusen, Germany), 50 µl of 10% sodium dodecyl sulfate (SDS), and 0.5 ml of acidic phenol (pH 3.75). Cells were disrupted using three treatments of 45 s in a Fastprep (Qbiogene) interspersed by a 90-s incubation on ice. After disruption the samples were centrifuged at 4°C for 10 min at 12,000 rpm. The aqueous phase was used for RNA isolation by extraction with chloroform-phenol-isoamylalcohol (1:1:1) twice and with chloroform-isoamylalcohol (1:24) once. This was followed by a cleaning step using RNeasy columns (Qiagen Sciences, MD), according to the RNeasy cleaning protocol, and treatment with DNase (Roche Diagnostics, Mannheim, Germany) on the column. After elution from the column, RNA concentrations were measured using the Nanodrop ND-1000 (Nanodrop Technologies, Wilmington, DE). Additionally, quality and concentration of the RNA were measured using an Experion Automated Electrophoresis System (Bio-Rad) with RNA StdSens microfluidic chips. Total RNA samples with a 16S/23S rRNA ratio of at least 2.0 were used for microarray experiments.
![]() View larger version (11K): [in a new window] |
FIG. 1. Growth (log10 CFU/ml) (A) and pH curves (B) for B. longum in glucose medium ( ), human breast milk ( ), and formula milk ( ). An arrow indicates the harvest point for RNA isolation.
|
![]() View larger version (12K): [in a new window] |
FIG. 2. Growth (OD600) (A) and pH (B) of B. longum on glucose ( ), GOS ( ), and no carbon source (x). An arrow indicates the harvest point for RNA isolation. Error bars indicate standard deviations.
|
Fluorescent labeling and hybridization to bifidobacterial microarrays.
The microarray slides were provided by TNO Nutrition and Food Research, The Netherlands (9). The array contained 2,000 spots of a B. longum LMG 13197 clone library containing an average insert size of approximately 1.5 kb covering 85% of the genome. Slides were incubated in preheated prehybridization solution (1% bovine serum albumin, 5x SSC [1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate], 0.1% SDS; filtered by 0.45-µm-pore-size filters) at 42°C for 45 min under rotation, washed, and dried under N2 flow. cDNA was dissolved in 4.5 µl of 0.1 M sodium carbonate buffer, pH 9.0, for 10 min at room temperature. An aliquot of 4.5 µl of the appropriate N-hydroxysuccinimide ester Cy dye was added and incubated in the dark at room temperature for 1 h. After incubation, unincorporated dyes were removed using Autoseq G50 columns (AmershamBiosciences, Freiburg, Germany). The Cy5- and Cy3-labeled probes were combined and dried using vacuum centrifugation. The samples were dissolved in 45 µl of Easyhyb buffer at 42°C for 15 min. The preparations were mixed well and spun down briefly. The subsequent steps were performed at 42°C. The prehybridized microarray slides were placed on hybridization chambers (Corning) moistened with MilliQ water. Hybridization mixtures were pipetted onto the slides and covered with glass lids. The chambers were immediately sealed, and the hybridization was then allowed to take place overnight at 42°C. Afterwards, the slides were washed by placing them in 50-ml tubes containing the following wash solutions at 37°C and agitated for 10 s: 1x SSC-0.2% SDS, 0.5x SSC, 0.2x SSC; finally the slides were washed twice in 0.2x SSC with agitation on a rotation plateau for 10 min. The slides were dried under N2 flow and scanned using ScanArray (Perkin Elmer) using two wavelengths corresponding to the excitation wavelengths of the Cy3 and Cy5 dyes. Scanned images of the slides were further processed with ImaGene, version 5.6, software.
Analysis of microarray data.
The microarray experimental design for the various samples is presented in Fig. 3. Samples for each condition were hybridized in duplicate. Human milk samples were compared to the formula milk and to samples of semisynthetic medium with glucose, which were also compared to GOS medium samples. For data analysis, flagged spots were eliminated, followed by background subtraction, log transformation, and Lowess fit normalization, using the Arraynorm program (http://genome.tugraz.at). Positive signals in biological replicates were analyzed for significant differential expression, using significance analysis of microarray with the Tigr Multiarray viewer program (37). Raw microarray data are available in the supplemental material.
![]() View larger version (18K): [in a new window] |
FIG. 3. Hybridization scheme for the different growth media tested. The arrowheads correspond to the Cy5-labeled samples, and the dots at the ends of the arrows correspond to the Cy3-labeled samples.
|
|
|
|---|
![]() View larger version (12K): [in a new window] |
FIG. 4. Venn diagram showing the number of differential and commonly hybridized clones in human breast milk (BM) compared to formula milk (FM) and human milk compared to glucose medium (Gl).
|
|
View this table: [in a new window] |
TABLE 1. B. longum genes upregulated in human milk compared to formula milk and glucose medium
|
(i) Carbohydrate metabolism.
Despite the addition of oligosaccharides to the formula milk, there were only a few genes implicated in carbohydrate metabolism that were similarly upregulated during growth in both formula and human milk. These included two putative glycosyl hydrolases (BL0529 and BL0544) (Tables 1 and 2) and an ATP-binding cassette (ABC)-type transporter (BL1164-BL1165), which were reported to be induced by FOS and lactose in semisynthetic medium (33), and a solute-binding protein of ABC transporter for sugars (BL1638) (Table 1).
|
View this table: [in a new window] |
TABLE 2. B. longum genes highly upregulated in formula compared to human milk
|
(iii) Colonization factors.
Two genes for cell surface proteins with high similarity to the type 2 glycoprotein-binding FimA fimbrial subunit and sortase-like fimbria-associated protein of the oral inhabitant Actinomyces naeslundii, in clones BL0675 and BL0676, respectively (Tables 1 and 2), were highly upregulated during growth of B. longum in human milk and even more in formula milk (Table 2). These proteins are involved in the adhesion of A. naeslundii to microbial oral biofilms and interaction with epithelial surfaces (29, 32). The putative FimA fimbrial subunit clone BL0675 contained other motifs: a putative immunoglobulin domain found in other cell-surface proteins and a Cna protein B-type domain found in the Staphylococcus aureus collagen-binding surface protein, which forms a repetitive β-sandwich structure that presents the ligand-binding domain away from the bacterial cell surface. The putative sortase-like fimbria-associated surface protein BL0676 contained the sortase LPXTG motif, which can function as an anchor on the cell wall. There is some tentative evidence that bifidobacteria have fimbria-like structures (38); however, further characterization is required to determine if they contribute to attachment or retention in the human intestine.
Genes upregulated specifically during growth in human milk.
Complex neutral oligosaccharides have been identified as the most likely prebiotic factors in human milk that stimulate the growth of bifidobacteria in the infant gut. This was supported by the high percentage of upregulated genes with predicted functions in the transport and metabolism of oligosaccharides in human compared to formula milk (Table 1). Significantly upregulated genes included those for putative N-acetylglucosamine-6-P deacetylase (BL1344) that catalyzes the deacetylation of N-acetyl-D-glucosamine 6-P and glucosamine-6-phosphate isomerase 1 (BL1343), which converts D-glucosamine-6-P to D-fructose-6-P and ammonia. N-Acetylglucosamine-containing oligosaccharides are regarded as the most bifidogenic oligosaccharides in human milk (24). Consequently B. longum would be able to use the N-acetylglucosamine-6-P in the human milk as a source of carbon and nitrogen. Other significant genes involved in carbohydrate metabolism upregulated specifically in human milk included a putative desulfatase (BL1642) and sugar permease subunit (BL1639) (Table 1). The latter forms part of a novel putative operon for metabolism of mucin sugars and human milk oligosaccharides in B. longum JCM1217. The complete operon encodes an ABC-type sugar transporter, lactose-N-biose phosphorylase, mucin desulfatase, galactose-1-phosphate uridylyltransferase, and UDP-glucose 4-epimerase. It has been speculated that this operon plays a critical role in the bifidogenic effect of human milk as the lacto-N-biose structure appears to be unique to human milk oligosaccharides (20); moreover, intestinal colonization by bifidobacteria would be enhanced by metabolism of mucin sugars. The putative (galacto)oligosaccharide transporter genes (BL1638 and BL1639) were previously found to be induced in vitro by FOS and expressed but not induced by lactose, raffinose, and maltose (33). Another gene involved in galactose metabolism, a galactose-1-phosphate uridylyltransferase gene (BL1211), was also specifically upregulated in breast milk.
As expected, since lactose is the predominant sugar in milk, a putative galactoside transporter gene for lactose transport, lacS (in clone BL0976), was upregulated in human milk although rather surprisingly not in formula milk. In fact, the amount of lactose reaching the newborn colon from breast as well as formula milk is low since it is absorbed in the ileum (11). However, the expression of lacS in B. longum NCC2705 was shown to be induced by several carbohydrates besides lactose including raffinose, maltose, and FOS (33, 34), so perhaps alternative sugar sources such as oligosaccharides specific to human milk promoted the upregulation. Significant upregulation of a putative β-galactosidase gene (lacZ, or BL0978) was also detected during growth in human milk and again not in formula milk; lacZ may be involved in hydrolysis of lactose to glucose and galactose, and/or alternatively in hydrolysis of terminal nonreducing β-D-galactose residues in milk oligosaccharides (8). A tentative lactaldehyde reductase (BL1673) gene that anaerobically reduces L-lactaldehyde (a product of both the L-fucose and L-rhamnose catabolic pathways to L-1,2-propanediol) was highly upregulated (30) only in human milk. It is noteworthy that L-fucose can be derived from human milk oligosaccharides (24) or mucin sugars. Putative genes for transketolase and transaldolase, involved in the pentose phosphate pathway (BL0715, BL0716, and BL0966), were upregulated to a great extent in the human milk compared to glucose medium. The pentose phosphate pathway, besides energy metabolism, is used for reducing agents, nucleotide formation, and reducing equivalents for synthetic purposes (41), suggesting an enhanced requirement for the latter in human milk. The bifidobacterial transaldolase protein was previously identified in vivo in fresh infant feces as well as in batch cultures in MRS medium via a proteomics approach (21). This supports the assumption that the genes expressed in these in vitro studies in milk can also be expected in vivo.
Highly upregulated genes in formula milk versus human milk.
There were 54 genes significantly upregulated in formula milk compared to human milk (and glucose medium), and the sequences of 17 clones have been determined (Table 2). Among the genes tentatively identified were those annotated for carbohydrate utilization, including a fragment of lacZ and a glucosyltransferase (BL1775 and BL0672). Overall, the number of differentially expressed genes involved in carbohydrate metabolism in formula milk was much lower than the number upregulated in human milk. One could speculate that this is due to a lower carbohydrate metabolism by B. longum in formula milk with GOS-lcFOS than in human milk, as the complexity of the carbohydrates in breast milk is much higher and is likely to induce a more complex carbohydrate metabolism than the formula milk. However, an insufficient number of clones were sequenced and identified to confirm this.
A few of the genes were also upregulated during growth of B. longum in human milk but more strongly in formula milk (Table 2). Specifically, the putative genes for regulation of glutamine synthase activity (BL0433 and BL0434) were more highly upregulated in formula milk than in human milk. The same was observed for the potential fimbria genes discussed above. Other formula-specific upregulated genes were predicted to be involved in transport such as a possible oxalate/formate transporter (BL0613), lipid metabolism (BL1537), and several genes of unknown function. A nucleotide hydrolyzing gene (BL1148) was found to be upregulated, as well as a ribonucleotide-2P reductase gene (BL0671), possibly due to the influence of the high level of nucleotides in formula compared to human milk. A putative bioY gene (BL1534) encoding part of a group of BioMNY proteins, which are considered to constitute tripartite biotin ABC-transporters in prokaryotes, was significantly upregulated in formula milk (17). This might be due to the presence of a high concentration of biotin in the formula milk (17 µg/liter) compared to the human milk used in this study; the level in the latter can vary from 5 to 11 µg/liter (4). Enzymes needed for synthesis of biotin appeared to be missing from the B. longum genome sequence (38).
Growth and transcriptomics of B. longum in GOS medium.
The transcriptome of B. longum growing on GOS medium compared with glucose medium was studied at mid-exponential growth (Fig. 2).
B. longum showed a relatively lower specific growth rate on GOS (µ = 0.24 ± 0.008 h–1) than on glucose (µ = 0.27 ± 0.044 h–1) with a concomitant slower production of acids. No growth was observed when the bacteria were grown on the semisynthetic medium without a carbon source (Fig. 2). Comparison of the transcriptome of B. longum grown in GOS medium with glucose showed 91 clones with significantly induced signals (Table 3). Sequencing of 36 clones revealed 25 involved in carbohydrate metabolism and transport, 1 in amino acid biosynthesis, 1 in ribosomal genes, 1 in genetic information, 1 in cellular processes, and 17 with unknown functions. The Venn diagram (Fig. 5) indicates that many more genes are upregulated with growth in GOS than with growth in glucose, probably due to the complex nature of GOS. The degradation of GOS requires a wide range of enzymes and transporters for uptake into the cell, and these are discussed below.
|
View this table: [in a new window] |
TABLE 3. B. longum induced genes by growth on GOS compared to glucose medium
|
![]() View larger version (14K): [in a new window] |
FIG. 5. Venn diagram showing the number of differential and commonly hybridized clones of GOS compared to glucose (Gl).
|
Shared upregulated genes in GOS medium with human or formula milks.
There was substantial overlap in the upregulated carbohydrate utilization genes between the GOS medium and human milk. Other putative genes for β-galactosidase (BL0978), galactoside transport (BL0976), galactose-1-phosphate uridylyltransferase (BL1211), glycosyltransferase (BL1674), and possible lactaldehyde reductase (BL1673) were detected and were also upregulated in human milk (Tables 1 and 3), indicating specifically induced carbohydrate metabolism. Genes BL1638 and BL1639, which form part of a predicted transporter in a novel gene cluster for metabolism of mucin sugars and human milk oligosaccharides as described above (Table 3) (20), and other genes for transport of sugars (BL1164-BL1165) were significantly upregulated in GOS medium and human and formula milks compared to glucose medium.
Conclusions.
The gene expression of B. longum, a common intestinal bacterium, was studied to gain insight into the effects of human breast milk versus formula milk on its functionality and associated health effects within the human intestine. Transcriptomes for human milk and GOS medium showed substantial overlap in the upregulation of carbohydrate utilization genes, presumably due to the large amount of GOS and galactose in these substrates. In contrast, there were fewer genes for sugar utilization shared between the human and formula milks among the clones sequenced despite the GOS-lcFOS (1:9) oligosaccharides present in the formula milk. Moreover, genes involved in carbohydrate metabolism formed the most dominant group specifically upregulated in breast milk. The latter included putative genes for N-acetylglucosamine degradation (BL1642) and metabolism of mucin and human milk oligosaccharides via the novel galactose/lacto-N-biose operon with associated sugar transporters (BL1638 and BL1639), of which only the transporters were also induced by GOS and formula milk with GOS-lcFOS. This supports the reported findings that the bifidogenic effect of human milk is to a great extent based on its unique and complex oligosaccharides. It was interesting that the putative genes for cell surface type 2 glycoprotein-binding fimbriae that are implicated in attachment and colonization in the intestine appeared to be upregulated only in the milks. Importantly, the present study generates leads for further investigation of genes involved in the metabolism and colonization of B. longum in the infant human intestine.
We are very grateful to volunteers for providing milk. We thank Michiel Wells for assistance with database searching and are grateful to Friesland Foods Corporate Research (Deventer, The Netherlands) for their kind gift of purified GOS.
Published ahead of print on 6 June 2008. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
R.G. and E.S.K. contributed equally to this work. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»