MALDI-TOF-MS Platform for Integrated Proteomic and Peptidomic

Jun 22, 2015 - Adulteration of ovine, caprine, and buffalo milks with more common bovine material occurs for economic reasons and seasonal availabilit...
0 downloads 0 Views 1MB Size
Subscriber access provided by NEW YORK UNIV

Article

A MALDI-TOF-MS platform for integrated proteomic and peptidomic profiling of milk samples allows rapid detection of food adulterations Mauro Sassi, Simona Arena, and Andrea Scaloni J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.5b02384 • Publication Date (Web): 22 Jun 2015 Downloaded from http://pubs.acs.org on June 24, 2015

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 free 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 accessible to all readers and 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.

Journal of Agricultural and Food Chemistry 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 35

Journal of Agricultural and Food Chemistry

A MALDI-TOF-MS platform for integrated proteomic and peptidomic profiling of milk samples allows rapid detection of food adulterations

Mauro Sassi *, Simona Arena * and Andrea Scaloni

Proteomics & Mass Spectrometry Laboratory, ISPAAM, National Research Council, Naples, Italy

* These authors gave an equal contribution to this work

Keywords: MALDI-TOF-MS; protein profiling, peptide profiling, milk analysis, speciation, thermal treatment.

Corresponding author:

Andrea Scaloni Proteomics & Mass Spectrometry Laboratory ISPAAM, National Research Council via Argine 1085, Naples, I-80147, Italy Email: [email protected] Phone: +39 (081) 5966006 Fax: +39 (081) 5965291

1   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 2 of 35

Abstract Adulteration of ovine, caprine and buffalo milks with more common bovine material occurs for economic reasons and seasonal availability. Frauds are also associated with the use of powdered milk instead of declared, fresh material. In this context, various analytical methods have been adapted to dairy science applications with the aim to evaluate adulteration of milk samples, although time-consuming, suitable only for speciation or thermal treatment analysis, or useful for a specific fraud type. An integrated MALDI-TOF-MS platform for the combined peptidomic and proteomic profiling of milk samples is here presented, which allows rapid detection of illegal adulterations due to the addition of either non-declared bovine material to water buffalo, goat and ovine milks, or of powdered bovine milk to the fresh counterpart. Peptide and protein markers of each animal milk were identified after direct analysis of a large number of diluted skimmed and/or enriched diluted skimmed filtrate samples. In parallel, markers of thermal treatment were characterized in different types of commercial milks. Principal components scores of ad hoc prepared species- or thermal treatment-associated adulterated milk samples were subjected to partial least squares regression, permitting a fast accurate estimate of the fraud extents in test samples either at protein and peptide level. With respect to previous reports on MALDI-TOF-MS protein profiling methodologies for milk speciation, this study extends this approach to the analysis of the thermal treatment and introduces and independent, complementary peptide profiling measurement, which integrates protein data with additional information on peptides, validating final results and ultimately broadening the method applicability.

2   

ACS Paragon Plus Environment

Page 3 of 35

Journal of Agricultural and Food Chemistry

Introduction Raw milk is one of the main constituent of the human diet, being an essential source of nutrients for the newborns and an important food for adult individuals. It is also used to prepare a multitude of dairy products that, depending on the making dairy animal species and transformation procedures, are characterized by different organoleptic/nutritional properties and economic values 1,2. In this context, the exclusive use of a specific, pure raw milk is mandatory in the manufacturing of traditional high-grade European Protected Designation of Origin cheeses, such as the water buffalo Mozzarella, Feta, Roquefort, Idiazabal, Pecorino Romano and Fiore Sardo. However, it has been estimated that buffalo, goat and sheep milk only represent 11, 2 and 1.4% of the raw material produced worldwide and used for dairy productions, compared to the bovine counterpart (about 85%) 3. The adulteration of ovine, caprine and buffalo milk with common bovine material occurs often because of the higher prices of the first ones, and the fluctuation in their availability due to the corresponding seasonal shortages. Fraudulent activities are also associated with the addition of powdered/frozen material to fresh milk in the period of low lactation 4,5. Adulteration practices mentioned above determined an ultimate need of rapid, sensitive, and reliable methods to determine the nature of raw milk (and of the corresponding dairy products). To this purpose, various analytical techniques and deriving methods were used and/or adapted to dairy science applications with the aim to evaluate the nature of milk biomolecules. In this context, near-infrared (NIR), infrared (IR) and nuclear magnetic resonance (NMR) spectroscopic, chromatographic, electrophoretic, immunoenzymatic, sensoristic, LC-electrospray (ESI)-mass spectrometric (MS) and polymerase chain reaction-based assays were used for milk speciation and thermal treatment (TT) determination (see Suppl. Tab. S1 for application details and related literature). Above-mentioned approaches have the disadvantage of being sometimes time-consuming as result of different preparation steps or long instrument acquisition times, suitable only for a specific milk assay (speciation or TT analysis), useful for a specific adulteration type and/or eventually not adaptable for global screening of milk samples. For example, the EU reference methodology to detect bovine proteins in dairy products is based on gel isoelectric focusing of -caseins after plasminolysis (EC Regulation No. EC 273/2008) 6. However, this method suffers of interpretation difficulties due to overlapping of species-specific bands and may request integration with laborious immunoblotting steps 7. Recently, matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS)-based approaches have been developed for biological fluid profiling, showing their potential in the quali-quantitative

3   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 4 of 35

analysis of milk proteins and peptides 8. MALDI-TOF-MS has the advantage of requiring small amounts of even heterogeneous samples, permitting their analysis in a very short time. In particular, detection of -lactalbumin, lactoglobulin and 2/3-casein isoforms during protein profiling (PR-P) of whey and skimmed milk samples has been used to reveal definite species-associated milk frauds

9-13

. Linear regression analysis of the resulting mass spectra

quantified non-declared milk material even in low amounts 13. Recently, MALDI-TOF-MS PR-P of skimmed samples combined with multivariate calculation techniques allowed an accurate estimate of the levels of adulteration in bovine-goat-sheep milk mixtures14. No similar applications were developed to assess fraudulent addition of powdered material to fresh milk. Direct serum or skimmed milk PR-P applications did not find a counterpart at peptide level. In fact, peptide profiling (PE-P) approaches were developed only when contemplating additional experimental steps consisting in the enzymatic digestion (with trypsin or plasmin) of the whole milk samples or of selected proteins preventively resolved by electrophoresis

4,5,15-20

. This experimental set up severely affected the time necessary for

entire analysis and highly complicated the nature of the resulting MALDI-TOF spectra, generally determining the need of additional chromatographic steps for an accurate quantification of the adulteration percentages. In this study, an integrated platform for the combined peptidomic and proteomic profiling of milk samples has been developed to rapidly reveal illegal adulterations due to the addition of either non-declared bovine material to water buffalo, goat and ovine milks, or of powdered bovine material to the fresh counterpart. To this purpose, specific peptide and protein markers of fresh bovine milk (BM), water buffalo milk (WM), ovine milk (OM) and goat milk (GM) were identified after direct analysis of skimmed samples. Markers of TT were characterized in fresh bovine milk (BM), pasteurized milk (PM), UHT milk (UM) and powdered milk (IM) samples. Principal components (PCs) scores of had hoc prepared species- or TT-associated adulterated milk samples were subjected to partial least squares (PLS) regression, allowing an accurate estimate of the fraud extents in test samples either at protein and peptide level.

Experimental Section Sample preparation Fresh raw BM, WM, OM and GM samples were collected from local farms after morning milking, stored at 4 °C, and assayed within 5 h. Bovine PM, UM and IM samples from different (20 in number) commercial brands were

