Rapid Quantification of Melamine in Different Brands/Types of Milk Powders Using Standard Addition Net Analyte Signal and Near-Infrared Spectroscopy

Abstract

Multivariate calibration (MVC) and near-infrared (NIR) spectroscopy have demonstrated potential for rapid analysis of melamine in various dairy products. However, the practical application of ordinary MVC can be largely restricted because the prediction of a new sample from an uncalibrated batch would be subject to a significant bias due to matrix effect. In this study, the feasibility of using NIR spectroscopy and the standard addition (SA) net analyte signal (NAS) method (SANAS) for rapid quantification of melamine in different brands/types of milk powders was investigated. In SANAS, the NAS vector of melamine in an unknown sample as well as in a series of samples added with melamine standards was calculated and then the Euclidean norms of series standards were used to build a straightforward univariate regression model. The analysis results of 10 different brands/types of milk powders with melamine levels 0~0.12% (w/w) indicate that SANAS obtained accurate results with the root mean squared error of prediction (RMSEP) values ranging from 0.0012 to 0.0029. An additional advantage of NAS is to visualize and control the possible unwanted variations during standard addition. The proposed method will provide a practically useful tool for rapid and nondestructive quantification of melamine in different brands/types of milk powders.

Bang-Cheng Tang 1 and Chen-Bo Cai 2 and Wei Shi 1 and Lu Xu 1,3

Academic Editor:Josep Esteve-Romero

1, Institute of Applied Chemistry, College of Material and Chemical Engineering, Tongren University, Tongren, Guizhou 554300, China

2, College of Chemistry and Life Science, Chuxiong Normal University, Chuxiong 675000, China

3, The Modernization Engineering Technology Research Center of Ethnic Minority Medicine of Hubei Province, College of Pharmacy, South-Central University for Nationalities, Wuhan 430074, China

Received 7 April 2016; Revised 31 May 2016; Accepted 7 June 2016

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Dairy products are essential components of a healthy diet for human and are very popular in all age groups owing to their high nutritional value and pleasurable flavor [1-3]. Recent years have seen greatly increasing production and consumption of dairy products in the world [4, 5]. An important factor that influences the textures and flavors of dairy products is the protein content, which has been adopted as a quality index by many industries. Unfortunately, because the traditional Kjeldahl method for analysis of total nitrogen content as an indication of protein levels is insufficient to distinguish between organic nitrogen and protein and nonprotein sources, unethical manufacturers deliberately added some illegal exogenous materials to dairy products to obtain an incorrectly higher readout of apparent protein content [6-8].

One of the most notorious exogenous adulterants used in dairy products is melamine, chemically known as 2,4,6-triamino-1,3,5-triazine [9-11]. Melamine is a nitrogen-rich (about 66.6%) heterocyclic triazine produced on a large scale (1.2 million tonnes in 2007) [12, 13]. It is primarily used in the synthesis of melamine formaldehyde resins for the production of paper finishers, flame retardant, commercial filters, moulding compounds, wrinkle-free textile, and many other materials [14-16]. As a very cheap and readily available industrial material, melamine was added to various food and food-related products, including milk, infant formula, frozen yogurt, biscuits, candy, coffee drinks, pet food, and feed [17-19], to increase the nitrogen level and to reduce costs. Although melamine has low toxicity by itself, when combined with cyanuric acid and uric acid, it can form insoluble crystals, which may lead to kidney stones, eventual renal failure, and ultimately death [20]. The most severe outbreak of melamine contamination occurred in 2008 in China, which had caused kidney stones in thousands of people and at least six deaths of young children [21].

