Next Article in Journal
Determination of Selected Isoquinoline Alkaloids from Mahonia aquifolia; Meconopsis cambrica; Corydalis lutea; Dicentra spectabilis; Fumaria officinalis; Macleaya cordata Extracts by HPLC-DAD and Comparison of Their Cytotoxic Activity
Previous Article in Journal
Development of Ultra-Performance Liquid Chromatography with Post-Column Fluorescent Derivatization for the Rapid Detection of Saxitoxin Analogues and Analysis of Bivalve Monitoring Samples
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Occurrence and Seasonal Variations of Aflatoxin M1 in Milk from Punjab, Pakistan

1
Department of Food Science & Human Nutrition, University of Veterinary and Animal Sciences, Lahore 54000, Pakistan
2
Department of Epidemiology & Public Health, University of Veterinary and Animal Sciences, Lahore 54000, Pakistan
3
Department of Dairy Technology, University of Veterinary and Animal Sciences, Lahore 54000, Pakistan
4
University Institute of Diet and Nutritional Sciences, Faculty of Allied Health Sciences, The University of Lahore, Lahore 54000, Pakistan
5
School of Exercise and Nutrition, Faculty of Health, Deakin University, Canberra 2600, Australia
6
Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, Jouf University, Sakaka 75471, Saudi Arabia
7
Student Research Committee, School of Medicine, Bam University of Medical Sciences, Bam 44340847, Iran
8
Zabol Medicinal Plants Research Center, Zabol University of Medical Sciences, Zabol 61615–585, Iran
9
Department of Nutrition and Dietetics, Faculty of Pharmacy, University of Concepción, Concepcion 4030000, Chile
10
Department of Clinical Oncology, Queen Elizabeth Hospital, 30 Gascoigne Road, Kowloon, Hong Kong 999077, China
*
Authors to whom correspondence should be addressed.
Submission received: 1 August 2019 / Revised: 16 September 2019 / Accepted: 28 September 2019 / Published: 2 October 2019

Abstract

:
The manifestation of aflatoxins in feed and food is a major issue in the world as its presence leads to some health problems. This study investigates the incidence of aflatoxin M1 (AFM1) contamination in raw milk samples which were collected from Punjab, Pakistan. The Cluster Random Sampling technique was used to collect 960 milk samples from five different regions, and samples were collected every month. The AFM1 level in raw milk was analyzed by the ELISA technique. The findings demonstrate that 70% of samples exceeded the United States permissible maximum residue limits (MRL 0.50 µg/L), with an overall AFM1 level that ranged from 0.3 to 1.0 µg/L. AFM1 contamination varied with the season: The highest average contamination was detected in winter (0.875 µg/L), followed by autumn (0.751 µg/L), spring (0.654 µg/L), and summer (0.455 µg/L). The Eastern region exhibited the highest average AFM1 contamination (0.705 µg/L). Milk samples from the Northern region were found to be widely contaminated, as 86.9% samples exceeded the US MRL, followed by the Eastern region, with 72.3% samples being contaminated with >0.5 µg/L AFM1. The study indicated that the raw milk supply chain was heavily contaminated. Recommendations and remedial measures need to be developed by regulatory authorities to improve the raw milk quality.
Key Contribution: The present study represents the first large-scale investigation of AFM1 contamination of the milk supply chain that covers the major areas of Punjab, Pakistan. This study also reports the seasonal changes in a year that might affect the milk AFM1 contamination.

1. Introduction

Aflatoxins belong to a widely studied group of mycotoxins. The major aflatoxins producing fungi are Aspergillus flavus and A. parasiticus [1]. Aflatoxin contamination in food, feed, and agricultural produce is a matter of colossal concern around the world because of their carcinogenic, metabolic, mutagenic, immunosuppressive, and teratogenic effects [2,3]. Four major types of aflatoxins (B1, B2, G1, and G2) can contaminate food and feed, posing serious health complications for human beings, as well as for animals [4]. Aflatoxin M1 (AFM1) is the monohydroxylated derivative of aflatoxin B1 (AFB1), developed in the liver of lactating animals during metabolism and further excreted into the raw milk of cattle usually fed with AFB1-contaminated feed [5]. The International Agency for Research on Cancer (IARC) classified AFM1 as a 2B carcinogenic group because it can damage DNA and may cause various types of cancers [6]. Chromosomal anomalies, cell transformation in mammals, and gene mutation are also a few of the known outcomes of AFM1 exposure [7].
Milk is one of the nutritious sources required for the better growth of infants and children, and at the same time, it is a rich source of nutrition for all age groups [8]. However, with an increasing demand for milk, it becomes challenging for the dairy sector to maintain a uniform and standardized quality in under developed countries. This situation is posing serious health threats to consumers. Research studies have verified that pasteurization, heat processing, and a few other techniques are ineffective for controlling AFM1 in raw milk [9,10]. Once aflatoxins contaminate the milk supply chain, it becomes impossible to completely purify raw milk [11,12]. The consumption of AFM1-contaminated dairy products would likely harm human health [10,13]. Recent studies have highlighted the alarming threats to health associated with the use of milk contaminated with aflatoxins. This is the reason why the detection of aflatoxins in agricultural commodities and development of effective strategies for their control are key research areas in the world [14]. The lethal nature and harmful impacts of aflatoxins on the health of humans and animals have increased the need for effective management. Furthermore, it is also evident that the occurrence of aflatoxins in the food supply chain is also affected by the season, weather, and contaminated feed ingredients [11,15,16].
Many countries have set a maximum residue level (MRL) in milk to ensure food safety. The MRL of AFM1 varies worldwide as agricultural practices, milk collection systems, and the food supply chain are different in various parts or regions of the world. The European Union (EU) established a principle of “As Low As Reasonably Achievable” (ALARA) regarding aflatoxin levels in food products. The ALARA implies that the maximum permissible level of AFM1 in milk in the EU is 0.05 µg/L [17], which is one tenth of the restrictions implemented in the United States (US) [18], Brazil [19], and China [20], and is also lower than Syria, where the limit is confined to 0.2 µg/L and 0.05 µg/kg for raw milk and powdered milk, respectively [21]. According to the report published in 2002, about 100 countries have set limits of AFM1 [2], while there are still many countries in the world which do not have regulations for AFM1 in raw milk and milk-based products [22].
Pakistan is the third largest milk producer in the world, with 40 billion liters of milk production annually [23]. During 2016, about 43,818 metric tons of milk was consumed nationwide [24], while the milk supply and demand gap is still wide. Despite the immense potential, a sizeable bulk is untapped in Pakistan due to ineffective and informal milk marketing. The major chunk of milk demand (94%) is being fulfilled through the casual supply chain (milk-man), while the packed milk industry bridges only 6%. This informal supply chain is likely to bear a high contamination of AFM1. As several studies have reported, AFB1-contaminated feed is the source of milk AFM1 contamination [25]. The consumption of these contaminated feed sources is more common in an informal supply chain compared to corporate farming [26]. Recently, Pakistan imposed 5 µg/L legal limits of the AFM1 level in milk for processing and 0.5 µg/L in milk for consumers, but the implementation is faulty due to the lack of strict monitoring facilities. This study was aimed at investigating the occurrence of AFM1 levels in widespread milk supply chains of Punjab regions by examining the milk throughout the year. Several efforts have been made to investigate the AFM1 contamination in milk of Punjab, Pakistan, but all of these studies were limited due to the sample size [27], area [28], season [15], milking time [29], and point of sample collection [28]. However, this study was conducted with the screening of the raw milk supply chain in a wider area, which constituted the major portion of milk supply network, and included all seasons of Punjab, Pakistan.

