N-glycan profile as a tool in qualitative and quantitative analysis of

antigen epitope when enters into human body43. These structures are all the key factors. 248 responsible for the change of N-glycan UPLC profile when ...
0 downloads 0 Views 1MB Size
Subscriber access provided by Macquarie University

Omics Technologies Applied to Agriculture and Food

N-glycan profile as a tool in qualitative and quantitative analysis of meat adulteration Zihe Shi, Binru Yin, Yuquan Li, Guanghong Zhou, Chunbao Li, Xinglian Xu, Xin Luo, Xibin Zhang, Jun Qi, Josef Voglmeir, and Li Liu J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.9b03756 • Publication Date (Web): 29 Aug 2019 Downloaded from pubs.acs.org on August 30, 2019

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 31

Journal of Agricultural and Food Chemistry

Graphic for table of contents 142x103mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

N-glycan profile as a tool in qualitative and quantitative analysis of meat adulteration

Zihe Shia, Binru Yina, Yuquan Lia, Guanghong Zhoub, Chunbao Li b, Xinglian Xu b, Xin Luoc, Xibin Zhangc,d, Jun Qie, Josef Voglmeir*,a,b, Li Liu*,a,b

a

Glycomics and Glycan Bioengineering Research Center (GGBRC), College of Food

Science and Technology, Nanjing Agricultural University, Nanjing, Jiangsu 210095, China b

Jiangsu Colleborative Innovation Center of Meat Production, Processing and Quality

Control, College of Food Science and Technology, Nanjing Agricultural University, Nanjing, Jiangsu 210095, China c

Lab of Beef Processing and Quality Control, College of Food Science and Engineering,

Shandong Agricultural University, Taian, Shandong 271018, China d

New Hope Liuhe Co. Ltd., Beijing 100102, China

e

Anhui Engineering Laboratory for Agro-products Processing, Anhui Agricultural

University, Hefei, Anhui 230036, China *Correspondence

should be addressed to:

E-mail: [email protected]: Fax: +86 25 84399553 Tel: +86 25 84399511 or E-mail: [email protected], Fax: +86 25 84399553 Tel: +86 25 84399512

ACS Paragon Plus Environment

Page 2 of 31

Page 3 of 31

Journal of Agricultural and Food Chemistry

Keywords: adulteration; N-glycan profiles; principal component analysis; PLS regression 1

Abstract: Adulteration of meat and meat products causes a concerning threat for

2

consumers. It is necessary to develop novel robust and sensitive methods which can

3

authenticate the origin of meat by qualitative and quantitative means to compensate

4

the drawbacks of the existing methods. This study has shown that the protein N-

5

glycosylation profiles of different meats are species specific and thus can be used for

6

meat authentication. Based on N-glycan pattern, the investigated five meat species (beef,

7

chicken, pork, duck and mutton) can be distinguished by principal component analysis

8

(PCA), and partial least square (PLS) regression was performed to build a calibration

9

and validation model for prediction of the adulteration ratio. Using this method, beef

10

samples adulterated with the lower value duck meat could be detected down to the

11

addition ratio as low as 2.2%. The most distinguishing N-glycans from beef and duck

12

were elucidated for the detailed structures.

13 14 15 16 17

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 4 of 31

18 19

Introduction

20

Meat and meat products are highly nutritious components of the diet and are widely

21

appreciated for their flavor1. More costly meats such as beef and mutton are reported to

22

be fraudulently substituted by cheaper meats such as chicken, pork and duck2-4.

23

Therefore, authentication of meat is essential to safeguard consumer rights, religious

24

beliefs, immune specific dietary requirements and even to some extent to protect

25

wildlife5-8.

26

Historically, the identification of meat species is challenging due to the change of

27

external morphological features after meat pre-processing such as deboning, mincing,

28

chopping, emulsification and other pre-processing procedures9. Different products and

29

consumption patterns require different adulteration identification techniques. At

30

present, many approaches are established for meat authentification including

31

biochemical methods

32

spectroscopic methods

33

the last 25 years was achieved by using DNA-based methods such as the polymerase

34

chain reaction (PCR). DNA shows much better stability and is generally universally

35

applicable, as all tissue samples contain DNA24. However, one drawback of DNA-based

36

methods is that it is challenging to quantify the extent of meat adulteration, as even

37

minute amounts of contaminants (i.e. during the meat processing) will be detected1, 25-26.

38

The comprehensive study on complex carbohydrates in tissues and cells (so-called

39

glycomics) has rapidly developed in recent years. Currently, glycomic techniques are

40

primarily applied in biomedical research for the evaluation of aberrant glycosylation

10-12, 20-23.

immunological methods13-14, molecular methods

15-19and

The greatest progress made in species identification over

ACS Paragon Plus Environment

Page 5 of 31

Journal of Agricultural and Food Chemistry

41

pattern as biomarkers27-28. The similar methodologies have also been proven as

42

powerful tools in food safety and food quality control29-30. For example, Shim et al.

43

described a method to characterize edible bird nest (which is made from solidified

44

swiftlet saliva) adulteration based on their glycan composition31. Another recent

45

example by Wang et al. showed that protein glycosylation of edible Gingko seeds varies

46

based on their geographic origin32. However, no study was yet conducted using the

47

analysis of protein glycosylation for the determination of meat species.

48

In this study, five meat species (beef, chicken, pork, duck and mutton) were

49

investigated based on their N-glycan profiles, which were analyzed using hydrophilic

50

interaction ultra-performance liquid chromatography (HILIC-UPLC). The obtained

51

UPLC data is evaluated with PCA (principal component analysis) for discrimination of

52

samples and PLS (partial least squares) regression analysis for determination of

53

adulteration ratio. This method was then applied to discriminate duck meat, which is

54