To prevent further contamination and frauds, maximum limits have been established for melamine in infant formula and other foods by many countries. The emergent need for regulation of melamine has promoted extensive and intensive laboratory efforts to develop rapid, widely available, and cost-effective methods for analysis of melamine in various samples [1, 12, 18, 22-24], including capillary electrophoresis [25], high-performance liquid chromatography (HPLC) [26, 27], LC with mass spectrometry (LC-MS) [28, 29], gas chromatography with MS (GC-MS) [30, 31], micellar LC [32, 33] matrix-assisted laser desorption/ionization MS (MALDI-MS) [34], nuclear magnetic resonance spectroscopy [35], vibrational spectroscopy [1] and imaging [36, 37], chemiluminescence analysis [38, 39], electrochemical analysis [40], and immunoassay [41]. Among various methods, near-infrared (NIR) spectroscopy, as a rapid and nondestructive analytical method widely used in food analysis [42, 43], is promising for high-throughput screening and detection of melamine [44]. Although NIR has a comparatively lower sensitivity and higher detection limit compared with other methods, it was demonstrated to be sufficient in detecting excessive use of melamine in dairy products and could provide a convenient tool to rapidly screen and quantify melamine in Chinese markets [17, 45]. For quantitative analysis, chemometric methods for multivariate calibration (MVC), such as partial least squares (PLS), support vector machines (SVM), and artificial neural network (ANN), have also been successfully used to calibrate melamine levels in adulterated samples [1].

Although MVC combined with NIR spectroscopy has shown good accuracy and precision in analysis of melamine for some specific samples, the practical application of ordinary MVC to different brands/types of samples can be largely limited because the prediction of a new sample from an uncalibrated group would be subject to a significant bias due to matrix effect. This problem can be solved by performing calibration transfer [46] or standard addition [47]. While calibration transfer is a complicated procedure with high requirements of the practitioners' expertise and its performance can be influenced by many factors, standard addition method is a relatively easy-to-use tool for analysis of samples with complex compositions and from diverse origins.

Net analyte signal (NAS) theory [48] has been proven to be useful in the development of new multivariate calibration methods, evaluation of the figures of merit of multivariate calibration, variable selection, outlier diagnosis, and data preprocessing. NAS is based on an intuitive idea, namely, separating a part of the signal that is directly related to the concentration of the analyte of interest from that of interfering components. Using the Euclidean norm of the computed NAS vector, the multivariate signal can be represented as a univariate scalar. In this way, the multivariate calibration can be reduced to a simple univariate linear regression, which is especially useful for model validation and prediction. A multivariate standard addition method using NAS (SANAS) [49, 50] has been suggested and demonstrated to be effective in overcoming the matrix effect or indirect interferences.

The objective of this work was to study the feasibility of using multivariate standard addition method and NIR spectroscopy for rapid quantification of melamine in milk powders of different brands/types. To control the possible variations in preparation and measurement of added samples, the SANAS method was used to analyze different brands/types of milk powder samples.

2. Materials and Methods

2.1. Preparation of Melamine-Adulterated Milk Powder Samples and Standard Addition

Ten different brands/types of milk powder samples were collected from the quality inspection departments of producers as shown in Table 1. The shelf lives of all the samples were equal or over six months and, to the date of analysis, no sample has expired. The milk powder samples were kept in a dark, cool, and dry area at about 25°C with complete packaging before preparation and analysis. A melamine (Sinopharm Chemical Reagent Co., Ltd., Shanghai, China) gradient consisting of 5 levels, namely, 0.01, 0.02, 0.04, 0.08, and 0.12 percent (w/w%), together with the unadulterated samples, was prepared for each of the 10 batches. Therefore, 60 adulterated milk powder samples were obtained and used for analysis by SANAS. For standard addition, a designed gradient of 4 melamine standards was added to each level of the above 60 samples as listed in Table 2.

Table 1: The original milk powder samples of 10 different batches.

Number	Codes of brands/types	Types a	Production date
1	Q1	Skimmed, regular	May 3, 2015
2	Y1	Semiskimmed, high-calcium	Apr. 13, 2015
3	Q2	Skimmed, regular	May 22, 2015
4	M1	Semiskimmed, regular	Jun. 7, 2015
5	A1	Whole, sweet	May 19, 2015
6	Y2	Whole, sweet, high-calcium	Apr. 15, 2015
7	M2	Whole, high-calcium	Jun. 7, 2015
8	Y3	Skimmed, regular	May 12, 2015
9	M3	Skimmed, high-calcium	Jun. 9, 2015
10	Q3	Whole, regular	Apr. 27, 2015

^a The term "regular" means no deliberate addition of calcium or other elements to the milk; the term "high-calcium" denotes addition of calcium to the samples by the producers; the term "sweet" denotes addition of sucrose.

Table 2: Composition of serial samples spiked with different amounts of melamine standard.