2. Results

The occurrence of AFM1 in raw milk samples (n = 960) collected from five regions of Punjab is presented in Figure 1, and the total number of samples per month in the regions exceeding the US maximum residual limit are shown in Table 1. The proportion of AFM1-contaminated milk samples was higher than 50% in all regions (Figure 2). The average value of AFM1 (0.705 ± 0.211) in milk samples was found to be higher in the Eastern region, which was significantly (P < 0.05) different from the average value of AFM1 in the Southern region (0.577 ± 0.219) and on par with the average of AFM1 in the Western (0.655 ± 0.193), Northern (0.639 ± 0.150), and Central (0.577 ± 0.219) regions. The value of AFM1 was found to be statistically non-significant for the Western, Northern, and Central regions.
In the current study, the findings were segregated in terms of regions, and the highest concentration of AFM1 was observed in the Eastern region, while the lowest was found in the Southern region (Table 1). The results of the Eastern region showed that the maximum range (0.52–1.67 μg/L) of AFM1 contamination in milk samples occurred during the month of November, and the lowest (0.16–0.63 μg/L) was seen in the month of May. Similarly, the maximum range of AFM1 contamination in other regions, i.e., the Western (0.13–1.56 μg/L), Northern (0.54–1.5 μg/L), and Central (0.36–1.63 μg/L) regions, occurred in the months of Nov–Dec. The contents of AFM1 in milk samples varied from 0.30 to 1.0 μg/L (Figure 1). The milk samples’ AFM1 contamination in all regions swiftly decreased from March to August. It was observed that AFM1 trends in all regions were almost similar to each other; however, after mid-April to August, the concentration remained below the US MRL (<0.5 μg/L). Immediately after August, a drastic increase in the AFM1 concentration in all regional milk samples was observed. The maximum contaminated samples (93.7%–100%) were in the month of February and the lowest (12.5%–43.7%) were in July.
One-way analysis of variance (ANOVA) was applied; to check for significant differences among the means of clusters and seasons, the Tukey HSD test was used. The results were considered statistically significant at p < 0.05. The clusters with the same superscript letters are statistically non-significant.
This study showed that, on average, 69% of milk samples displayed AFM1 contamination that was more than 0.5 μg/L US MRL. The highest contaminated samples (87%) were in the Northern region and the lowest (58%) were in the Southern region (Figure 2).
The extent of detection of AFM1 was at the maximum during winter (0.875 µg/L), followed by autumn (0.751 µg/L), spring (0.654 µg/L), and summer (0.455 µg/L). The results given in Table 2 showed that irrespective of regions, 54%–88% of milk samples were contaminated with AFM1 in the range of 0.5–1.00 µg/L for the winter, autumn, and spring season. Additionally, in the summer season, 54%–78% milk samples in all regions were below the US MRL, i.e., <0.5 µg/L.

3. Discussion