often used for the adulteration of beef due to its intense flavor and dark meat color33-34,

55

from beef samples. Multivariate data analysis revealed the quantity of added duck meat

56

based on the glycosylation pattern present in the analyzed meat samples.

57 58

Materials and Methods

59

Materials. HPLC-grade acetonitrile (ACN) was purchased from Merck (Darmstadt,

60

Germany). Supelclean ENVI-Carb solid-phase extraction (SPE) columns, trifluoroacetic

61

acid (TFA), and dextran oligomers (molecular range 300 - 3000 Da) were purchased

62

from Sigma (St. Louis, USA). 2-aminobenzamide (2-AB) was obtained from J&K

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

63

Chemicals (Beijing, China). PNGase F (50 mU) was provided by Qlyco Ltd. (Nanjing,

64

China).

65 66

Sample preparation. Fifteen meat samples for each species were purchased from

67

local supermarkets or slaughterhouses. After removal of connective tissue and visible

68

fat, approximately 200 g of each meat sample was minced to homogeneity. To mimic the

69

actual meat adulteration, duck meat in various concentrations (5%, 10%, 20%, 30%,

70

40% and 50% w/w) were added into beef for qualitative and quantitative binary

71

adulteration analysis, and chicken meat, pork and chicken + pork (in 10% and 20%

72

concentrations, respectively) were added into beef for trinary species adulteration

73

detection. All different types of samples were stored at -18 °C prior to sample analysis.

74 75

Preparation of meat N-glycans. Approximately 100 mg of minced meat was first

76

treated with 660 µL of chloroform and 330 µL of methanol, and the sample was then

77

centrifuged at 13,000 g for 5 min after vigorous mixing. The interface was recovered by

78

the careful removal the top (methanol) and bottom (chloroform) phase, and re-

79

suspended in 100 µL of 40% trichloroacetic acid (TCA, v/v). After centrifugation at

80

13,000 g for 10 min at 4°C, the supernatant was removed. The obtained (glyco-) protein

81

pellet was then washed with 1.8 mL of deionized water and centrifuged at 13,000 g for

82

5 min at 4°C. After removal of the supernatant, the pellet was then re-solubilized in 100

83

µL of 6 M urea. After the addition of 46 µL of sodium-phosphate solution (500 mM, pH

84

7.5), 23 µL of sodium dodecyl sulfate solution (2% w/v SDS in 1 M β-mercaptoethanol)

85

and 56 µL of distilled water, the mixtures were heated and boiled at 100 °C for 5 min.

ACS Paragon Plus Environment

Page 6 of 31

Page 7 of 31

Journal of Agricultural and Food Chemistry

86

The solutions were cooled down and 38 µL of Triton-X100 solution (10% w/v) together

87

with 200 µL of PNGase F enzyme solution (50 mU) were added, prior to the incubation

88

at 37°C for 16 h. The reaction mixture was centrifuged at 13,000 g for 10 min at 4 °C to

89

remove debris and the cleared supernatant was purified using solid phase extraction

90

(SPE) columns (which were pre-treated by using 3 mL of an aqueous 80% acetonitrile

91

solution containing 0.1% TFA, and then equilibrated with 3 mL of deionized water). The

92

enzymatically released carbohydrates were then eluted from the SPE column using 1.5

93

mL of aqueous 40% acetonitrile solution containing 0.1% TFA (v/v).

94 95

N-glycan labeling and analysis. N-glycans derived from the meat samples were

96

fluorescently labeled with 10 µL of 2-AB labeling reagent (48 mg 2-AB and 64 mg of

97

NaCNBH3 dissolved in 7 mL dimethylsulfoxide and 3 mL acetic acid), and the reaction

98

mixtures were incubated at 65°C for 2 h. 10 µL of the 2-AB labeled sample were diluted

99

with 30 μL deionized water and 35 μL of acetonitrile prior to HILIC-UPLC analysis

100

(Nexera, Shimadzu Corporation, Kyoto, Japan), and separated by an Acquity UPLC BEH

101

Glycan column (2.1×150 mm, 1.7 mm particle size; Waters, Ireland) at a column

102

temperature of 60°C. Fluorescence detection was set at the excitation/emission

103

wavelengths of 330 and 420 nm, respectively. The separation of N-glycans was carried

104

out using a previously developed linear gradient of aqueous ammonium formate (50

105

mM, pH 4.5) and acetonitrile35 (Supplementary Table S1). Elution profiles were

106

compared with a 2AB-labeled dextran oligomers (2–20 glucose units), and the absolute

107

retention time of each N-glycan peak was converted to relative retention times of the

108

dextran oligomers (GU values).

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

109 110

MALDI-TOF MS/MS analysis. Peak fractions of the 2AB-labeled UPLC fractions were

111

manually collected, dried by centrifugal evaporation, and analyzed using matrix-

112

assisted laser desorption/ionization time of-flight mass spectrometry (MALDI-TOF MS).

113

These analyses were performed on a Bruker Autoflex Speed instrument (equipped with

114

a 1000 Hz Smartbeam-II laser) in a positive ion mode using with 2, 5-dihydroxybenzoic

115

acid (10 mg/ml in acetone) as matrix. Mass fragmentation experiments (MS/MS) were

116

performed by laser induced dissociation. The acquired data were processed using the

117

Bruker Flex Analysis software (version 3.3), and the mass peaks derived from the MS-

118

and MS/MS-spectra were evaluated using GlycoWorkbench (version 1.1)36 .

119 120

Data analysis. N-glycan peaks exceeding a relative peak area over 1% were

121

considered for statistical analysis. To identify the various types of meat and the ratios of

122

adulterated meat blends, multivariate data analysis was carried out. PCA and PLS

123