Sample levels	Melamine content(%, w/w)	Gradient of melamine standard added(%, w/w)
1	0	0.01, 0.02, 0.04, and 0.08
2	0.01	0.01, 0.03, 0.07, and 0.11
3	0.02	0.02, 0.06, 0.10, and 0.14
4	0.04	0.04, 0.08, 0.12, and 0.20
5	0.08	0.04, 0.08, 0.16, and 0.24
6	0.12	0.04, 0.12, 0.20, and 0.36

2.2. NIR Spectral Measurements

The NIR diffuse reflectance spectra of milk powder samples were measured in the spectral range from 4000 to 10000 cm^-1 on an Antaris II Fourier transform-NIR spectrometer (Thermo Electron Co., Waltham, Massachusetts, USA) using the RESTLT 3.0 software. All samples were measured with a PbS detector and an internal gold background as the reference. The resolution was 8 cm^-1 and the scanning interval was 3.857 cm^-1 . Therefore, each spectrum had 1557 individual data points for chemometric analysis. For each object, 32 scans were performed and more scans did not enhance the spectral signals significantly.

2.3. Computation of Net Analyte Signal (NAS) Vector

For a detailed description of the net analyte signal (NAS) theory, one can refer to [48]. In this work, only a brief introduction is presented. For bilinear spectral data like NIR, the NAS of a chemical component is defined as the part of its pure spectrum that is orthogonal to the spectral space spanned by all the interfering components. The NAS vector of the k th component in a multicomponent mixture is defined as [figure omitted; refer to PDF] where _{R - k} is a matrix whose columns contain the pure spectrum of each component in the mixture except the k th component; _{s k} is the pure spectrum of the k th component; and I is a unit matrix and the superscript "+" represents the Moore-Penrose pseudo inverse of a matrix. The spectral variations caused by instrumental and environmental disturbances are also included in _{R - k} . Using the data matrix of calibration mixtures reconstructed by principal components analysis (PCA) or partial least squares (PLS), a rank annihilation procedure is adopted to compute _{R - k} : [figure omitted; refer to PDF] where _{R r e c} is the PCA- or PLS-reconstructed data matrix of interfering components using A significant components; _{c ^ k} is the analyte concentration vector _{c k} explained by the A significant components; and r is a vector including the spectrum of the k th analyte. Although r can be the pure spectrum of analyte k , _{s k} , or the spectrum of a mixture including analyte k , the pure spectrum _{s k} is the best choice because it contains maximal information on the k th analyte. α is a scalar factor which can be computed as [figure omitted; refer to PDF]

2.4. Standard Addition Method Using NAS

To do standard addition, a set of n standard solutions of the k th analyte were added to each of the samples. The concentration data of standards were collected in the vector _{c s} . The NIR data of the serial spiked samples were collected in the matrix, _{R s a} ( n × p ), where p is the number of NIR channels. The spectrum of a sample without standard addition is collected in a column vector, _{r u n} . By subtracting _{r u n} from each column of _{R s a} , one can obtain a matrix, _{R s} ( n × p ), which contains the spectra of the standards in the presence of the matrix effect. _{R s} is subject to PCA and the reconstructed _{R s} is used to rebuild _{R s a} and _{R s a , r e c} is obtained, and then the matrix of interfering components, _{R - k} , can be readily computed according to (2) and (3).

Subsequently, the NAS vectors of the k th analyte for the standard addition samples can be computed as [figure omitted; refer to PDF]

By plotting the Euclidean norm of the row vectors of N A _{S s a , k} ( n × p ) against _{c s} , the standard addition curve as for the univariate SA can be obtained to derive the concentration of k th analyte in the unknown sample.

2.5. Software

All the data analysis was performed using MATLAB 7.10.0 (R2010a) platform (MathWorks, USA). The data preprocessing and SANAS algorithms were performed based on in-house computational coded scripts written by authors in MATLAB.

3. Results and Discussion