AFM1 contamination of food and feed is a global issue, especially in developing countries. The Punjab province of Pakistan is the major cash crop-producing and livestock-keeping area [30]. AFM1 contamination screening in the milk of this area depicts that the average AFM1 contamination in all regions of Punjab was above the recommended US regulations (0.5 μg/L). These findings are similar to those of Jawaid, et al. [31], who determined the AFM1 contamination in milk from the Sind province of Pakistan using ELISA and found 94% contaminated samples, such that 70% of samples surpassed the established standards of US regulations. Asghar, et al. [32] collected 156 fresh milk samples from local markets of Karachi, Pakistan, during 2016–2017. The results of these samples were analyzed by the ELISA technique and the mean contamination of AFM1 in milk was reported as 0.34 µg/L, with 32.7% samples being above the US allowable limits. Hussain, Anwar, Asi, Munawar and Kashif [27] performed HPLC with a fluorescent detection test to probe AFM1 contamination of the milk of cow, buffalo, camel, goat, and sheep. The contamination ranged from 16.7% to 34.5% for the five species. The severity of AFM1 contamination in Pakistani milk was further reported by Iqbal, Asi and Jinap [15], who analyzed 221 samples by HPLC and reported 45% contamination during the winter season and 32% contamination during the summer. Another study conducted in Central areas of Punjab Pakistan, found that 53% and 43% of buffalo and cow milk samples were contaminated, with an average AFM1 of 0.03 and 0.04 μg/L, respectively [22]. The findings further inferred that the AFM1 concentration was higher in milk samples obtained from urban or semi-urban areas compared to those from rural localities [22].
Many contemporary studies conducted across the world, especially in Asia, agree with the findings of the current study. Studies from border-sharing neighbor countries like India, having a similar geography, have reported ample incidences of AFM1 contamination in milk samples. Rastogi, et al. [33] carried out a competitive ELISA technique to analyze 87 milk samples and infant milks collected from the Indian market and found that 99% of samples were contaminated with AFM1 levels higher than the EU standard limits, whereas 9% of samples surpassed the US regulatory levels. The results of the current study are supportive of the reviews of Nile, et al. [34], where 46% of samples from different milking animals were contaminated. The mean value of AFM1 in buffalo, cow, goat, and sheep milk was 0.026, 0.018, 0.014, and 0.017 μg/L, respectively. In comparison to Pakistan, India has lower milk contamination issues; however, the difference in weather conditions, geography, and socio-economic status of both countries is not huge.
The climatic changes in Pakistan are the major reasons for AFM1 presence in the food supply chain. Usually, under climatic variations, farmers and farm owners face feed shortage and use alternatives, which become a source of AFM1 contamination. Cotton is the main crop of the Southern region, whereas corn is that of the Northern region. Corn dried under the sun in the open fields of Punjab, Pakistan, and stored in conditions that are favorable for fungal growth, results in the higher occurrence of AFB1. When the AFB1-contaminated sources are fed to the animals, they are secreted in the form of AFM1 in milk. In general, the weather conditions of the Southern and Northern regions are distinctly different due to variation in altitude. Hence, the Southern region is at the height of 400 ft, with a hot and dry climate, whereas the Northern region is around 2000 ft above the sea level and generally ambient and moist. The level of AFM1 was highly correlated with the topography [35] of the country, as indicated by the findings that contamination (%) was the highest in the Northern region and kept on decreasing, moving through to the Southern region (Figure 2).
The results were further compared with studies from another border-sharing countries, i.e., China. Raw milk and other dairy products (135 samples) from Northeast China revealed that 41% of samples were contaminated, with an average AFM1 concentration ranging from 0.32 to 0.50 μg/L [36]. In other research, an attempt was made to examine the 72 raw milk samples from China, which reported that 59.7% of samples were contaminated, with an AFM1 concentration ranging from 0.01 to 0.42 μg/L [16]. Further findings identified that the concentration of AFM1 was higher during the winter season (October–February). It is significant to compare the findings from the present study to those of China as both nations are neighbors and share more or less the same weather conditions. The results of the present study were further compared with Iran, which shares a border and weather conditions as well. Kamkar [37] found the presence of AFM1 in 76.5% of raw milk samples, with concentrations between 0.02 and 0.28 μg/L. The study reported that 40% of samples exceeded the EU maximum tolerance limit. An increased concentration of AFM1 was found during winter and it was shown to be inversely proportional to the temperature [15]. In another study performed by Hashemi [38], it was shown that the temperature had no role in the contamination of AFM1, but the study by Kamkar [37] reported a positive correlation between these two variables. Our results are in agreement with those of Kamkar [37]. Another study conducted in Pakistan found a higher concentration of AFM1 in the winter season [39] compared to the other seasons and these findings are in agreement with the present and previous studies [37]. It is worth mentioning that no study related to milk contamination with AFM1 was found in Afghanistan.
Numerous contemporary research studies across the globe have reported the contamination of AFM1 in milk. There are many reports on AFM1 contamination in milk, including results from Slovenia [40], North Africa [41], Turkey [42,43], Brazil [44], Portugal [6,45], Morocco [46], Syria [21], Sudan [47], Serbia [48], Spain [49], Korea [50], and Croatia [51].
In the current study, AFM1 contamination in all regions portrayed a similar trend, which indicated that it was being influenced by common contributors. It was clear from Table 2 that the amount of AFM1 contamination in milk was significantly higher in the winter season (0.875 ± 0.406 µg/L) compared to the autumn (0.751 ± 0.148 µg/L), spring (0.654 ± 0.037 µg/L), and summer (0.455 ± 0.052 µg/L) seasons. A lower amount of AFM1 contamination was found in the summer season, which was significantly different from all other seasons. The current study revealed that the type of animal feed and its handling are two major contributors to raw milk contamination. For example, the handling of cotton seed cake and corn, which are traditionally and nutritionally essential feed ingredients and valued throughout the year. The harvesting, drying, and storage of these feed ingredients in humid weather increases the concentration of toxins in feed and ultimately become the cause of higher AFM1 in milk [11,16,52]. The current findings are in line with the study of Hussain, Anwar, Munawar and Asi [22], which demonstrated that an inadequate availability of green fodder and excessive usage of concentrated feed and some of its ingredients, such as corn, wheat straw, cottonseed cake, paddy straw, and wheat bran, increase the AFM1 contamination in milk.
Previous studies [39,53] have pointed out that the contents of AFM1 are at the maximum during winter and autumn compared to other seasons. Higher AFM1 contamination in milk during winter is similar to the milk samples studied in Iran during winter (average 0.121 μg/L) [54]. Ghajarbeygi, et al. [55] reported similar findings of higher contamination during winter in Iran to those presented in Croatia [51,52]. In another study conducted in Pakistan by Asi, Iqbal, Ariño and Hussain [39], higher levels of AFM1 were documented in the winter season compared to the summer season, with significant difference (p < 0.05). The present study results are in exact agreement with previous studies conducted in Pakistan [15,22,27,31,35,56], where they reported higher contamination of milk with AFM1 in the winter season. The present study found high levels of AFM1 contamination in milk samples collected in the autumn season (0.751 μg/L) after winter (0.875 μg/L); hence, the results are in agreement with those of Aslam, et al. [57], who reported significantly higher AFM1 levels in milk during the autumn season (2.60 μg/L). Similar studies conducted by Hojjatollah and Sepideh [58] also reported an increased level of AFM1 contamination in milk during the winter and autumn seasons, with an average value of 54.24 and 34.92 ng/L, respectively. Winter and autumn are both critical seasons indicating the maximum presence of AFM1 (Table 2). At the same time, the majority of milk samples had higher levels of AFM1 contamination (than US permissible limits), regardless of seasons and regions (Table 1 and Table 2).
In a recent study, Miocinovic, et al. [59] concluded that AFM1 contamination in milk is greater in the autumn season. The results infer that the contamination of AFM1 in milk is influenced by seasonal variations, as reported in previous studies [39,60]. Jalili and Scotter [61] assessed the association between AFM1 contamination and feed consumed by animals in different seasons. The feeding pattern of animals more often alters from season to season. For instance, a number of research studies [22,62] have pointed out that the contamination of milk with AFM1 was influenced with a seasonal impact. The level of contamination in milk was observed to be greater in the cold season compared to the hot season in some other places [22,37,63].
Owing to the insufficient availability of feed during winter, farmers have to supplement animal feed, which is composed of various feed ingredients. The compound feed as a practice is prepared with discarded sources of grains and hence is susceptible to increases in aflatoxin levels in storage. In Pakistan, Asi, Iqbal, Ariño and Hussain [39] reported a higher AFM1 concentration in the milk of animals usually fed with compound feed compared to lactating animals fed on grazing or fresh green feed. It is a reality that animals have fractional access to grassland or pastures. Diet requirements are met with feed supplements. In addition to these factors, a higher milk yield in winter is a motivator for farmers to feed maximum inputs to animals that are in the form of corn, cotton seeds or cotton seed cake, raw rice bran/rice polish, wheat bran, and gluten. These ingredients are most vulnerable to fungal attack [22,36] and give the maximum output in terms of milk yield in winter. Another factor linked to the higher contamination in winter is the crop harvesting time. In Punjab, Pakistan, the corn-harvesting time is October, while the cotton-harvesting time is August to September, so they become a cheaper source of feed in these months. Dairy farmers grow these crops and feed animals without knowing their aflatoxin contamination. The commercial dairy farmers and farm farmers do not have the concept of standardized and hygienic feed usage. However, there are quite of few research studies which have shown no difference in AFM1 levels in samples obtained during winter and summer seasons [15,39,55,58]. At the same time, the results of these studies confirmed that the contamination of milk with AFM1 might be influenced by seasonal variations in different geographical locations. Despite the above factors, the weather conditions in Pakistan are unpredictable and variations in weather throughout the year may bring about abrupt changes in the contamination of AFM1.

4. Conclusions