4   

ACS Paragon Plus Environment

Page 5 of 35

Journal of Agricultural and Food Chemistry

purchased from supermarkets and used before their expiration date (Suppl. Tab. S2). In all cases, 25 different samples were analyzed for each type of milk. For adulteration studies, binary mixtures of BM in WM, OM and GM counterparts, or IM in BM were prepared at 1, 2, 5, 10, 15, 20, 30, 40 and 50% v/v values. Three independent samples were prepared for each specific binary mixture. All milk samples or their binary mixtures were defatted by centrifugation at 3000 rpm, for 30 min, at 4 °C. A part of the corresponding skimmed material (20 µl) was directly analyzed by MALDI-TOF-MS for the generation of the corresponding protein profiles. The remaining material (200 µl) was diluted 1:1 v/v with water, transferred onto Amicon Ultra devices (3 kDa cut-off) (Millipore) and filtrated by centrifugation at 10,000 rpm, at 4 °C. Two technical replicates were prepared for each filtrate. Peptides from filtrates were directly analyzed by MALDI-TOF-MS or extracted with C18ZipTip devices (Millipore) and then subjected to MS analysis.

MALDI-TOF-MS profiling measurements Before MALDI-TOF-MS PR-P, skimmed milk samples were diluted 1:100 with water. Each sample (0.5 μl) was mixed with 0.5 μl of a solution of sinapinic acid (10 mg/ml) (Bruker Daltonics, Bremen, Germany) in 0.1% v/v trifluoroacetic acid (TFA), 30% v/v acetonitrile (ACN), placed onto a Ground Steel Target (Bruker Daltonics) and dried at room temperature. Before MALDI-TOF-MS PE-P, peptide samples (0.5 μl) were mixed with 0.5 μl of a solution of α-cyano-4-hydroxycinnamic acid (25 mg/ml) (Bruker Daltonics) in 0.1% v/v TFA, 30% v/v ACN, placed onto the instrument target reported above and dried at room temperature. Protein and peptide samples were analyzed by MALDI-TOF-MS using an UltraflexExtreme mass spectrometer (Bruker Daltonics) equipped with the FlexControl software package (version 3.4, Bruker Daltonics)

21

. In all cases,

spectral acquisition methods were developed to maximize the number of signals present in the mass spectra and the corresponding signal to noise ratios. In the case of PR-P measurements, spectra were recorded in the positive linear mode (LM) (laser frequency, 1000 Hz; ion source 1 voltage, 25.19 kV; ion source 2 voltage, 23.59 kV; lens voltage, 7.50 kV; sample rate, 0.31; mass range, 5000-20,000 Da). Five independent spectra (1000 shots at random positions on the same target place, for spectrum) were automatically collected, calibrated externally by using the Protein Calibration Standard 1 (Bruker Daltonics) and subsequently analyzed. In the case of PE-P measurements, spectra were

5   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 6 of 35

recorded both in the positive LM (laser frequency, 1000 Hz; ion source 1 voltage, 25.19 kV; ion source 2 voltage, 23.94 kV; lens voltage, 6.50 kV; sample rate, 0.16; mass range, 500-5000 Da), or in the positive reflector mode (RM) (laser frequency, 1000 Hz; ion source 1 voltage, 25.19 kV; ion source 2 voltage, 22.50 kV; lens voltage, 8.50 kV; reflector, 26.83 kV; reflector 2, 13.75 kV; sample rate, 0.63; mass range, 500 to 5000 Da). Ten independent spectra (500 shots at random positions on the same target place, for spectrum) for each peptide mixture were automatically collected, externally calibrated by using the Peptide Calibration Standard 2 (Bruker Daltonics) and subsequently analyzed.

Bioinformatics and statistics FlexAnalysis (version 3.4) and ClinProt Tools (version 2.2) software packages (Bruker Daltonics) were used for the analysis of all MALDI-TOF-MS data, which included spectral mass adjustment (compression by a factor of 10 in the total mass range), optional smoothing (using the Savitsky–Golay algorithm with a frame size of 25 Da), spectral baseline subtraction, normalization, internal peak alignment and peak picking. Pretreated data were then subjected to visualization and statistical analysis. Peptides or proteins showing a statistically significant difference in signal intensity or mass value were determined by means of Wilcoxon- (PWKW), Anderson-Darling- (PAD) and t- (PTTA) test. Class prediction model was set up by Genetic Alghorithms (GA). Discriminant peaks were considered those presenting at least a p-value < 0.000001 among the PWKW, PAD and PTTA ones. Finally, it was performed a principal component analysis (PCA) of the spectra, which was carried out by an external MATLAB software tool integrated into ClinPro Tools software. Regarding adulterated milk, corresponding principal components (PCs) scores were subjected to PLS regression by XLSTAT software (Microsoft). Adulteration prediction at two points (5 and 15% v/v) was obtained by the forecast function of Excel software (Microsoft), using PLS results as matrix reference data.

Results and Discussion In order to define an optimal method for the suitable discrimination of milk samples by MALDI-TOF-MS profiling, dedicated procedures were set up for the combined analysis of the corresponding protein and peptide components. Due to the composition complexity of each milk type, various experimental protocols were developed to this purpose,

6   

ACS Paragon Plus Environment

Page 7 of 35

Journal of Agricultural and Food Chemistry

which were compared a posteriori for their discrimination performances. PR-P assays were performed directly on diluted skimmed milk samples

9-14

. On the other hand, PE-P measurements were realized on diluted low-molecular