The NIR spectra of the original milk powder samples and pure melamine as well as the melamine-adulterated samples are demonstrated in Figure 1. For ease of peak attribution, chemical bonds are denoted as atom-atom, where an atom can be carbon (C), hydrogen (H), oxygen (O), and nitrogen (N). For the spectra of milk powder, the peak around 4258 cm^-1 is the combination absorbance of C-H symmetric stretching and C-H bending, and those at 4335 cm^-1 can be attributed to the combination absorbance of C-H antisymmetric stretching and C-H bending. Other peak assignments are as follows: 4750 cm^-1 , combination of the basebands of N-H stretching and bending; 5157 cm^-1 , combination of the basebands of O-H stretching and bending; ~5700 cm^-1 , the first overtones of C-H stretching in various groups; ~6500 cm^-1 , the first overtone of N-H stretching; ~6900 cm^-1 , the first overtone of O-H stretching; and ~8300 cm^-1 , the second overtones of C-H stretching in various groups. By comparing the spectra of milk powder and melamine, the most significant difference is the intensive peak of melamine at 6811 cm^-1 , which could be attributed to the baseband of N-H antisymmetric stretching. Figure 2 demonstrates the second-order derivative (D2) spectra of the original and melamine-adulterated milk powder samples. However, because the added concentrations of melamine were very low (0.01~0.12%, w/w), no significant changes were found from the spectra of the adulterated milk powder by the naked eye. Multivariate calibration models were developed on the raw and D2 spectra using SANAS.

Figure 1: The NIR spectra of the original milk powder samples, pure melamine, and adulterated milk powder samples.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

Figure 2: The second-order derivative (D2) spectra of the original and melamine-adulterated milk powder samples.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

For each unknown sample, in order to compute the NAS vector, the spectrum of the original sample was subtracted from each of the melamine-spiked samples to obtain the matrix _{R s} . The resultant spectra in _{R s} for the sample Q1 in Table 1 are plotted in Figure 3. Seen in Figure 3, the high similarity between the raw spectra in _{R s} and that of pure melamine indicates that the computation of _{R s} using the method proposed in SANAS was very effective and helpful for accurate estimation of the NAS vectors. To rebuild _{R s} , _{R s a , r e c} was obtained using the primary principal components (PCs) of _{R s} . Because there were only 4 columns in _{R s} and its data structure was simple, the number of primary PCs was estimated according to an intuitive criterion that 95 percent of the total variances of _{R s} should be explained. The NAS vectors of the 4 standard-added samples were computed according to (2)-(4). Figure 4 demonstrates the computed NAS vectors of the 4 standard addition samples for the milk powder Q1.

Figure 3: The _{R s} matrix computed for milk powder sample Q1.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

Figure 4: The NAS vectors of milk powder sample Q1.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

Theoretically, the melamine level in an unknown sample could be estimated by plotting the Euclidean norm of the computed NAS vectors, ( N A _{S s a , k} ) , against the concentrations of added standards, _{c s} . However, because the number (60) of unknown samples to be analyzed was large, least squares regression (LSR) was performed between ( N A _{S s a , k} ) and _{c s} : [figure omitted; refer to PDF]

By setting the value of N A _{S s a , k} to be zero, the melamine level of the unknown samples can be computed as _{c u n k n o w n} = b / a . For all the 60 unknown samples, the correlation coefficients ( ^{r 2} ) of the above LSR were over 0.953, indicating that the computation of NAS vectors was very accurate. The prediction results of the 60 melamine-adulterated milk powder samples are summarized in Table 3. Seen in Table 3, the root mean squared error of prediction (RMSEP) of melamine ranged from 0.0012 to 0.0029 with different levels of melamine to be analyzed. The prediction accuracy of SANAS was not significantly influenced by the content of melamine. Moreover, the prediction performance by SANAS using the raw spectra and D2 spectra was similar although the results with D2 spectra seemed to be slightly better than those by the raw spectra, indicating that the unwanted spectral variations caused by baseline shifts were well controlled during standard addition. The predicted melamine levels were plotted against the reference values as shown in Figure 5, also indicating that the prediction errors were low and approximately uniform for different melamine levels to be analyzed.

Table 3: The prediction results of melamine levels in 10 different batches of milk powder samples by SANAS using the raw and D2 spectra.

Melamine level(w/w, %)	Raw spectra	D2 spectra
RSD a (%)	RMSEP b	RSD (%)	RMSEP
0	-	0.0019	-	0.0014
0.01	19.5	0.0019	17.3	0.0016
0.02	9.8	0.0019	9.2	0.0017
0.04	4.8	0.0018	3.1	0.0012
0.08	3.8	0.0029	2.6	0.0020
0.12	1.3	0.0014	1.2	0.0014

^a RSD: relative standard deviation.

^b RMSEP: root mean squared error of prediction.