analysis were performed using SPSS statistics package (version 19) and Matlab (version

124

2016b). The limit of detection (LOD) in the multivariate domain was calculated

125

according to latest development in PLS calibration37-38. The RSD of the relative peak

126

areas of each peak from 15 replicates was calculated using Microsoft Excel 2010.

127 128

Results

129

Meat N-glycan analysis and species classification.

ACS Paragon Plus Environment

Page 8 of 31

Page 9 of 31

Journal of Agricultural and Food Chemistry

130

N-glycans were released by PNGase F and separated using a 55 min HILIC-UPLC

131

program with the N-glycans being eluted between 15 and 40 min. For comparing UPLC

132

chromatograms, the retention time of glycan peaks was standardized into GU (glucose

133

unit) value, based on the retention times of dextran oligomer standards, using a linear

134

regression model (Fig. S1). Each of the five tested meat samples showed a distinct

135

overall N-glycan profile (Fig. 1A). The N-glycan profile of chicken sample bore the

136

largest number of peaks, whereas the beef sample showed the least number of peaks.

137

Duck and chicken meat N-glycans contained more relatively smaller structures

138

compared to others. Generally, although the majority of the peaks were commonly

139

shared among the tested samples, the presence of some species specific peaks and the

140

quantitative difference in relative peak area percentage of each peak made the N-

141

glycome profile of each species meat sample different from others.

142

In order to distinguish different meat species based on the protein glycosylation, the

143

N-glycan UPLC patterns of all samples were compared in terms of the relative level of

144

each target peak (percentage of the UPLC area of the target peak over the areas of all

145

peaks). Due to the difference of relative level of each peak in different meats, the overall

146

UPLC profile enables the discrimination of different meats. The identification of

147

characteristic peaks in each species from a total of 15 replicates was a critical step, and

148

may greatly affect the accuracy of the species prediction22. The UPLC N-glycan profiles

149

of 15 individuals of each species were shown in Fig. S2. The N-glycan peak with a

150

relative peak area percentage above 1% and presents in all samples was selected for

151

the PAC dataset building up. Variability of RPA of each peak within one species was

152

examined by the relative standard deviations (RSD) value of the relative peak areas

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

153

(RPA). For each species dataset, the RSD values of most peaks (account for over 60% of

154

the total peak area) were under 10%, with the maximum value being less than 17%,

155

indicating the low variability of the dataset (Table 1).

156

As shown in Figure 1B, the dataset obtained was transformed into 3 principle

157

components named PC1, PC2 and PC3. PC1 explained 45.82% of variance, and PC2 and

158

PC3 explained 33.98% and 10.76% of variance, respectively. The percentage of the total

159

variances which are able to be explained by the 3 PCs was over 90%, and therefore the

160

PCA results was regarded as reliable (>85%)39. The resulting score plot, based on the

161

PCA analysis of N-glycan profiles of the meat samples, allowed the obvious

162

discrimination of the 5 tested meat species.

163 164

Quantitative and qualitive identification of binary adulterated meat samples.

165

After the analysis of the meat samples of individual species, samples of beef blended

166

with various levels (between 0% and 50%) of added duck meat were examined (Fig.

167

2A). The overall N-glycan UPLC profile of the analyzed samples gradually changed with

168

increasing amounts of added duck meat. When compared with the N-glycan profile

169

derived from the pure beef sample, the peak area of each of the four peaks (peak 1

170

eluted at GU 6.12, peak 2 eluted at GU 6.89, peak 3 eluted at GU 7.30, and peak 4 eluted

171

at GU 7.67) over all peak areas increased, whereas the relative area of peak 5 (eluted at

172

GU 8.92) and peak 6 (eluted at GU 9.12) decreased. The relative area of peaks 1-6 were

173

further visualized using boxplot analysis (Fig. 2B) and it was shown that the peak area

174

percentage of the target peaks obviously changed with the added level of duck meat,

ACS Paragon Plus Environment

Page 10 of 31

Page 11 of 31

Journal of Agricultural and Food Chemistry

175

indicating the possible ability of discriminating the meat adulteration based on the

176

change of relative amount of each peak.

177

Principal component analysis was further applied in qualitative discrimination of

178

mixed meat samples. Beef sample was mixed with lower value duck meat at various

179

percentages including 5%, 10%, 20%, 30%, 40% and 50%. First, three principal

180

components were chosen and the resulted scatter plot showed (Fig. 3A) that the mixed

181

beef-duck samples can be separated from pure beef and other pure lower value meat

182

species. Importantly, cluster of mixed beef-duck samples were obviously placed closer

183

to beef and duck, while far from pork and chicken, indicating the ability of identifying

184

the spiked meat from the mixed sample using meat N-glycan profile. Subsequently, 2D

185

principal component analysis of beef, duck and mixed beef-duck samples exhibited

186

better ability of distinguishing all mixed samples from pure beef and pure duck (Fig. 3B).

187

Furthermore, when comparing with beef only, 2D principal components analysis of

188

mixed beef-duck samples containing various levels of duck meat showed that the mixed

189

sample with each different level of spiked duck can be isolated from others, indicating

190

the ability of discriminating the addition level of adulterated meat using this method

191

(Fig. 3C). Generally, the built PCA models based on meat N-glycan UPLC profile allowed

192

to effectively tell whether the known meat samples are pure or adulterated.

193

Furthermore, it can also identify which meat species and to which level this meat

194

species is used in adulteration.

195

In order to quantify the amount of duck meat spiked in beef, a PLS (partial least

196

squares) regression model, which is an effective and convenient quantification method

197

for binary mixtures, was employed (Fig. 4). Beef spiked with seven different levels of

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

198

duck meat (0%, 5%, 10%, 20%, 30%, 40% and 50%) were investigated for the

199

quantification of the adulteration. 114 of mixed duck/beef meat samples were used for