mass skimmed milk filtrates that were then straightly analyzed by MALDI-TOF-MS or extracted with concentration devices before MS analysis. Attempts to assay peptides directly in diluted skimmed milk samples were unsuccessful in terms of recognition performances (data not shown). Different instrument acquisition ranges were assessed for protein and peptide analysis; in the latter case, the effect of a spectral acquisition in LM and RM was also evaluated (see below). Optimized procedures were used either for milk speciation, for the characterization of commercial bovine milk samples subjected to various TTs and finally for the detection of adulterated test samples as result of the addition of non-declared material. In the whole, a dataset of 3500 MALDI-TOF mass spectra was recorded for the speciation of raw milk from four animal species. In particular, 500 mass spectra were acquired for the corresponding PR-P (5 spectral replicas for 25 milk samples from totally different individuals of each animal species); a representative example is reported in Fig. 1A-D. A thousand mass spectra were acquired for the corresponding PE-P of milk filtrates (10 spectral replicas for 25 milk samples from totally different individuals of each animal species) either in LM and RM; representative examples are reported in Suppl. Fig. S1A-D and Suppl. Fig. S2A-D. Five hundred mass spectra were acquired for the corresponding PE-P of milk filtrates subjected to C18ZipTip enrichment (5 spectral replicas for 25 milk samples from totally different individuals of each animal species) either in LM and RM; representative examples are reported in Suppl. Fig. S3A-D and Fig. 2A-D. A dataset analogous to that reported above was recorded in the case of BM and commercial thermally-treated PM, UM and IM samples. A representative example of the corresponding PR-P spectra is reported in Fig. 3A-D. Similarly, exemplificative PE-P spectra of milk filtrates acquired either in LM and RM are shown in Suppl. Fig. S4A-D and Suppl. Fig. S5A-D, respectively. Peptide MALDI-TOF spectra of milk filtrates subjected to C18ZipTip enrichment recorded either in LM and RM are presented in Suppl. Fig. S6A-D and Fig. 4A-D, respectively. Finally, a dataset of 3780 MALDI-TOF mass spectra was recorded for the detection of adulterated milk samples as result of the addition of non-declared material. Experiments were performed either to identify adulteration from raw BM added to WM, OM, or GM, and from IM added to BM. In particular, 540 mass spectra were acquired for the corresponding PR-P (5 spectral replicas of the 9 percentage values assayed for each kind of adulteration; experiments

7   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 8 of 35

were performed on triplicate samples) (data not shown); 1080 mass spectra were acquired for the corresponding PE-P of milk filtrates (10 spectral replicas of the 9 percentage values assayed for each kind of adulteration; experiments were performed on triplicate samples) either in LM and RM (data not shown); 540 mass spectra were acquired for the corresponding PE-P of milk filtrates subjected to C18ZipTip enrichment (5 spectral replicas of the 9 percentage values assayed for each kind of adulteration; experiments were performed on triplicate samples) either in LM and RM (data not shown). In each case, statistical analysis of mass spectra of pure milk samples allowed to identify discriminant signals based on either intensity or mass values. Class prediction model was set up by GA; marker proteins or peptides were identified by means of their PWKW, PAD and and PTTA p-values. Corresponding data for milk speciation are shown in Suppl. Tab. S3 (PR-P - 23 markers), Suppl. Tab. S4 (PE-P in LM - 65 markers), Suppl. Tab. S5 (PE-P in RM - 106 markers), Suppl. Tab. S6 (PE-P in LM after C18ZipTip enrichment - 63 markers) and Suppl. Tab. S7 (PE-P in RM after C18ZipTip enrichment - 95 markers), together with additional statistical parameters. Corresponding data for the characterization of thermally-treated commercial bovine milks are shown in Suppl. Tab. S8 (PR-P - 10 markers), Suppl. Tab. S9 (PE-P in LM - 27 markers), Suppl. Tab. S10 (PE-P in RM - 47 markers), Suppl. Tab. S11 (PE-P in LM after C18ZipTip enrichment - 46 markers) and Suppl. Tab. S12 (PE-P in RM after C18ZipTip enrichment - 53 markers), together with additional statistical parameters. As expected, a number of common peptide markers occurred between analyses performed in LM and RM. In particular, 43, 39, 18 and 23 common peptide markers were observed during PE-P on milk filtrates for speciation, on C18ZipTip-enriched milk samples for speciation, on milk filtrates for TT characterization, on C18ZipTip-enriched milk samples for TT characterization, respectively. Importantly, data from PWKW, PAD and PTTA analysis defined the recognition capability of each specific assay; corresponding results are reported in Tab. 1 and 2. They demonstrate that milk speciation (100% recognition capability) was guaranteed in the case of MALDI-TOF-MS profiling measurements performed either on proteins from milk filtrates, on peptides from milk filtrates analyzed in RM, and on peptides from C18ZipTip-enriched milk filtrates analyzed both in LM and RM, notwithstanding the heterogeneous polypeptide composition of this biological fluid

16,18,22

(Tab. 1). Conversely, the

characterization of the TT of commercial bovine milks was ensured only in the case of peptides from C18ZipTipenriched milk filtrates analyzed in RM (Tab. 2). The latter result is probably dependent on the higher complexity of the protein and peptide mixtures present in commercial milk samples 5,15,17,19,20,23-26; it also underlines that the method

8   

ACS Paragon Plus Environment

Page 9 of 35

Journal of Agricultural and Food Chemistry

used for spectral acquisition can highly affect the recognition capability of the whole platform. On this basis, particular care is necessary for the initial set up of the PR-P and PE-P experimental conditions to guarantee the maximum number of signals present in the mass spectra and of the corresponding signal to noise ratios. These instrumental setting have to be maintained during the whole analysis of the different milk samples with the aim of not compromising the recognition capability of the whole platform. As expected, PCA of the data reported in Suppl. Tab. S3-12 were in line with the recognition capability values reported in Tab. 1 and 2. For example, Fig. 5A and 5B show the PCA 3D scatter plots of results from BM, WM, GM and OM proteins analyzed in LM and from BM, WM, GM and OM C18ZipTip-enriched peptides in RM, respectively, which were used for milk speciation. A good separation of the colored spots was evident in both cases. Similarly, Fig. 6A and 6B show the PCA 3D scatter plots of results from BM, PM, UM and IM proteins analyzed in LM, and from BM, PM, UM and IM C18ZipTip-enriched peptides in RM, respectively, which were used for the recognition of the material thermal treatment. Although more appreciable in a three-dimensional representation, well-resolved colored spots were evident also in these cases. Best dot separation occurred in the case of peptides. Less resolved PCA 3D scatter plots were obtained in the case of the remaining methods used for milk speciation (Suppl. Fig. S7A-C) and for the recognition of milk TT (Suppl. Fig. S8A-C). On the basis of the data reported in Tab. 1-2, Fig. 5-6, Suppl. Fig. S7-S8 A-C and the capability of this MS platform of performing independent analyses on different molecules (proteins and peptides), we conclude that complementary information associated with the highest recognition capability values for milk speciation and for thermal treatment characterization are obtained when diluted skimmed milk samples and C18ZipTip-enriched diluted skimmed milk filtrates are analyzed for the respective protein and peptide content in LM and RM, respectively.

Marker assignment Marker protein signals recognized for milk speciation (Suppl. Tab. S3 and Fig. 1) and for TT recognition (Suppl. Tab. S8 and Fig. 3) were associated with specific components based on literature data

9-14,17

and mass calculations. As

expected, most of the signals were assigned to -Lactalbumin, -Lactoglobulin, as well as to γ2/3-caseins and proteoso peptones resulting from the activity of milk plasmin on β-Casein

17,27-29

. In particular, γ2-caseins and 3-

caseins were identified as marker signals in BM, WM, GM and OM; -Lactalbumins were assigned markers in BM,

9   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 10 of 35

WM and GM; β-CN fragments (98-207) were discriminant signals in GM and OM; -Lactoglobulin, proteoso peptones p.p.8.I., and β-Casein fragments (1-68) and (69-209) were discriminant molecules in GM, BM and WM, respectively (Tab. 3). The latter complementary peptides originated after specific proteolysis at K68; this lysine is exclusively present within the water buffalo β-Casein sequence 4,17. When compared to data reported above, PR-P of thermally-treated milk samples identified a reduced number of protein markers, i.e. proteoso peptones p.p.8.I. (BM, PM and UM), -Lactoglobulin (BM and PM), -Lactalbumin adduct with lactose (m = +324) (IF), γ2-casein (UM) and its adduct with lactose (IM) (Tab. 4). Various non-assigned marker species were specific for UM and IM; the corresponding hypothetical lactose adducts were present only in IM. Finally, marker peptide signals from C18ZipTip-enriched diluted skimmed milk filtrates recognized for speciation (Suppl. Tab. S7 and Fig. 2) and for TT recognition (Suppl. Tab. S12 and Fig. 4) were associated with specific components based on MALDI-TOF/TOF and nanoLC-ESI-LIT-MS/MS analyses (see Supporting Information for experimental details). Fragmentation data were obtained from peptide mixtures either loaded on the MALDI plates also used for PE-P experiments and directly analyzed by nanoLC-ESI-LIT-MS/MS analysis; in the latter case, a specific data searching for peptides having a mass value coincident with that of the marker species was realized. This combined approach identified 38 marker peptides for milk speciation as deriving from -Casein (BM, WM, GM and OM), -S1-Casein (BM and OM), -S2-Casein (GM), Glycosylation-dependent cell adhesion molecule 1 (BM and WM), -Casein (WM, GM and OM), Cytoglobin isoform X1 (WM), Serum amyloid A protein (GM and OM), Butyrophilin subfamily 1 member A1 (GM), Lactoperoxidase (GM) and Cohesin subunit SA-3 (OM), among which some bearing specific modifications (Tab. 5). Similarly, 28 marker peptides for TT recognition were identified as deriving from -Casein (BM, PM, UM and IM), -S1-Casein (BM, PM, UM and IM), -S2-Casein (PM, UM and IM), Glycosylation-dependent cell adhesion molecule 1 (BM) and Polymeric immunoglobulin receptor (PM) (Tab. 6). Either in the case of milk speciation and thermal treatment recognition, the slight difference in the nature of the marker molecules identified with PR-P and PE-P experiments should be associated with the molecular mass value of intact proteins, which may exceed the instrument mass range we used, and the resistance of specific proteins to proteases that may hamper the production of the corresponding peptides.

