The establishment of the intestinal microbiota during and after birth is critical for the maturation and development of important body systems, including the digestive, immune and central nervous system(Reference Ratsika, Codagnone and O’mahony1) Disruptions in early life microbial development have been linked to a variety of diseases including increased risk of obesity, asthma and neurodevelopmental disorders(Reference Sarkar, Yoo and Dutra2). Following birth, the appearance of a diverse community of microbes takes place in a process called ecological succession(Reference Bäckhed, Roswall and Peng3). Early colonisers are deemed especially influential to this succession and determine the success of subsequent colonisers(Reference Laforest-Lapointe and Arrieta4) Environmental factors, particularly diet, are key determinants of the early colonisers and subsequent intestinal ecosystem.
Breast milk is the ideal source of nutrients, energy and bioactive compounds for newborns and is vital for postnatal growth and development(Reference Kelly and Coutts5–Reference Butte, Garza and Smith7). Infant formula is the alternative to provide or complement optimal nutrition for newborns when breast-feeding is not available or is inadequate. Human milk oligosaccharides (HMO) are the third most abundant solid component (> 10 g/l) after lactose (> 60 g/l) and lipids (> 30 g/l) in breast milk(Reference Zivkovic, German and Lebrilla8–Reference Lawson, O’Neill and Kujawska11). While non-digestible by host enzymes, HMO can be utilised by certain human Bacteroides, Akkermansia muciniphila and Bifidobacterium strains and rodent Enterococcus gallinarum strains, which are considered principal drivers of diet–microbe interactions in early life(Reference Zivkovic, German and Lebrilla8,Reference Lawson, O’Neill and Kujawska11–Reference Lewis, Totten and Smilowitz14) . HMO also partially incorporate into the systemic circulation, commonly detected in the plasma and urine of breast/HMO-fortified formula-fed infants and suckling rats(Reference Vazquez, Santos-Fandila and Buck15–Reference Marriage, Buck and Goehring19).
HMO profiles in breast milk, mainly fucosylated HMO, depend on the mother’s glycosyltransferase FUT2 ( Se) and FUT3 (Le) genotypes(Reference Totten, Zivkovic and Wu20). Mothers with a functional FUT2 gene allele, known as secretors, are able to synthesise α1–2 fucosylated HMO, and those with functional FUT3 gene allele, known as Lewis positives, secrete high content of α1–3/-4 fucosylated HMO(Reference Totten, Zivkovic and Wu20–Reference Difilippo, Pan and Logtenberg23). In multiple birth cohorts, the fucosylated HMO in secretors constitute up to 45–60 % of total HMO and significantly differ in concentration and composition from that of non-secretors(Reference Charbonneau, O’Donnell and Blanton10,Reference Lewis, Totten and Smilowitz14,Reference Liu, Yan and Wang24–Reference Siziba, Mank and Stahl27) . The overall intestinal microbiota, faecal glycan profiles, and growth and clinical characteristics of the infants from these cohorts, however, did not vary greatly based on their mother’s secretor status(Reference Lewis, Totten and Smilowitz14,Reference Liu, Yan and Wang24,Reference Moossavi, Atakora and Miliku26,Reference Sprenger, Lee and De Castro28) . The extensive coexistence of multiple Bifidobacterium species in infants’ microbiota, which is generally considered as a symbiotic relationship between host and the gut microbiota obtained through long-term coevolution, likely reflects their resilience to the variable HMO profiles in mothers of varying secretor status(Reference Lawson, O’Neill and Kujawska11,Reference Lewis, Totten and Smilowitz14,Reference Moossavi, Atakora and Miliku26,Reference Schell, Karmirantzou and Snel29,Reference Duranti, Lugli and Milani30) .
Beyond the fucosylated HMO, the sialylated HMO count for 15–25 % of total HMO and have been shown to be significantly less abundant in the milk of Malawian mothers with severely stunted infants(Reference Charbonneau, O’Donnell and Blanton10,Reference Liu, Yan and Wang24) . This association with growth was further verified by the administration of purified sialylated bovine milk oligosaccharides to germ-free mice and pigs colonised with a bacterial strain consortium from an undernourished Malawian infant. Supplementation was able to produce significant growth promotion and metabolic changes in liver, muscle and brain(Reference Charbonneau, O’Donnell and Blanton10). Similarly, adding purified sialylated bovine milk oligosaccharides improved bone growth in germ-free mice colonised with gut microbiota from a 6-month-old stunted infant by decreasing osteoclastgenesis while sparing osteoblast activity, actions which the fucosylated HMO 2′FL failed to elicit(Reference Cowardin, Ahern and Kung31).
The Malawian study highlights the structural specificity of milk oligosaccharides and the potential implications for infant growth and development. From the perspective of long-term evolution, the fucosylated milk oligosaccharides are less likely central to human lactation due to their significant variations across populations, while the unusual presence of relatively consistent sialylation of milk oligosaccharides may be more influential in the fundamental development of newborns(Reference Tao, Wu and Kim32,Reference Hobbs, Jahan and Ghorashi33) . It is also worth noting that the growth-promoting effects of sialylated HMO have been observed and validated in animal models of undernutrition(Reference Charbonneau, O’Donnell and Blanton10,Reference Román Riechmann, Moreno Villares and Domínguez Ortega34,Reference Chouraqui35) . Although the addition of these components has the potential to bring infant formula closer to human breast milk, due to the diversity of structures and complexity of purification, before the industrialised manufacture of several HMO in 2015, no single HMO was used in commercial infant formula(Reference Bych, Mikš and Johanson9,Reference Zhou, Jiang and Wang36) . In most clinical HMO studies, the two prominent HMO, 2′FL and lacto-N-neotetraose (LNnT), have been supplemented into infant formula and resulted in similar growth (body and head circumference) to that of normal formula/breastfed infants(Reference Marriage, Buck and Goehring19,Reference Román Riechmann, Moreno Villares and Domínguez Ortega34,Reference Puccio, Alliet and Cajozzo37) .
To our knowledge, in vivo studies comparing the function of fucosylated and sialylated HMO are limited. To test our hypothesis that fucosylated and sialylated HMO exert different impacts on neonatal growth and microbial development, we used an artificial rearing system called the pup-in-a-cup model(Reference Puiman and Stoll38–Reference Beierle, Chen and Hartwich40) to examine the effects of two HMO alone or in combination in rats from postnatal days 4 to 21. Two commercially available HMO, 2′-O-fucosyllactose (2′FL) and 3′-O-sialyllactose (3′SL) were selected as representative fucosylated and sialylated HMO.
Materials and methods
Preparation of rat milk substitutes
Basal rat milk substitute was prepared with ingredients listed in Table 1 and served as the control diet (CTR). Three HMO interventions were prepared with basal milk substitute and the addition of 1·2 g/l 2′-O-fucosyllactose (96·1 w/w%, 2′FL), or 1·2 g/l 3′-O-sialyllactose (97·5 w/w%, 3′SL), or 0·6 g/l 2′-O-fucosyllactose and 0·6 g/l 3′-O-sialyllactose (2’FL + 3’SL). The basal milk substitute served as an appropriate control formula that was based on mature cows’ milk which has an extremely low oligosaccharide concentration (0·1 g/l) compared with mature human milk (5–15 g/l)(Reference Quinn, O’Callaghan and Tobin41). The concentration of HMO in the present study was selected based on the average intake of total fucosylated or sialylated lactose (∼500 mg/kg per d, including 2′FL; 3′FL; 3′SL; 6′SL; DFLac, difucosyl-lactose) in infants at 4–12 months of age(Reference Plows, Berger and Jones42,Reference Larsson, Lind and Laursen43) . Ingredients were homogenised by immersion dispersers (Kinematica AG) at a speed of 15,500 rpm/s for 10 min × 3 times to avoid delamination. The homogenised rat milk substitutes were placed in 90 ml sterile containers and stored at 4°C until used. Rat milk substitutes were prepared every 3 d following the same protocol.
Ingredients found in the rat milk substitute
* Carnation evaporated milk provided (per 100 g) 134 kcal; saturated fat 4·6 g, polyunsaturated fat 0·2 g, monounsaturated fat 2·3 g; cholesterol 29 mg; Na 106 mg; K 303 mg; total carbohydrate 10 g and protein 7 g, 3′SL < 20 mg, 6′SL < 2 mg, 2′FL < 1 mg. SL, sialyllactose, FL, fucosyllactose.
† Mazola maize oil provided (per 14 g) 120 kcal; saturated fat 2 g, polyunsaturated fat 7 g, monounsaturated fat 4 g and vitamin E 2 mg.
‡ Vitamin mixture AIN-93VX provided (per 1 g) thiamin HCl 6 mg, riboflavin 6 mg, pyridoxine HCl 7 mg, niacin 180 mg, Ca 16 mg, folic acid 2 mg, biotin 0·2 g, cyanocobalamin (vitamin B12) 25 g, vitamin A palmitate 4000 μg, vitamin E acetate 75 μg, vitamin D3 1000 μg and vitamin K1 0·75 mg.
Animals and establishment of pup-in-a-cup model
Ethical approval was granted by the University of Calgary Animal Care Committee (no. AC19-0104) and followed the guidelines of the Canadian Council on Animal Care. Examination of both male and female rats is important in nutrition research but based on the complexity of establishing the pup-in-a-cup model, the limited number of ‘cups’ that can be placed in the heated water bath at one time and in view of recent human clinical studies showing no sex differences in growth, head circumference and gastrointestinal tolerance(Reference Marriage, Buck and Goehring19,Reference Román Riechmann, Moreno Villares and Domínguez Ortega34,Reference Puccio, Alliet and Cajozzo37) with HMO supplementation, we examined only one sex (males) in the present study. Fifty-four male Sprague–Dawley rats at 4 d of age with an average body weight (BW) of 10–12 g were obtained from the University of Calgary Life & Environmental Science Animal Care facility and n 13–14 pups/group randomised to CTR, 2′FL, 3′SL and 2′FL + 3′SL. The entire artificial rearing system was placed in a humidity-controlled room maintained at 30°C with a 12-h reverse light–dark cycle.
The procedure for establishing the pup-in-a-cup model was as follows: (1) Prior to surgery, hypothermic anaesthesia was applied to pre-weighed pups; (2) A cannula made of polyethylene (PE)10 tubing was pierced through the lining of the cheek and reinforced in place (Fig. 1(a)); (3) After cannulation, all pups remained on a heating blanket until warmed to normal body temperature and were then individually placed in a pre-heated and labelled foam cup (11 cm diameter × 15 cm deep). The temperature inside the cups was maintained at 34–37°C during the first week and then gradually reduced by 2–3°C/week; (4) The cheek cannula was connected to a multi-channel syringe pump (NE-1200 Twelve Channel Syringe Pump, Bio-Lynx Scientific Equipment Inc.) programmed for the delivery of the milk substitutes. Cannulated pups were randomly allocated to one of four rat milk substitutes (CTR, 2′FL, 3′SL, 2′FL + 3′SL) that were delivered with the same flow speed adjusted based on age and daily average BW as flow speed = (0·35 + 0·02(Age-4)) × BW/feeding hours (Fig. 1(b)) with a 15-min stop flow interval programmed for every 2 h of feeding; (5) Milk substitutes were replaced twice a day with syringe replacement and a flush of the tubing.
Establishment of the pup-in-a-cup model. Pups at postnatal day 4 with cheek cannula in place (a), milk flow rate (% increase) (b), growth curve (c) and survival rate (d) of pups during 18 d of artificial rearing. No significant differences in daily body weight were observed between HMO interventions. Milk flow rate (% increase) = (flow speed at age (n) − flow speed at age (n–1))/ flow speed at age (n–1) × 100. Survival rate (%) = survival pup numbers/total pup number per diet group × 100. HMO, human milk oligosaccharide.
Daily monitoring, sample collection and measurements
Pups were weighed twice a day throughout the 18 d of artificial rearing (Fig. 1(c)). Bowel movements and food intake behaviours were monitored daily. A dose of 0·1–0·2 ml deoxycholate solution (0·5 g/100 ml) was given to the pups if abdominal bloating occurred. The pups were euthanised if the bloating had not subsided after 48 h of treatment. Three faecal samples (week (W)1, W2 and W3) for each pup were collected by homogenising the faeces collected at 5–7, 12–14, 19–21 d of age. Pups underwent an oral glucose tolerance test at 18 d of age, as previously described(Reference Chleilat, Klancic and Ma44). The fat mass and lean mass of pups at weaning were measured using dual-energy X-ray absorptiometry (Hologic ODR 4500; Hologic Inc.) under light anaesthetisation with isoflurane. Following a 3 h fast, animals were euthanised by over anaesthetisation and rapid decapitation. Blood, intestinal tissues and caecal digesta samples were collected, and the brain, liver and caecum weighed.
DNA extraction and gut bacterial community profiling
To evaluate the effects of HMO interventions on the development of faecal microbiota of pups throughout the 18 d of artificial rearing, total bacterial DNA was isolated from the W1, W2 and W3 faecal samples of thirty-six weanling pups (n 108). Prior to DNA extraction using FastDNA spin kit (MP Biomedicals), faecal samples were pre-treated with bead-beating (MP Biomedicals) for 40 s × 3 times. Purified DNA was quantified using Quant-iT™ PicoGreen™ dsDNA Assay Kit (Invitrogen) and diluted to 20 ng/µl for use.
The V3–V4 region of the16S rRNA gene was sequenced on 2 × 300 bp Illumina MiSeq platform (Centre for Health Genomics and Informatics) and analysed in QIIME2 platform using DADA2 for sequences quality control and denoising(Reference Bolyen, Rideout and Dillon45,Reference Callahan, McMurdie and Rosen46) (QIIME2 2020·11). Only the amplicon sequence variants with a frequency greater than 10 were retained for downstream analysis. The taxonomy of amplicon sequence variants was assigned by aligning to Silva 138 reference database. Amplicon sequence variants classified as genus Lactobacillus were filtered out and additionally aligned to Genome Taxonomy Database as references release 95 (GTDB, https://gtdb.ecogenomic.org/) to reflect the current taxonomy of Lactobacillaceae (Reference Zheng, Wittouck and Salvetti47). Additional quantitative PCR assays were conducted to quantify the abundance of total bacteria using a universal bacterial primer (forward: 5′-CGGCAACGAGCGCAACCC-3′ and reverse: 5′-CCATTGTAGCACGTGTGTAGCC-3′) targeting the conservation region of 16S rRNA gene(Reference Denman and McSweeney48,Reference Bahl, Bergström and Licht49) and E. gallinarum species using a primer (forward: 5′-TTACTTGCTGATTTTGATTCG-3′ and reverse: 5′-TGAATTCTTCTTTGAAATCAG-3′) targeting the E. gallinarum-specific region of gene sodA, as previously described(Reference Akazawa, Tsujikawa and Fukuda12,Reference Jackson, Fedorka-Cray and Barrett50) . Amplicon sequence was validated on Sanger Sequencing Economy platform (Centre for Health Genomics and Informatics).
Statistical analyses
Outcomes with multiple measurements over time were assessed using repeated-measures ANOVA in R version 4.0.0 (2020-05-18). Body composition, total bacteria abundance and organ weight at weaning and the AUC values from the 2-h oral glucose tolerance test were analysed using linear mixed-effects models in R version 4.0.0 (2020-05-18). Dietary treatment was treated as a fixed factor; rats were considered as experimental units and their random effect was removed. Diet and/or age effects on the relative abundance of 16S rRNA gene sequence variants and weighted UniFrac distances of faecal microbiota were analysed using the Kruskal–Wallis rank-sum test in R version 4.0.0 (2020-05-18). All data are presented as mean values with their standard error of the mean. P values with Bonferroni-adjustment < 0·05 were considered significant. To determine the adequate sample size to identify significant differences in the relative abundance of bacterial genera in the faecal matter of the pups, a sample size calculation utilising a power of 0·80 and α probability of 0·05 was performed using https://www.stat.ubc.ca/∼rollin/stats/ssize/n2.html. The effect size was estimated using the results from a previous study(Reference Chleilat, Klancic and Ma44). From the power analysis, it was determined that thirty-six rats (n 9/group) were needed to complete the study.
Results
Establishment of the pup-in-a-cup model
To formulate a rat milk substitute that is close to rat’s milk in macronutrients and energy, the basal rat milk substitutes used in the present study contained 12 g fat, 5·25 g protein and 7·5 g total carbohydrates. The calculated energetic density was 6359·68 kJ/l (152 kcal/100 ml), which is comparable with that of actual rat milk and the previously described formula for neonate rats 6485–8452 kJ/l (155–202 kcal/100 ml)(Reference Charbonneau, O’Donnell and Blanton10,Reference Autran, Schoterman and Jantscher-Krenn51,Reference Nicholas and Hartmann52) . The energetic content of the other three rat milk substitutes was considered the same as the basal milk substitute, as the 1·2 g/l HMO-related energy difference is as low as 19·6–20·4 kJ/l(Reference Flint53). The total volume of milk substitutes delivered increased from 4·75 (sem 0·12) ml/d at 4 d of age to 25·25 (sem 0·51) ml/d at weaning, corresponding to an increase in energetic intake from 30·17 (sem 0·75) kJ/d (7·21 (sem 0·18) kcal/d) to 160·58 (sem 3·26) kJ/d (38·38 (sem 0·78) kcal/d). A total of thirty-nine of fifty-four (72·2 %) pups survived until weaning (Fig. 1(d)) which did not differ according to treatment (P > 0·05). The average BW of pups increased from 13·56 (sem 0·34) g to 37·14 (sem 0·75) g (Fig. 1(b)). All the surviving pups opened their eyes at 12–14 d of age. The highest mortality was observed at 8–12 d of age (Fig. 1(d)). No dietary differences were observed in the age at which eyes opened or mortality (P > 0·05; Fig. 1(d)).
Effect of human milk oligosaccharides on physical outcomes of artificially reared suckling pups
To determine whether HMO differentially affected the growth and physiological development of artificially reared newborn rats, Lee index (Fig. 2(a1)), body composition (Fig. 2(a2,a3)) and organ weights (Fig. 2(b1–3)) at weaning, average daily weight gain (Fig. 2(c1)), as well as glucose tolerance (Fig. 2(c2,c3)) at 18 d of age were measured. The Lee index (weight 0·33/naso-anal length × 1000) of weanling pups(Reference Stephens54) did not differ between groups (P > 0·05). Mean weanling fat mass (%) and lean mass (%) reached 9·71 (sem 2·75) % and 72·49 (sem 6·38) %, respectively, and did not differ between groups (P > 0·05). At weaning, the brain, the liver and the caecum weighed 1·31 (sem 0·13) g, 1·49 (sem 0·23) g and 0·18 (sem 0·05) g, respectively, and did not differ by group (P > 0·05). The pups gained 0·37 (sem 0·03) g/d during the first 3 d (postnatal days 4–7) of artificial rearing and the rate of growth increased to 1·00 (sem 0·20) g/d and 2·4 (sem 0·28) g/d in days 8–14 and 15–21, respectively. Blood glucose concentrations were significantly affected by time (Fig. 2(c2), P = 0·001) but not affected by HMO interventions (P = 0·99), or the interaction of time × HMO interventions (P = 0·48). No dietary difference was observed in the total AUC (Fig. 2(c3)) for the 2-h oral glucose tolerance test, growth performance and physical outcomes described above (P > 0·05).
Effect of HMO interventions on physical outcomes in artificially reared suckling rats. Difference in Lee index (a1), fat mass % (a2), lean mass % (a3), brain weight (b1), liver weight (b2), caecum weight (b3), average daily weight gain (c1), glucose tolerance (c2) and glucose AUC (c3) in suckling rats reared with CTR: basal rat milk substitute; 2′FL: CTR + 1·2 g/l 2′-fucosyllactose; 3′SL: CTR + 1·2 g/l 3′-sialyllactose; 2’FL + 3′SL: CTR + 0·6 g/l 2′-fucosyllactose + 0·6 g/L 3′-sialyllactose. Data with different superscripts represent significant difference (P < 0·05). Lee index = (weight 0·33/length) × 1000 (g 0·33/cm). HMO, human milk oligosaccharide.
Stepwise development of gut microbiota in artificially reared pups
The inclusion of 1·2 g/l HMO in the rat milk substitute did not influence the structure of faecal bacterial communities (online Supplementary Table S1, Fig. 3(a), R = –0·01, P = 0·479). However, across the 18 d of artificial rearing with rat milk substitutes, the development of pups’ faecal microbiota was governed by age and demonstrated a stepwise pattern (online Supplementary Table S1, Fig. 3(b), R = 0·33, P = 0·001, Fig. 3(d)). The weighted UniFrac distances of W1, W2 and W3 faecal samples, which depict the dissimilarities of microbial community structure, were significantly different from each other (online Supplementary Table S1, Fig. 3(b); P > 0·05). The richness (observed features) and evenness did not differ between W1 and W2 faecal samples but were significantly increased in W3 faecal samples (P = 0·0001; online Supplementary Table S1). The total bacteria, quantified as Log10 copy number of 16S rRNA gene ng−1 genomic DNA using qPCR, were observed to be present at a consistent level across groups and did not differ significantly between time points (P = 0·284) or dietary treatments (P = 0·156) (mean overall = CTR, 7·79 (sem 0·09); 2′FL, 7·79 (sem 0·08); 3′SL, 7·91 (sem 0·11); 2′FL + 3′SL, 7·56 (sem 0·09)). E. gallinarum, recently shown to be a degrader of 3′SL in suckling rats(Reference Akazawa, Tsujikawa and Fukuda12,Reference Jackson, Fedorka-Cray and Barrett50) , was only detected in four out of eighty faecal samples assessed by qPCR (data not shown).
Influence of HMO interventions and age on faecal microbiota of artificially reared suckling rats. Principal coordinate analysis (PCoA) (a), (b) and ternary plots (c), (d) demonstrating the structural and compositional differences in faecal microbiota of suckling rats fed with different HMO (a), (c) or at different ages (b), (d). CTR, basal rat milk substitute; 2′FL, CTR + 1·2 g/l 2′-fucosyllactose; 3′SL, CTR + 1·2 g/l 3′-sialyllactose; 2′FL + 3′SL, CTR + 0·6 g/l 2′-fucosyllactose + 0·6 g/l 3′-sialyllactose. HMO, human milk oligosaccharide.
The inclusion of 1·2 g/l 2′FL did not influence the composition of faecal microbiota in artificially reared suckling rats (Fig. 3(c), P > 0·05). However, the presence of 3′SL (3′SL or 2′FL + 3′SL) in rat milk substitutes significantly reduced the relative abundance of genus Proteus in W1 faecal samples (Fig. 4(a), online Supplementary Table S2, P = 0·009) and genus Terrisporobacter in W3 faecal samples (Fig. 4(b), online Supplementary Table S2, P = 0·008). The relative abundance of twenty out of the thirty most abundant faecal bacterial genera of artificially reared suckling rats significantly changed (P < 0·05) over time (online Supplementary Table S2). The ternary plots (Fig. 3(d)) and boxplots (Fig. 5(a1)–(c3)) depicted the abundant changes of the representative genera. Specifically, in W3 samples, the abundance of Enterococcus (Fig. 5(a1)), which was the most abundant genus in W1 and W2 samples (> 50 %), decreased significantly (P = 0·0001). Conversely, the relative abundance of Escherichia-Shigella, Ligilactobacillus, Blautia, Lachnoclostridium, Clostridium, Clostridium sensu stricto 1, Terrisporobacter and Proteus in W3 samples was significantly (P < 0·05) higher than that in W1 and W2 faecal samples (Fig. 5(a2)–(c3), online Supplementary Table S2).
Inclusion of 3′SL reduced select bacterial genera in the faecal microbiota of artificially reared suckling rats. Relative abundance (%) of Proteus in week 1 (W1) samples (a) and Terrisporobacter in week 3 (W3) samples (b) in faecal microbiota of artificially reared suckling rats. Data with asterisk (*) are significantly different (P < 0·05) between diet groups. CTR, basal rat milk substitute; 2′FL, CTR + 1·2 g/l 2′-fucosyllactose; 3′SL, CTR + 1·2 g/l 3′-sialyllactose; 2′FL + 3′SL, CTR + 0·6 g/l 2′-fucosyllactose + 0·6 g/l 3′-sialyllactose.
Bacterial genera showing significant differences in faecal microbiota of artificially reared suckling rats. Data with asterisk (*) are significantly different (P < 0·05) between age groups.
Discussion
In the present study, we established an artificial rearing system to examine the effects of select HMO during the suckling period. We selected the most abundant fucosylated HMO, 2′FL and one of the abundant sialyllactose that typically remains stable during lactation, 3′SL to determine if supplementation with physiological doses of these HMO affected growth and microbial outcomes in rat pups over an 18-d investigative period(Reference Chaturvedi, Warren and Altaye55–Reference Michaelsen, Larsen and Thomsen58). With the exception of the potential pathogen inhibitory effects of 3′SL, no other distinct differences in growth performance and physiological outcomes were observed in the rats. Longitudinal changes in the faecal microbiota of artificially reared suckling rats were primarily governed by age and not affected by the presence of 2′-FL and/or 3′-SL in rat milk substitutes.
The pup-in-a-cup model allowed us to assess whether the addition of two HMO to the control milk substitute, meant to imitate human infants solely fed with infant formula, had any impacts on developmental parameters. Although it is challenging to assess the average milk intake of dam-fed suckling rats, we formulated a milk substitute that is close to rat milk in macronutrient content and energy density and increased the flow rate with age and BW until weaning(Reference Charbonneau, O’Donnell and Blanton10,Reference Autran, Schoterman and Jantscher-Krenn51,Reference Nicholas and Hartmann52) . Compared with inducing overnutrition or undernutrition in suckling rats by adjusting litter size at birth, the established artificial rearing system herein is capable of precisely controlling the nutritional supply and thus the nutritional status of newborns(Reference Seidler, Bell and Slotkin59,Reference Plagemann, Harder and Rake60) . The cheek cannula we adopted is less invasive and movement restricting than the intragastric cannula used in a previous pup-in-a-cup rodent model(Reference Penttila, Flesch and McCue61). However, feeding with a uniform flow rate via the cheek cannula is also a leading cause of death during the age of 6–13 d, the critical stage for structural and functional changes in lung development of newborn rats(Reference Roberts, Weesner and Bucher62,Reference Bucher and Roberts63) . Rat pups with poor adaptability to increasing flow rate (needed to maintain proper weight gain and growth) were prone to die from choking on the milk and suffocation when the flow rate increased sharply after the age of 8 d. In addition to selecting pups with similar BW, further improvement is needed in the model to fulfil the nutritional requirement of pups while reducing choking-associated mortality. More importantly, single-cup feeding enables the random allocation of dietary treatments within and between litters and eliminates the interference of litter effects on neonatal gut microbiota(Reference Azagra-Boronat, Massot-Cladera and Mayneris-Perxachs64,Reference Dai, Yang and Yuan65) ; therefore, the pups instead of the litters can be treated as experimental units, which greatly reduces the use of animals.
Neonatal rats are born with an immature gut with respect to the stage of functional and immunological development, which develops slowly during early and mid-lactation and matures rapidly around weaning(Reference Puiman and Stoll38,Reference Sangild66) . The general stepwise pattern of gut microbiota observed in normally suckled and now our artificially reared suckling rats consistently demonstrates these developmental milestones(Reference Akazawa, Tsujikawa and Fukuda12,Reference Azagra-Boronat, Massot-Cladera and Mayneris-Perxachs64,Reference Dai, Yang and Yuan65,Reference Matsui, Akazawa and Tsujikawa67) . The shared nature of the microbial shifts between breastfed and our artificially reared rats also points towards a similar developmental physiology of the intestine, particularly, the dramatic shift in the luminal environment towards anaerobic conditions around weaning. In normal suckled pups, the abundance of E. gallinarum, the major SL-degrading species recently identified in rats, markedly increases from 7 to 12 d and decreases thereafter, indicating a critical period of colonisation with this species(Reference Akazawa, Tsujikawa and Fukuda12,Reference Matsui, Akazawa and Tsujikawa67) . E. gallinarum, however, was absent in most of our pups, which were separated from their mothers at 4 d of age, indicating a possible interruption of the maternal–offspring microbial transmission. In addition to the distinct milk oligosaccharide-degrading bacterial consortium in infants (Bacteroides, Akkermansia muciniphila and Bifidobacterium) and suckling rats (E. gallinarum), the increase in strict anaerobes in exclusively breastfed infants, represented by HMO-degrading bifidobacteria, starts from the first days of life and dominates the gut microbiota from the first weeks of life(Reference Akazawa, Tsujikawa and Fukuda12,Reference Matsui, Akazawa and Tsujikawa67–Reference Korpela and de Vos69) . This high bifidobacteria abundance was not seen in our artificially reared rats, revealing a divergent developmental manner in intestinal physiology(Reference Akazawa, Tsujikawa and Fukuda12,Reference Matsui, Akazawa and Tsujikawa67–Reference Korpela and de Vos69) . These newfound differences identified by our pup-in-a-cup model underscore the need to further probe the limitations of the rat model for human-directed HMO microbiome research.
Despite differences in dominant HMO-degrading bacteria in rats and humans, the findings of a series of necrotising enterocolitis (NEC) studies in preterm infants and neonatal rats may shed some additional insights into the benefits of HMO(Reference Autran, Schoterman and Jantscher-Krenn51,Reference Schanler, Lau and Hurst70–Reference Masi, Embleton and Lamb75) . A lower risk of NEC has consistently been observed in breastfed preterm infants compared with their formula-fed counterparts(Reference Jantscher-Krenn, Zherebtsov and Nissan73–Reference Masi, Embleton and Lamb75). A neonatal rat NEC model has consequently been used to explore the reasons for this differential risk, which identified HMO, specifically, disialyllacto-N-tetraose as the determinant of the NEC-protective effects of mother’s milk(Reference Jantscher-Krenn, Zherebtsov and Nissan73). Two NEC-disialyllacto-N-tetraose clinical studies subsequently confirmed disialyllacto-N-tetraose deficits in the milk of mothers of infants who developed NEC and determined that a threshold level of 241 nmol/ml disialyllacto-N-tetraose in mother’s milk could be a potential biomarker to predict the risk of NEC development in preterm infants(Reference Autran, Kellman and Kim74,Reference Masi, Embleton and Lamb75) . Notably, associations between the relative abundance of Bifidobacterium longum and NEC were also shown in one of the two NEC cohort studies(Reference Masi, Embleton and Lamb75). Even though the mechanism of the NEC-protective effects of HMO is not yet fully understood, these results collectively imply a possible existence of a mechanism in humans and rats that is not mediated by a common HMO-degrading microbiome as a core factor. Instead, the presence of most HMO structures in serum and their excretion in urine could suggest a possible association between the systemic absorption of HMO and their general protective effects in infants and suckling rats(Reference Vazquez, Santos-Fandila and Buck15–Reference Marriage, Buck and Goehring19). Importantly, the reduced abundance of the potential pathogens, Proteus and Terrisporobacter, following 3′SL ingestion by our rats could be seen as protective given that Proteus has been associated with Crohn’s disease and inflammation(Reference Zhang, Hoedt and Liu76) and Terrisporobacter has been linked to oxidative stress and inflammation in preterm infants fed formula v. human milk(Reference Cai, Zhang and Morales77). Even with the lower dose of 3′SL in the combination treatment (0·6 g/l compared with 1·2 g/l in the 3′SL alone group), the inhibitory effect on Proteus and Terrisporobacter was still present. This might be attributed to the previously shown ability of 3′SL exposure to cause Caco-2 cells to change their surface glycan profile(Reference Bode78,Reference Bode79) . It is possible that even the lower dose in our combination treatment was able to modify the glycan content on the surface of the rats’ epithelial cells and receptor sites for these potential pathogens.
Compared with human and rat milk, mature bovine milk contains almost no non-fucosylated oligosaccharides and a significantly lower content of 3′SL (0·04–0·12 g/l)(Reference Hobbs, Jahan and Ghorashi33,Reference Tao, DePeters and German80,Reference Nakamura, Kawase and Kimura81) . The basal rat milk substitute we used was formulated based on evaporated bovine milk, which would consequently mimic the 2′FL and 3′SL deficits in our artificially reared pups. Overall, when supplementing with these two HMO, we observed the potential pathogen inhibitory effect of 3′SL but did not identify any differences in growth performance, body composition and organ weights in pups fed with 2′FL or 3′SL milk replacer. Although this lack of difference in growth is similar to human studies with HMO supplementation(Reference Marriage, Buck and Goehring19,Reference Román Riechmann, Moreno Villares and Domínguez Ortega34,Reference Puccio, Alliet and Cajozzo37) that included both male and female infants, we only examined male rats and it would be important to repeat this work with female rat pups. Given previous research(Reference Autran, Kellman and Kim74), it is possible that the HMO could have affected immune response without affecting growth and the other parameters we examined; however, we did not assess any immune parameters. Furthermore, the reduction in potential pathobionts such as Terrisporobacter, while important in healthy normally developing infants, could be even more relevant in preterm or undernourished infants. It would be important in future studies to examine a broader scope of metabolic and immunologic outcomes(Reference Charbonneau, O’Donnell and Blanton10). Given that LNnT, difucosyllactose (DFL), lacto-N-tetraose (LNT) and 6′SL can now be produced on a commercial large scale(Reference Bych, Mikš and Johanson9), using additional HMO structures in future studies would be helpful in identifying whether there are more structure-specific attributes of HMO. In addition, given the relevance between the protective effects of sialylated HMO and host abnormalities (premature birth, extremely low BW or malnutrition)(Reference Charbonneau, O’Donnell and Blanton10,Reference Jantscher-Krenn, Zherebtsov and Nissan73–Reference Masi, Embleton and Lamb75) , further investigation into the protective effects of specific HMO on the host under abnormal conditions may have greater clinical significance in neonatal research.
Acknowledgements
The authors wish to thank Shelly Wegener and Dr. Richard Pon, Centre for Health Genomics and Informatics at the University of Calgary, for support with the 16S rRNA sequencing.
This work was supported by a research grant from the Natural Sciences and Engineering Research Council of Canada (NSERC RGPIN/03773-2016). W. W. is supported by a University of Calgary Eyes High Postdoctoral Fellowship.
W. W. designed research, performed animal experiments, prepared samples, performed bioinformatics analysis and wrote paper. C. M. and J. S. helped with the pup-in-a-cup model and provided resources. N. A. C., E. W. N. T., D. E. L., F. C., K. M. S. and K. S. assisted with animal experiments. R. A. R. designed research, assisted with interpretation, obtained funding and had primary responsibility for final content. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Supplementary material
For supplementary material referred to in this article, please visit https://doi.org/10.1017/S0007114521005146