200

calibration and 28 mixed meat samples were used for validation. The correlation

201

coefficients of calibration (R2C >0.99), the validation (R2V>0.99), the lower root mean

202

square errors of calibration (RMSEC, 1.35%) and the prediction (RMSEP, 1.58%) of the

203

PLS model confirmed its effectiveness and accuracy40. In addition, the residual

204

predictive deviation (RPD) values calculated to be 13.06 (>8) for both calibration and

205

validation, indicating the suitability of this method for analytical tasks41-42. This result

206

showed that the built statistical model was robust enough to predict the addition of

207

duck meat in beef samples with a predicted limit of detection (LOD) of 2.2%. However,

208

given that adulteration of meat samples in industry or on the market usually contains

209

higher amounts of contaminant, we suggest that this method is feasible for the

210

detection of duck meat in beef samples in practice.

211 212

Analysis of trinary adulterated meat samples.

213

In order to further verify the applicability of the N-glycosylation based method in

214

multi-species meat authentification, trinary-species mixed meat samples were also

215

examined. Beef samples were blended with the same concentration of both pork and

216

chicken meat (10% and 20% were used) and compared with pure beef and the binary

217

beef mixtures containing either pork or chicken meat. As shown in Figure 5, the overall

218

N-glycan UPLC profile of all the samples showed differences in relative amount of some

219

peaks. For instance, the peak area of the five highlighted peaks (peak 1 eluted at GU

220

5.20, peak 2 eluted at GU 5.96, peak 3 eluted at GU 6.96, and peak 4 eluted at GU 7.67,

ACS Paragon Plus Environment

Page 12 of 31

Page 13 of 31

Journal of Agricultural and Food Chemistry

221

peak 6 eluted at GU 9.98) over all peaks areas increased when compared with the N-

222

glycan profile derived from the pure beef sample, whereas the relative area of peak 5

223

(eluted at GU 8.12) decreased.

224

To further analyze the difference of the samples more accurately, PCA was performed

225

to discriminate all the mixed meat samples. Beef-pork-chicken samples presented to be

226

closer to beef, pork and chicken meat samples, whereas to be farther to duck meat

227

samples, indicating the ability of preliminarily identifying the species of trinary-species

228

meat samples (Fig. 6A). Meanwhile, when compared with binary mixed meat samples

229

such as beef-pork and beef-chicken, clusters of beef-pork-chicken meat samples were

230

able to be distinguished from others (Fig. 6B). Furthermore, Figure 6C displayed that

231

both mixed binary and ternary samples containing different concentrations of spiked

232

meats can all be distinguished from each other, indicating the ability of quantitatively

233

discriminating the trinary-species adulterated meat using this method.

234 235

N-glycan structural identification.

236

In order to know the N-glycan structures which are different between beef and duck

237

samples, the 6 N-glycan peaks, experienced the most significant change during the

238

adulteration of beef with duck, were further analyzed by MALDI-TOF-MS and MS/MS

239

for detailed structures. As shown in Table 2, the mass spectra of these 6 peculiar peaks

240

could be identified by the obtained m/z values (Fig. 7), and the detailed structures were

241

proposed based on the MS/MS analysis of the fragmentation of the peaks (Fig.S3). It can

242

be seen from Figure 2A that peak 2 and peak 4, featured with tri-antenna bearing core

243

fucosylation and terminal galactosylation in β1,4 linkage, are the most abundant

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

244

structures specific for duck meat. Peak 5 and 6 were the most abundant structures in

245

beef, with the former being a hybrid type N-glycan which is not widely present in

246

nature and the latter containing α1,3 linked terminal galactose which is known to be the

247

antigen epitope when enters into human body43. These structures are all the key factors

248

responsible for the change of N-glycan UPLC profile when meat samples are mixed. It

249

can be seen from these detailed characteristic structures that there are significant

250

compositional differences between different species. The structural specificity is

251

determined by the genes encoding the glycoenzymes which are involved in the

252

synthesis of N-glycans. Therefore, the significance of N-glycan structural difference

253

between meats is stable and reliable for species discrimination. Furthermore, by

254

knowing the detailed structures of these peaks, it will enable the detection of beef

255

adulteration with duck using mass spectrometry to directly detect the characteristic

256

structures in duck, instead of using UPLC which requires fluorescent labeling and time-

257

cost pretreatment procedure.

258 259 260

Discussion

261

The protein N-glycosylation of biological materials is generally species-specific and

262

therefore extensively studied on its potential for acting as biomarkers in various areas,

263

especially in medical biology44-46. In first studies, identification of protein N-

264

glycosylation has been reported to be a reliable tool in food safety and food quality

265

control research30,

266

structural difference was found between human milk and bovine milk49. These findings

47-48.

Our previous study also showed that significant N-glycan

ACS Paragon Plus Environment

Page 14 of 31

Page 15 of 31

Journal of Agricultural and Food Chemistry

267

suggest the idea that N-glycan profiles may be of interest for meat species specific

268

biomarker and therefore can be used as a tool in meat authentication. As shown in this

269

study, the meat protein N-glycan UPLC profile is shown to be species specific and thus a

270

reliable method for meat species identification and for the detection of meat

271

adulteration which normally occurs in high value meat such as beef. Additionally, it can

272

also tell the species of the meat which is spiked in the major carrier and even the ratio

273

of the adulteration, for both binary and multi-species adulteration. This is for the first

274

time, to the authors’ knowledge, that protein N-glycosylation is explored in meats and

275

used to establish a novel robust and versatile method for qualitative and quantitative

276

analysis of meat authentication detection. On the other hand, only 5 meat species were

277

employed in this study and it certainly needs to extend the range to validate the

278

applicability of this method when used for meat species other than the ones listed in

279

this study.

280