10   

ACS Paragon Plus Environment

Page 11 of 35

Journal of Agricultural and Food Chemistry

Adulteration Platform capability to recognize milk adulteration due to the addition of either non-declared BM to WM, GM, and OM, or of IM to fresh BM was verified with PR-P and PE-P by analyzing nine impure milk samples, which were prepared in triplicate following the addition of different BM volumes into WM, GM and OM to generate 1, 2, 5, 10, 15, 20, 30, 40, and 50% v/v binary mixtures. Similarly, adulterated samples having identical mixed concentrations were prepared when IM was added to BM. After spectral processing according to the procedures already set up for pure materials, MALDI-TOF-MS data of these binary mixtures were subjected to PCA analysis. Resulting PCA data also included variance, PC score and coefficient values. These information were then used for PLS regression calculations. Fig. 7 shows the PLS regression plot of each binary mixture, as deriving from PR-P and PE-P analysis. An optimal linearity of the data was observed in each case. Normalized values for two adulterated test samples analyzed in triplicate were then used to get a prediction of the corresponding percentages of adulteration. Plots of the predicted adulteration levels, as estimated from the above-mentioned PLS model based on the corresponding PR-P and PE-P data, are reported in Fig. 8. Data from PR-P and PE-P analyses matched each other and were in good agreement with real adulteration percentages (5 and 15%) within the experimental error. These results demonstrated the good capability of this platform to predict illegal adulteration levels of unknown milk samples, when they are processed together with proper adulteration standards. They also demonstrated the importance of having analytical system that provides simultaneously two independent adulteration measurements realized on diluted skimmed milk samples and C18ZipTip-enriched diluted skimmed milk filtrates, which are analyzed for the respective protein and peptide content, respectively. As example, Fig. 9 shows the simultaneous detection of bovine protein and peptide markers in the PR-P and PE-P spectra of an adulterated GM sample containing 5% v/v BM.

Conclusions The development of a versatile and fast method for the detection of fraudulent adulterations in milk is a very important issue, due to the economic impact that these frauds have in dairy productions 6. As mentioned before, a number of analytical methods have been developed to this purpose, which have the disadvantage of being sometimes time-consuming, suitable only for a specific milk assay (speciation or TT analysis), useful for a unique animal species to animal species milk discrimination and/or eventually not adaptable for global screening of milk samples (Suppl.

11   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 12 of 35

Tab. S1). To overcome these limits and taking advantage of previous applications of the ClinProt technology for biomarker recognition and disease prediction in humans

30-31

, a MALDI-TOF-MS platform for combined PR-P and

PE-P of milk samples has been here developed, which recognizes markers either for speciation in BM, WM, GM and OM, as well as for thermal treatment in BM, PM, UM and IM, quickly assigning the nature of a pure unknown material. This platform was developed taking advantage of other applicative studies in the biomedical research field 32-34

. Platform capability to recognize milk adulteration due to the addition of small volumes of either non-declared

BM to WM, GM, and OM, or of IM to fresh BM was verified by analyzing different ad hoc prepared binary mixtures and specific test samples. In the latter case, concentration values were predicted according to independent PR-P and PE-P analyses and verified a posteriori to be coincident with real data. With respect to previous reports describing the use of PR-P analysis on diluted skimmed samples for milk speciation 9-14

, this study extends this method to the analysis of TT and introduces and independent, complementary

measurement, based on an ad hoc developed PE-P approach that has been set up to integrate protein data with additional information on other molecules, i.e. peptides, validating final results. When coupled to a robust statistical evaluation of the resulting data, this platform ultimately permits the typing of milk and the correlation of MALDI spectra with quantitative adulteration values. By introducing PE-P analysis in the sample processing pipeline, this integrated profiling approach also permits the detection of adulterations due to the addition of powdered material into fresh milk, broadening the applicative scenario of the MALDI profiling techniques in dairy science. The highthroughput character of the MALDI ionization technique and its versatility in answering to a number of issues related to the certification of raw materials used in dairy industry suggest the massive use of this platform in milk analysis with the aim to preserve product tradition and value.

Acknowledgments This work was supported by funds from the Italian Ministry of Economy and Finance to CNR for the project "Innovazione e sviluppo del Mezzogiorno - Conoscenze Integrate per Sostenibilità ed Innovazione del Made in Italy Agroalimentare" - Legge n.191/2009, and from Regione Campania for the Project “QUARC - Qualità delle produzioni tipiche campane ed il suo territorio: approcci innovativi ed integrati per rafforzare la competitività del sistema Agroalimentare” (A.G.C.6 - settore 1), bando CAMPUS.

12   

ACS Paragon Plus Environment

Page 13 of 35

Journal of Agricultural and Food Chemistry

