New Insights into Porcine Milk - ACS Publications - American

Dec 21, 2018 - New Insights into Porcine Milk N‑Glycome and the Potential. Relation ... KEYWORDS: breed, lactation time, N-glycome, porcine milk N-g...
0 downloads 0 Views 2MB Size
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