The most commonly used approaches for meat species identification at present

281

include PCR, a DNA-based molecular method, and spectroscopic methods. However,

282

sample pollution and DNA degradation during meat processing may greatly affect the

283

PCR results50. The N-glycan profiles of meats were shown to stay stable upon various

284

processing procedures such as boiling, roasting, frying, microwave and high-pressure

285

treatments (data now shown), indicating that our method probably bears advantage

286

over other methods such as PCR when used for processed meat adulteration detection.

287

Spectroscopic methods perform a rapid analysis of meat samples with minimum pre-

288

processing requirements, while the data processing is complicated. For example, each

289

type of product matrix (species and processing) needs a separate calibration which is

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

290

time consuming. In this study, the established protein glycosylation based method

291

presents to be relatively sensitive, specific and repeatable. The length of workflow

292

reported here was used to establish reliability of the analysis, but this procedure can be

293

shortened based on our previous study, according to which a rapid sample preparation

294

can be done within 4 h35.

295

The limit of detection (LOD) level for duck meat adulterated beef sample was found

296

to be 2.2% using the N-glycan profile based method in this study, which is low enough

297

for meat authentication on market as the adulteration ratio in reality is generally over

298

10% for the producers to make profit. However, this LOD value is not as low as those of

299

other meat adulteration methods and optimization of this method, such as the

300

optimization of meat sample preparation, N-glycan purification and enrichment, novel

301

derivatization and even the UPLC detection, needs to be performed to realize the low

302

LOD comparable with that of other methods.

303

The structures of the six peculiar peaks from beef and duck meat N-glycans were also

304

identified in details based on MALDI-TOF-MS and MS/MS analysis. The structural

305

differences between beef and duck are relatively big. It is well known that

306

glycoenzymes are responsible for the glycan synthesis and their abundance specificity

307

varies with species, organs and even tissues. It therefore can be deducted from the

308

structures that the level of each of glycosyltransferases differs between animals which

309

eventually resulted in the significant structural difference. In addition to understanding

310

the working principle of the novel UPLC-based method established in this study, the

311

detailed structure identification also reveals that the N-glycan based approach for meat

312

adulteration is not only limited to UPLC detection, but can also be cooperated with

ACS Paragon Plus Environment

Page 16 of 31

Page 17 of 31

Journal of Agricultural and Food Chemistry

313

MALDI-TOF mass spectrometry analysis. In other words, this finding will enable the

314

detection of meat adulteration using mass spectrometry through directly detecting the

315

characteristic structures species without the pretreatments required by UPLC. The MS-

316

based method will thus be more rapid, easy and robust than the UPLC-based method.

317

Given all the above-stated reliabilities, this proof-of-concept work may also be of

318

interest for the analysis of other adulterated meat blends, or even other type of foods, in

319

future.

320 321

Abbreviations

322

HILIC-UPLC, hydrophobic interaction ultra-performance liquid chromatography; PCA,

323

principal component analysis; PLS, partial least square; RMSEC, root mean square

324

errors of calibration; RMSEP, root mean square errors of prediction; RPD, residual

325

predictive deviation; LOD, limit of detection; PCR, polymerase chain reaction; ACN,

326

acetonitrile;

327

aminobenzamide; TCA, trichloroacetic acid; MALDI-TOF MS, matrix-assisted laser

328

desorption/ionization time of-flight mass spectrometry; GU, glucose unit; PC, principle

329

component; R2C, correlation coefficients of calibration; R2V, correlation coefficients of

330

validation.

SPE,

supelclean

ENVI-Carb

solid-phase

331 332

Funding Sources

ACS Paragon Plus Environment

extraction;

2-AB,

2-

Journal of Agricultural and Food Chemistry

333

This work was financially supported by the National Natural Science Foundation of

334

China (NSFC, grant number 31871793 to J.V.) and National Key Research and

335

Development Plan (grant number 2018YFC1602804 to C.B. L.).

336 337 338 339 340

Supporting Information The Supporting Information is available.

341 342 343

References

344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360

1. 2. 3. 4. 5. 6. 7.

Rahmati, S.; Julkapli, N. M.; Yehye, W. A.; Basirun, W. J., Identification of meat origin in food products–A review. Food Control 2016, 68, 379-390. Henchion, M.; McCarthy, M.; Resconi, V. C.; Troy, D., Meat consumption: trends and quality matters. Meat Sci. 2014, 98 (3), 561-568. Kumar, A.; Kumar, R. R.; Sharma, B. D.; Gokulakrishnan, P.; Mendiratta, S. K.; Sharma, D., Identification of Species Origin of Meat and Meat Products on the DNA Basis: A Review. Crit. Rev. Food Sci. Nutr. 2013, 55 (10), 1340-1351. Ballin, N. Z.; Vogensen, F. K.; Karlsson, A. H., Species determination - Can we detect and quantify meat adulteration? Meat Sci. 2009, 83 (2), 165-174. Fajardo, V.; González, I.; Rojas, M.; García, T.; Martín, R., A review of current PCR-based methodologies for the authentication of meats from game animal species. Trends Food Sci. Technol. 2010, 21 (8), 408-421. von Bargen, C.; Brockmeyer, J.; Humpf, H.-U., Meat Authentication: A New HPLC–MS/MS Based Method for the Fast and Sensitive Detection of Horse and Pork in Highly Processed Food. J. Agric. Food Chem. 2014, 62 (39), 9428-9435. Hossain, M. A.; Ali, M. E.; Abd Hamid, S. B.; Asing; Mustafa, S.; Mohd Desa, M. N.; Zaidul, I. S., Double Gene Targeting Multiplex Polymerase Chain Reaction-Restriction Fragment Length

ACS Paragon Plus Environment

Page 18 of 31

Page 19 of 31

361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408