References 1. Bertoni, G.; Calamari, L.; Maianti, M.G. Producing specific milks for speciality cheeses. Proc. Nutr. Soc. 2001, 60, 231-246 2. Johnson, M.E.; Lucey, J.A. Major technological advances and trends in cheese. J. Dairy Sci. 2006, 89, 11741178 3. Gerosa, S.; Skoet, J. Milk a vailability: trends in production and demand and medium-term outlook. FAO of the United Nations ESA Working paper 2012, 12-01 4. Di Luccia, A.; Picariello, G.; Trani, A.; Alviti, G.; Loizzo, P.; Faccia, M.; Addeo, F. Occurrence of betacasein fragments in cold-stored and curdled river buffalo (Bubalus bubalis L.) milk. J. Dairy Sci. 2009, 92, 1319-1329 5. Calvano, C.D.; Monopoli, A.; Loizzo, P, Faccia, M., Zambonin, C. Proteomic approach based on MALDITOF MS to detect powdered milk in fresh cow's milk. J. Agric. Food Chem. 2013, 61, 1609-1617 6. Commission Regulation (EU) No. 273/2008, Methods for the analysis and the evaluation of the quality milk and dairy products. Off. J. Eur. Commun. 2008, L88, 1-115 7. Addeo, F.; Garro, G.; Intorcia, N.; Pellegrino, L., Resmini, P., Chianese, L. Gel electrophoresis and immunoblotting for the detection of casein proteolysis in cheese. J. Dairy Res. 1995, 62, 297-309 8. Villar-Garea, A.; Griese, M.; Imhof, A. Biomarker discovery from body fluids using mass spectrometry. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 2007, 849, 105-114 9. Angeletti, R.; Gioacchini, A.M.; Seraglia, R.; Piro, R.; Traldi, P. The potential of matrix-assisted laser desorption/ionization mass spectrometry in the quality control of water buffalo mozzarella cheese. J. Mass Spectrom. 1998, 33, 525-553 10. Fanton, C.; Delogu, G.; Maccioni, E.; Podda, G.; Seraglia, R.; Traldi, P. Matrix-assisted laser desorption/ionization mass spectrometry in the dairy industry 2. The protein fingerprint of ewe cheese and its application to detection of adulteration by bovine milk. Rapid Commun. Mass Spectrom. 1998, 12, 1569– 1573 11. Cozzolino, R.; Passalacqua, S.; Salemi, S.; Malvagna, P.; Spina, E.; Garozzo, D. Identification of adulteration in milk by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. J. Mass Spectrom. 2001, 36, 1031–1037 12. Cozzolino, R.; Passalacqua, S.; Salemi, S.; Garozzo, D. Identification of adulteration in water buffalo mozzarella and in ewe cheese by using whey proteins as biomarkers and matrix-assisted laser desorption/ionization mass spectrometry. J. Mass Spectrom. 2002, 37, 985–991 13. Cunsolo, V.; Muccilli, V.; Saletti, R.; Foti, S. MALDI-TOF mass spectrometry for the monitoring of shedonkey’s milk contamination or adulteration. J. Mass Spectrom. 2013, 48, 148–153 14. Nicolaou, N.; Xu, Y.; Goodacre, R. MALDI-MS and multivariate analysis for the detection and quantification of different milk species. Anal. Bioanal. Chem. 2011, 399, 3491-3502 15. Calvano, C.D.; de Ceglie, C.; Monopoli, A.; Zambonin, C.G. Detection of sheep and goat milk adulterations by direct MALDI-TOF MS analysis of milk tryptic digests. J. Mass. Spectrom. 2012, 47, 1141–1149 16. Baum, F.; Fedorova, M.; Ebner, J.; Hoffmann, R.; Pischetsrieder, M. A sensitive and effective proteomic approach to identify she-donkey's and goat's milk adulterations by MALDI-TOF MS fingerprinting. J. Proteome Res. 2013, 12, 5447-5462 17. Somma, A.; Ferranti, P.; Addeo, F.; Mauriello, R.; Chianese, L. Peptidomic approach based on combined capillary isoelectric focusing and mass spectrometry for the characterization of the plasmin primary products from bovine and water buffalo beta-casein. J. Chromatogr. A. 2008, 1192, 294-300 18. Cuollo, M.; Caira, S.; Fierro, O.; Pinto, G.; Picariello, G.; Addeo, F. Toward milk speciation through the monitoring of casein proteotypic peptides. Rapid Commun. Mass Spectrom. 2010, 24, 1687-1696 19. De Simone, C.; Picariello, G.; Mamone, G.; Stiuso, P.; Dicitore, A.; Vanacore, D.; Chianese, L.; Addeo, F.; Ferranti, P. Characterisation and cytomodulatory properties of peptides from Mozzarella di Bufala Campana cheese whey. J. Pept. Sci. 2009, 15, 251-258 20. Holland, J.W.; Gupta, R., Deeth, H.C.; Alewood, P.F. Proteomic analysis of temperature-dependent changes in stored UHT milk. J. Agric. Food Chem. 2011, 59, 1837-1846 21. Mastrogiacomo, R.; D’Ambrosio, C.; Niccolini, A.; Serra, A.; Gazzano, A.; Scaloni, A.; Pelosi, P. An odorant-binding protein is abundantly expressed in the nose and in the seminal fluid of the rabbit. PLoS One 2014, 9, e111932

13   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 14 of 35

22. Dallas, D.C., Guerrero, A., Parker, E.A., Garay, L.A., Bhandari, A., Lebrilla, C.B., Barile, D., German, J.B. Peptidomic Profile of Milk of Holstein Cows at Peak Lactation. J. Agric. Food Chem. 2014, 62, 58-65 23. Arena, S.; Renzone, G.; Novi, G.; Paffetti, A.; Bernardini, G.; Santucci, A.; Scaloni, A. Modern proteomic methodologies for the characterization of lactosylation protein targets in milk. Proteomics 2010, 10, 34143434 24. Arena, S.; Renzone, G.; Novi, G.; Scaloni, A. Redox proteomics of fat globules unveils broad protein lactosylation and compositional changes in milk samples subjected to various technological procedures. J. Proteomics 2011, 74, 2453-2475 25. Renzone, G.; Arena, S.; Scaloni, A. Proteomic characterization of intermediate and advanced glycation endproducts in commercial milk samples. J. Proteomics 2015, 117, 12-23 26. Meltretter, J.; Wust, J.; Pischetsrieder, M. Modified peptides as indicators for thermal and nonthermal reactions in processed milk. J. Agric. Food Chem. 2014, 62, 10903-10915 27. Eigel, W.N. Identification of proteose-peptone component 5 as a plasmin-derived fragment of bovine betacasein. Int. J. Biochem. 1981, 13, 1081-1086 28. Andrews, A.T. Proteinases in normal bovine milk and their action on caseins. J. Dairy Res. 1983, 50, 45-55 29. Eigel, W.N.; Hofmann, C.J.; Chibber, B.A.K.; Tomik, J.M.; Keenen, T.W.; Mertz, E.T. Plasmin-mediated proteolysis of casein in bovine milk. Proc. Natl. Acad. Sci. USA 1979, 76, 2244-2248 30. Solassol, J.; Jacot, W.; Lhermitte, L.; Boulle, N.; Maudelonde, T.; Mangé, A. Clinical proteomics and mass spectrometry profiling for cancer detection. Expert Rev. Proteomics 2006, 3, 311-320 31. Magni, F.; Van Der Burgt, Y.E.; Chinello, C.; Mainini, V.; Gianazza, E.; Squeo, V.; Deelder, A.M.; Kienle, M.G. Biomarkers discovery by peptide and protein profiling in biological fluids based on functionalized magnetic beads purification and mass spectrometry. Blood Transfus. 2010, 8, Suppl 3:s92-97 32. Terracciano, R.; Preianò, M.; Palladino, G.P.; Carpagnano, G.E.; Barbaro, M.P.; Pelaia, G.; Savino, R.; Maselli, R. Peptidome profiling of induced sputum by mesoporous silica beads and MALDI-TOF MS for non-invasive biomarker discovery of chronic inflammatory lung diseases. Proteomics 2011, 11, 3402-3414 33. Preianò, M.; Pasqua, L.; Gallelli, L.; Galasso, O.; Gasparini, G., Savino, R.; Terracciano, R. Simultaneous extraction and rapid visualization of peptidomic and lipidomic body fluids fingerprints using mesoporous aluminosilicate and MALDI-TOF MS. Proteomics 2012, 12, 3286-3294 34. Preianò, M.; Falcone, D.; Maggisano, G.; Montalcini, T.; Navarra, M.; Paduano, S.; Savino, R.; Terracciano, R. Assessment of pre-analytical and analytical variables affecting peptidome profiling of gingival crevicular fluid by MALDI-TOF mass spectrometry. Clin. Chim. Acta 2014, 437, 120-128 35. Chianese, L.; Caira, S.; Lilla, S.; Pizzolongo, F.; Ferranti, P.; Pugliano, G.; Addeo, F. Primary structure of water buffalo alpha-lactalbumin variants A and B. J. Dairy Res. 2004, 71, 14-19

14   

ACS Paragon Plus Environment

Page 15 of 35

Journal of Agricultural and Food Chemistry

Table 1. ClinProt model generation report for the speciation of milk samples from different animal species (bovine, buffalo, goat and ovine).

  Milk speciation Molecular targets/Acquisition method

Cross validation (%)

Recognition capability (%)

Proteins from diluted samples/LM

95.11

100.00

Peptides from C18ZipTip-enriched diluted milk filtrates/RM

89.43

100.00

Peptides from diluted milk filtrates/RM

90.10

100.00

Peptides from C18ZipTip-enriched diluted milk filtrates/LM

89.11

100.00

Peptides from diluted milk filtrates/LM

77.47

97.39

 

15   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 16 of 35

