Subscriber access provided by University of Winnipeg Library
Article
New insights into porcine milk N-glycome and the potential relation with offspring gut microbiome Chunlong Mu, Zhipeng Cai, Gaorui Bian, Yamin Du, Shouqing Ma, Yong Su, Li Liu, Josef Volgmeir, Ruihua Huang, and Weiyun Zhu J. Proteome Res., Just Accepted Manuscript • DOI: 10.1021/acs.jproteome.8b00789 • Publication Date (Web): 21 Dec 2018 Downloaded from http://pubs.acs.org on December 23, 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 41 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Journal of Proteome Research 1
1
New insights into porcine milk N-glycome and the potential relation with
2
offspring gut microbiome
3 4
Chunlong Mu †‡§⊥ , Zhipeng Cai ‖⊥ , Gaorui Bian †‡§ , Yamin Du ‖ , Shouqing Ma †‡§ , Yong
5
Su†‡§, Li Liu‖, Josef Volgmeir‖, Ruihua Huang†, Weiyun Zhu*†‡§
6
† Laboratory
7
‡Jiangsu
8
for International Research on Animal Gut Nutrition, ‖Glycomics and Glycan Bioengineering
9
Research Center (GGBRC), College of Food Science and Technology, Nanjing Agricultural
of Gastrointestinal Microbiology, College of Animal Science and Technology,
Key Laboratory of Gastrointestinal Nutrition and Animal Health, §National Center
10
University, Nanjing 210095, China
11
⊥These
authors contributed equally to this work.
12
ACS Paragon Plus Environment
Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 2 of 41 2
14
ABSTRACT: N-glycans are an important source of milk oligosaccharides. In addition to
15
free oligosaccharides found in milk, N-glycans can also be utilized by gut microbes. A
16
potential for milk N-glycans to act as gut microbe regulators in suckling animals has
17
attracted considerable attention; however, sow milk N-glycans and their potential effects
18
upon the piglet’s gut microbes in vivo remain unknown. In the present study, we profiled
19
the milk N-glycans of Meishan and Yorkshire sows during lactation using UPLC and a mass
20
spectrometry-based glycome method, and explored the correlations between milk N-
21
glycans and offspring gut microbiota. Twenty-two N-glycan structures were identified in
22
sow milk, among which 36% (8 out of 22) were fucosylated, 41% (9 out of 22) were
23
sialylated, and 14% (3 out of 22) were high mannosylated. An N-glycan with a NeuGc
24
structure was identified in sow milk for the first time. No compositional differences between
25
the two breeds or between different lactation times were found in porcine milk N-linked
26
oligosaccharides (PNOs); however, the abundances of different structures within this class
27
did vary. The relative abundances of fucosylated PNO3 (GlcNAc4-Man3-Fuc) and sialylated
28
PNO18 (GlcNAc4-Man3-Gal2-NeuAc) increased during lactation, and Meishan sows
29
demonstrated a higher (P < 0.05) abundance of mannosylated PNO10 (GlcNAc2-Man6)
30
and sialylated PNO17 (GlcNAc5-Man3-Gal-NeuAc) than Yorkshire sows. Apparent
31
correlations between milk N-glycans and offspring gut microbial populations were found;
32
for example, mannosylated PNO21 (GlcNAc2-Man9) was positively correlated with OTU706
33
(Lactobacillus amylovorus) and OTU1380 (Bacteroides uniformis). Overall, our results
34
indicate that the sow milk N-glycome of Meishan and Yorkshire sows differs in N-glycome
35
characteristics, and that this is correlated to abundances of certain piglet gut microbes.
ACS Paragon Plus Environment
Page 3 of 41 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Journal of Proteome Research 3
36
These findings provide a reference for future elucidation of the involvement of gut microbes
37
in milk N-glycan metabolism, which is important to the health both of large domestic
38
animals and humans.
39
KEYWORDS: breed, lactation time, N-glycome, porcine milk N-glycan oligosaccharides
ACS Paragon Plus Environment
Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 4 of 41 4
41
■ INTRODUCTION
42
Milk oligosaccharides are essential nutrients for neonates, and also stimulate the growth
43
of beneficial microbes and participate in pathogen inhibition in the gut.1 Milk
44
oligosaccharides include both free and conjugated oligosaccharides, and free
45
oligosaccharides have already been shown to be utilized by intestinal microbes.1 In
46
addition to free oligosaccharides, N-glycans, one of the major conjugated sugars, are also
47
an important component of milk oligosaccharides and have demonstrated pathogen-
48
inhibiting activity.2 N-glycans are oligosaccharides attached to the asparagine residues of
49
a protein via N-acetylglucosamine linkages. The compositions and structures of free
50
oligosaccharides in sow milk have been well-elucidated3, 4; however, the N-glycan
51
component in sow milk has yet to be explored.
52
Breed differences affect milk composition in sows and thus nutrients for neonates. For
53
example, Meishan (a typical obese Chinese breed) sows have a higher concentration of
54
milk fats and calcium than Large White (Yorkshire, a lean type breed) sows.5, 6 Interestingly,
55
Meishan and Yorkshire pigs differ in genes encoding fucosyltransferases7 which mediate
56
the transfer of fucose onto N-glycan chains. In a survey of the genotype frequency of the
57
porcine fucosyltransferase 1 gene, at the 307 bp loci, Meishan pigs were found to have a
58
higher frequency of GG motifs, while Yorkshire pigs presented with both GG and AG
59
motifs.7 Glycoprotein profiles in the ileum mucosa of Yorkshire pigs showed that those with
60
the AG motif have a lower abundance of sialylated glycans than those with an AA motif.8
61
These results suggest potential important differences in glycan patterns among pig breeds.
62
However, whether breed affects milk N-glycan composition has yet to be determined.
ACS Paragon Plus Environment
Page 5 of 41 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Journal of Proteome Research 5
63
Recently, milk glycans have been shown to be an important substrate for gut
64
microbiota. Both glycoproteins and their glycan components can be utilized by specific
65
microbes. Bifidobacterium infantis use endo-β-N-acetylglucosaminidase to release
66
complex N-glycans from human milk glycoproteins.9 Bacteroides fragilis can also
67
deglycosylate glycoprotein by producing endo-β-N-acetylglucosaminidase which in turn
68
generates glycans that can be used for bacterial growth.10 Bacteroides thetaiotaomicron
69
can
70
glucosamine phosphorylase for degradation of mannosylated N-glycans.11 The human milk
71
N-glycome can promote the growth of Bifidobacterium longum subsp. infantis, which
72
utilizes N-glycans with a preference for small structures.12 These results suggest that
73
variations among preferred milk glycan substrates by disparate microbial species could
74
affect gut microbial composition. In our previous study, we cross-fostered newborn piglets
75
between genetically obese Meishan and lean Yorkshire sows, and found that milk
76
appeared to affect the fecal microbiota composition of the pre-weaned piglets.13 However,
77
little is known about the in vivo effects of sow milk N-glycan composition on gut microbial
78
colonization in offspring.
secrete
endo-β-N-acetylhexosaminidase
and
β-1,4-D-mannosyl-N-acetyl-D-
79
In the present study, we further analyzed the milk N-glycan compositions of Meishan
80
and Yorkshire sows during lactation. We also delineated the microbial population
81
succession in piglets to explore potential correlations between the piglet gut microbiota and
82
milk N-glycans. The interest of profiling sow milk and the pig microbiome is to provide a
83
rational scientific basis for future studies on the impact of glycans on gut microbiota in this
84
and other large domestic species. The pattern of early colonization and acquisition of gut
ACS Paragon Plus Environment
Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 6 of 41 6
85
microbiota of the piglet is similar to that in human infants during breastfeeding in many
86
aspects, such as increasing α diversity, gradually increasing Bacteroidetes spp.
87
abundance, and decreasing Proteobacteria abundance.13, 14 In the present study, we use
88
the piglet as a model to elucidate the relationship between milk N-glycans and offspring
89
gut microbiota composition. We hypothesize that milk N-glycans may be related to the
90
alteration of early-life gut microbiota colonization in offspring. By profiling the composition
91
of milk N-glycans and examining the in vivo correlation of these glycans to gut microbiota,
92
we intend to provide the basis for future research into the involvement of specific microbes
93
in milk N-glycan metabolism, as well as underscore its importance in animal and human
94
health. These data may also provide a reference for regulating milk nutrition and seeking
95
strategies for assisting infants subject to milk deficiency.
96
■ MATERIALS AND METHODS
97
Ethics
98
Animals were managed throughout the study in accordance with requirements for the
99
Experimental Animal Care and Use guidelines of Chinese Science and Technology
100
Committee, 1998.
101
Experimental design
102
All Meishan and Yorkshire pigs were raised under the same conditions on a commercial
103
farm in Jiangsu Province, China. A cross-fostering approach was adopted with a
104
randomized complete block design [2 treatments (breed) *2 block (nursing) with 10
105
replicates per group] as described previously.13 Briefly, after vaginal delivery, half of the
106
piglets in a litter of a Meishan sow were fostered onto a Yorkshire sow before the piglets
ACS Paragon Plus Environment
Page 7 of 41 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Journal of Proteome Research 7
107
suckled their mothers’ colostrum. The other half of each litter remained with the birth
108
mother. In total, 10 sows of each breed with litters of 10 or 12 piglets were used for the
109
cross-fostering operation in delivery-matched pairs. The gestational length of all sows was
110
113 –114 d. Diet composition was the same for Yorkshire and Meishan sows, as shown in
111
Table S1. Finally, four groups of piglets were generated as follows: Meishan piglets
112
fostered by their birth mother (Mm), Yorkshire piglets fostered by Meishan sows (My),
113
Meishan piglets fostered by Yorkshire sows (Ym), and Yorkshire piglets fostered by their
114
birth mother (Yy). From day 14 after birth, all suckling piglets were offered creep feed ad
115
libitum and had free access to water.
116
Milk samples from the sows of both breeds were collected when the piglets were 1, 3,
117
7, 14 and 21 days old (n = 6 for each sows at sample days). The sample replicates were
118
chosen in order to balance the cost of replicates and the number of replicates necessary
119
for statistical analysis. The fecal samples from offspring piglets of the chosen sows, from
120
which milk samples were collected, were used for microbiome analysis. The first fecal and
121
milk samples were collected within 12 h after birth, while other samples were collected
122
between 7:00 and 9:00 AM on the sampling day. A total of 60 samples were used for glycan
123
analysis.
124
Glycoprofiling of Porcine Milk N-glycan Oligosaccharides (PNOs)
125
N-glycans preparation. The stored frozen milk samples were completely thawed to room
126
temperature and protein was subsequently precipitated using chloroform and methanol.15
127
The N-glycans were obtained by N-glycanase F (PNGase F) digestion followed by solid-
ACS Paragon Plus Environment
Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 8 of 41 8
128
phase extraction with non-porous graphitized carbon cartridges to remove the buffer salts
129
prior to 2-AB labelling.
130
aminobenzamide (2-AB) prior to ultra performance liquid chromatography (UPLC) and
131
exoglycosidase assays.17
16
Then, the purified N-glycans were fluorescently labeled with 2-
132
Hydrophilic interaction ultra performance liquid chromatography (HILIC-UPLC) for
133
profiling of milk N-glycome. 2-AB-labeled N-glycan samples were profiled on an Acquity
134
UPLC BEH Glycan Column (2.1×150 mm, 1.7 µm particle size; Waters, Ireland) at a
135
column temperature of 60°C, using the Nexera UPLC system (Shimadzu Corporation,
136
Kyoto, Japan) consisting of a LC-30AT pump system and a RF-20Axs fluorescence
137
detector (excitation 330 nm, emission 420 nm). Solvent A and B were aqueous ammonium
138
formate buffer (50 mM, pH 4.5) and acetonitrile respectively. The separation was
139
performed using a linear gradient of 5–12% of solvent A from 0-6 min, solvent B was then
140
increased to 44.1% over 39 min followed by a further increase to 100% over 3 min and
141
held at 100% for 1 min. The solvent A (aqueous ammonium formate buffer, 50 mM, pH
142
4.5) was then decreased to 5% in 7 min and the column was equilibrated with the initial
143
conditions for 7 min. The flow rate was 0.5 ml/min from 0-44.5 min, reduced to 0.25 ml/min
144
from 44.5-55 min and increased back to 0.5 ml/min from 55-57 min. Different flow rates
145
wee used to reduce the backpressure of BEH Glycan Column. The method of
146
chromatographic separation is based on the manual of BEH Glycan Column (2.1×150 mm,
147
1.7 μm particle size). The optimal flow rate for the elution of N-glycans is 0.5mL/min. After
148
that, it is necessary to wash the column several minutes in 100% solvent A (aqueous
149
ammonium formate buffer, 50 mM, pH 4.5) before next sample. However, the backpressure
ACS Paragon Plus Environment
Page 9 of 41 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Journal of Proteome Research 9
150
will be very high in this condition. Therefore, we reduced the flow rate to 0.25mL/min during
151
the chromatographic separation. This procedure can benefit the column by avoiding the
152
high pressure and maximizing the lifetime of column and system. The column was
153
calibrated with a 2AB-labeled dextran standard (2-20 glucose units), and the absolute
154
retention time of each N-glycan peak was converted to relative retention units based on
155
the dextran-ladder-derived glucose units (GU value) to avoid experimental error. UPLC
156
fractions with selected peaks were manually collected for mass spectrometry assay.
157
Exoglycosidase assay. In order to elucidate the glycan structure, N-glycans were
158
subjected to exoglycosidase assays. All the enzymes employed in this study were obtained
159
from Prozyme (San Leandro, CA, USA) and used according to the manufacturer’s
160
instruction. After Jack bean α-mannosidase digestion, the sample was divided into five
161
parts. One part was directly analyzed by HILIC-UPLC, whereas the others were subjected
162
to further sequential glycosidase treatment including: sialidase A, bovine kidney α-
163
fucosidase, bovine testis β-galactosidase, alpha-galactosidase, and jack bean β-N-
164
acetylhexosaminidase. Samples were then analyzed by UPLC.
165
Mass Spectrometry (MS). 2AB-labelled N-glycans from UPLC were analyzed by
166
matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (Autoflex;
167
Bruker Daltonics, Bremen, Germany). The MS spectra were acquired in monoisotopic
168
resolution using a Bruker Autoflex Speed (equipped with a 2000 Hz Smartbeam™‐II laser)
169
instrument (Bruker Daltonics, Bremen, Germany). 6-aza-2-thiothymine was used as the
170
matrix (3 mg/ml in 70% aqueous acetonitrile solution). MS/MS was performed by laser-
171
induced dissociation. Spectra were processed with the manufacturer’s software (Bruker
ACS Paragon Plus Environment
Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 10 of 41 10
172
Flex Analysis software version 3.3.80) using the SNAP algorithm with a signal/noise
173
threshold of 6 for MS (unsmoothed) and 3 for MS/MS (four times smoothed). The MS and
174
MS/MS spectra were interpreted using GlycoWorkbench version 1.1 (database: CFG,
175
Carbbank and Glycosciences),18 with the attached labeling reagent 2AB being taken into
176
account. Permethylation of N-glycan was carried out according to the method of Harvey
177
and coworkers to protect the sialylated glycans from the desialylation during the MS
178
detection.19
179
Analysis of 16S rRNA pyrosequencing data
180
In our previous study,13 fecal samples were collected from piglets at day 1, 3, 7, and 14 of
181
age. The total genomic DNA was isolated from fecal samples by the bead-beating method
182
as described by Zoetendal and colleagues.20 Amplicon pyrosequencing was performed
183
from the A-end using a 454/Roche A sequencing primer kit on a Roche Genome
184
Sequencer GS-FLX Titanium platform. The processes of amplification, amplicon
185
purification, check of amplicon quantity and quality have been described previously.13 The
186
microbial diversity and composition at the phylum and genus level have been reported
187
previously.13 To link PNOs with fecal bacteria species, we further analyzed the operational
188
taxonomic units (OTUs) using the microbiota dataset.
189
The sequence filtering and chimeric trimming were conduced using QIIME software
190
package.21 Using the average neighbor algorithm with a cutoff of 97% similarity, sequences
191
were clustered to OTUs. Representative sequences from each OTU were taxonomically
192
classified with a 90% confidence level using Ribosomal Database Project classifier.22 The
193
heatmap of microbial composition was visualized using the Heatmap Illustrator version
ACS Paragon Plus Environment
Page 11 of 41 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Journal of Proteome Research 11
194
1.0.23
195
Statistical analysis
196
The data of porcine milk N-glycan oligosaccharides at different time points were analyzed
197
using two-way repeated analysis of variance procedures as implemented in MetaboAnalyst
198
3.0.24 Comparisons of N-glycans at each time point were conducted using Student’s t test
199
(SPSS version 20.0; SPSS Inc., Chicago, Illinois). P-values < 0.05 were regarded as
200
statistically significant. The relative abundances of bacterial taxa at the OTU level were
201
analyzed using the general lineal mode procedures as implemented in SPSS (SPSS
202
version 20.0; SPSS Inc., Chicago, Illinois). In case of multiple comparisons, P-values were
203
adjusted with Benjamini-Hochberg multiple testing correction, limiting the overall false
204
discovery rate to 5% (q < 0.05). 0.05< q