Journal of Agricultural and Food Chemistry

8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21. 22. 23. 24. 25.

Polymorphism Assay Discriminates Beef, Buffalo, and Pork Substitution in Frankfurter Products. J. Agric. Food Chem. 2016, 64 (32), 6343-6354. Kesmen, Z.; Celebi, Y.; Güllüce, A.; Yetim, H., Detection of seagull meat in meat mixtures using real-time PCR analysis. Food Control 2013, 34 (1), 47-49. Ballin, N. Z., Authentication of meat and meat products. Meat Sci. 2010, 86 (3), 577-587. Li, Y.; Zhang, Y.; Li, H.; Zhao, W.; Guo, W.; Wang, S., Simultaneous determination of heat stable peptides for eight animal and plant species in meat products using UPLC-MS/MS method. Food Chem. 2018, 245, 125-131. Von Bargen, C.; Brockmeyer, J.; Humpf, H. U., Meat authentication: a new HPLC-MS/MS based method for the fast and sensitive detection of horse and pork in highly processed food. J. Agric. Food Chem. 2014, 62 (39), 9428-9435. Watson, A. D.; Gunning, Y.; Rigby, N. M.; Philo, M.; Kemsley, E. K., Meat Authentication via Multiple Reaction Monitoring Mass Spectrometry of Myoglobin Peptides. Anal. Chem. 2015, 87 (20), 10315-10322. Djurdievic, N.; Sheu, S. C.; Hsieh, Y. H. P., Quantitative detection of poultry in cooked meat products. J. Food Sci. 2005, 70 (9), C586-C593. Mandli, J.; El Fatimi, I.; Seddaoui, N.; Amine, A., Enzyme immunoassay (ELISA/immunosensor) for a sensitive detection of pork adulteration in meat. Food Chem. 2018, 255, 380-389. Ali, M. E.; Ahamad, M. N. U.; Asing; Hossain, M. A. M.; Sultana, S., Multiplex polymerase chain reaction-restriction fragment length polymorphism assay discriminates of rabbit, rat and squirrel meat in frankfurter products. Food Control 2018, 84, 148-158. Kumar, D.; Singh, S. P.; Karabasanavar, N. S.; Singh, R.; Umapathi, V., Authentication of beef, carabeef, chevon, mutton and pork by a PCR-RFLP assay of mitochondrial cytb gene. J. Food Sci. Technol. 2014, 51 (11), 3458-3463. Lubis, H.; Salihah, N. T.; Hossain, M. M.; Ahmed, M. U., Development of fast and sensitive realtime qPCR assay based on a novel probe for detection of porcine DNA in food sample. Lwt-Food Sci.Technol. 2017, 84, 686-692. Floren, C.; Wiedemann, I.; Brenig, B.; Schutz, E.; Beck, J., Species identification and quantification in meat and meat products using droplet digital PCR (ddPCR). Food Chem. 2015, 173, 1054-1058. Fang, X.; Zhang, C., Detection of adulterated murine components in meat products by TaqMan(c) real-time PCR. Food Chem. 2016, 192, 485-490. Bilge, G.; Velioglu, H. M.; Sezer, B.; Eseller, K. E.; Boyaci, I. H., Identification of meat species by using laser-induced breakdown spectroscopy. Meat Sci. 2016, 119, 118-122. Boyaci, I. H.; Temiz, H. T.; Uysal, R. S.; Velioglu, H. M.; Yadegari, R. J.; Rishkan, M. M., A novel method for discrimination of beef and horsemeat using Raman spectroscopy. Food Chem. 2014, 148, 37-41. Xiong, Z.; Sun, D.-W.; Pu, H.; Zhu, Z.; Luo, M., Combination of spectra and texture data of hyperspectral imaging for differentiating between free-range and broiler chicken meats. LWT Food Sci.Technol. 2015, 60 (2), 649-655. Hu, Y.; Zou, L.; Huang, X.; Lu, X., Detection and quantification of offal content in ground beef meat using vibrational spectroscopic-based chemometric analysis. Sci. Rep. 2017, 7 (1), 15162. Rodriguez-Ramirez, R.; Gonzalez-Cordova, A. F.; Vallejo-Cordoba, B., Review: Authentication and traceability of foods from animal origin by polymerase chain reaction-based capillary electrophoresis. Anal. Chim. Acta 2011, 685 (2), 120-126. Kumar, A.; Kumar, R. R.; Sharma, B. D.; Gokulakrishnan, P.; Mendiratta, S. K.; Sharma, D., Identification of species origin of meat and meat products on the DNA basis: a review. Crit. Rev. Food Sci. Nutr. 2015, 55 (10), 1340-1351.

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455

