Validation of the Analytical Method
Method validation was performed by evaluating the following analytical performance parameters: selectivity, matrix-effect, linearity, trueness (recovery), precision (intra-day and inter-day repeatability), limit of detection (LOD) and limit of quantification (LOQ).
The selectivity of the method was evaluated by analyzing blank matrix samples of whole powdered milk and different infant formulas types. For all the evaluated matrices, the absence of interference signals eluting at the same retention time of the selected mycotoxins was verified. Figure 1 displays a chromatogram of a follow-up formula spiked with the target compounds at 0.15 µg kg-1 (concentration of infant formulas reconstituted with water). Mycotoxin identification was performed by comparing the analyte retention time in the samples with those of the standard solutions. Confirmation was performed by comparing the signal intensity ratios of the two ion transitions of each analyte in the sample with those of the standard solutions (Table 1).
To investigate matrix effects, whole powdered milk and infant formulas from different producers available in the local market (starter milk-based infant formula, n = 4; follow-up milk-based infant formula, n = 3; lactose-free infant formula, n = 2) were selected. Calibration curves for the target compounds in matrix extracts (matrix-matched calibration) and in methanol/water (1:1 v/v) were prepared at four concentration levels, ranging from 0.1 to 1.5 ng mL-1; and each solution was then analyzed in triplicate. The calibration curve slopes were compared by an analysis of covariance (ANCOVA), considering a significance level of 5% (García-Campaña et al. 1997). Prior to performing the ANCOVA, the homogeneity of the residual variances of all the calibration curves was confirmed (p-values > 0.05) using a modified Levene test (Brown and Forsythe 1974). A strong suppression of the analytical signal was observed for the target mycotoxins in soy-based infant formula matrix. Therefore, the analytical method was considered inadequate for the determination of target mycotoxins in soy-based infant formula samples.
For ochratoxin A and aflatoxin M1, no significant differences between the slopes of the calibration curves for all evaluated matrices and the slope of the calibration curve in solvent was observed (p-values > 0.05). An absence of matrix effects for these compounds in all studied matrices was, thus, observed. However, significant differences between the slopes of the calibration curves prepared in solvent and matrices for at least two of the target aflatoxins in each studied matrix were observed (p-values < 0.05). Matrix effects for the target compounds in different matrices are displayed in Table 2, where positive values indicate signal enhancement and negative values indicate signal suppression.
Table 2
Matrix effect evaluation for infant formulaS and powdered milk
| | Milk-based infant formulas (started) n = 4 | | Milk-based infant formulas (follow-up) n = 3 | | Milk-based infant formulas (lactose free) n = 2 | | Powdered Milk (whole) n = 1 | | |
Mycotoxins | | a(p-value) | b(p-value) | Matrix Effect (%) | | a(p-value) | b(p-value) | Matrix Effect (%) | | a(p-value) | b(p-value) | Matrix effect (%) | | b(p-value) | Matrix effect (%) | | c(p-value) |
Aflatoxin M2 | | 0.07 | 0.00 | 39 | | 0.09 | 0.00 | 75 | | 3.03 | 0.00 | 29 | | 0.00 | 74 | | 0.00 |
Aflatoxin M1 | | 0.72 | 0.07 | 8 | | 0.84 | 0.29 | 11 | | 2.03 | 0.26 | 6 | | 0.22 | 4 | | 0.89 |
Aflatoxin G2 | | 0.32 | 0.12 | -6 | | 0.73 | 0.08 | 9 | | 3.85 | 0.98 | 0 | | 0.01 | -12 | | 0.00 |
Aflatoxin G1 | | 0.07 | 0.00 | -17 | | 0.18 | 0.07 | -5 | | 0.19 | 0.02 | -8 | | 0.00 | -29 | | 0.00 |
Aflatoxin B2 | | 0.10 | 0.00 | -15 | | 0.79 | 0.38 | -5 | | 0.25 | 0.00 | -11 | | 0.00 | -24 | | 0.00 |
Aflatoxin B1 | | 0.16 | 0.00 | -11 | | 0.06 | 0.00 | -10 | | 0.13 | 0.01 | -6 | | 0.00 | -19 | | 0.00 |
Ochratoxin A | | 0.18 | 0.20 | 6 | | 0.33 | 0.51 | -2 | | 1.34 | 0.25 | 3 | | 0.71 | 5 | | 0.72 |
aObtained by analysis of covariance comparing the slopes of the calibration curves prepared for each matrix (different producers); b obtained by analysis of covariance comparing the slopes of the matched calibration matrix with the slopes of the calibration curve in solvent; c obtained by analysis of covariance comparing all the slopes of the calibration curves prepared in the matrix. At p-Value > 0.05, there is no significant difference between the slopes of the calibration curves, evaluated with a level of confidence of 95% |
The linearity of the calibration curves used in routine analyses was evaluated using matrix-matched calibration curves at four concentration levels, between 0.1 and 1.5 ng mL-1 and the solutions were analyzed in triplicate. The ordinary least squares method was applied for the elaboration of the calibration curves. The resulting correlation coefficients were always higher than 0.99. Linearity was assessed according to the procedure described by Souza and Junqueira (2005). As assumptions for regression analysis, the homocedasticity, independency and normality of the regression residuals were verified. Initially, the outliers were successively investigated by the Jacknife standardized residuals test (Belsley et al. 1980), until non-detection or a maximum exclusion of 22.2% in the number of original results (Horwitz 1995). The homocedasticity of residuals, verified by a modified Levene test (Brown and Forsythe 1974), was confirmed (p-values > 0.05). The independency of residuals, verified by the Durbin-Watson statistic (Durbin and Watson 1951), was confirmed (p-values > 0.05). The normality of residuals, verified by Ryan-Joiner test (Ryan and Joiner 1976), was also confirmed (p-values > 0.05). The regression significance and the lack-of-fit were performed by an analysis of variance (ANOVA) (Draper and Smith 1998). A high regression significance (p-values < 0.001) and non-significant lack-of-fit (p-values > 0.05) were observed, attesting the linearity of the evaluated calibration curves.
Trueness and intra-day precision were evaluated by recovery studies using follow-up infant formula spiked with the investigated mycotoxins at three concentration levels (0.05, 0.1 and 0.15 µg kg-1), with four replicates for each level. The resulting values are displayed in Table 3, where the precision is expressed by the relative standard deviation (RSD) and trueness by the recovery values. The results were evaluated according to the European Commission Decision 2002/657/EC (EC 2002). The recovery values ranged from 75 to 111%, with RSD lower than 12% for all investigated mycotoxins. The RSD for intermediate precision (inter day) was also lower than 12%. Hence, the method showed good trueness and precision for the evaluated concentration levels.
Table 3
validation Parameters of the method for infant formulas
| | | | 0.05 µg kg− 1 | | 0.1 µg kg− 1 | | 0.15 µg kg− 1 |
Mycotoxins | LOD | LOQ | | Rec | aRSD | bRSD | | Rec | aRSD | | Rec | aRSD |
Aflatoxin M2 | 0.003 | 0.011 | | 89.0 | 3.1 | 4.6 | | 87.3 | 2.0 | | 95.4 | 7.4 |
Aflatoxin M1 | 0.005 | 0.016 | | 88.5 | 8.2 | 7.9 | | 110.6 | 3.9 | | 99.1 | 10.1 |
Aflatoxin G2 | 0.008 | 0.028 | | 103.1 | 5.0 | 13.5 | | 85.3 | 7.3 | | 88.0 | 6.9 |
Aflatoxin G1 | 0.002 | 0.008 | | 83.5 | 7.4 | 10.8 | | 79.6 | 2.5 | | 79.8 | 7.5 |
Aflatoxin B2 | 0.004 | 0.012 | | 92.1 | 5.9 | 4.88 | | 87.6 | 4.9 | | 91.4 | 4.9 |
Aflatoxin B1 | 0.004 | 0.013 | | 104.0 | 6.2 | 12.1 | | 93.5 | 2.9 | | 98.3 | 2.3 |
Ochratoxin A | 0.018 | 0.059 | | 75.0 | 12.2 | 10.3 | | 107.5 | 4.3 | | 97.3 | 11.7 |
LOD, Limit of detection (µg kg − 1); LOQ, Limit of Quantification (µg kg − 1); Concentration of reconstituted infant formulas; Rec (%), recovery; aRSD (%) relative standard deviation (intra-day, n = 4); bRSD (%) relative standard deviation (inter-day, n = 4 days). |
Infant formulas spiked with all the target mycotoxins at 0.05 µg kg-1 were used to calculate the limits of detection (LOD) and limits of quantification (LOQ), considering signal-to-noise ratios of 3 and 10, respectively. The summary results are displayed in Table 3. The sensitivity of the validated method was considered adequate for routine analysis. The LOQ of the method for aflatoxins B1, B2, G1 and G2 were much lower than the maximum permissible Brazilian limit (1 µg kg-1) for total aflatoxins (sum of the aflatoxins B1, B2, G1 and G2) in infant formulas. The LOQ for aflatoxin M1 (0.016 µg kg-1) was lower than the restriction limit set by the EU regulation (0.025 µg kg-1) for infant formulas (EC 2006) and comparable to the 0.010 and 0.0125 µg kg-1 values reported respectively by Zhang et al. (2013) and Tonon et al. (2018), who also employed the LC-MS/MS technique.
Sample Analyses
In order to protect public health, regulations for the control of aflatoxin M1 contamination levels in milk and infant formula have been established in several countries (FAO 2004; EC 2006). The current regulation in Brazil, IN nº 88 / 2021, set the maximum limit of 1 µg kg-1 for total aflatoxins (sum of G2, G1, B1 and B2) in infant formula. The limit regulating the presence of aflatoxin M1 in infant formula, however, has not yet been established in Brazil (Brasil 2021). For evaluation of sample contamination, the reconstitution stated on the formula containers and the limits according to EU regulation were considered.
The incidence and concentrations of aflatoxin M1 in the analyzed samples are displayed in Table 4. Aflatoxin M1 was found in 27 (22.0%) of the analyzed samples. In Brazil, incidences of aflatoxin M1 in infant formulas of 0, 23.1 and 43.8% have been reported respectively by Iha et al. (2013), Tonon at al. (2018) and Ishikawa et al. (2016). Incidences of aflatoxin M1 in infant formulas ranging from 1–100% have been reported in other countries (Kim et al. 2000; Baydar et al. 2007; Alvito et al. 2010; Gómez-Arranz and Navarro-Blasco 2010; Meucci et al. 2010; Beltran et al. 2011; Zhang et al. 2013, Omar, 2016 and Akhtar et al., 2017).
Table 4
occurrence OF AFLATOXIN M1 IN the analyzed INFANT FORMULA SAMPLES
Infant formula | No. of analyzed samples | No. of positive samples, (%) | Frequency distribution, µg kg− 1 (%) |
<LOD | LOD < LOQ | LOQ < 0.03 | 0.03 < 0.06 |
Follow-up | 52 | 15 (28.8) | 37 (71.2) | 11 (21.2) | 02 (3.8) | 02 (3.8) |
Starter | 49 | 10 (20.4) | 39 (79.6) | 07 (14.3) | 02 (4.1) | 01 (2.0) |
Specials | 17 | 01 (5.9) | 16 (94.1) | - | 01 (5.9) | - |
Toddler | 05 | 01 (20.0) | 04 (80.0) | 01 (20.0) | - | - |
Total | 123 | 27 (22.0) | 96 (78.0) | 19 (15.4) | 05 (4.1) | 03 (2.4) |
LOD: limit of detection, 0.005 µg kg− 1 |
LOQ: limit of quantification, 0.016 µg kg− 1 |
Aflatoxin M1 was quantified (≥ LOQ) in eight of the analyzed samples (6.5%), at concentration levels ranging from 0.016 to 0.056 µg kg-1, averaging 0.031 µg kg-1.
The maximum aflatoxin M1 concentration found herein (0.056 µg kg-1) was similar to the value previously reported by Ishikawa et al. (2016), of 0.046 µg kg-1, and about three-fold higher than the value reported by Tonon et al. (2018) (0.017 µg kg-1) in Brazil. Concerning studies carried out in other countries, the maximum value observed in the present study was below the values reported by Quevedo-Garza et al. (2020), Akhtar et al. (2017) and Omar et al. (2016), of 0.450 µg kg-1, 0.108 µg kg-1 and 0.154 µg kg-1, respectively, and higher than the values reported by Alvito et al. (2010), Gómez-Arranz and Navarro-Blasco (2010) and Beltran et al. (2011), of 0.041, 0.012 and 0.006 µg kg-1 reported, respectively.
The average aflatoxin M1 concentration observed in this study (0.031 µg kg-1, calculated using concentrations ≥ LOQ), was two-fold higher than the value of 0.015 µg kg-1 reported by Tonon et al. (2018), and similar to the concentration of 0.024 µg kg-1 reported by Ishikawa et al. (2016) in Brazil. Compared to results for other countries, the mean value reported herein is 10-fold higher that that reported by Gómez-Arranz and Navarro-Blasco (2010), of 0.0031 µg kg-1, similar to the values reported respectively by Quevedo-Garza et al. (2020) and Akhtar et al. (2017), of 0.040 µg kg-1 and 0.037 µg kg-1, and half the value of 0.06 µg kg-1 reported by both Kim et al. (2000) and Baydar et al. (2007).
The concentration levels (0.040, 0.044 and 0.057 µg kg-1) found for aflatoxin M1 in three (2.4%) of the analyzed samples exceeded the maximum permissible level set by EU regulations for infant formulas, of 0.025 µg kg-1. In other studies carried out in Brazil, this level was exceeded in three (18.8%) of the samples analyzed by Ishikawa et al. (2016), while all samples investigated by Iha et al. (2013) and Tonon et al. (2018) exceeded the established level. In studies carried out in other countries by Quevedo-Garza et al. (2020), Akhtar et al. (2017) and Omar et al. (2016), the maximum allowed level was exceeded respectively by 20.0, 30.8 and 85.0% of the analyzed samples, while none of the samples analyzed by Gómez-Arranz and Navarro-Blasco (2010) exceeded the permitted level.
Figure 2 displays the chromatograms for one of the samples naturally contaminated by aflatoxin M1 (0.044 µg kg− 1). The other investigated mycotoxins were not detected in the analyzed samples.
Estimated Daily Intake (Edi) Of Aflatoxin M1
The EDI for aflatoxin M1 for infants was calculated based on the approach applied by Ishikawa et al. (2016), considering the exclusive use of infant formulas. Reconstituted infant formula consumption of 590, 642, 560 and 452 g, respectively, for infant ages one week, one month, six months and 12 months were considered (Brasil 2005). The average weights considered for male and female were, respectively, 3.3 and 3.2 kg for infant ages one week, 4.5 and 4.2 kg for one month, 7.9 and 7.3 kg for six months, and 9.6 and 8.9 kg for twelve months (WHO 2006). The International Programme on Chemical Safety/Global Environment Monitoring System criteria (replacing all non-detectable results to LOD and all non-quantified results to LOQ), used when more than 80% of the samples present results are lower than the LOD, was adopted to calculate the average aflatoxin M1 concentration (0.0083 µg kg-1) in reconstituted, ready-to-use, infant formulas (IPCS/GEMS 1995). The results are displayed in the Table 5.
Table 5
Estimated daily intake of aflatoxin m1 for infants through infant formula coNsumption
Age | Infant formula consumption (g) | Average weight (Kg) | Estimated Daily Intake (ng kg− 1 body weigh day− 1) |
Male | Female | Male | Female |
1 week | 590 | 3.3 | 3.2 | 1.48 | 1.53 |
1 month | 642 | 4.5 | 4.2 | 1.18 | 1.26 |
6 months | 560 | 7.9 | 7.3 | 0.59 | 0.63 |
12 months | 452 | 9.6 | 8.9 | 0.39 | 0.42 |
Average of aflatoxin M1 (0.0083 µg kg− 1 or ng g− 1) in reconstituted infant formulas calculated according to the International Programme on Chemical Safety/Global Environment Monitoring System criteria (IPCS/GEMS, 1995). |
The EDI values found for aflatoxin M1 in the present study (0.39 to 1.53 ng kg-1 bw day-1) were higher than the values reported by Ishikawa et al. (2016) for infant formula intake by infants up to 12 months old in Brazil (0.078 to 0.306 ng kg-1 bw day-1).
The EDI results obtained in our study were higher than the values of 0.14 ng kg-1 bw day-1 for infant formulas and 0.014 ng kg-1 bw day-1 for pre-term feeding previously reported in Spain (Gómez-Arranz and Navarro-Blasco, 2010). In Mexico, the EDI for aflatoxin M1 reported by Quevedo-Garza et al., (2020) ranged from 1.56 to 14 ng kg-1 bw day-1, representing the values estimated for one year-old infants when fed starter or follow-on formulas, respectively