Table 2. ClinProt model generation report for the recognition of thermally-treated commercial milk samples (raw, pasteurized, UHT and powdered).   Recognition of milk thermal treatment Molecular targets/Acquisition method

Cross validation (%)

Recognition capability (%)

Proteins from diluted samples/LM

87.82

97.19

Peptides from C18ZipTip-enriched diluted milk filtrates/RM

89.69

100.00

Peptides from diluted milk filtrates/RM

90.07

96.46

Peptides from C18ZipTip-enriched diluted milk filtrates/LM

86.94

98.04

Peptides from diluted milk filtrates/LM

82.44

98.00

16   

ACS Paragon Plus Environment

Page 17 of 35

Journal of Agricultural and Food Chemistry

Table 3. Marker species recognized during MALDI-TOF-MS protein profiling for milk speciation. MS experiments were performed in linear mode as reported in the experimental section. Difference value between the maximal and the minimal average peak/area intensity of all classes (DAve) (see ClinProt Tools software for details), protein/polypeptide name, and experimental/theoretical MH+ (average) value are reported. Theoretical MH+ values were calculated according to the amino acid sequence entries present in the UniProt Knowledge Database and the corresponding post-translational modifications. P, phosphate; f, fragment. Lactalbumin variant B mass value was calculated according to dedicated studies 35.

Bovine milk DAve

Description

Exp. MH+

Theor. MH+ (UniProt Code)

36.26

-Lactalbumin

14181

14179 (P00711)

29.46

2-Casein/-Casein f(106–209)

11827

11825 (P02666)

17.55

3-Casein/-Casein f(108–209)

11560

11560 (P02666)

15.97

Unknown

12006

10.29

Proteoso peptone p.p.8.I/-Casein f(29–105)

8640

8638 (P02666)

33.66

-Casein f(69–209)

15749

15750 (Q9TSI0)

30.90

Unknown

11991

30.39

3-Casein/-Casein f(108–209)

11521

11520 (Q9TSI0)

27.62

-Casein f(1–68)P5

8251

8252 (Q9TSI0)

20.56

2-Casein/-Casein f(106–209)

11786

11785 (Q9TSI0)

16.63

Lactalbumin variant B

14236

14236

36.26

-Lactalbumin variant A

14186

14187 (P00712)

29.02

Unknown

11882

21.33

3-Casein/-Casein f(108–207)

11288

16.48

Unknown

11922

15.17

Unknown

5891

14.06

2-Casein/-Casein f(106–207)

11554

11553 (Q95L76)

9.48

-Casein f(98–207)

12466

12466 (Q95L76)

8.66

-Lactoglobulin

18185

18187 (P02756)

11287 (Q95L76)

Ovine milk 36.25

2-Casein/-Casein f(106–207)

11536

11535 (P11839)

31.26

3-casein/-Casein f(108–207)

11271

11270 (P11839)

22.59

Unknown

8573

14.40

-Casein f(98–207)

12449

9.48

Unknown

9712

12448 (P11839)

17   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 18 of 35

Table 4. Marker species recognized during MALDI-TOF-MS protein profiling of thermally-treated commercial milk samples. MS experiments were performed in linear mode as reported in the experimental section. Difference value between the maximal and the minimal average peak/area intensity of all classes (DAve) (see ClinProt Tools software for details), protein/polypeptide name, and experimental/theoretical MH+ (average) value are reported. Theoretical MH+ values were calculated according to the bovine amino acid sequence entries present in the UniProt Knowledge Database and the corresponding post-translational modifications. The Amadori adduct is indicated with + Lactose; f, fragment.

Raw milk DAve

Description

Exp. MH+

Theor. MH+ (UniProt Code)

10.88

Proteoso peptone p.p.8.I/-Casein f(29–105)

8639

8638 (P02666)

4.46

-Lactoglobulin variant A

18362

18364 (P02754)

Pasteurized milk 10.88

Proteoso peptone p.p.8.I/-Casein f(29–105)

8639

8638 (P02666)

4.46

-Lactoglobulin variant A

18362

18364 (P02754)

UHT milk 38.27

Unknown

8007

24.64

Unknown

5326

10.88

Proteoso peptone p.p.8.I/-Casein f(29–105)

8639

8638 (P02666)

6.19

2-Casein + Lactose/-Casein f(106–209) + Lactose

12150

12149 (P02666)

Powdered milk 40.19

Unknown

6787

24.82

Unknown

6754

22.77

Unknown at m/z 6787 + Lactose

7111

16.40

Unknown at m/z 6754 + Lactose

7079

7.51

-Lactalbumin + Lactose

14502

14503 (P00711)

6.19

2-Casein/-Casein f(106–209) + Lactose

12150

12149 (P02666)

18   

ACS Paragon Plus Environment

Page 19 of 35

Journal of Agricultural and Food Chemistry

Table 5. Marker species recognized during MALDI-TOF-MS peptide profiling for milk speciation. MS experiments were performed in reflectron mode as reported in the experimental section. Most intense MH+ signals reported in Supplementary Tables S7 were subjected to fragmentation analysis for structural assignment (see Supporting Information for experimental details). Difference value between the maximal and the minimal average peak/area intensity of all classes (DAve) (see ClinProt Tools software for details), UniProt/NCBI sequence code (the latter highlighted in italics), protein name, experimental/theoretical MH+ (monoisotopic) value, modifications(s), mass error (in Da), peptide sequence, MASCOT/SEQUEST score (the latter highlighted in italics), and the mass spectrometry method used for peptide assignment (LIFT or CID) are reported. Theoretical MH+ values were calculated according to the indicated UniProt/NCBI sequence entries. N.d., not determined.

Bovine milk DAve

Accession code

Description

Exp. MH+

Theor. MH+

709.31

P02666

-Casein

1589.97

476.00

P02666

-Casein

1264.74

Unknown

2132.00

-S1-Casein

1250.67

Unknown

1347.78

57.52 56.66

P02662

56.16

Modification(s)

Error (Da)

Sequence

Score

LIFT

1589.93

-0.04

EPVLGPVRGPFPIIV

54



1264.70

-0.04

EPVLGPVRGPFP

55

N.d. 1250.60

-0.07

QKEDVPSERY

✓ ✓

37

N.d.

✓ ✓

44.51

P80195

Glycosylation-dependent cell adhesion molecule 1

2723.33

2723.36

0.03

ILNKPEDETHLEAQPTDASAQFIR

126

30.65

P02666

-Casein

1414.76

1414.71

-0.05

VPYPQRDMPIQA

47

Unknown

2528.43

Description

Exp. MH+

827.39

Unknown

2541.46

N.d.



564.57

Unknown

2145.91

N.d.



Glycosylation-dependent cell adhesion molecule 1

1693.92

Unknown

2555.35

-Casein

2026.95

Unknown

2379.29

Glycosylation-dependent cell adhesion molecule 1

2724.35

Unknown

1731.78

11.41

CID

N.d.

✓ ✓ ✓

Buffalo milk DAve

313.56

Accession code

594080212

168.81 119.65

P11840

97.33 95.92 56.19

594080212

Theor. MH+

Modification(s)

1693.89

Error (Da)

-0.03

SSRQPQNQNPKLPLS

Score

75

N.d. 2026.89

Q->pyroE (Nterm); Q->E (at Q3, Q5 or Q7)

-0.06

QEQNQEQPIRCEKEER

2724.30

-0.05

ILNEPEDETHLEAQPTDASAQFIR N.d.

ACS Paragon Plus Environment

LIFT

✓ ✓

49

N.d.

19   

Sequence

✓ ✓

135

✓ ✓

CID

Journal of Agricultural and Food Chemistry