26. Sentandreu, M. A.; Fraser, P. D.; Halket, J.; Patel, R.; Bramley, P. M., A Proteomic-Based Approach for Detection of Chicken in Meat. J. Proteome Res. 2010, 9 (7), 3374-3383. 27. Gavrilov, Y.; Shental-Bechor, D.; Greenblatt, H. M.; Levy, Y., Glycosylation May Reduce Protein Thermodynamic Stability by Inducing a Conformational Distortion. J. Phys. Chem. Lett. 2015, 6 (18), 3572-3577. 28. Vanhooren, V.; Liu, X. E.; Franceschi, C.; Gao, C. F.; Libert, C.; Contreras, R.; Chen, C., N-glycan profiles as tools in diagnosis of hepatocellular carcinoma and prediction of healthy human ageing. Mech. Ageing. Dev. 2009, 130 (1-2), 92-97. 29. Mullins, R. J.; James, H.; Platts-Mills, T. A.; Commins, S., Relationship between red meat allergy and sensitization to gelatin and galactose-alpha-1,3-galactose. J. Allergy Clin. Immunol. 2012, 129 (5), 1334-1342 e1. 30. Meli, V. S.; Ghosh, S.; Prabha, T. N.; Chakraborty, N.; Chakraborty, S.; Datta, A., Enhancement of fruit shelf life by suppressing N-glycan processing enzymes. Proc. Natl. Acad. Sci. U. S. A. 2010, 107 (6), 2413-2418. 31. Shim, E. K.; Chandra, G. F.; Pedireddy, S.; Lee, S. Y., Characterization of swiftlet edible bird nest, a mucin glycoprotein, and its adulterants by Raman microspectroscopy. J. food sci. technol. 2016, 53 (9), 3602-3608. 32. Wang, T.; Hu, X. C.; Cai, Z. P.; Voglmeir, J.; Liu, L., Qualitative and Quantitative Analysis of Carbohydrate Modification on Glycoproteins from Seeds of Ginkgo biloba. J. Agric. Food Chem. 2017, 65 (35), 7669-7679. 33. Qiao, D.; Xu, J.; Qin, P.; Yao, L.; Lu, J.; Eremin, S.; Chen, W., Highly Simple and Sensitive Molecular Amplification-Integrated Fluorescence Anisotropy for Rapid and On-Site Identification of Adulterated Beef. Anal. Chem. 2018, 90 (12), 7171-7175. 34. Zheng, X.; Li, Y.; Wei, W.; Peng, Y., Detection of adulteration with duck meat in minced lamb meat by using visible near-infrared hyperspectral imaging. Meat Sci. 2019, 149, 55-62. 35. Du, Y. M.; Xia, T.; Gu, X. Q.; Wang, T.; Ma, H. Y.; Voglmeir, J.; Liu, L., Rapid Sample Preparation Methodology for Plant N-Glycan Analysis Using Acid-Stable PNGase H+. J. Agric. Food Chem. 2015, 63 (48), 10550-10555. 36. Ceroni, A.; Maass, K.; Geyer, H.; Geyer, R.; Dell, A.; Haslam, S. M., GlycoWorkbench: A Tool for the Computer-Assisted Annotation of Mass Spectra of Glycans+. J. Proteome Res. 2008, 7 (4), 1650-1659. 37. Shen, T.; Kong, W.; Liu, F.; Chen, Z.; Yao, J.; Wang, W.; Peng, J.; Chen, H.; He, Y., Rapid Determination of Cadmium Contamination in Lettuce Using Laser-Induced Breakdown Spectroscopy. Molecules 2018, 23 (11). 38. VALDERRAMA, P.; B, J. W.; BRAGA; POPPI, R. J., Variable Selection, Outlier Detection, and Figures of Merit Estimation in a Partial Least-Squares Regression Multivariate Calibration Model. A Case Study for the Determination of Quality Parameters in the Alcohol Industry by NearInfrared Spectroscopy. J. Agric. Food Chem. 2007, 55 (21), 8331-8338. 39. He, Y.; Feng, S.; Deng, X.; Li, X., Study on lossless discrimination of varieties of yogurt using the Visible/NIR-spectroscopy. Food Res. Int. 2006, 39 (6), 645-650. 40. Qiao, L.; Tang, X.; Dong, J., A feasibility quantification study of total volatile basic nitrogen (TVBN) content in duck meat for freshness evaluation. Food Chem. 2017, 237, 1179-1185. 41. Kartakoullis, A.; Comaposada, J.; Cruz-Carrión, A.; Serra, X.; Gou, P., Feasibility study of smartphone-based Near Infrared Spectroscopy (NIRS) for salted minced meat composition diagnostics at different temperatures. Food Chem. 2019, 278, 314-321. 42. Conzen, J., Multivariate Calibration: A practical guide for developing methods in the quantitative analytical chemistry. Bruker Optik GmbH: Ettlingen, Germany 2006.

ACS Paragon Plus Environment

Page 20 of 31

Page 21 of 31

456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481

Journal of Agricultural and Food Chemistry

43. Commins, S. P.; James, H. R.; Kelly, L. A.; Pochan, S. L.; Workman, L. J.; Perzanowski, M. S.; Kocan, K. M.; Fahy, J. V.; Nganga, L. W.; Ronmark, E.; Cooper, P. J.; Platts-Mills, T. A., The relevance of tick bites to the production of IgE antibodies to the mammalian oligosaccharide galactose-alpha-1,3-galactose. J. Allergy Clin. Immunol. 2011, 127 (5), 1286-1293. 44. Li, D.; Zhang, P.; Li, F.; Chi, L.; Zhu, D.; Zhang, Q., Recognition of N-glycoforms in human chorionic gonadotropin by monoclonal antibodies and their interaction motifs. J. Biol. Chem. 2015, 290 (37), 22715-22723. 45. Ito, E.; Oka, R.; Ishii, T.; Korekane, H.; Kurimoto, A.; Kizuka, Y.; Kitazume, S.; Ariki, S.; Takahashi, M.; Kuroki, Y.; Kida, K.; Taniguchi, N., Fucosylated surfactant protein-D is a biomarker candidate for the development of chronic obstructive pulmonary disease. J. Proteomics 2015, 127 (Pt B), 386-394. 46. Pan, S.; Chen, R.; Tamura, Y.; Crispin, D. A.; Lai, L. A.; May, D. H.; McIntosh, M. W.; Goodlett, D. R.; Brentnall, T. A., Quantitative glycoproteomics analysis reveals changes in N-glycosylation level associated with pancreatic ductal adenocarcinoma. J. Proteome Res. 2014, 13 (3), 12931306. 47. Cao, X.; Song, D.; Yang, M.; Yang, N.; Ye, Q.; Tao, D.; Liu, B.; Wu, R.; Yue, X., Comparative Analysis of Whey N-Glycoproteins in Human Colostrum and Mature Milk Using Quantitative Glycoproteomics. J. Agric. Food Chem. 2017, 65 (47), 10360-10367. 48. Xu, H.; Zheng, L.; Xie, Y.; Zeng, H.; Fan, Q.; Zheng, B.; Zhang, Y., Identification and determination of glycoprotein of edible brid's nest by nanocomposites based lateral flow immunoassay. Food Control 2019, 102, 214-220. 49. Wang, W. L.; Du, Y. M.; Wang, W.; Conway, L. P.; Cai, Z. P.; Voglmeir, J.; Liu, L., Comparison of the bifidogenic activity of human and bovine milk N-glycome. J. Funct. Food. 2017, 33, 40-51. 50. Di Pinto, A.; Bottaro, M.; Bonerba, E.; Bozzo, G.; Ceci, E.; Marchetti, P.; Mottola, A.; Tantillo, G., Occurrence of mislabeling in meat products using DNA-based assay. J. Food Sci. Technol. 2015, 52 (4), 2479-2484.