The strength of this study lies with the results being obtained from a considerable number of samples (n = 960), which were collected during four seasons. In the current investigation, raw milk samples were collected from January to December. This research-based investigation confirms that the level of AFM1 in raw milk in Pakistan is due to erratic weather conditions and elevated levels of AFB1 contamination in animal feed. Regarding the seasonal impact on AFM1 in all regions, high levels were observed in the winter season, particularly from November to February. These observations are in agreement with previously published studies of elevated levels of AFM1 in winter. The increase in AFM1 level was higher compared to previous studies conducted in Pakistan, which shows a mounting concern of AFM1 in Pakistan. This study urges food authorities and livestock departments of Pakistan to mobilize food safety laws in the milk supply chain and initiate official regulations for the control of AFB1 in animal feed sources. Government agencies should introduce subsidy programs to introduce good storage practices, pre- and post-harvest control measures, effective grain drying techniques, and support to protect animal feed ingredients from the weather. The introduction of contamination-free forage sources, e.g., high-yield green fodder varieties, grasses, and silage-making techniques, could restrict AFM1 in the milk supply chain. Linkages between grain producers, feed manufacturers, dairy farmers, food, and livestock departments represent a dire need to combat increasing levels of contamination.

5. Materials and Methods

5.1. Sampling Plan

The Punjab province consists of thirty-six districts, which were divided into five regions, as follows: North (Attock, Chakwal, Rawalpindi, Sargodha, Khushab, Mianwali, Mandi Bhaudin, Gujrat, and Jhelum), South (Bahwalpur, Bhawalnagar, Multan, Khanewal, Lodhran, Muzaffargarh, Rahimyar Khan, Rajanpur, and Vehari), East (Gujranwala, Hafizabad, Pakpattan, Okara, Sahiwal, Kasur, Sailkot, and Norrowal), West (Dera Ghazi Khan, Bhakkar, and Layyah), and Central (Lahore, Faisalabad, Jhang, Toba Tek Singh, Chiniot, Sheikhupura, and Nankana Sahib), as shown in Figure 3 [64]. Sampling selection consisted of three stages. In the first stage, a district was selected from each region using Cluster Random Sampling. At the second stage, only one tehsil (an administrative area) was selected from each district. In the third stage, a milk shop was selected from each tehsil using random sampling. Milk samples were procured from local milk shops every month for one year.

5.2. Sampling

In this study, during the year 2015, a total of 960 milk samples were collected from 16 shops at random from each of the clusters (North, South, East, West, and Central) of Punjab, Pakistan. The plastic bottles, with a volume capacity of 200 mL, were filled with milk samples and placed in sampling coolers with icepacks to ensure that the temperature of the milk samples stayed below 10 °C. In most of the cases, samples were analyzed immediately; however, where the delay was unavoidable, samples were stored in a refrigerator.

5.3. Aflatoxin Determination

The competitive enzyme immunoassay using ELISA kit Ridascreen® Fast Aflatoxin M1, R5812 (R-Biopharm AG, Darmstadt Germany) was used for the quantitative analysis of AFM1 in the milk samples. Sample preparation was done by following the instructions given in the manual. Fifty (50 μL) samples or standards were added to a 96-well plate. Conjugate and anti-aflatoxin antibody solutions were added (50 μL) to each well. The solutions were mixed gently without splashing reagents from the wells by back and forth sliding on a flat surface for 10–20 seconds and incubated for 10 minutes at room temperature (20–25 °C). Liquid contents of the antibody wells were poured out by tapping the microwell holder upside down vigorously (three times in a row) against absorbent paper to ensure the complete removal of liquid from the wells. Again, each well was filled with 250 μL of washing buffer using an 8-channel pipette and poured out, and this step was repeated two times. After washing, 100 μL of substrate was added to the wells and gently mixed, as done previously, and incubated at room temperature in the dark for 5 minutes. Finally, the reaction was stopped by adding 100 μL of stop solution to each well and mixed gently, avoiding the splashing of reagents. Absorbance at 450 nm was measured by an ELISA plate reader (ELX800, Bio-Tek Instruments, USA).

5.4. Method Validation

The method was validated according to the criteria set by Commission Regulations 401/2006 [65]. The validation parameters were calculated in terms of linearity, limit of detection (LOD), limit of quantification (LOQ), linearity, repeatability (precision), and intermediate reproducibility (intra laboratory reproducibility), expressed as the relative standard deviation (RSD%) and trueness as recovery (Table 3 and Table 4).
Linearity was tested based on the six-point calibration curve acquired for an intermediate reproducibility (iR basis). Calibration curves were constructed as follows: the logit value for each standard was plotted against the decadic logarithm (log10) of the corresponding concentration. The logit value for each standard (log it STDX) was calculated using the following equation:
log i t S T D X = ln O D S T D X O D S T D 0 O D S T D X
where ODSTDX is the average optical density (duplicate measurements) of the standards recorded at 450 nm, and ODSTD 0 is the average optical density (duplicate measurements) of the standard containing no AFM1 (0.00 µg/L) recorded at 450 nm.
A standard curve was constructed each time using series of standards of AFM1 solutions provided with the test kit. The sample concentration was calculated based on the standard curve (Figure 4). The plot of the regression line showed a linear response (r2 > 0.99) within the working range studied (0.125 to 2.00 µg/L).
In order to calculate the precision, milk was analyzed without any purification step and after spiking at three different AFM1 concentration levels (0.00, 0.25, 0.50, and 0.75 µg/L), on six different days, and each time in duplicate. Consequently, the number of AFM1 determinations at each level was 12, which gave 48 ELISA measurements for this matrix. All repeatability and intermediate reproducibility values (i.e., CV(r) and CV (iR), in %) were within 7.0% to 9.8% and 5.4% to 24.6%, respectively. Recoveries were calculated using the results obtained from the repeatability and intermediate reproducibility. All recovery values were within 79.3% to 86.0%, with relative standard deviation (RSD%) from 9.6% to 13.1% and the limit of detection (LOD) of 0.10 µg/L.

5.5. Statistical Analysis

Numerical variable data were presented as frequencies, percentages, and the average ± standard deviation. The bar chart was used to express the total % of AFM1 samples exceeding US MRL <0.5 µg/L in selected clusters, while the line chart indicated levels of averages of AFM1 (µg/L) in different months of the year 2015. One-way analysis of variance (ANOVA) was applied; to check for significant differences among means of clusters and seasons, the Tukey HSD test was used. The results were considered statistically significant at p < 0.05. Correlation and regression analyses were applied to calculate R2 using the statistical package Q Stat.net.

Author Contributions

N.A., M.N. conceived the experimental design and collected milk samples and data from field. N.N., M.-u.-D.A. and S.I. helped with laboratory analyses and drafted the manuscript. A.R., M.I., T.A.G., M.A., B.S. performed statistical analyses and helped in writing results. J.S.-R., M.M. and W.C.C. designing experiment and writing the paper. All authors read and approved the final version of the document.