Page 20 of 35

55.09

Unknown

2237.23

N.d.



52.43

Unknown

2291.17

N.d.



50.79

Q3HW31

-Casein

887.49

887.46

50.35

E1BMU2

Cytoglobin isoform X1

1989.83

1989.96

43.21

594080212

Glycosylation-dependent cell adhesion molecule 1

1606.89

25.18

594080212

Glycosylation-dependent cell adhesion molecule 1

1519.85

24.69

Unknown

2131.43

N.d.



23.15

Unknown

1603.67

N.d.



-0.03

YPFPGPIP

21



0.13

KVPGEMEIERRERSEE

21



1606.86

-0.03

SRQPQNQNPKLPLS

23



1519.83

-0.02

RQPQNQNPKLPLS

33



M6Ox

19.93

Q3HW31

-Casein

1314.77

1314.74

-0.03

SLVYPFPGPIPK

67

18.00

594080212

Glycosylation-dependent cell adhesion molecule 1

1319.73

1319.71

-0.02

RQPQNQNPKLP

24

17.13

Q9TSI0

-Casein

1227.74

1227.71

-0.03

LVYPFPGPIPK

41

10.62

Unknown

2907.30

N.d.



8.79

Unknown

2305.35

N.d.





✓ ✓ ✓

Goat milk DAve

Accession code

Description

Exp. MH+

Theor. MH+

Modification(s)

Error (Da)

Sequence

Score

993.43

P02670

-casein

1957.82

1957.79

Q1->pyroE

-0.03

QEQNQEQPICCEKDER

3.75

858.47

W5PJE0

Serum amyloid A protein

1899.05

1898.94

-0.11

SGKDPNHFRPAGLPDKY

39



576.54

A5JST2

Serum amyloid A protein

1834.90

1834.85

-0.05

QGWGTFLREAGQGAKDM

110



336.00

P33048

-Casein

1115.61

1115.56

-0.05

REQEELNVV

32



221.76

A5JST2

Serum amyloid A protein

1703.84

1703.81

Q1->pyroE

-0.03

QGWGTFLREAGQGAKD

92



171.67

P02670

-casein

2922.27

2922.30

Q1->pyroE

0.03

QEQNQEQPICCEKDERFFDDKIAK

4.15



131.71

P02670

-casein

2794.11

2794.20

Q1->pyroE

0.09

QEQNQEQPICCEKDERFFDDKIA

4.33



126.32

P02670

-casein

2251.83

2251.93

Q1->pyroE

0.10

QEQNQEQPICCEKDERFF

3.64



23.52

P02670

-casein

2610.15

22.48

A5JST2

Serum amyloid A protein

1388.69

17.91

Unknown

947.48

N.d.



16.79

Unknown

1048.54

N.d.



Q1->pyroE

N.d. 1388.66

-0.03

20   

ACS Paragon Plus Environment

AYRDMKEANYK

LIFT

CID ✓

✓ 48



Page 21 of 35

Journal of Agricultural and Food Chemistry

15.96

A3EY52

Butyrophilin member A1

15.35

A3F9D6

15.04 11.37

subfamily

1

1460.72

1461.62

0.90

MGEDSASGDIETLH

77



Lactoperoxidase

1344.76

1344.73

-0.03

GQVWEESLKRL

60



P33049

-S2-casein

1209.68

1209.66

-0.02

TNAIPYVRYL

A5JST2

Serum amyloid A protein

1850.82

1850.84

0.02

QGWGTFLREAGQGAKDM

Unknown

1189.62

11.10

Q1->pyroE; M17Ox

✓ 102

N.d.

✓ ✓

Ovine milk DAve

Accession code

Description

Exp. MH+

Theor. MH+

Modification(s)

Error (Da)

Sequence

810.55

P02669

-casein

2016.85

2016.84

Q1->pyroE

-0.01

QEQNQEQRICCEKDER

563.21

W5PJR5

Serum amyloid A protein

1761.85

1761.85

Q1->pyroE

0.00

QGLGTFLREAGQGAKDM

102



455.43

B0LRN2

Serum amyloid A protein

1625.88

1625.81

D->N D12)

-0.07

DPNHFRPAGLPDKY

60



302.85

P42819

Serum amyloid A

1652.80

1652.81

0.01

DPNHFRPPGLPDKY

57



174.48

P11839

-casein

1630.83

1630.79

-0.04

REQEELNVVGETVE

39

118.53

Unknown

2311.86

N.d.



98.44

Unknown

2853.27

N.d.