482

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 22 of 31

Tables Table 1 The relative peak area (RPA) of characteristic peaks in five meat species samples (n = 15).

Peak

Beef

Chicken

Pork

Duck meat

Mutton

Average

RSD

Average

RSD

Average

RSD

Average

RSD

Average

RSD

of RPA (%)

(%)

of RPA (%)

(%)

of RPA (%)

(%)

of RPA (%)

(%)

of RPA (%)

(%)

1

2.18

13.27

5.75

14.80

2.03

14.64

5.34

15.44

1.24

12.33

2

1.91

13.54

2.82

10.10

4.01

12.08

1.69

6.79

1.44

15.19

3

3.54

8.34

9.79

7.26

4.75

11.42

11.84

6.40

4.40

9.75

4

2.85

4.05

9.54

5.37

4.65

9.22

4.31

13.99

3.45

9.36

5

5.07

11.04

8.73

7.84

2.06

16.61

8.42

6.05

4.71

8.38

6

7.98

8.63

7.89

8.59

4.53

14.30

7.86

8.81

1.84

15.44

7

9.27

6.06

8.60

7.10

5.59

12.61

2.62

15.54

5.24

9.78

8

10.12

6.45

5.53

12.81

13.75

2.86

2.24

7.94

3.62

13.15

9

12.21

8.76

9.51

9.69

13.52

8.03

13.45

9.20

15.31

9.79

10

12.33

6.51

6.50

8.52

4.10

9.54

5.73

8.14

12.16

9.84

11

10.47

9.72

3.54

13.97

4.53

9.02

5.58

9.84

14.49

8.10

12

6.00

9.59

6.54

9.10

22.15

7.62

2.77

12.73

4.87

14.21

13

3.23

11.63

7.15

8.12

3.18

9.64

8.02

9.45

8.76

8.30

14

0.00

0.00

3.18

12.82

0.00

0.00

0.00

0.00

6.46

6.82

15

0.00

0.00

0.00

0.00

0.00

0.00

0.00

0.00

4.15

8.47

No.

ACS Paragon Plus Environment

Page 23 of 31

Journal of Agricultural and Food Chemistry

Table 2 Structures of beef and duck N-glycans based on the results of MALDI-TOF-MS/MS. Number Structure Name GU M/Z Source Beef ,

M5

6.12

1377.92

FA3G1

6.89

1970.58

Duck

A3G2

7.30

1987.48

Duck

FA3G2

7.67

2132.83

Duck

FM5A1G1S1

8.92

2359.47

Beef

FA2G(4)2G(3)2

9.12

2254.20

Beef

1

duck

2

3

4

5

6 M: mannose; A: N-acetyl-hexosamine; G: galactose; F: fucose; S: sialic acid.

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figures

Figure 1 Identification of five meat species. A: UPLC- profiles of N-glycomes derived from five different species of meat samples. B: PCA score plot of UPLC N-glycome profiles of five meat species. (One column)

ACS Paragon Plus Environment

Page 24 of 31

Page 25 of 31

Journal of Agricultural and Food Chemistry

Figure 2 UPLC-based N-glycosylation pattern analysis of beef samples adulterated with different levels of duck. A: UPLC- profiles of N-glycome derived from mixed meat samples. B: boxplot of relative peak areas of 6 selected N-glycans. In each box, the central mark indicates the median, and the bottom and top edges of the box indicate the 25th and 75th percentiles, respectively. (Two column)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 3 PCA score plot of UPLC-profiles of different mixed meat samples. A: classification of mixed meat samples, beef, pork, chicken meat and duck meat. B: separation of mixed meat samples from beef and duck meat. C: identification of seven different levels of mixed meat samples. (Two column)

ACS Paragon Plus Environment

Page 26 of 31

Page 27 of 31

Journal of Agricultural and Food Chemistry

Figure 4 PLS calibration and validation models for quantitative determination of duck adulterated beef samples. (One column)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 5 UPLC- profiles of N-glycome derived from mixed meat samples. (One column)

ACS Paragon Plus Environment

Page 28 of 31

Page 29 of 31

Journal of Agricultural and Food Chemistry

Figure 6 PCA score plot of UPLC N-glycome profiles of mixed meat samples. A: classification of trinary mixed meat samples, beef, pork, chicken meat and duck meat. B: separation of mixed meat samples from each pure species beef, pork and chicken meat. C: identification of mixed meat samples with different concentrations of adulterated meat. (Two columns)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 7 Structural analysis of 2-AB labeled N-glycans derived from beef and duck by MALDI-TOFMS. (Two Columns)

ACS Paragon Plus Environment

Page 30 of 31

Page 31 of 31

Journal of Agricultural and Food Chemistry

Graphic for table of contents

ACS Paragon Plus Environment