Figure 5: The prediction results and reference values of melamine levels by SANAS.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

To further evaluate the figures of merit of the method, the selectivity, sensitivity in terms of limit of detection (LOD), linearity (Pearson's ^{r 2} ), and the accuracy and precision in terms of mean relative standard deviation (RSD) were computed and listed in Table 4. In China and the US, the maximum residue levels (MRL) for infant formula are 1.0 mg/kg and 2.5 mg/kg for milk and other dairy products, respectively [51]. Although the LOD (0.0025%) of this method was much worse than the regulation standards, this method provides the potential of rapid analysis and screening of the frauds in Chinese markets, where the practical melamine contents were much higher (up to 2563 mg/kg) [52].

Table 4: The figures of merit of melamine analysis in 10 different batches of milk powder samples by D2 spectra.

Selectivity	Sensitivity (LOD b , w/w%)	Linearity ( r 2 )	Accuracy (RSD a )	Precision (RSD)
0.614	0.0025	0.953	5.80%	3.74%

^a RSD: relative standard deviation.

^b LOD: limit of detection.

4. Conclusions

The feasibility of using NIR and SA for rapid quantification of melamine in different brands/types of milk powders was investigated. The analysis results for the 10 batches of melamine-adulterated milk powder samples demonstrate that SANAS is an effective method for SA multivariate calibration, which can visualize and control the spectral variations caused during SA in univariate regression. Moreover, the calibration accuracy was not significantly influenced by melamine levels to be analyzed. Compared with traditional multivariate calibration, combination of NIR and SANAS will provide a more practically applicable method for analysis of melamine in different brands/types of milk powder without requiring complex calibration transfer procedures.

Acknowledgments

The authors are grateful to the financial support from the Research Projects of Guizhou Science and Technology (no. QKHJZLKT[ [figure omitted; refer to PDF] ]15), the Open Research Program (no. GCTKF2014007) of State Key Laboratory Breeding Base of Green Chemistry Synthesis Technology (Zhejiang University of Technology), the Research Fund for the Doctoral Program of Tongren University (no. trxyDH1501), the Open Research Program (no. 2015ZY006) from the Modernization Engineering Technology Research Center of Ethnic Minority Medicine of Hubei province (South-Central University for Nationalities), and the research funds from the Education Department of Guizhou Province (no. QJHKYZ [ [figure omitted; refer to PDF] ]498).

References

[1] E. Domingo, A. A. Tirelli, C. A. Nunes, M. C. Guerreiro, S. M. Pinto, "Melamine detection in milk using vibrational spectroscopy and chemometrics analysis: a review," Food Research International , vol. 60, pp. 131-139, 2014.

[2] A. Hilding-Ohlsson, J. A. Fauerbach, N. J. Sacco, M. C. Bonetto, E. Corton, "Voltamperometric discrimination of urea and melamine adulterated skimmed milk powder," Sensors , vol. 12, no. 9, pp. 12220-12234, 2012.

[3] G. B. Sanvido, J. S. Garcia, Y. E. Corilo, B. G. Vaz, J. J. Zacca, R. G. Cosso, M. N. Eberlin, M. G. Peter, "Fast screening and secure confirmation of milk powder adulteration with maltodextrin via electrospray ionization-mass spectrometry [ESI(+)-MS] and selective enzymatic hydrolysis," Journal of Agricultural and Food Chemistry , vol. 58, no. 17, pp. 9407-9412, 2010.

[4] E. Lu, "Radical shifts in China's milk market," 2009, http://www.china.org.cn/business/news/2009-01/14/content_17105973.htm

[5] X. Lu, H. Zong, "The problems and countermeasures after China's dairy enters the adjustment period," Agricultural Economic Problems , vol. 7, pp. 5-50, 2008.

[6] Z. C. Y. Chan, W.-F. Lai, "Revisiting the melamine contamination event in China: implications for ethics in food technology," Trends in Food Science & Technology , vol. 20, no. 8, pp. 366-373, 2009.

[7] X. Pei, A. Tandon, A. Alldrick, L. Giorgi, W. Huang, R. Yang, "The China melamine milk scandal and its implications for food safety regulation," Food Policy , vol. 36, no. 3, pp. 412-420, 2011.

[8] C. Xiu, K. K. Klein, "Melamine in milk products in China: examining the factors that led to deliberate use of the contaminant," Food Policy , vol. 35, no. 5, pp. 463-470, 2010.

[9] H. Miao, S. Fan, Y.-N. Wu, L. Zhang, P.-P. Zhou, J.-G. Li, H.-J. Chen, Y.-F. Zhao, "Simultaneous determination of melamine, ammelide, ammeline, and cyanuric acid in milk and milk products by gas chromatography-tandem mass spectrometry," Biomedical and Environmental Sciences , vol. 22, no. 2, pp. 87-94, 2009.

[10] H. Xin, R. Stone, "Tainted milk scandal. Chinese probe unmasks high-tech adulteration with melamine," Science , vol. 322, no. 5906, pp. 1310-1311, 2008.

[11] N. Yan, L. Zhou, Z. Zhu, X. Chen, "Determination of melamine in dairy products, fish feed, and fish by capillary zone electrophoresis with diode array detection," Journal of Agricultural and Food Chemistry , vol. 57, no. 3, pp. 807-811, 2009.

[12] P. Lutter, M.-C. Savoy-Perroud, E. Campos-Gimenez, L. Meyer, T. Goldmann, M.-C. Bertholet, P. Mottier, A. Desmarchelier, F. Monard, C. Perrin, F. Robert, T. Delatour, "Screening and confirmatory methods for the determination of melamine in cow's milk and milk-based powdered infant formula: validation and proficiency-tests of ELISA, HPLC-UV, GC-MS and LC-MS/MS," Food Control , vol. 22, no. 6, pp. 903-913, 2011.

[13] World Health Organization, Food and Agriculture Organization of the United Nations Expert Meeting to Review Toxicological Aspects of Melamine and Cyanuric Acid , World Health Organization, Geneva, Switzerland, 2008.

[14] K. H. Lund, J. H. Petersen, "Migration of formaldehyde and melamine monomers from kitchen- and tableware made of melamine plastic," Food Additives and Contaminants , vol. 23, no. 9, pp. 948-955, 2006.

[15] R. E. Martin, C. B. Hizo, A. M. Ong, O. M. Alba, H. Ishiwata, "Release of formaldehyde and melamine from tableware manufactured in the Philippines," Journal of Food Protection , vol. 55, no. 8, pp. 632-635, 1992.

[16] L. C. Mecker, K. M. Tyner, J. F. Kauffman, S. Arzhantsev, D. J. Mans, C. M. Gryniewicz-Ruzicka, "Selective melamine detection in multiple sample matrices with a portable Raman instrument using surface enhanced Raman spectroscopy-active gold nanoparticles," Analytica Chimica Acta , vol. 733, pp. 48-55, 2012.

[17] L. J. Mauer, A. A. Chernyshova, A. Hiatt, A. Deering, R. Davis, "Melamine detection in infant formula powder using near- and mid-infrared spectroscopy," Journal of Agricultural and Food Chemistry , vol. 57, no. 10, pp. 3974-3980, 2009.

[18] F. Sun, W. Ma, L. Xu, Y. Zhu, L. Liu, C. Peng, L. Wang, H. Kuang, C. Xu, "Analytical methods and recent developments in the detection of melamine," Trends in Analytical Chemistry , vol. 29, no. 11, pp. 1239-1249, 2010.

[19] S. Yang, J. Ding, J. Zheng, B. Hu, J. Li, H. Chen, Z. Zhou, X. Qiao, "Detection of melamine in milk products by surface desorption atmospheric pressure chemical ionization mass spectrometry," Analytical Chemistry , vol. 81, no. 7, pp. 2426-2436, 2009.

[20] C. A. Brown, K.-S. Jeong, R. H. Poppenga, B. Puschner, D. M. Miller, A. E. Ellis, K.-I. Kang, S. Sum, A. M. Cistola, S. A. Brown, "Outbreaks of renal failure associated with melamine and cyanuric acid in dogs and cats in 2004 and 2007," Journal of Veterinary Diagnostic Investigation , vol. 19, no. 5, pp. 525-531, 2007.

[21] E. Chan, S. Griffiths, C. Chan, "Public-health risks of melamine in milk products," The Lancet , vol. 372, no. 9648, pp. 1444-1445, 2008.

[22] Y. Liu, E. E. D. Todd, Q. Zhang, J.-R. Shi, X.-J. Liu, "Recent developments in the detection of melamine," Journal of Zhejiang University SCIENCE B , vol. 13, no. 7, pp. 525-532, 2012.

[23] K. Rovina, S. Siddiquee, "A review of recent advances in melamine detection techniques," Journal of Food Composition and Analysis , vol. 43, pp. 25-38, 2015.

[24] Y. Wu, Y. Zhang, "Analytical chemistry, toxicology, epidemiology and health impact assessment of melamine in infant formula: recent progress and developments," Food and Chemical Toxicology , vol. 56, no. 2, pp. 325-335, 2013.

[25] J. Xia, N. Zhou, Y. Liu, B. Chen, Y. Wu, S. Yao, "Simultaneous determination of melamine and related compounds by capillary zone electrophoresis," Food Control , vol. 21, no. 6, pp. 912-918, 2010.

[26] X.-T. Peng, Z.-G. Shi, Y.-Q. Feng, "Rapid and high-throughput determination of melamine in milk products and eggs by full automatic on-line polymer monolith microextraction coupled to high-performance liquid chromatography," Food Analytical Methods , vol. 4, no. 3, pp. 381-388, 2011.

[27] N. Poorjafari, A. Zamani, M. Mohseni, A. Parizanganeh, "Assessment of residue melamine in dairy products exhibited in Zanjan market, Iran by high-performance liquid chromatography method," International Journal of Environmental Science & Technology , vol. 12, no. 3, pp. 1003-1010, 2015.

[28] D. He, X. Zhang, B. Gao, L. Wang, Q. Zhao, H. Chen, H. Wang, C. Zhao, "Preparation of magnetic molecularly imprinted polymer for the extraction of melamine from milk followed by liquid chromatography-tandem mass spectrometry," Food Control , vol. 36, no. 1, pp. 36-41, 2014.

[29] Z. Meng, Z. Shi, S. Liang, X. Dong, Y. Lv, H. Sun, "Rapid screening and quantification of cyromazine, melamine, ammelide, ammeline, cyanuric acid, and dicyandiamide in infant formula by ultra-performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry and triple quadrupole mass spectrometry," Food Control , vol. 55, pp. 158-165, 2015.

[30] Y. Fan, X. Ma, Z. Li, M. Chen, "Fast derivatization followed by gas chromatography-mass spectrometry for simultaneous detection of melamine, ammeline, ammelide, and cyanuric acid in fish and shrimp," Food Analytical Methods , vol. 9, no. 1, pp. 16-22, 2016.

[31] Y. Zhang, X. Ma, Y. Fan, "A rapid and sensitive method for determination of melamine in fish, shrimp, clam, and winkle by gas chromatography-mass spectrometry with microwave-assisted derivatization," Food Analytical Methods , vol. 7, no. 9, pp. 1763-1769, 2014.

[32] B. Beltran-Martinavarro, J. Peris-Vicente, S. Carda-Broch, J. Esteve-Romero, "Development and validation of a micellar liquid chromatography-based method to quantify melamine in swine kidney," Food Control , vol. 46, pp. 168-173, 2014.

[33] J. Peris-Vicente, J. Albiol-Chiva, P. Roca-Genoves, J. Esteve-Romero, "Advances on melamine determination by micellar liquid chromatography: a review," Journal of Liquid Chromatography & Related Technologies , vol. 39, no. 7, pp. 325-338, 2016.

[34] H.-W. Tang, K.-M. Ng, S. S.-Y. Chui, C.-M. Che, C.-W. Lam, K.-Y. Yuen, T.-S. Siu, L. C.-L. Lan, X. Che, "Analysis of melamine cyanurate in urine using matrix-assisted laser desorption/ionization mass spectrometry," Analytical Chemistry , vol. 81, no. 9, pp. 3676-3682, 2009.

[35] B. Tomita, C. Y. Hse, "Analysis of cocondensation of melamine and urea through formaldehyde with carbon 13 nuclear magnetic resonance spectroscopy," Japan Wood Research Society , vol. 41, no. 3, pp. 349-354, 1995.

[36] X. Fu, M. S. Kim, K. Chao, J. Qin, J. Lim, H. Lee, A. Garrido-Varo, D. Perez-Marin, Y. Ying, "Detection of melamine in milk powders based on NIR hyperspectral imaging and spectral similarity analyses," Journal of Food Engineering , vol. 124, no. 1, pp. 97-104, 2014.

[37] Y. Huang, K. Tian, S. Min, Y. Xiong, G. Du, "Distribution assessment and quantification of counterfeit melamine in powdered milk by NIR imaging methods," Food Chemistry , vol. 177, pp. 174-181, 2015.

[38] V. De Lourdes Mendes Finete, M. M. Gouvêa, F. F. De Carvalho Marques, A. D. P. Netto, "Is it possible to screen for milk or whey protein adulteration with melamine, urea and ammonium sulphate, combining Kjeldahl and classical spectrophotometric methods?," Food Chemistry , vol. 141, no. 4, pp. 3649-3655, 2013.

[39] J. Zhang, M. Wu, D. Chen, Z. Song, "Ultrasensitive determination of melamine in milk products and biological fluids by luminol-hydrogen peroxide chemiluminescence," Journal of Food Composition and Analysis , vol. 24, no. 7, pp. 1038-1042, 2011.

[40] K. Rovina, S. Siddiquee, "Electrochemical sensor based rapid determination of melamine using ionic liquid/zinc oxide nanoparticles/chitosan/gold electrode," Food Control , vol. 59, pp. 801-808, 2016.

[41] Y. Gong, M. Zhang, M. Wang, Z. Chen, X. Xi, "Development of immuno-based methods for detection of melamine," Arabian Journal for Science and Engineering , vol. 39, no. 7, pp. 5315-5324, 2014.

[42] M. Garcia-Alvarez, J. F. Huidobro, M. Hermida, J. L. Rodriguez-Otero, "Major components of honey analysis by near-infrared transflectance spectroscopy," Journal of Agricultural and Food Chemistry , vol. 48, no. 11, pp. 5154-5158, 2000.

[43] N. Sinelli, L. Cerretani, V. D. Egidio, A. Bendini, E. Casiraghi, "Application of near (NIR) infrared and mid (MIR) infrared spectroscopy as a rapid tool to classify extra virgin olive oil on the basis of fruity attribute intensity," Food Research International , vol. 43, no. 1, pp. 369-375, 2010.

[44] S. A. Haughey, S. F. Graham, E. Cancouët, C. T. Elliott, "The application of Near-Infrared Reflectance Spectroscopy (NIRS) to detect melamine adulteration of soya bean meal," Food Chemistry , vol. 136, no. 3-4, pp. 1557-1561, 2013.

[45] A. K. Ghazi-Tehrani, H. N. Pontell, "Corporate crime and state legitimacy: the 2008 Chinese melamine milk scandal," Crime, Law and Social Change , vol. 63, no. 5, pp. 247-267, 2015.

[46] M. Kompany-Zareh, F. van den Berg, "Multi-way based calibration transfer between two Raman spectrometers," Analyst , vol. 135, no. 6, pp. 1382-1388, 2010.

[47] H. Zou, K. Xu, Y. Feng, B. Liang, "Application of first order derivative UV spectrophotometry coupled with H-point standard addition to the simultaneous determination of melamine and dicyandiamide in milk," Food Analytical Methods , vol. 8, no. 3, pp. 740-748, 2015.

[48] R. Hajian, E. Mousavi, N. Shams, "Net analyte signal standard addition method for simultaneous determination of sulphadiazine and trimethoprim in bovine milk and veterinary medicines," Food Chemistry , vol. 138, no. 2-3, pp. 745-749, 2013.

[49] K. Asadpour-Zeynali, P. Najafi-Marandi, "Bismuth modified disposable pencil-lead electrode for simultaneous determination of 2-nitrophenol and 4-nitrophenol by net analyte signal standard addition method," Electroanalysis , vol. 23, no. 9, pp. 2241-2247, 2011.

[50] B. Hemmateenejad, S. Yousefinejad, "Multivariate standard addition method solved by net analyte signal calculation and rank annihilation factor analysis," Analytical and Bioanalytical Chemistry , vol. 394, no. 7, pp. 1965-1975, 2009.

[51] World Health Organization Toxicological and Health Aspects of Melamine and Cyanuric Acid , WHO, 2009.

[52] General Administration of Quality Supervision, Inspection and Quarantine the People's Republic of China (AQSIQ), Report on the results of special inspection of infant formula milk powder in China, 2008