(D1

LIFT



W5Q6Z7

Cohesin subunit SA-3

2981.27

2981.24

71.12

B0LRN2

Serum amyloid A protein

953.49

46.23

P04653

-S1-casein

1529.90

35.47

Unknown

2104.08

N.d.



35.33

Unknown

2002.30

N.d.



25.50

Unknown

1616.68

N.d.



-S1-casein

1430.79

P04653

-0.03

QQELEELLQSSFLDEDEVYSLAATLK

27



953.45

-0.04

DPNHFRPA

40



1529.82

-0.08

SPEVLNENLLRFV

65



1430.76

-0.03

21   

ACS Paragon Plus Environment

SPEVLNENLLRF

CID ✓

89.56

23.68

Q1->pyroE

or

Score

46



Journal of Agricultural and Food Chemistry

Page 22 of 35

Table 6. Marker species recognized during MALDI-TOF-MS peptide profiling of thermally-treated commercial milk samples. MS experiments were performed in reflectron mode as reported in the experimental section. Most intense MH+ signals reported in Supplementary Tables S12 were subjected to fragmentation analysis for structural assignment (see Supporting Information for experimental details). Difference value between the maximal and the minimal average peak/area intensity of all classes (DAve) (see ClinProt Tools software for details), UniProt/NCBI sequence code (the latter highlighted in italics), protein name, experimental/theoretical MH+ (monoisotopic) value, modifications(s), mass error (in Da), peptide sequence, MASCOT/SEQUEST score (the latter highlighted in italics), and the mass spectrometry method used for peptide assignment (LIFT or CID) are reported. Theoretical MH+ values were calculated according to the indicated UniProt/NCBI sequence entries. N.d., not determined.

Raw milk DAve

Accession code

Description

Exp. MH+

Theor. MH+

144.16

P02666

-Casein

1589.97

1589.94

32.55

Unknown

1490.91

N.d.



18.21

Unknown

1979.13

N.d.



-Casein

1151.73

Unknown

2131.84

15.94

P02666

14.19

Modification(s)

1151.69

Error (Da)

Sequence

Score

LIFT

-0.03

EPVLGPVRGPFPIIV

54



-0.04

GPVRGPFPIIV

26

N.d.

CID

✓ ✓

11.05

P80195

Glycosylation-dependent cell adhesion molecule 1

2723.31

2723.36

0.05

ILNKPEDETHLEAQPTDASAQFIR

126

10.03

P02662

-S1-Casein

1250.65

1250.60

-0.05

QKEDVPSERY

37

Error (Da)

Sequence

Score

LIFT

CID

✓ ✓

Pasteurized milk DAve

Accession code

Description

Exp. MH+

Theor. MH+

129.77

P02662

-S1-Casein

1337.71

1337.68

-0.03

HIQKEDVPSER

23





26.06

P02663

-S2-Casein

923.57

923.53

-0.04

IPYVRYL

46





18.96

P81265

Polymeric receptor

768.52

769.42

0.90

AAPAGAAIQ

29



15.94

P02666

-Casein

1151.73

1151.69

-0.04

GPVRGPFPIIV

28



Error (Da)

Sequence

Score

immunoglobulin

Modification(s)

UHT milk DAve

Accession code

Description

Exp. MH+

Theor. MH+

262.34

P02666

-Casein

1284.68

1284.64

-0.04

RELEELNVPGE

59

120.64

P02666

-Casein

1994.16

1994.15

-0.01

LYQEPVLGPVRGPFPIIV

45

76.91

P02663

-S2-Casein

1251.66

1251.75

0.09

TKVIPYVRYL

29

Modification(s)

22   

ACS Paragon Plus Environment

LIFT

CID ✓



✓ ✓

Page 23 of 35

41.20

Journal of Agricultural and Food Chemistry

P02662

-S1-Casein

1689.83

1689.82

-0.01

PIGSENSEKTTMPLW

112



Error (Da)

Sequence

Score

LIFT

-0.03

HIQKEDVPSER

53



-0.03

NENLLRF

39



-0.05

RELEELNVPGE

58

Powdered milk DAve

Accession code

Description

Exp. MH+

Theor. MH+

129.77

P02662

-S1-Casein

1337.71

1337.68

97.04

P02662

-S1-Casein

905.51

905.48

73.78

P02666

72.10

-Casein

1283.70

Unknown

791.45

Modification(s)

1283.65

E->Q (E2, E4, E5 or E11) E->Q (E2, E4, E5, E11 or E14)

N.d.

63.23

P02666

-Casein

1624.87

1624.85

54.19

P02754

-Lactoglobulin

1246.70

1245.58

53.94

P02666

-Casein

1396.77

1396.74

33.77

P02666

-Casein

1068.61

28.44

P02666

-Casein

25.01

P02663

-S2-Casein



CID





-0.02

RELEELNVPGEIVE

51



-1.12

TPEVDDEALEK

77



-0.03

RELEELNVPGEI

30



1068.56

-0.05

YPVEPFTES

26



1311.70

1311.71

0.01

YQEPVLGPVRGP

58



1595.85

1594.71

S7Phospho

-1.14

KTVDMESTEVFTK

98



E->Q (E2, E4, E5 or E11)

-0.04

RELEELNVPGEIV

34



E->Q (E2, E4, E5 or E11)

21.07

P02666

-Casein

1495.84

1495.80

20.13

P02662

-S1-Casein

1099.64

1099.61

-0.03

EVLNENLLR

46



17.41

P02666

-Casein

788.43

788.41

-0.02

RELEEL

27



16.39

P02666

-Casein

1097.63

1097.59

-0.04

RELEELNVP

43



12.48

P02666

-Casein

1154.66

1154.61

-0.05

LTLTDVENLH

67



8.86

P02666

-Casein

1125.63

1125.61

-0.02

VVPPFLQPEV

43



5.69

P02662

-S1-Casein

2012.42

2014.00

1.57

IPNPIGSENSEKTTMPLW

46



E->Q (E2, E4 or E5) Q7->E

23   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 24 of 35

 

Legend to Figures

Fig. 1 MALDI-TOF-MS PR-P in LM of diluted milk samples. Spectra refer to BM (A), WM (B), GM (C) and OM (D). MH+ signals corresponding to identified marker molecules are highlighted by reporting the associated mass values.

Fig. 2 MALDI-TOF-MS PE-P in RM of diluted milk filtrates subjected to C18ZipTip enrichment as resulting from BM (A), WM (B), GM (C) and OM (D). MH+ signals corresponding to identified marker molecules are highlighted by reporting the associated mass values.

Fig. 3 MALDI-TOF-MS PR-P in LM of diluted thermally-treated commercial milk samples. Spectra refer to fresh BM (A), PM (B), UM (C) and IM (D). MH+ signals corresponding to identified marker molecules are highlighted by reporting the associated mass values.

Fig. 4 MALDI-TOF-MS PE-P in RM of milk filtrates subjected to C18ZipTip enrichment from thermally-treated commercial milk samples. Spectra refer to fresh BM (A), PM (B), UM (C) and IM (D). MH+ signals corresponding to identified marker molecules are highlighted by reporting the associated mass values.

Fig. 5 PCA scatter plot of protein (A) and peptide (B) MALDI-TOF-MS data in milk samples having different animal origin. Data referring to BM, WM, GM and OM are highlighted in red, green, blue and light blue, respectively. MALDI-TOF-MS PR-P experiments were performed in LM on diluted milk samples. MALDI-TOF-MS PE-P experiments were performed in RM on diluted milk filtrates subjected to C18ZipTip enrichment.

Fig. 6 PCA scatter plot of protein (A) and peptide (B) MALDI-TOF-MS data in milk samples subjected to different thermal treatment. Data referring to fresh BM, PM, UM and IM are highlighted in red, green, blue and light blue, respectively. MALDI-TOF-MS PR-P experiments were performed in LM on diluted milk samples.

24   

ACS Paragon Plus Environment

Page 25 of 35

Journal of Agricultural and Food Chemistry

MALDI-TOF-MS PE-P experiments were performed in RM on diluted milk filtrates subjected to C18ZipTip enrichment.

Fig. 7 PLS Plot of MALDI-TOF-MS data deriving from protein (A, C, E and G) and peptide (B, D, F and H) profiling experiments on adulterated milk samples resulting from the addition of BM to WM (A, B), GM (C, D), and OM (E, F), or of fresh BM to IM (G, H). MALDI-TOF-MS PR-P experiments were performed in LM on diluted milk samples. MALDI-TOF-MS PE-P experiments were performed in RM on diluted milk filtrates subjected to C18ZipTip enrichment. The different dots show the concentration of BM in the mixture in relation to the normalized PC scores.

Fig. 8 Plots showing the predicted levels of milk adulteration after MALDI-TOF-MS protein (A, C, E and G) and peptide (B, D, F and H) profiling experiments. The different panels represent data for adulterated samples resulting from the addition of BM to WM (A, B), GM (C, D), and OM (E, F). Data for adulterated samples resulting from the addition of fresh BM to IM are also shown (G, H). MALDI-TOF-MS PR-P experiments were performed in LM on diluted milk samples. MALDI-TOF-MS PE-P experiments were performed in RM on diluted milk filtrates subjected to C18ZipTip enrichment. Blue squares represent the training data set, while the red rhombus indicate the test set. The predicted levels of milk adulteration are reported ± S.D.

  Fig. 9 Cropped regions of the MALDI-TOF-MS PR-P and PE-P spectra of a GM sample containing 5% v/v BM. MALDI-TOF-MS PR-P was performed in LM on diluted milk samples. MALDI-TOF-MS PE-P was performed in RM on diluted milk filtrates subjected to C18ZipTip enrichment. MH+ signals corresponding to GM and BM marker molecules are highlighted by reporting the associated mass values; those associated with BM adulteration are marked with an asterisk.

 

25   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

 

Figure 1 26   

ACS Paragon Plus Environment

Page 26 of 35

Page 27 of 35

Journal of Agricultural and Food Chemistry

 

Figure 2 27   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 28 of 35

 

Figure 3 28   

ACS Paragon Plus Environment

Page 29 of 35

Journal of Agricultural and Food Chemistry

 

Figure 4

29   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 30 of 35

 

Figure 5

30   

ACS Paragon Plus Environment

Page 31 of 35

Journal of Agricultural and Food Chemistry

 

Figure 6

31   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 7 

32   

ACS Paragon Plus Environment

Page 32 of 35

Page 33 of 35

Journal of Agricultural and Food Chemistry

 

Figure 8   33   

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 34 of 35

   

Figure 9   34   

ACS Paragon Plus Environment

Page 35 of 35

Journal of Agricultural and Food Chemistry

TOC Graphics 266x151mm (96 x 96 DPI)

ACS Paragon Plus Environment