Funding

Funding to complete research was provided by Higher Education Commission, Islamabad-Pakistan in the form of scholarship [Pin #. 112-24484-2-VA1-275 (50021115)] to Naveed Akbar.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kumar, P.; Mahato, D.K.; Kamle, M.; Mohanta, T.K.; Kang, S.G. Aflatoxins: A global concern for food safety, human health and their management. Front. Microbiol. 2017, 7, 2170. [Google Scholar] [CrossRef] [PubMed]
  2. Van Egmond, H.P.; Jonker, M.A. Worldwide regulations on aflatoxins—The situation in 2002. J. Toxicol. Toxin Rev. 2004, 23, 273–293. [Google Scholar] [CrossRef]
  3. Sarma, U.P.; Bhetaria, P.J.; Devi, P.; Varma, A. Aflatoxins: Implications on health. Indian J. Clin. Biochem. 2017, 32, 124–133. [Google Scholar] [CrossRef] [PubMed]
  4. Kowalska, A.; Walkiewicz, K.; Koziel, P.; Muc-Wierzgon, M. Aflatoxins: Characteristics and impact on human health. Postepy Hig. Med. Dosw. 2017, 71, 315–327. [Google Scholar] [CrossRef]
  5. Prandini, A.; Tansini, G.; Sigolo, S.; Filippi, L.; Laporta, M.; Piva, G. On the occurrence of aflatoxin M1 in milk and dairy products. Food Chem. Toxicol. 2009, 47, 984–991. [Google Scholar] [CrossRef] [PubMed]
  6. Martins, M.L.; Martins, H.M. Aflatoxin M1 in yoghurts in portugal. Int. J. Food Microbiol. 2004, 91, 315–317. [Google Scholar] [CrossRef]
  7. WHO. Iarc Monographs on the Evaluation of Carcinogenic Risks to Humans; WHO International Agency for Research on Cancer: Lyon, France, 2006. [Google Scholar]
  8. Lucey, J.A. Raw milk consumption: Risks and benefits. Nutr. Today 2015, 50, 189–193. [Google Scholar] [CrossRef] [PubMed]
  9. Iqbal, S.Z.; Paterson, R.R.M.; Bhatti, I.A.; Asi, M.R. Survey of aflatoxins in chillies from pakistan produced in rural, semi-rural and urban environments. Food Addit. Contam. B 2010, 3, 268–274. [Google Scholar] [CrossRef]
  10. Rahmani, J.; Alipour, S.; Miri, A.; Fakhri, Y.; Riahi, S.M.; Keramati, H.; Moradi, M.; Amanidaz, N.; Pouya, R.H.; Bahmani, Z.; et al. The prevalence of aflatoxin M1 in milk of middle east region: A systematic review, meta-analysis and probabilistic health risk assessment. Food Chem. Toxicol. 2018, 118, 653–666. [Google Scholar] [CrossRef]
  11. Fallah, A.A.; Rahnama, M.; Jafari, T.; Saei-Dehkordi, S.S. Seasonal variation of aflatoxin M1 contamination in industrial and traditional iranian dairy products. Food Control 2011, 22, 1653–1656. [Google Scholar] [CrossRef]
  12. Ismail, A.; Goncalves, B.L.; de Neeff, D.V.; Ponzilacqua, B.; Coppa, C.; Hintzsche, H.; Sajid, M.; Cruz, A.G.; Corassin, C.H.; Oliveira, C.A.F. Aflatoxin in foodstuffs: Occurrence and recent advances in decontamination. Food Res. Int. 2018, 113, 74–85. [Google Scholar] [CrossRef] [PubMed]
  13. Duarte, S.C.; Almeida, A.M.; Teixeira, A.S.; Pereira, A.L.; Falcão, A.C.; Pena, A.; Lino, C.M. Aflatoxin M1 in marketed milk in portugal: Assessment of human and animal exposure. Food Control 2013, 30, 411–417. [Google Scholar] [CrossRef]
  14. Razzaghi-Abyaneh, M.; Chang, P.-K.; Shams-Ghahfarokhi, M.; Rai, M. Global health issues of aflatoxins in food and agriculture: Challenges and opportunities. Front. Microbiol. 2014, 5, 420. [Google Scholar] [CrossRef] [PubMed]
  15. Iqbal, S.Z.; Asi, M.R.; Jinap, S. Variation of aflatoxin M1 contamination in milk and milk products collected during winter and summer seasons. Food Control 2013, 34, 714–718. [Google Scholar] [CrossRef]
  16. Xiong, J.L.; Wang, Y.M.; Ma, M.R.; Liu, J.X. Seasonal variation of aflatoxin M1 in raw milk from the Yangtze River Delta region of China. Food Control 2013, 34, 703–706. [Google Scholar] [CrossRef]
  17. EFSA. Opinion of the scientific panel on contaminants in the food chain [contam] related to aflatoxin b1 as undesirable substance in animal feed. EFSA J. 2004, 2, 39. [Google Scholar] [CrossRef]
  18. FDA. Sec. 527.400 Whole Milk, Low Fat Milk, Skim Milk e Aflatoxin M1 (cpg 7106.10) fda/ora Compliance Policy Guides; FDA USA: Silver Spring, MD, USA, 2005.
  19. Oliveira, C.A.F.; Franco, R.C.; Rosim, R.E.; Fernandes, A.M. Survey of aflatoxin M1 in cheese from the north-east region of São Paulo, Brazil. Food Addit. Contam. B 2011, 4, 57–60. [Google Scholar] [CrossRef]
  20. Zheng, N.; Sun, P.; Wang, J.Q.; Zhen, Y.P.; Han, R.W.; Xu, X.M. Occurrence of aflatoxin M1 in uht milk and pasteurized milk in China market. Food Control 2013, 29, 198–201. [Google Scholar] [CrossRef]
  21. Ghanem, I.; Orfi, M. Aflatoxin M1 in raw, pasteurized and powdered milk available in the Syrian market. Food Control 2009, 20, 603–605. [Google Scholar] [CrossRef]
  22. Hussain, I.; Anwar, J.; Munawar, M.A.; Asi, M.R. Variation of levels of aflatoxin M1 in raw milk from different localities in the central areas of Punjab, Pakistan. Food Control 2008, 19, 1126–1129. [Google Scholar] [CrossRef]
  23. Zia, U.-E. Analysis of milk marketing chain—Pakistan. Ital. J. Anim. Sci. 2007, 6, 1384–1386. [Google Scholar] [CrossRef]
  24. G.O.P. Ministry of Finance. Economic Survery of Pakistan; G.O.P. Ministry of Finance: Islamabad, Pakistan, 2016–2017.
  25. Mahmoudi, R.; Norian, R. Aflatoxin b1 and M1 contamination in cow feeds and milk from Iran. Food Agric. Immunol. 2015, 26, 131–137. [Google Scholar] [CrossRef]
  26. Lanyasunya, T.; Wamae, L.; Musa, H.; Olowofeso, O.; Lokwaleput, I. The risk of mycotoxins contamination of dairy feed and milk on smallholder dairy farms in Kenya. Pak. J. Nutr. 2005, 4, 162–169. [Google Scholar]
  27. Hussain, I.; Anwar, J.; Asi, M.R.; Munawar, M.A.; Kashif, M. Aflatoxin M1 contamination in milk from five dairy species in Pakistan. Food Control 2010, 21, 122–124. [Google Scholar] [CrossRef]
  28. Muhammad, K.; Tipu, M.Y.; Abbas, M.; Khan, A.M.; Anjum, A.A. Monitoring of aflatoxin M1 in market raw milk in Lahore City, Pakistan. Pak. J. Zool. 2010, 42, 697–700. [Google Scholar]
  29. Iqbal, S.Z.; Asi, M.R.; Selamat, J. Aflatoxin M1 in milk from urban and rural farmhouses of Punjab, Pakistan. Food Addit. Contam. B 2014, 7, 17–20. [Google Scholar] [CrossRef] [PubMed]
  30. Younas, M.; Yaqoob, M. Feed resources of livestock in the Punjab, Pakistan. Livest. Res. Rural. Dev. 2005, 17, 2005. [Google Scholar]
  31. Jawaid, S.; Talpur, F.N.; Nizamani, S.M.; Afridi, H.I. Contamination profile of aflatoxin M1 residues in milk supply chain of Sindh, Pakistan. Toxicol. Rep. 2015, 2, 1418–1422. [Google Scholar] [CrossRef] [Green Version]
  32. Asghar, M.A.; Ahmed, A.; Asghar, M.A. Aflatoxin M1 in fresh milk collected from local markets of Karachi, Pakistan. Food Addit. Contam. B 2018, 11, 167–174. [Google Scholar] [CrossRef]
  33. Rastogi, S.; Dwivedi, P.D.; Khanna, S.K.; Das, M. Detection of aflatoxin M1 contamination in milk and infant milk products from Indian markets by ELISA. Food Control 2004, 15, 287–290. [Google Scholar] [CrossRef]
  34. Nile, S.H.; Park, S.W.; Khobragade, C.N. Occurrence and analysis of aflatoxin M1 in milk produced by Indian dairy species. Food Agric. Immunol. 2016, 27, 358–366. [Google Scholar] [CrossRef]
  35. Iram, W.; Anjum, T.; Abbas, M.; Khan, A.M. Aflatoxins and ochratoxin a in maize of punjab, pakistan. Food Addit. Contam. B 2014, 7, 57–62. [Google Scholar] [CrossRef] [PubMed]
  36. Pei, S.C.; Zhang, Y.Y.; Eremin, S.A.; Lee, W.J. Detection of aflatoxin M1 in milk products from China by ELISA using monoclonal antibodies. Food Control 2009, 20, 1080–1085. [Google Scholar] [CrossRef]
  37. Kamkar, A. A study on the occurrence of aflatoxin M1 in raw milk produced in Sarab city of Iran. Food Control 2005, 16, 593–599. [Google Scholar] [CrossRef] [Green Version]
  38. Hashemi, M. A survey of aflatoxin M1 in cow milk in Southern Iran. J. Food Drug Anal. 2016, 24, 888–893. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Asi, M.R.; Iqbal, S.Z.; Ariño, A.; Hussain, A. Effect of seasonal variations and lactation times on aflatoxin M1 contamination in milk of different species from Punjab, Pakistan. Food Control 2012, 25, 34–38. [Google Scholar] [CrossRef]
  40. Torkar, K.G.; Vengušt, A. The presence of yeasts, moulds and aflatoxin M1 in raw milk and cheese in Slovenia. Food Control 2008, 19, 570–577. [Google Scholar] [CrossRef]
  41. Elgerbi, A.M.; Aidoo, K.E.; Candlish, A.A.; Tester, R.F. Occurrence of aflatoxin M1 in randomly selected North African milk and cheese samples. Food Addit. Contam. 2004, 21, 592–597. [Google Scholar] [CrossRef]
  42. Kabak, B.; Ozbey, F. Aflatoxin M1 in uht milk consumed in turkey and first assessment of its bioaccessibility using an in vitro digestion model. Food Control 2012, 28, 338–344. [Google Scholar] [CrossRef]
  43. Tekinşen, K.K.; Eken, H.S. Aflatoxin M1 levels in uht milk and kashar cheese consumed in Turkey. Food Chem. Toxicol. 2008, 46, 3287–3289. [Google Scholar] [CrossRef]
  44. Garrido, N.S.; Iha, M.H.; Santos Ortolani, M.R.; Duarte Fávaro, R.M. Occurrence of aflatoxins M1 and M2 in milk commercialized in Ribeirão Preto-SP, Brazil. Food Addit. Contam. 2003, 20, 70–73. [Google Scholar] [CrossRef] [PubMed]
  45. Martins, M.L.; Martins, H.M. Aflatoxin M1 in raw and ultra high temperature-treated milk commercialized in Portugal. Food Addit. Contam. 2000, 17, 871–874. [Google Scholar] [CrossRef] [PubMed]
  46. Zinedine, A.; González-Osnaya, L.; Soriano, J.M.; Moltó, J.C.; Idrissi, L.; Mañes, J. Presence of aflatoxin M1 in pasteurized milk from Morocco. Int. J. Food Microbiol. 2007, 114, 25–29. [Google Scholar] [CrossRef] [PubMed]
  47. Elzupir, A.O.; Elhussein, A.M. Determination of aflatoxin M1 in dairy cattle milk in Khartoum State, Sudan. Food Control 2010, 21, 945–946. [Google Scholar] [CrossRef]
  48. Kos, J.; Lević, J.; Đuragić, O.; Kokić, B.; Miladinović, I. Occurrence and estimation of aflatoxin M1 exposure in milk in Serbia. Food Control 2014, 38, 41–46. [Google Scholar] [CrossRef]
  49. Rodriguez Velasco, M.L.; Calonge Delso, M.M.; Ordonez Escudero, D. Elisa and hplc determination of the occurrence of aflatoxin M(1) in raw cow’s milk. Food Addit. Contam. 2003, 20, 276–280. [Google Scholar] [CrossRef]
  50. Kim, E.K.; Shon, D.H.; Ryu, D.; Park, J.W.; Hwang, H.J.; Kim, Y.B. Occurrence of aflatoxin M1 in Korean dairy products determined by ELISA and HPLC. Food Addit. Contam. 2000, 17, 59–64. [Google Scholar] [CrossRef] [PubMed]
  51. Bilandžić, N.; Božić, Đ.; Đokić, M.; Sedak, M.; Kolanović, B.S.; Varenina, I.; Cvetnić, Ž. Assessment of aflatoxin M1 contamination in the milk of four dairy species in Croatia. Food Control 2014, 43, 18–21. [Google Scholar] [CrossRef]
  52. Bilandžić, N.; Božić, Đ.; Đokić, M.; Sedak, M.; Kolanović, B.S.; Varenina, I.; Tanković, S.; Cvetnić, Ž. Seasonal effect on aflatoxin M1 contamination in raw and uht milk from Croatia. Food Control 2014, 40, 260–264. [Google Scholar] [CrossRef]
  53. Hussain, I.; Anwar, J. A study on contamination of aflatoxin M1 in raw milk in the punjab province of Pakistan. Food Control 2008, 19, 393–395. [Google Scholar] [CrossRef]
  54. Rahimirad, A.; Maalekinejad, H.; Ostadi, A.; Yeganeh, S.; Fahimi, S. Aflatoxin M1 Concentration in Various Dairy Products: Evidence for Biologically Reduced Amount of AFM1 in Yoghurt. Iran. J. Public Health 2014, 43, 1139–1144. [Google Scholar] [PubMed]
  55. Ghajarbeygi, P.; Palizban, M.; Mahmoudi, R.; Jahed Khaniki, G.; Pakbin, B. Aflatoxin M1 contamination of cow’s raw milk in different seasons from Qazvin Province, Iran. J. Biol. Today World 2016, 5, 173–176. [Google Scholar] [CrossRef]
  56. Yunus, A.W.; Imtiaz, N.; Khan, H.; Ibrahim, M.N.M.; Zafar, Y. Aflatoxin contamination of milk marketed in Pakistan: A longitudinal study. Toxins 2019, 11, 110. [Google Scholar] [CrossRef] [PubMed]
  57. Aslam, N.; Tipu, M.Y.; Ishaq, M.; Cowling, A.; McGill, D.; Warriach, H.M.; Wynn, P. Higher levels of aflatoxin M1 contamination and poorer composition of milk supplied by informal milk marketing chains in Pakistan. Toxins 2016, 8, 347. [Google Scholar] [CrossRef] [PubMed]
  58. Hojjatollah, S.; Sepideh, T. The effect of milk composition, yeast-mould numbers and seasons on aflatoxin M1 amounts in camel milk. J. Food Saf. 2017, 37, e12300. [Google Scholar]
  59. Miocinovic, J.; Keskic, T.; Miloradovic, Z.; Kos, A.; Tomasevic, I.; Pudja, P. The aflatoxin M1 crisis in the serbian dairy sector: The year after. Food Addit. Contam. B 2017, 10, 1–4. [Google Scholar] [CrossRef] [PubMed]
  60. Tomašević, I.; Petrović, J.; Jovetić, M.; Raičević, S.; Milojević, M.; Miočinović, J. Two year survey on the occurrence and seasonal variation of aflatoxin M1 in milk and milk products in Serbia. Food Control 2015, 56, 64–70. [Google Scholar] [CrossRef]
  61. Jalili, M.; Scotter, M. A review of aflatoxin M1 in liquid milk. IJHSE 2015, 2, 283–295. [Google Scholar]
  62. Riazipour, M.; Tavakkoli, H.R.; Razzaghi Abyane, M.; Rafati, H.; Momtaz, S.S. Measuring the amount of M1 aflatoxin in pasteurized milks. Kowsar Med. J. 2010, 15, 89–93. [Google Scholar]
  63. Ruangwises, S.; Ruangwises, N. Occurrence of aflatoxin M1 in pasteurized milk of the school milk project in Thailand. J. Food Prot. 2009, 72, 1761–1763. [Google Scholar] [CrossRef]
  64. Khan, H.A.; Siddique, M.J.; Arshad, M.J.; Khan, Q.M.; Rehman, S.U. Sero-prevalence of peste despetits ruminants (PPR) goats in Punjab Province of Pakistan. Pak. Vet. J. 2007, 27, 109–112. [Google Scholar]
  65. Commission, E. Commission Regulation 401/2006 of 23 February 2006 laying down the methods of sampling and analysis for the official control of the levels of mycotoxins in foodstuffs. OJEU 2006, 70, 34. [Google Scholar]
Figure 1. Occurrence of aflatoxin M1 (AFM1) (µg/L) in raw milk among five regions of Punjab, Pakistan, during the year 2015 (n = 960).
Figure 1. Occurrence of aflatoxin M1 (AFM1) (µg/L) in raw milk among five regions of Punjab, Pakistan, during the year 2015 (n = 960).
Toxins 11 00574 g001
Figure 2. Percentage of milk samples exceeding the US AFM1 maximum residues level (<0.5 µg/L) throughout the year, in five regions of Punjab, Pakistan.
Figure 2. Percentage of milk samples exceeding the US AFM1 maximum residues level (<0.5 µg/L) throughout the year, in five regions of Punjab, Pakistan.
Toxins 11 00574 g002
Figure 3. Geographical locations of milk sampling sites in different colors from five selected regions of Punjab, Pakistan. Data Source ArcGIS Map 10.3.
Figure 3. Geographical locations of milk sampling sites in different colors from five selected regions of Punjab, Pakistan. Data Source ArcGIS Map 10.3.
Toxins 11 00574 g003
Figure 4. The plot of the regression line shows a linear response (r2 > 0.99) within the working range studied (0.125 µg/L to 2.00 µg/L).
Figure 4. The plot of the regression line shows a linear response (r2 > 0.99) within the working range studied (0.125 µg/L to 2.00 µg/L).
Toxins 11 00574 g004
Table 1. Aflatoxin M1 range (µg/L) and percentage of milk samples exceeding the US maximum residues level (MRL <0.5 µg/L) of aflatoxin M1 in five clusters of Punjab, Pakistan, during 2015 (n = 960).
Table 1. Aflatoxin M1 range (µg/L) and percentage of milk samples exceeding the US maximum residues level (MRL <0.5 µg/L) of aflatoxin M1 in five clusters of Punjab, Pakistan, during 2015 (n = 960).
ClusterTotal SamplesMonthsAverage AFM1 ± SD
JanuaryFebruaryMarchAprilMayJuneJulyAugustSeptemberOctoberNovemberDecember(µg/L)
EasternN = 192% samples exceeding US MRL93.793.77568.756.22543.743.77510010093.70.705 ± 0.211a
AFM1 range
(µg/L)
0.35–1.640.45–1.650.34–1.350.17–0.870.16–0.630.18–0.780.24–0.670.15–0.650.25–1.630.54–1.560.52–1.670.46–1.65
WesternN = 192% samples exceeding US MRL10010068.77543.737.531.212.562.581.287.5750.655 ± 0.193a,b
AFM1 range
(µg/L)
0.55–1.460.68–1.060.33–1.230.32–0.920.23–0.920.29–0.740.22–0.660.19–0.560.19–0.990.17–1.440.13–1.560.32–1.46
NorthernN = 192% samples exceeding US MRL10093.787.587.562.368.768.768.781.287.5751000.639 ± 0.150a,b
AFM1 range
(µg/L)
0.5–1.350.46–1.560.45–1.330.45–0.950.19–0.750.21–0.780.24–0.760.15–0.940.24–1.190.43–0.750.22–0.960.54–1.50
CentralN = 192% samples exceeding US MRL10010087.562.5012.512.56.281.210010087.50.637 ± 0.249a,b
AFM1 range
(µg/L)
0.55–1.440.56–1.430.46–0.850.35–1.230.16–0.470.12–0.640.16–0.660.04–0.890.34–1.590.64–1.000.65–1.560.36–1.63
SouthernN = 192% samples exceeding US MRL87.51007581.231.218.7018.756.2757581.30.577 ±
0.219b
AFM1 range
(µg/L)
0.47–1.140.56–1.30.43–0.880.33–0.990.23–0.750.3–0.650.19–0.490.13–1.130.19–0.840.26–1.30.13–0.870.40–1.50
Region wise, 16 samples were analyzed monthly in the year 2015. In total, 192 samples/region/year were analyzed. The total number of samples was 960. Underlined values are the ranges of the lowest and highest AFM1 (µg/L).
Table 2. Incidence of AFM1 in raw milk samples from five regions of Punjab, Pakistan, with respect to season.
Table 2. Incidence of AFM1 in raw milk samples from five regions of Punjab, Pakistan, with respect to season.
Range of AFM1 Concentration (µg/L)SouthernNorthernWesternEasternCentralAverage AFM1 + SD
WinterAverage AFM1 (µg/L)0.808 ± 0.2730.874 ± 0.2730.888 ± 0.2840.918 ± 0.3440.888 ± 0.2750.875 ± 0.406 a
n = 240Numbers of samples
(% frequency distribution)
<0.505(10%)2(4%)4(8%)3(6%)2(4%)
0.51–1.0030(63%)32(67%)27(56%)26(54%)34(71%)
>1.0013(27%)14(29%)17(35%)19(40%)12(25%)
AutumnAverage AFM1 (µg/L)0.605 ± 0.2060.594 ± 0.1550.767 ± 0.3440.902 ± 0.3320.890 ± 0.1880.751 ± 0.148 b
n = 160Numbers of samples
(% frequency distribution)
<0.508(25%)8(25%)6(19%)0(0%)0(0%)
0.51–1.0023(72%)24(75%)22(69%)23(72%)28(88%)
>1.001(3%)0(0%)4(13%)9(28%)4(03%)
SpringAverage AFM1 (µg/L)0.606 ± 0.1480.664 ± 0.1900.688 ± 0.2470.687 ± 0.2870.625 ± 0.1670.654 ± 0.037 c
n = 160Numbers of samples
(% frequency distribution)
<0.507(22%)7(22%)10(31%)6(19%)8(25%)
0.51–1.0025(78%)24(75%)19(59%)21(66%)23(72%)
>1.000(0%)1(3%)3(9%)5(16%)1(8%)
SummerAverage AFM1 (µg/L)0.416 ± 0.1850.506 ± 0.2020.458 ± 0.1800.507 ± 0.2410.391 ± 0.2510.455 ± 0.052 d
n = 400Numbers of samples
(% frequency distribution)
<0.5060(75%)40(50%)51(64%)43(54%)62(78%)
0.51–1.0019(24%)38(48%)29(36%)32(40%)16(20%)
>1.001(1%)2(3%)0(0%)5(6%)2(3%)
Overalln = 960Average AFM1 (µg/L)0.58 ± 0.220.64 ± 0.150.66 ± 0.190.71 ± 0.210.64 ± 0.250.640 ± 0.053
Statistical analysis: One-way analysis of variance (ANOVA) was applied; to check for significant differences among means of regions and seasons, the Tukey HSD test was used. The results were considered statistically significant at p < 0.05. Means with superscript letters a, b, c and d show significant differences between seasons.
Table 3. Calibration data for the coefficient of determination obtained by the linear model.
Table 3. Calibration data for the coefficient of determination obtained by the linear model.
AnalyteUnitConcentration RangeSlopeInterceptCoefficient of Determination
R2
Standard Deviation of Residuals
MinMaxCentral ValueSlope = 0
(Y/N)
Central ValueIntercept = 0 (Y/N)
AFM1µg/L0.1252.00.93N0.04N0.9960.044
Table 4. Validation parameters of the ELISA method used for the quantification of AFM1.
Table 4. Validation parameters of the ELISA method used for the quantification of AFM1.
Spiking Levels (µg/L)No. of Days x No. of ReplicatesMedianRepeatabilityIntermediate ReproducibilityLOD
(µg/L)
LOQ
(µg/L)
Recovery
(%)
RSD(Rec) Corrected
(%)
SD(r)CV(r)
(%)
RSD (iR)CV (iR)
(%)
iR
0.256 × 20.220.0219.80.0580.05324.60.1460.100.2586.09.6
0.506 × 20.390.0329.50.1020.04411.50.123--77.015.4
0.756 × 20.600.0437.00.1160.0325.40.089--79.313.1
SD, standard deviation; (r), repeatability; CV, coefficient of variation; iR, intermediate reproducibility; RSD (rec), corrected relative standard deviation of the recovery.

Share and Cite

MDPI and ACS Style

Akbar, N.; Nasir, M.; Naeem, N.; Ahmad, M.-u.-D.; Iqbal, S.; Rashid, A.; Imran, M.; Aslam Gondal, T.; Atif, M.; Salehi, B.; et al. Occurrence and Seasonal Variations of Aflatoxin M1 in Milk from Punjab, Pakistan. Toxins 2019, 11, 574. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins11100574

AMA Style

Akbar N, Nasir M, Naeem N, Ahmad M-u-D, Iqbal S, Rashid A, Imran M, Aslam Gondal T, Atif M, Salehi B, et al. Occurrence and Seasonal Variations of Aflatoxin M1 in Milk from Punjab, Pakistan. Toxins. 2019; 11(10):574. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins11100574

Chicago/Turabian Style

Akbar, Naveed, Muhammad Nasir, Naureen Naeem, Mansur-ud-Din Ahmad, Sanaullah Iqbal, Anjum Rashid, Muhammad Imran, Tanweer Aslam Gondal, Muhammad Atif, Bahare Salehi, and et al. 2019. "Occurrence and Seasonal Variations of Aflatoxin M1 in Milk from Punjab, Pakistan" Toxins 11, no. 10: 574. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins11100574

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop