Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES
Food and Beverage Chemistry/Biochemistry
Capillary electrophoresis analysis of bovine milk oligosaccharides permits an assessment of the influence of diet and the discovery of nine abundant sulfated analogues Sara Vicaretti, Nina Mohtarudin, Alexander Garner, and Wesley F. Zandberg J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b01041 • Publication Date (Web): 10 May 2018 Downloaded from http://pubs.acs.org on May 13, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 38
Journal of Agricultural and Food Chemistry
Capillary Electrophoresis Analysis of Bovine Milk Oligosaccharides Permits an Assessment of the Influence of Diet and the Discovery of Nine Abundant Sulfated Analogues. Sara D Vicaretti1, Nina A Mohtarudin2, Alexander M Garner2, and Wesley F Zandberg1* 1
Department of Chemistry and 2Biochemistry, The University of British Columbia, Okanagan
* Corresponding author contact details: Wesley F Zandberg Tel: 250-807-9821; Fax: 250-807-8001; Email:
[email protected] 1 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 2 of 38
1
Abstract.
2
Bovine milk oligosaccharides (BMOs), like their analogues in human milk, have important
3
prebiotic functions. Environmental factors have previously been linked to variation in BMO
4
structures, thus to test the hypothesis that the bovine diet may lead to these changes in relative
5
BMO abundances, a rapid capillary electrophoresis (CE)-based work flow was developed to
6
profile the BMOs extracted from the milk of cows fed distinctly different diets. Over the first
7
week of lactation, few significant differences were observed between the different diet groups,
8
with the dominant changes being clearly linked to lactation period. CE analyses indicated the
9
presence of ten unusually anionic BMOs which were predicted to be phosphorylated and sulfated
10
species. Nine unique sulfated BMOs were detected by high-resolution accurate mass
11
spectrometry, none of which have been previously described in bovine milk. The biosynthesis of
12
these was in direct competition with 3’-sialyllactose, the most abundant BMO in bovine milk.
13 14
Keywords: milk; sulfated oligosaccharides; capillary electrophoresis; mass spectrometry; diet.
2 ACS Paragon Plus Environment
Page 3 of 38
15
Journal of Agricultural and Food Chemistry
Introduction
16
Oligosaccharides are among the most abundant components of human milk, being found
17
in both their free forms, typically with the disaccharide lactose found at the reducing end, and
18
conjugated to other biomolecules such as milk proteins and lipids. The abundance of free,
19
unconjugated human milk oligosaccharides (HMOs) is especially notable in light of the fact that
20
these are indigestible to neonates. Though not directly used to meet the newborns’ caloric needs,
21
HMOs are important bioactive molecules with numerous health-promoting roles.1,2 For example,
22
HMOs function as prebiotics that encourage the colonization of the infant gastrointestinal (GI)
23
tract with a protective microbiome.3,4 In addition, HMOs often bear glycoepitopes that are also
24
observed on the protein-linked oligosaccharides comprising the mucus lining the gut, permitting
25
them to serve as soluble decoys capable of blocking the attachment of pathogens to the intestinal
26
epithelium. Recently, it has been demonstrated that HMOs may bind to lectins expressed on the
27
surface of dendritic cells5–7 in the intestinal epithelium, and in this capacity they may directly
28
influence inflammatory processes8 critical to both the short- and long-term health of neonates.
29
The ability of HMOs to engage cell surface lectins, coupled with evidence of systemic
30
absorption from the GI tract,9 is also consistent with observations that these compounds have
31
important functions beyond their roles in the GI tract.10
32
The numerous biological functions of HMOs are clearly dependent on their specific
33
chemical structures. To date, over 250 HMOs have been detected (although many still await
34
structural identification),1,2 most of which are based on a lactose core structure (Gal-β1,4-Glc,
35
where Gal and Glc are the monosaccharides D-galactose and D-glucose, respectively) that has
36
been further elongated with Gal, N-acetyl-D-glucosamine (GlcNAc), L-fucose and 5-N-acetyl-D-
37
neuraminic acid (Neu5Ac) residues (Figure S1). Like HMOs, bovine milk oligosaccharides
3 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 4 of 38
38
(BMOs) are composed of the same five monosaccharides, in addition to the non-human,
39
mammal-specific monosaccharide 5-N-glycolyl-D-neuraminic acid (Neu5Gc; Figure S2).
40
Relative to HMOs, BMOs are found at roughly 20-fold lower absolute concentrations, and have
41
higher relative levels of acidic structures containing Neu5Ac/Neu5Gc in preference over the
42
fucosylated analogues that are more prevalent in human milk. To date, over 50 distinct BMOs
43
have been detected,11,12 a number of which are structurally identical to HMOs. Although BMOs
44
have not received as much research attention, a clearer understanding of their structures and the
45
factors affecting their biosynthesis is warranted given that they may confer HMO-like health
46
benefits to milk and dairy product consumers13 and have been established as viable additives to
47
infant formulas, that are otherwise nearly devoid of HMOs/BMOs, prompting efforts to recover
48
BMOs from dairy processing streams.14
49
The dependence on the health-promoting effects of milk oligosaccharides on their
50
structures, and the notable variation in structures observed between women, not to mention the
51
differences between human and bovine milk, has prompted research into the biochemical basis
52
for milk oligosaccharide variation and the subsequent impact on infant health.15,16 It is known
53
that genetic factors are responsible for much of the observed variation in the relative HMO
54
concentrations in the milks collected from healthy mothers; these are thought to have a
55
significant impact on infant.4,13,16,17 Genetic sources of variation in milk oligosaccharides from
56
other mammals is less clear. Sundekilde et al. have identified differences in the BMO pools in
57
milk from Jersey and Holstein-Friesian cows, with the former producing higher relative amounts
58
of Fuc- and Neu5Ac-containing structures.18 Meanwhile, the non-genetic factors influencing
59
milk oligosaccharide composition in either humans or cows still await a more rigorous
60
characterization. It is known that the levels of both HMOs15 and BMOs19,20 decrease during the
4 ACS Paragon Plus Environment
Page 5 of 38
Journal of Agricultural and Food Chemistry
61
course of lactation. Liu et al. have recently demonstrated that BMO levels exhibit a seasonal
62
variability, with generally higher levels observed in the autumn.21 A similar trend has also been
63
observed in human milk, with total HMO levels significantly decreasing in the milk of mothers
64
when food resources are more scarce.16 These studies, paired with evidence of clear variations
65
among HMO levels observed between ethnically (and presumably genetically) similar women
66
living in differing locations,22 suggests that the environment may play an important role in
67
regulating the biosynthesis of milk oligosaccharides. We hypothesize that the maternal diet may
68
represent a significant non-genetic factor impacting both HMO and BMO biosynthesis, thus
69
explaining the seasonal and geographical trends noted to date. Such a dietary link would also be
70
consistent with reports indicating that over nourishment/obesity,1 malnourishment13 and
71
hyperglycemia1 all impact the structures and/or relative concentrations of HMOs.
72
To test the hypothesis that the maternal diet impacts milk oligosaccharide biosynthesis
73
we compared the BMOs produced by dairy cows fed different diets over the first week of
74
lactation. It was reasoned that any correlations would be more readily detectible in dairy cows
75
since these lack the obvious genetic differences observed amongst HMOs, and since, unlike
76
people, dairy cattle within a single farm are all fed essentially identical diets. A recent study
77
comparing pasture-only-fed cows on organic or conventional diets has indicated that the bovine
78
diet may have had a slight impact on the levels two (of 11 analyzed) BMOs 23; we reasoned that
79
BMO differences would be more readily apparent by comparing cows fed grass alone to those
80
whose diets also included corn and grain. A second aim of this research was to develop a rapid
81
analytical work flow permitting the accurate investigation of the biosynthetic relationships
82
between BMOs (or HMOs). In this respect, it should be noted that while mass spectrometry
83
(MS) is one of the most sensitive and chemically-informative techniques available for
5 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 6 of 38
84
BMO/HMO analysis, the vastly different response factors for different oligosaccharides—for
85
example, neutral versus acidic (i.e., Neu5Ac/Neu5Gc-containing) species—complicates a clear
86
assessment of the biosynthetic relationships between BMOs present within a single sample,
87
especially in the absence of suitable stable isotopic analogues.24 A solution to this challenge is to
88
analyze these compounds after labelling their reducing ends with fluorophores (or
89
chromophores) to permit their detection by optical techniques that are insensitive to the BMO
90
structures4,11,22. Thus, a capillary electrophoresis (CE) method employing laser-induced
91
fluorescence (LIF) detection was used here to rapidly (within 8 min) establish the relative
92
quantities of 34 BMOs shared in common with all the individual milk samples collected; these
93
BMOs were also analyzed in pooled samples collected from three dairy farms and a range of
94
processed liquid milk samples sold by local grocery stores. High-performance liquid
95
chromatography (HPLC) coupled with high-resolution accurate mass-MS was also utilized to
96
characterize unique BMOs identified by CE-LIF, at least ten of which have not yet been
97
previously described in either bovine or human milk.
98
Materials and Methods
99
Chemical and General Details. The following chemicals were purchased from Sigma Aldrich
100
(Saint Louis, MO, USA) and used as received: HPLC- or analytical grade acetonitrile (ACN),
101
methanol (MeOH), chloroform (CHCl3), trifluoroacetic acid (TFA), sodium acetate (NaOAc),
102
sodium cyanoborohydride (NaBH3CN), dimethylsulfoxide (DMSO), acetic acid (AcOH), and
103
4,5-dimethyl-1,2-benzyldiamine (DMBA). The following chemicals were purchased from Alfa
104
Aeser or J.T. Baker (Thermo Fisher Scientfic; Tewksbury, USA; Ontario, Canada) and used as
105
received: Sodium hydroxide (NaOH), and ethanol (EtOH). All BMO and/or HMO standards
106
(Figure S1) were purchased from Dextra Laboratories Inc. (Reading, UK) or Carbosynth Ltd.
6 ACS Paragon Plus Environment
Page 7 of 38
Journal of Agricultural and Food Chemistry
107
(Compton, UK). Unless otherwise noted, 18 MΩ water was provided by a Barnstead E-Pure
108
water purification system (Thermo Fisher Scientific; Waltham, MA, USA) for preparing all
109
aqueous solutions, including uHPLC mobile phases.
110
(APTS) was synthesized exactly according to a published procedure25 with the exception that the
111
total reaction time was reduced from 30 h to 3 h. APTS was further purified over a graphite
112
column26 and quantified by fluorescence, using commercially-available APTS (SCIEX;
113
Washington, DC, USA) as a calibration standard. APTS was stored as a 100 mM stock in 0.9 M
114
citric acid at -20 °C. A DMBA stock solution (240 mM in 2M AcOH) was prepared freshly
115
before use. The following enzymes were purchased from New England Biolabs (Ipswich, MA,
116
USA) and used exactly as described in the manufacturer’s instructions: calf intestinal alkaline
117
phosphatase,
118
perfringens) was obtained from Sigma-Aldrich.
119
Milk Sample Collection and BMO Preparation. Bovine milk samples were collected at Grass
120
Roots Dairies (Salmon Arm, B.C., Canada) and Hutley Acres Dairy (Armstrong, B.C., Canada)
121
in 50 mL plastic milk testing tubes; samples were frozen immediately after collection,
122
transported to the lab on ice and stored at -20 °C prior to thawing at 4 °C before extraction. Cows
123
at Grass Roots were fed exclusively grass diets; those at Hutley Acres were fed a diet consisting
124
of corn and alfalfa silage, earlage, and grain. Milk from three animals per farm was collected:
125
two Holsteins and a Holstein-Jersey cross (Hutley) and two milking Shorthorn and a Shorthorn-
126
Holstein cross (Grass Roots). Three time points were sampled: as soon as possible post-calving
127
(day 0), after 1 day (i.e., 24 h post-calving) and one week. At each farm, milk was collected from
128
sequentially freshening cows to ensure that the bovine diets were as similar as possible. A single
129
sample (labelled “tank”) was also collected from the storage tank at each farm; these pooled
β-1,3-galactosidase,
and
8-aminopyrene-1,3,6-trisulfonic acid
β-hexosaminidasef.
Neuraminidase
(Clostridium
7 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 8 of 38
130
samples were compared with another unprocessed milk sample (Riverbreeze Dairy, Armstrong,
131
B.C., Canada) in addition to processed milk samples obtained from a local (Kelowna, B.C.,
132
Canada) grocery store. All samples were collected in the late spring/early summer.
133
BMO-rich fractions were prepared from whole milk by liquid-liquid extractions as
134
described elsewhere.12,18,19,27 A 2 mL aliquot of milk was mixed with 400 µL of 18 MΩ water
135
and centrifuged for 30 min at 15,100 × g at room temperature. The fat (top layer) was removed,
136
leaving a protein- and oligosaccharide-rich bottom layer that was mixed with 2:1
137
chloroform:methanol (v/v) at a 4:1 solvent:sample ratio. This mixture was centrifuged for 30 min
138
at 3,270 × g at 4 °C and the top, oligosaccharide-rich layer was removed and mixed with ethanol
139
at a 2:1 (ethanol:sample) ratio. This mixture was stored overnight at -42 °C to facilitate the
140
precipitation of any remaining proteins, which were pelleted by centrifugation (30 min, 3,270 ×
141
g, 4 °C). The supernatant was dried in vacuo using a SpeedVac Concentrator (Thermo-Fisher
142
Scientific) at ambient temperature. Further BMO purification, including the elimination of the
143
majority of the lactose present (Table S1), was performed using pre-conditioned graphitized
144
carbon SPE cartridges26 (Supelco, ENVICarb, 200 mg). Acidic- (50% ACN/0.1% TFA) and
145
neutral- (20% ACN) BMO-containing SPE fractions were combined and dried in vacuo.
146
Desalted BMOs were re-dissolved in 18 MΩ H2O and lyophilized in five equivalent fractions for
147
(1) determining total reducing sugar levels, (2) neutral monosaccharide analysis, (3)
148
quantification of total Neu5Ac/Neu5Gc levels, (4) BMO profiling by CE-LIF and (5) targeted
149
HPLC-MS analysis.
150
BMO Analyses. (1) Reducing sugars were quantified by absorbance according to the method of
151
Mopper and Grindler.28 (2) BMOs equivalent to 10 nmol reducing ends were fortified with 5
152
nmol L-rhamnose (Rha) and D-N-acetylgalactosamine (GalNAc), as internal standards, 8 ACS Paragon Plus Environment
Page 9 of 38
Journal of Agricultural and Food Chemistry
153
lyophilized in 200 µL tubes, and subsequently hydrolyzed to their constituent monosaccharides
154
using 100 µL 2 M TFA (100 °C, 2 h)29 after which they were snap frozen in liquid nitrogen and
155
concentrated using a SpeedVac. Samples were re-N-acetylated29 and labelled with APTS,25 at 60
156
°C for 2 h in the dark (Figure S3). Monosaccharides were separated and quantified by CE-LIF
157
(ProteomeLab PA800; Beckman-Coulter) under normal polarity (+30 kV) using 240 mM
158
NaBO4, plus 0.1% (w/v) polyethylene glycol, pH 9.5 as the background electrolyte. All
159
monosaccharide peaks were identified by means of external standards and normalized to the Rha
160
peak. The ratios of the monosaccharides Gal, Fuc, and GlcNAc were all determined relative to
161
the single reducing end Glc residue of BMOs. (3) BMO-bound Neu5Ac and Neu5Gc were
162
quantified exactly as described in the accompanying manuscript.30 (4) BMOs were labelled with
163
APTS (16 h, room temperature) exactly as previously described (Figure S3).25,26 Labelling
164
reactions were quenched by adding 100 µL H2O before all samples were analyzed by CE-LIF
165
using reversed polarity (-30 kV) and a proprietary N-CHO separation buffer (SCIEX) as the
166
background electrolyte. All peaks not attributable to the APTS reagents were manually
167
integrated using 32 Karat (version 8.0) software. Some APTS-labelled samples treated with
168
glycosidases. In these assays 17 µL of sample were mixed with 2 µL 10x reaction buffer
169
(provided with each enzyme) and 1 µL enzyme, and incubated overnight at 37 °C after which
170
they were analyzed, without work up, by CE. (5) BMO samples were prepared for uHPLC-MS
171
analysis by first reducing them (1 M NaBH4 in 50 mM NH4OH; 2 h; 65 °C) to avoid the
172
resolution of α- and β-anomers;12 after neutralizing these reactions with AcOH, BMOs were
173
desalted by graphite SPE as described above. HPLC-MS parameters are described in Table S2.
174
Data Acquisition and Statistical Analysis. The HPLC-MS data acquisition and processing was
175
performed using the MassHunter Workstation software suite (Agilent Technologies), with the 9 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 10 of 38
176
following version numbers: Data Acquisition Workstation (v B.08.00). HMO structures were
177
characterized using the Agilent MassHunter Qualitative Analysis (v B.07.00, Service Pack 2)
178
find-by-formula algorithm using a ±10 ppm mass window. Data reduction and statistical
179
calculations were performed using Microsoft Excel (Microsoft Corporation; Redmond, USA).
180
CE peak areas were divided by their respective migration times and expressed as a
181
fraction of the total integrated area of all 34 peaks detectible in individual electropherograms,
182
excluding lactose. In several instances where peaks could not be detected, due to poor resolution
183
or low overall BMO abundance, the relevant region of baseline was manually integrated. Median
184
BMO levels at differing stages of lactation (n = 6 cows) were deemed statistically significant if
185
the two-tailed Mann-Whitney U values were lower than Ucrit at P < 0.05 or 0.01. Differences
186
between mean BMO levels in the two diet groups (n = 3 cows/group) were assessed using a two-
187
tailed Student’s t-test for heteroscedastic samples, with p < 0.05 considered significant.
188
Associations between BMOs were determined by calculating the two-tailed Spearman rank
189
correlations for all biological replicates (N = 18; 6 cows x 3 time points) with ρcrit values above
190
0.472 or 0.600 considered significant at p = 0.05 and 0.01, respectively.
191
Results and Discussion
192
Monosaccharide Analysis. To initially test the hypothesis that the bovine diet impacts BMO
193
structures and levels, the monosaccharide content of BMOs extracted from milk obtained from
194
cows fed exclusively grass (gr) or those fed a more typical diet of corn silage, earlage and grain
195
(c/g) were compared (Figure 1). As previously noted,19 a consistent drop in the total levels of
196
reducing sugars (Figure 1A) was observed over the first week of lactation, although after one
197
week the average BMO levels still exceeded those detected in the bulk storage tank at each farm.
198
However, bovine diet did not appear to influence total levels of BMOs observable at any time 10 ACS Paragon Plus Environment
Page 11 of 38
Journal of Agricultural and Food Chemistry
199
point. BMOs were completely hydrolyzed to their constituent monosaccharides which were
200
fluorescently labeled with APTS and analyzed by CE (Figure 1B) permitting their levels to be
201
determined relative to the reducing end Glc; note, however, that traces or residual lactose would
202
lead to an underestimation of the relative levels of these monosaccharides. The relative levels of
203
the neutral monosaccharides Fuc, Gal, and GlcNAc (Figure 1C) remained constant during the
204
first seven days of lactation and did not vary with bovine diet. Likewise, the acidic
205
monosaccharides (Neu5Ac and Neu5Gc), which were quantified by HPLC-MS, did not vary
206
significantly between the two diet groups, although they both rapidly decreased during the first
207
week of lactation (Figure 1D). It is, however, unclear based on these data only, if the declining
208
Neu5Ac/Neu5Gc levels were due to the general decrease in total BMO levels (Figure 1A) or a
209
more specific change in BMO biosynthetic pathways, such as a reduction in the levels or
210
activities of the required Neu5Ac transferases or, conversely, an increase in a competitive
211
biosynthetic process. Nevertheless, these data indicate that grass- versus corn/grain-containing
212
bovine diets do not noticeably affect the relative levels of monosaccharides comprising BMOs.
213
Characterization of Bovine Milk Oligosaccharides by Capillary Electrophoresis. CE-LIF
214
has proven to be an ideal separation technique for highly polar and/or charged analytes like
215
oligosaccharides, having previously been used to characterize HMOs,31 and porcine32 or equine33
216
milk oligosaccharides. Prior to CE analysis, oligosaccharides are typically labelled at their
217
reducing ends with the highly fluorescent and charged compound APTS25,26,29 (Figure S3)
218
enabling the efficient separation of both neutral and acidic oligosaccharides, and ensuring
219
equivalent response factors across all oligosaccharide structures. Here, CE-LIF was used as a
220
screening tool to rapidly profile the intact BMOs (Figure 2). Peak areas for thirty four BMOs,
221
encompassing all those that could be reliably detected in most milk samples, were determined
11 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 12 of 38
222
(peaks are numbered sequentially in Figure 2A). Nine individual BMOs were identified based on
223
their co-migration with authentic standards and chemical characteristics of others could be
224
deduced on the basis of their sensitivity to well-characterized exo-glycosidases (Table S3).
225
Relative BMO levels (Figure 2B and Figure S4) were determined by expressing each individual
226
peak area as a percentage of the total integrated area, i.e., the sum of all 34 peaks with the
227
exception of peak 14, which was composed of the co-migrating pair 6’SNL and lactose, traces of
228
which remained after the BMO extraction process. Six sialylated (i.e., Neu5Ac-containing)
229
BMOs were detected in all samples (Table S3), the most abundant of which was 3’SL, averaging
230
an excess of 29% abundance across all time points (Figure 2B). 3’SL was, in fact, the most
231
abundant BMO detected in accordance with previous studies.11,12,19,27 However, with few
232
exceptions such as BMO3 and BMO23 (Figure 2B), no significant differences were observed
233
between the relative levels of BMOs between the grass and grain/corn diet groups (Figure S4),
234
although the divergence amongst pooled samples of raw milk from three different farms, and
235
processed milks from local grocery stores (Table S4) suggests that some environmental influence
236
on BMO composition exists. Nevertheless, the stage of lactation was clearly the dominant source
237
of variation in relative BMO levels with several significantly increasing (BMO 3, 4, 5, 20 and
238
23) to the decrement of others (BMO 12 and 13; Figure 4). Decreases in BMO13 (3’SL) and
239
BMO12, the latter containing Neu5Ac/Neu5Gc as demonstrated by its sensitivity to
240
neuraminidase, indicate that the decrease in absolute Neu5Ac/Neu5Gc levels observed over the
241
first week of lactation (Figure 1D) must be at least partially due to alterations in BMO
242
biosynthetic pathways. In summary, our data do not support the hypothesis that the bovine diet
243
affects BMO structures, although notable differences in the pooled samples from both processed
12 ACS Paragon Plus Environment
Page 13 of 38
Journal of Agricultural and Food Chemistry
244
and non-processed sources suggest that a larger sample size may permit a more definitive
245
conclusion.
246
Characterization of Unexpectedly Acidic BMOs. The significant decrease observed in 3’SL
247
(BMO13) levels, which dropped from a median level of 59% relative abundance to 29% during
248
the first week of lactation, mirrored a nearly identical (2 to 24%), and also significant, increase in
249
BMO3—which was the second most abundant BMO observed in mature milk—over the same
250
time span, thus making the identification of this compound imperative. Also intriguing, was the
251
observation that all BMOs from 1-10 possess unexpectedly high ionic mobilities (Figure 2B).
252
Note that in CE, analyte mobility is dictated by both size (i.e., hydrodynamic radius) and charge,
253
with the larger, neutral BMOs eluting later than the smaller, more anionic ones eluting earlier
254
under the reversed polarity and zero electroosmotic flow CE conditions used here. Thus, for
255
example, attaching a single α-1,3-Gal residue to lactose (BMO14) creates a later eluting
256
trisaccharide (BMO18) while linking a Neu5Ac residue to lactose, which also generates a
257
trisaccharide (3’SL; BMO13), results in a faster migrating species than lactose itself, due to the
258
increased negative charge imparted by the Neu5Ac moiety. Thus, BMOs 1-10 represent species
259
that are apparently either smaller or more highly charged than 3’SL, prompting a more detailed
260
investigation of these structures initially based on their sensitivity (or resistance) to both
261
enzymatic or controlled chemical hydrolysis (Figure 4). It was reasoned that APTS-labelled
262
BMOs would be suitable substrates for neuraminidase or weak acid hydrolysis, both of which
263
would only cleave Neu5Ac or Neu5Gc residues (Figure 4A). Four Neu5Ac-containing BMO
264
standards indeed proved to be neuraminidase substrates, being completely hydrolyzed to either
265
lactose or lactosamine, both of which migrated more slowly than the parent BMOs due to the
266
loss of negative charge (Figure 4B). In contrast, when the BMOs extracted from a colostrum
13 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 14 of 38
267
sample were subjected to the same conditions, six BMOs (11-13, 17, 21, and 25) were obviously
268
hydrolyzed, as evidenced by their complete loss from the electropherogram produced from the
269
neuraminidase-treated sample (Figure 4C), while the levels of BMOs 1-10 remained essentially
270
unchanged (only BMOs 7-10 are shown here; BMOs 1 – 11 are depicted in Figure S6). Note that
271
under these hydrolysis conditions, no detectible lactosamine was produced indicating that in
272
these samples the lactosamine-based BMOs 3’ and 6’SLN only made minor contributions to the
273
total BMO pool.
274
Neuraminidases are known to be sensitive to both the regiochemistry of the Neu5Ac
275
glycosidic bond (α2,3 vs. 2,6 vs. 2,8) in addition to modifications to the Neu5Ac moiety (i.e., O-
276
acetylation) or the nature of the underlying glycan substrate.34–36 To account for the possibility
277
that some Neu5Ac-bearing BMOs could be resistant to neuraminidases, a BMO sample was
278
subjected to weak acid hydrolysis, under conditions which we30 and others37,38 have shown to be
279
sufficient to cleave acid-labile Neu5Ac/Neu5Gc glycosidic bonds. Under these conditions all
280
neuraminidase-labile BMOs previously identified were completely eliminated from the resulting
281
electropherogram, while BMOs 1-4, and 6 remained resistant (Figure 4D). Careful inspection of
282
the data indicated that BMO5 was hydrolyzed, apparently yielding BMO6 which increased
283
considerably in the acid-treated sample. The relative order of migration (BMO5 preceding
284
BMO6) is consistent with BMO5 losing negative charge (i.e., an acidic sugar residue). The slight
285
increase in BMO3 relative to, for example BMO4, may indicate that the former was also
286
produced under the acidic conditions, although in this instance the precursor BMO is less clear.
287
The resistance of BMO7-10 to acid hydrolysis was less obvious (Figure S6) as several new peaks
288
appeared in this region after hydrolysis preventing clear peak assignments.
14 ACS Paragon Plus Environment
Page 15 of 38
Journal of Agricultural and Food Chemistry
289
Since some of the highly mobile BMOs (1-4, and 6) appeared to be resistant to both
290
neuraminidase and weak acid hydrolysis, we hypothesized that the cause of the increased anionic
291
character (relative to the known Neu5Ac-containing BMOs) was due to the presence of
292
phosphate or sulfate moieties, both of which are smaller than Neu5Ac yet possessing, at the pH
293
of the CE buffer, an equivalent negative charge. Note that several phosphorylated lactose11 and
294
lactosamine analogues,11,39 including a Neu5Ac-bearing phospho-lactosamine,40 have been
295
previously reported. To test for the presence of phosphorylated species, an APTS-labelled BMO
296
sample was treated with alkaline phosphatase (Figure 4E and Figure S7). Under these conditions
297
a single peak, BMO2, appeared to be obviously hydrolyzed; since BMO2 was resistant to both
298
acid-hydrolysis and neuraminidase it can be concluded that this was likely a phospho-
299
lactosamine or phosphor-lactose analogue lacking a Neu5Ac residue. It was reasoned that the
300
significantly increasing levels of BMO3, apparently at the expense of 3’SL (BMO13; Figure 3),
301
suggested that the former peak corresponded to 3’-sulfo-lactose, which would be consistent with
302
this peak’s resistance to phosphatase, neuraminidase and acid-catalyzed hydrolysis and also
303
account for their biosynthetic competition. 3’-Sulfo-lactose, which was commercially-available,
304
was indeed observed to co-migrate with BMO3. To the best of our knowledge, the presence of
305
3’-sulfo-lactose in bovine milk has never previously been reported (although it has been
306
observed in canine milk41 and larger, sulfated HMOs42 have also been characterized) which is
307
striking in light of the fact that after one week of lactation it was the second most abundant
308
BMO, accounting for nearly a quarter of the total relative levels (Figure 2B and Table S3), and
309
was also detected at high levels in pooled samples of fresh and store-bought bovine milk (Table
310
S4).
15 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 16 of 38
311
Characterization of sulfated and phosphorylated BMOs by mass spectrometry.
312
Electropherograms of APTS-labelled BMOs indicated the presence of at least 10 BMOs with
313
mobilities exceeding that of 6’SL (BMO11), and with resistance to chemical and neuraminidase-
314
catalyzed hydrolysis, suggesting that these were phosyphorylated and/or sulfated species, a
315
hypothesis proven in part by sensitivity to alkaline phosphatase (BMO2) and co-migration with a
316
standard (BMO3; 3’-sulfo-lactose). HPLC-MS was therefore used to further assess the extracted
317
BMO samples for the presence of phosphorylated and/or sulfated analogues (Figure 5 and Table
318
1). Consistent with the CE data, abundant amounts of 3’-sulfo-lactose (Figure 5, peak g) were
319
observed, matching the accurate mass, retention time and collision induced dissociation (CID)
320
spectrum of an authentic standard (Figure S5). Two phosphorylated disaccharides analogues,
321
composed of an N-acetylhexosamine and a hexose residue, (peaks a and e) were also detected,
322
one of which is likely 6-phospho-lactosamine—since 3’sialyl-6-phospho-lactosamine, although
323
not detected here, has previously been reported (with the phosphate borne on the reducing
324
GlcNAc). The other phosphylated disaccharide is hypothesized to be either 3’ or 6’-phospho-
325
lactosamine.40 In addition, masses consistent with eight sulfated BMOs were detected by HPLC-
326
MS, all of which contained a lactose moiety at the reducing end. These sulfated BMOs included
327
a disaccharide hypothesized to be 6’-sulfo-lactose (peak d), three different trisaccharides (peaks
328
h – j), two tetrasaccharides, one Neu5Ac-containing trisaccharide (peak f), hypothesized to be 6’-
329
sulfo-3’-sialyllactose, and a single trisaccharide bearing two sulfate groups (peak c). The MS/
330
MS spectra of both 3’-sulfo-lactose (peak g) and the putative 6’-sulfo-lactose (peak d) contained
331
a prominent peak at m/z = 241.0023, consistent with the neutral loss of the reducing end Glc
332
residue (Figure S5A-B). An identical fragment was observed for trisaccharide (i), indicating that
333
this compound also contained the sulfate on the non-reducing (rather than internal) hexose.
16 ACS Paragon Plus Environment
Page 17 of 38
Journal of Agricultural and Food Chemistry
334
Alternatively, trisaccharide (j) appeared to fragment via the neutral loss of the reducing end Glc
335
moiety, suggesting the presence of a sulfate on the internal Gal residue. Not detected by MS
336
were any compounds consistent with sulfo-lactosamine or phospho-lactose analogues.
337
Sulfated oligosaccharides have been abundantly observed across multiple tissue types, for
338
example in glycosaminoglycans such as heparin sulfate, and glycoproteins like the epithelial
339
cell-expressed sulfomucins43 (reviewed in reference
340
adhesion ligands like GlyCAM-1.45,46 Intriguingly, the latter two classes of glycoproteins have
341
each been demonstrated to bear both 3’- and 6’-sulfo-galactose residues, analogous to those
342
borne by BMOs like 3’- and the putative 6’-sulfo-lactose (Figure 5). GlyCAM-1 has also been
343
demonstrated to be decorated with glycans terminating in 6’-sulfo-sialyl-LewisX epitopes
344
wherein the same non-reducing galactose residue is elaborated with both Neu5Ac (at the 3’ OH)
345
and sulfate (at the 6’ OH) moieties;45 indeed, the only difference between this epitope and the
346
sulfated/sialylated BMO described herein (Figure 5, peak f) would be the substitution of a
347
reducing GlcNAc for the Glc moiety of the BMO. The existence of (at least) eight sulfated
348
BMOs, comprising a major fraction of the total pool, is likely significant to human health for
349
several reasons. First, in humans, the highly glycosylated mucus glycoproteins (i.e., sulfomucins)
350
protecting the gastrointestinal (GI) tract exhibit a gradient of increasing oligosaccharide
351
sulfation, with the highest levels observed in the large intestine and colon, locations in which
352
BMOs/HMOs are thought to exert their primary influence.1–4 Sulfatase-catalyzed uncapping of
353
sulfated GI mucins is considered to be the rate-limiting step in the subsequent degradation of
354
mucus-linked oligosaccharides by GI microbes,44 and such a loss of sulfation on colonic mucins
355
has been correlated with inflammatory GI disorders such as ulcerative colitis.36,47 Thus, the
356
exogenous supply of sulfated BMOs may limit the degradation of host tissues during periods of
44
) and endothelially-expressed cell-
17 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 18 of 38
357
GI microbial dysbiosis. Interestingly, Bifidobacterium breve,48 and B. bifidum,49 and Bacteroides
358
thetaiotaomicron,50 all prominent human GI microbes, have each been recently observed to
359
possess sulfatase encoding genes permitting them to metabolize sulfated monosaccharides.
360
Therefore, it is very probable that in addition to serving as pathogen decoys, sulfated BMOs
361
exert specific prebiotic functions in the GI tract. Finally, sulfate-containing glycans are known to
362
be ligands for lectins that play key roles in modulating the immune system. For example, 6’-
363
sulfo-lactose has been demonstrated to be an antagonist of L-selectin46 and thus may directly
364
interfere with leukocyte recruitment to the GI tract. Alternatively, 3’-sulfo-lactose has shown to
365
be recognized and bound by galectin-4 which is predominantly expressed in the GI tract where it,
366
paradoxically, may both exacerbate and/or reduce mucosal inflammation.51 Similarly, 3’-sulfo-
367
lactose has also been shown to bind to P-selectin when presented in the context of the
368
glycolipids,52 and thus may also function as a P-selectin antagonist
369
Interrelationship between BMO concentrations. Spearman rank correlations, which have
370
previously been employed to deduce the directionality and strength of the biosynthetic
371
relationships between HMOs,22,53 were likewise calculated here for BMOs (Figure 6). Recently,
372
a similar correlative analysis on the seasonal influences on BMO concentrations was reported by
373
Liu and co-workers,21 albeit with several notable differences from the approach used here.
374
Specifically, Liu et al. sampled milk from 15 cows (at eight time points over eight months),
375
following the relative abundance of 14 BMOs; the correlation analysis here was performed using
376
six cows (at three time points over one week) but it included 33 BMOs, nine which were
377
identified with standards (Table S3). Several clear correlations between the Neu5Ac and Neu5Gc
378
levels (Figure 1D) and the relative abundances of specific BMOs validated the use of this
379
correlative approach. For example, Neu5Ac and Neu5Gc levels were positively correlated (ρ =
18 ACS Paragon Plus Environment
Page 19 of 38
Journal of Agricultural and Food Chemistry
380
0.664) with each other, consistent with Neu5Gc being directly biosynthesized from Neu5Ac
381
(Figure S2). Similarly, the major neuraminidase-sensitive BMOs (12 and 13) also were
382
positively correlated with total Neu5Ac/Neu5Gc levels (ρ >0.65 in all four cases); based on this
383
evidence it is proposed that BMO12 corresponds to the Neu5Gc-analogue of 3’SL, this peak
384
being clearly lactose-derived but not co-migrating with any Neu5Ac-containing standard (Figure
385
4B and 4C). Interestingly, BMO11 (6’SL) was neither correlated with total Neu5Ac levels nor
386
3’SL as might be predicted if there was competition between different Neu5Ac transferases (e.g.,
387
α2,3 vs. α2,6) for a limited pool of the donor monosaccharide cytidine-5’-monophosphate
388
(CMP)-Neu5Ac. Neuraminidase-sensitive BMO17 exhibited a similar lack of correlation
389
between relative abundance and total Neu5Ac/Neu5Gc levels. These data suggest that
390
Neu5Ac/Neu5Gc donors were not supply-limited, at least during the first week of lactation, and
391
that any correlations detected between the Neu5Ac/Neu5Gc-containing BMOs and others, such
392
as the sulfated species, were likely due to changes in the expression levels of other glycosyl- or
393
sulfo-transferases, or enzymes involved in the activities of these. Consistent with this, is the
394
significant negative correlation (ρ = -0.618) that was detected between 3’-sulfo-lactose (BMO3)
395
and 3’-SL (BMO13), a trend that is in agreement with the significant reversal in the relative
396
abundances of these BMOs over the first week of lactation (Figure 3). The correlation analysis
397
revealed that 3’-sulfo-lactose was also significantly correlated with BMO5 (ρ = 0.833), BMO15
398
(ρ = 0.717) and BMO20 (ρ = 0.558), two of which are known to be sulfated: BMO5, based on its
399
electrophoretic mobility, and BMO20 based on the observation that β-galactosidase specifically
400
converted this oligosaccharide to BMO8, a sulfated compound as inferred from its mobility
401
(Table S3 and Figure S8). As observed for 3’-sulfo-lactose, putative sulfated BMOs 5, 15 and 20
402
were all significantly negatively correlated with both total Neu5Ac/Neu5Gc levels and 3’SL, the
19 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 20 of 38
403
most abundant Neu5Ac-containing BMO. Based on these correlations, it is hypothesized that
404
BMOs 20 and 15, corresponding to sulfated tetra- and trisaccharides, respectively, are related in
405
that they share a common 3’-sulfo-lactose reducing end. BMO5 is putatively assigned as the
406
sulfate- and Neu5Ac-containing compound identified by MS (Figure 5, peak f), an assignment
407
consistent with its sensitivity to weak acid hydrolysis (Figure 4D) to yield BMO6 which is
408
hypothesized to be 6’-sulfo-lactose. BMO4 is hypothesized to be the disulfated trisaccharide
409
based on its relative electrophoretic mobility and its significant negative correlation to BMOs 8
410
and 10 (ρ = -0.651 and -0.680, respectively) which are likely two of three mono-sulfated
411
trisaccharides, both of which can serve as direct precursors to BMO4. Several experiments are
412
currently in progress in order to unambiguously correlate the sulfated BMOs detected by MS
413
with their corresponding CE peaks.
414
In conclusion, the CE analyses described herein permitted the rapid characterization of
415
the relative abundances of 33 BMOs in milk collected from exclusively grass-fed or grain/corn-
416
fed cows at matched time points during the first week of lactation. Few significant differences
417
were observed between individual animals in these groups, although differences were observed
418
in pooled samples of mature milk. The unexpectedly high electrophoretic mobilities of ten
419
BMOs led us to propose the existence of sulfated species, ten of which were detected by MS and
420
tandem MS. Correlation analyses revealed that the biosynthesis of several sulfated BMOs was in
421
direct competition with that of Neu5Ac/Neu5Gc-containing analogues, 3’-sialyllactose in
422
particular. Although these sulfated BMOs have not been previously described, they are likely to
423
be important to human health given their high relative abundance in mature, processed milk and
424
given that they are close structural analogues to the highly sulfated oligosaccharides lining the
425
human GI tract.
20 ACS Paragon Plus Environment
Page 21 of 38
Journal of Agricultural and Food Chemistry
426
Acknowledgments. The authors thank Mike Broersma, Ray Vanderhorst, and Garry and Kathy
427
Wikkerink for graciously collecting and donating the milk samples used in this study.
428
Funding Sources. This research was supported by funding from the Natural Science and
429
Engineering Research Council of Canada (Discovery Grant, 2016-03929). Infrastructure was
430
obtained with the support of the Canada Foundation for Innovation (project number 35246) and
431
the British Columbia Knowledge Development Fund. SDV was supported by scholarships from
432
the British Columbia Proteomics Network and a UBC undergraduate research award; AMG was
433
supported by an NSERC undergraduate student research award (USRA).
434
Supporting Information. The Supporting Information is available free of charge on the ACS
435
publication website. Included are eight figures and four tables.
436
References
437
(1)
Smilowitz, J. T.; Lebrilla, C. B.; Mills, D. A.; German, J. B.; Freeman, S. L. Breast milk
438
oligosaccharides: Structure-function relationships in the neonate. Annu. Rev. Nutr. 2014,
439
34, 143–169.
440
(2)
glycobiome on the neonate gut microbiota. Annu. Rev. Anim. Biosci. 2015, 3, 419–445.
441 442
Pacheco, A. R.; Barile, D.; Underwood, M. A.; Mills, D. A. The impact of the milk
(3)
Sela, D. A.; Li, Y.; Lerno, L.; Wu, S.; Marcobal, A. M.; Bruce German, J.; Chen, X.;
443
Lebrilla, C. B.; Mills, D. A. An Infant-associated bacterial commensal utilizes breast milk
444
sialyloligosaccharides. J. Biol. Chem. 2011, 286, 11909–11918.
445
(4)
Matsuki, T.; Yahagi, K.; Mori, H.; Matsumoto, H.; Hara, T.; Tajima, S.; Ogawa, E.;
446
Kodama, H.; Yamamoto, K.; Yamada, T.; et al. A key genetic factor for fucosyllactose
447
utilization affects infant gut microbiota development. Nat. Commun. 2016, 7, 11939.
448
(5)
Noll, A. J.; Gourdine, J. P.; Yu, Y.; Lasanajak, Y.; Smith, D. F.; Cummings, R. D.
21 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Galectins are human milk glycan receptors. Glycobiology 2016, 26, 655–669.
449 450
(6)
Noll, A. J.; Yu, Y.; Lasanajak, Y.; Duska-McEwen, G.; Buck, R.; Smith, D. F.;
451
Cummings, R. D. Human DC-SIGN binds specific human milk glycans. Biochem. J.
452
2016, 473, 1343–1353.
453
Page 22 of 38
(7)
Shams-Ud-Doha, K.; Kitova, E. N.; Kitov, P. I.; St-Pierre, Y.; Klassen, J. S. Human milk
454
oligosaccharide specificities of human galectins. Comparison of electrospray ionization
455
mass spectrometry and glycan microarray screening results. Anal. Chem. 2017, 89, 4914–
456
4921.
457
(8)
Kurakevich, E.; Hennet, T.; Hausmann, M.; Rogler, G.; Borsig, L. Milk Oligosaccharide
458
Sialyl(α2,3)-lactose activates Intestinal CD11c+ Cells through TLR4. Proc. Natl. Acad.
459
Sci. 2013, 110, 17444–17449.
460
(9)
Goehring, K. C.; Kennedy, A. D.; Prieto, P. A.; Buck, R. H. Direct Evidence for the
461
presence of human milk oligosaccharides in the circulation of breastfed infants. PLoS One
462
2014, 9, e101692.
463
(10)
Victora, C. G.; Bahl, R.; Barros, A. J. D.; França, G. V. A.; Horton, S.; Krasevec, J.;
464
Murch, S.; Sankar, M. J.; Walker, N.; Rollins, N. C.; et al. Breastfeeding in the 21st
465
century: Epidemiology, mechanisms, and lifelong effect. The Lancet. 2016, 387, 475–490.
466
(11)
Albrecht, S.; Lane, J. A.; Mariño, K.; Al Busadah, K. A.; Carrington, S. D.; Hickey, R.
467
M.; Rudd, P. M. A comparative study of free oligosaccharides in the milk of domestic
468
animals. Br. J. Nutr. 2014, 111, 1313–1328.
469
(12)
Aldredge, D. L.; Geronimo, M. R.; Hua, S.; Nwosu, C. C.; Lebrilla, C. B.; Barile, D.
470
Annotation and structural elucidation of bovine milk oligosaccharides and determination
471
of novel fucosylated structures. Glycobiology 2013, 23, 664–676.
22 ACS Paragon Plus Environment
Page 23 of 38
472
Journal of Agricultural and Food Chemistry
(13)
Charbonneau, M. R.; O’Donnell, D.; Blanton, L. V.; Totten, S. M.; Davis, J. C. C.;
473
Barratt, M. J.; Cheng, J.; Guruge, J.; Talcott, M.; Bain, J. R.; et al. Sialylated milk
474
oligosaccharides promote microbiota-dependent growth in models of infant
475
undernutrition. Cell 2016, 164, 859–871.
476
(14)
Lee, H.; Cuthbertson, D. J.; Otter, D. E.; Barile, D. Rapid screening of bovine milk
477
oligosaccharides in a whey permeate product and domestic animal milks by accurate mass
478
database and tandem mass spectral library. J. Agric. Food Chem. 2016, 64, 6364–6374.
479
(15)
Thurl, S.; Munzert, M.; Henker, J.; Boehm, G.; Mller-Werner, B.; Jelinek, J.; Stahl, B.
480
Variation of human milk oligosaccharides in relation to milk groups and lactational
481
periods. Br. J. Nutr. 2010, 104 (9), 1261–1271.
482
(16)
Davis, J. C. C.; Lewis, Z. T.; Krishnan, S.; Bernstein, R. M.; Moore, S. E.; Prentice, A.
483
M.; Mills, D. A.; Lebrilla, C. B.; Zivkovic, A. M. Growth and morbidity of Gambian
484
infants are influenced by maternal milk oligosaccharides and infant gut microbiota. Sci.
485
Rep. 2017, 7, 40466.
486
(17)
Thurl, S.; Munzert, M.; Henker, J.; Boehm, G.; Müller-Werner, B.; Jelinek, J.; Stahl, B.;
487
Muller-Werner, B.; Jelinek, J.; Stahl, B.; et al. Variation of human milk oligosaccharides
488
in relation to milk groups and lactational periods. Br. J. Nutr. 2010, 104, 1261–1271.
489
(18)
Sundekilde, U. K.; Barile, D.; Meyrand, M.; Poulsen, N. A.; Larsen, L. B.; Lebrilla, C. B.;
490
German, J. B.; Bertram, H. C. Natural variability in bovine milk oligosaccharides from
491
Danish Jersey and Holstein-Friesian breeds. J. Agric. Food Chem. 2012, 60, 6188–6196.
492
(19)
Tao, N.; DePeters, E. J.; German, J. B.; Grimm, R.; Lebrilla, C. B. Variations in bovine
493
milk oligosaccharides during early and middle lactation stages analyzed by high-
494
performance liquid chromatography-chip/mass spectrometry. J. Dairy Sci. 2009, 92,
23 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
2991–3001.
495 496
(20)
Tao, N.; DePeters, E. J.; Freeman, S.; German, J. B.; Grimm, R.; Lebrilla, C. B. Bovine milk glycome. J. Dairy Sci. 2008, 91, 3768–3778.
497 498
Page 24 of 38
(21)
Liu, Z.; Auldist, M.; Wright, M.; Cocks, B.; Rochfort, S. Bovine milk oligosaccharide
499
contents show remarkable seasonal variation and intercow variation. J. Agric. Food Chem.
500
2017, 65, 1307–1313.
501
(22)
McGuire, M. K.; Meehan, C. L.; McGuire, M. A.; Williams, J. E.; Foster, J.; Sellen, D.
502
W.; Kamau-Mbuthia, E. W.; Kamundia, E. W.; Mbugua, S.; Moore, S. E.; et al. What’s
503
normal? Oligosaccharide concentrations and profiles in milk produced by healthy women
504
vary geographically. Am. J. Clin. Nutr. 2017, 105, 1086–1100.
505
(23)
Schwendel, B. H.; Wester, T. J.; Morel, P. C. H.; Fong, B.; Tavendale, M. H.; Deadman,
506
C.; Shadbolt, N. M.; Otter, D. E. Pasture feeding conventional cows removes differences
507
between organic and conventionally produced milk. Food Chem. 2017, 229, 805–813.
508
(24)
Hong, Q.; Ruhaak, L. R.; Totten, S. M.; Smilowitz, J. T.; German, J. B.; Lebrilla, C. B.
509
Label-free absolute quantitation of oligosaccharides using multiple reaction monitoring.
510
Anal. Chem. 2014, 86, 2640–2647.
511
(25)
Evangelista, R. A.; Liu, M. S.; Chen, F. T. A. Characterization of 9-aminopyrene-1,4,6-
512
trisulfonate derivatized sugars by capillary electrophoresis with laser-induced fluorescence
513
detection. Anal. Chem. 1995, 67, 2239–2245.
514
(26)
Danyluk, H. J.; Shum, L. K.; Zandberg, W. F. A rapid procedure for the purification of 8-
515
aminopyrene trisulfonate (APTS)-labelled glycans for capillary electrophoresis (CE)-
516
based enzyme assays. Methods Mol. Biol. 2017, 1588, 223–236.
517
(27)
Wu, S.; Tao, N.; German, J. B.; Grimm, R.; Lebrilla, C. B. Development of an annotated
24 ACS Paragon Plus Environment
Page 25 of 38
Journal of Agricultural and Food Chemistry
library of neutral human milk oligosaccharides. J. Proteome Res. 2010, 9, 4138–4151.
518 519
(28)
chromatography. Anal. Biochem. 1973, 56 (2), 440–442.
520 521
Mopper, K.; Melvin Gindler, E. A new noncorrosive dye reagent for automatic sugar
(29)
Chen, F. T. A.; Dobashi, T. S.; Evangelista, R. A. Quantitative analysis of sugar
522
constituents of glycoproteins by capillary electrophoresis. Glycobiology 1998, 8, 1045–
523
1052.
524
(30)
Wylie, A.; Zandberg, W. F. Quantitation of sialic acids in infant formulas by liquid
525
chromatography-mass spectrometry: An assessment of different protein sources and
526
discovery of new analogues. J. Agric. Food Chem. 2018. (accepted)
527
(31)
Albrecht, S.; Schols, H. A.; Van Den Heuvel, E. G. H. M.; Voragen, A. G. J.; Gruppen, H.
528
CE-LIF-MSn profiling of oligosaccharides in human milk and feces of breast-fed babies.
529
Electrophoresis 2010, 31, 1264–1273.
530
(32)
Difilippo, E.; Pan, F.; Logtenberg, M.; Willems, R.; Braber, S.; Fink-Gremmels, J.;
531
Schols, H. A.; Gruppen, H. Milk oligosaccharide variation in sow milk and milk
532
oligosaccharide fermentation in piglet intestine. J. Agric. Food Chem. 2016, 64, 2087–
533
2093.
534
(33)
Difilippo, E.; Willems, H. A. M.; Vendrig, J. C.; Fink-Gremmels, J.; Gruppen, H.; Schols,
535
H. A. Comparison of milk oligosaccharides pattern in colostrum of different horse breeds.
536
J. Agric. Food Chem. 2015, 63, 4805–4814.
537
(34)
Corfield, A. P.; Wagner, S. A.; O’Donnell, L. J.; Durdey, P.; Mountford, R. A.; Clamp, J.
538
R. The roles of enteric bacterial sialidase, sialate O-acetyl esterase and glycosulfatase in
539
the degradation of human colonic mucin. Glycoconj. J. 1993, 10, 72–81.
540
(35)
Hunter, C. D.; Khanna, H.; Richards, M. R.; Derestani, R. R.; Zou, C.; Klassen, J. S.;
25 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
541
Cairo, C. Human neuraminidase isoenzymes show variable activities for 9-O-Acetyl-
542
sialoside substrates. ACS Chem. Biol. 2018, 13, 922-932.
543
(36)
Corfield, A. P.; Wagner, S. A.; Clamp, J. R.; Kriaris, M. S.; Hoskins, L. C. Mucin
544
degradation in the human colon: Production of sialidase, sialate O-acetylesterase, N-
545
acetylneuraminate lyase, arylesterase, and glycosulfatase activities by strains of fecal
546
bacteria. Infect. Immun. 1992, 60, 3971–3978.
547
(37)
Page 26 of 38
Samraj, A. N.; Pearce, O. M. T.; Läubli, H.; Crittenden, A. N.; Bergfeld, A. K.; Banda, K.;
548
Gregg, C. J.; Bingman, A. E.; Secrest, P.; Diaz, S. L.; et al. A red meat-derived glycan
549
promotes inflammation and cancer progression. Proc. Natl. Acad. Sci. U. S. A. 2015, 112,
550
542–547.
551
(38)
Varki, A.; Diaz, S. The release and purification of sialic acids from glycoconjugates:
552
Methods to minimize the loss and migration of O-Acetyl Groups. Anal. Biochem. 1984,
553
137, 236–247.
554
(39)
colostrum. Br. J. Nutr. 2000, 84, S69-S74.
555 556
Gopal, P. K.; Gill, H. S. Oligosaccharides and glycoconjugates in bovine milk and
(40)
Parkkinen, J.; Finne, J. Occurrence of N-acetylglucosamine 6-phosphate in complex
557
carbohydrates. Characterization of a phosphorylated sialyl oligosaccharide from bovine
558
colostrum. J. Biol. Chem. 1985, 260, 10971–10975.
559
(41)
unusual lactose sulfate in dog milk. Carbohydr. Res. 1999, 318, 123–128.
560 561
Bubb, W. A.; Urashima, T.; Kohso, K.; Nakamura, T.; Arai, I.; Saito, T. Occurrence of an
(42)
Guérardel, Y.; Morelle, W.; Plancke, Y.; Lemoine, J.; Strecker, G. Structural analysis of
562
three sulfated oligosaccharides isolated from human milk. Carbohydr. Res. 1999, 320,
563
230–238.
26 ACS Paragon Plus Environment
Page 27 of 38
564
Journal of Agricultural and Food Chemistry
(43)
Thomsson, K. A.; Bäckström, M.; Holmén Larsson, J. M.; Hansson, G. C.; Karlsson, H.
565
Enhanced detection of sialylated and sulfated glycans with negative ion mode nanoliquid
566
chromatography/mass spectrometry at high pH. Anal. Chem. 2010, 82, 1470–1477.
567
(44)
the human body. Biol. Chem. 1998, 379, 1–18.
568 569
Nieuw Amerongen, A. V; Bolscher, J. G.; Bloemena, E.; Veerman, E. C. Sulfomucins in
(45)
Hemmerich, S.; Bertozzi, C. R.; Rosen, S. D.; Leffler, H. Identification of the sulfated
570
monosaccharides of GlyCAM-1, an endothelial-derived ligand for L-selectin.
571
Biochemistry 1994, 33, 4820–4829.
572
(46)
Bruehl, R. E.; Bertozzi, C. R.; Rosen, S. D. Minimal sulfated carbohydrates for
573
recognition by L-selectin and the MECA-79 antibody. J. Biol. Chem. 2000, 275, 32642–
574
32648.
575
(47)
Raouf, A. H.; Tsai, H. H.; Parker, N.; Hoffman, J.; Walker, R. J.; Rhodes, J. M. sulphation
576
of colonic and rectal mucin in inflammatory bowel disease: Reduced sulphation of rectal
577
mucus in ulcerative colitis. Clin. Sci. 1992, 83, 623–626.
578
(48)
Egan, M.; Jiang, H.; Motherway, M. O. C.; Oscarson, S.; Van Sinderen, D.
579
Glycosulfatase-encoding gene cluster in Bifidobacterium Breve UCC2003. Appl. Environ.
580
Microbiol. 2016, 82, 6611–6623.
581
(49)
Katoh, T.; Maeshibu, T.; Kikkawa, K. ichi; Gotoh, A.; Tomabechi, Y.; Nakamura, M.;
582
Liao, W. H.; Yamaguchi, M.; Ashida, H.; Yamamoto, K.; et al. Identification and
583
Characterization of a sulfoglycosidase from Bifidobacterium bifidum implicated in mucin
584
glycan utilization. Biosci. Biotechnol. Biochem. 2017, 81, 2018–2027.
585 586
(50)
Cartmell, A.; Lowe, E. C.; Baslé, A.; Firbank, S. J.; Ndeh, D. A.; Murray, H.; Terrapon, N.; Lombard, V.; Henrissat, B.; Turnbull, J. E.; et al. How members of the human gut
27 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
587
microbiota overcome the sulfation problem posed by glycosaminoglycans. Proc. Natl.
588
Acad. Sci. 2017, 114, 7037–7042.
589
(51)
Page 28 of 38
Bum-Erdene, K.; Leffler, H.; Nilsson, U. J.; Blanchard, H. Structural characterisation of
590
human galectin-4 N-terminal carbohydrate recognition domain in complex with glycerol,
591
lactose, 3′-sulfo-lactose, and 2′-fucosyllactose. Sci. Rep. 2016, 6, 20289.
592
(52)
binding to sulfatides on tumor cells. Glycobiology 2007, 17, 185–196.
593 594
Garcia, J.; Callewaert, N.; Borsig, L. P-selectin mediates metastatic progression through
(53)
Smilowitz, J. T.; O’Sullivan, A.; Barile, D.; German, J. B.; Lonnerdal, B.; Slupsky, C. M.
595
The human milk metabolome reveals diverse oligosaccharide profiles. J. Nutr. 2013, 143,
596
1709–1718.
597
Figure captions
598
Figure 1. Quantification of total BMO and monosaccharide levels in milk from corn/grain (c/g)-
599
or grass (gr)-fed cows. (A) Total reducing BMO levels. (B) Representative CE electropherogram
600
containing standards and a BMO hydrolysate used for (C) determining the levels of neutral
601
monosaccharides relative to Glc. (D) Absolute amounts of the acidic monosaccharides, Neu5Ac
602
and Neu5Gc. n = 3 biological replicates for the three time points in each diet group; pooled tank
603
samples were extracted and analyzed in triplicate. All data are reported ± 1 standard error of the
604
mean (SEM).
605
Figure 2. Relative BMO levels established by CE with fluorescence detection. (A)
606
Representative CE electropherogram of milk (colostrum) from a grass (only)-fed cow. Peak
607
areas for the first 34 BMOS, in order of migration, were determined for all samples. (B) Relative
608
levels of selected BMOs. Open circles represent each biological replicate; closed squares are
28 ACS Paragon Plus Environment
Page 29 of 38
Journal of Agricultural and Food Chemistry
609
mean values with the bars corresponding to the SEM. * indicates a significant difference
610
(Student’s t test; p < 0.05) between the two groups that time point.
611
Figure 3. Relative changes in major BMOs as a function of lactation time. BMOs are numbered
612
as shown if Figure 2. Each box-whisker plot contains six biological replicates, three animals
613
from each farm. Center lines show medians, box limits indicate the 25 and 75 percentiles,
614
whiskers extend to 1.5 times the interquartile range from the 25 and 75 percentiles, points
615
beyond these ranges were assigned as outliers. Statistical differences between groups were
616
assessed using a Mann-Whitney U test (* p < 0.05; ** p < 0.01).
617
Figure 4. Characterization of highly anionic BMOs based on their sensitivity or resistance to
618
weak acid- or enzyme-catalyzed hydrolysis. (A) Both mild acid (acetic acid = AcOH) and
619
neuraminidase were predicted to selectively hydrolyze Neu5Ac/Neu5Gc-bearing, APTS-labelled
620
BMOs which was proven (B) using commercially-available standards. (C) Neuraminidase failed
621
to hydrolyze BMOs of higher mobility than 6’SL (BMO11) but BMO11-13, 17, 21 and 25 were
622
completely digested (D) All neuraminidase-sensitive BMOs were also acid-labile under
623
conditions that failed to obviously hydrolyze BMO1-4, and 6, while relative levels of BMO3 and
624
6 increased. (E) Alkaline phosphatase obviously hydrolyzed BMO2. Note that different milk
625
samples were used to prepare Figure 2, and 4C- 4E; 2,4C and 4E were colostrum (day 0) samples
626
while 4D was a pooled sample of mature milk.
627
Figure 5. Sulfated and phosphorylated BMOs detected by MS. Extracted ion chromatograms of
628
ten putative BMOs (Table 1). All BMOs reported were from the same HPLC-MS injection. For
629
comparative purposes the intensity of peak g (3’-sulfo-lactose) and peaks h – j have been scaled
630
down 500- and 100-fold, respectively; peak c has been scaled up 10-fold.
th
th
th
th
29 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 30 of 38
631
Figure 6. Spearman rank correlation heatmap depicting the relationships between the 33 BMOs
632
detected in all individual milk samples. BMOs are numbed as depicted in Figure 2; structures are
633
indicated only where BMOs have been identified based on their co-migration with commercial
634
standards (see Figure S1 for full names and a list of all standards tested). Color indicates
635
directionality (red = positive; blue = negative) and intensity indicates the magnitude of each
636
association.
30 ACS Paragon Plus Environment
Page 31 of 38
Journal of Agricultural and Food Chemistry
Tables Table 1. Retention time (RT) and MS parameters for the HPLC-MS analysis of sulphated and phosphorylated bovine milk oligosaccharides. Peaka
RT (min)
Formula
Compositionb
Exact mass (Da)
δ (ppm)
a
5.0
C14H28NO14P-
PO3 + Hex + HexNAc
464.1160
+ 1.94
b
5.0
C24H42O24S
SO3+ 4Hex
746.1805
- 2.41
c
8.2
C18H34O22S2
2 SO3 + 3Hex
666.0983
-9.18
d
10.6
C12H24O14S
SO3 + 2Hex
424.0925
+ 8.92
e
12.2
C14H28NO14P- PO3 + Hex + HexNAc
464.1192
- 4.96
f
13.0
C23H41NO22S
SO3 + Neu5Ac+2Hex
715.1794
+ 6.57
g
15.0
C12H24O14S
3’-sulfo-lactose
424.0896
- 2.12
h
17.0
C18H34O19S
SO3 + 3Hex
586.1455
- 6.82
i
17.7
C18H34O19S
SO3 + 3Hex
586.1468
- 9.04
j
18.8
C18H34O19S
SO3 + 3Hex
586.1361
+ 9.21
k
22.0
C24H42O24S
SO3 + 4Hex
746.1800
- 1.74
a
Peaks are labelled as in Figure 5. bHex = hexose; HexNAc = N-acetylhexosamine (most likely GlcNAc); SO3 and PO3 = sulfate and phosphate, respectively.
31 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 32 of 38
Figure graphics
Figure 1.
32 ACS Paragon Plus Environment
Page 33 of 38
Journal of Agricultural and Food Chemistry
Figure 2.
33 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 34 of 38
Figure 3.
34 ACS Paragon Plus Environment
Page 35 of 38
Journal of Agricultural and Food Chemistry
Figure 4.
35 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 36 of 38
Figure 5
36 ACS Paragon Plus Environment
Page 37 of 38
Journal of Agricultural and Food Chemistry
Figure 6.
37 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 38 of 38
Graphic for table of contents
38 ACS Paragon Plus Environment