Co-occurrence of mycotoxins and other fungal metabolites in total mixed rations of cows from dairy farms in Punjab, Pakistan

After India and the USA, Pakistan is the third country leading in global dairy production, a sector of very high socioeconomic relevance in Asia. Mycotoxins can affect animal health, reproduction and productivity. This study analysed a broad range of co-occurring mycotoxins and fungal secondary metabolites derived from Alternaria, Aspergillus, Fusarium, Penicillium and other fungal species. To complete this, a validated multi-metabolite liquid chromatography/electrospray ionization-tandem mass spectrometric (LC/ESI–MS/MS) method was employed, detecting 96 of > 500 tested secondary fungal metabolites. This first preliminary study demonstrated that total mixed rations (TMRs) (n = 30) from big commercial dairy cattle farms (> 200 lactating cows) in Punjab, Pakistan, presented ubiquitous contamination with mixtures of mycotoxins. The mean of mycotoxins per sample was 14, ranging from 11 to 20 mycotoxins among all TMR samples. Metabolites derived from other fungi and Fusarium spp. showed the highest levels, frequency and diversity among the detected fungal compounds. Among the most prevalent mycotoxins were Fusarium toxins like fumonisins B1 (FB1) (93%), B2 (FB2) (100%) and B3 (FB3) (77%) and others. Aflatoxin B1 (AFB1) was evidenced in 40% of the samples, and 7% exceeded the EU maximum limit for feeding dairy cattle (5 µg/kg at 88% dry matter). No other mycotoxin exceeds the EU guidance values (GVs). Additionally, we found that dietary ingredients like corn grain, soybean meal and canola meal were related to increased contamination of some mycotoxins (like FB1, FB2 and FB3) in TMR from the province of Punjab, Pakistan. Among typical forage sources, the content of maize silage was ubiquitous. Individually, the detected mycotoxins represented relatively low levels. However, under a realistic scenario, long-term exposure to multiple mycotoxins and other fungal secondary metabolites can exert unpredictable effects on animal health, reproduction and productivity. Except for ergot alkaloids (73%), all the groups of metabolites (i.e. derived from Alternaria spp., Aspergillus spp., Fusarium spp., Penicillium spp. and other fungi) occurred in 100% of the TMR samples. At individual levels, no other mycotoxins than AFB1 represented a considerable risk; however, the high levels of co-occurrence with several mycotoxins/metabolites suggest that long-term exposure should be considered because of their potential toxicological interactions (additive or synergistic effects). Supplementary Information The online version contains supplementary material available at 10.1007/s12550-023-00502-5.


Introduction
Located in Asia, the continent with the highest milk production worldwide, Pakistan is the third major milk producer after India and the USA (Food and Agriculture Organization of the United Nations 2021).The dairy cow diets contain various ingredients, including roughages, cereal grains and agro-industrial by-products (FAO, IDF, IFCN 2014).Crops and feedstuffs are vulnerable to mould infection and colonization with successive mycotoxin contamination during the complete feed production chain (pre-and post-harvest) influenced by several biotic and abiotic factors.The livestock industry endures severe economic losses due to the adverse effects of contaminated feed on animal health and the final quality of the products (Bryden 2012).The climatic conditions of Pakistan typically favour mycotoxin contamination in agricultural commodities (Ashiq 2015).According to a survey, South Asia was, in the last decade, the world's region with the highest occurrence of aflatoxin B1 (AFB1) Extended author information available on the last page of the article (82%) in feed samples.South Asia, along with Sub-Saharan Africa, showed the highest median values of AFB1-positive feed samples (≥ 20 µg/kg) (Gruber-Dorninger et al. 2019).AFB1 is a public health concern because of its proven carcinogenic properties (Massey et al. 1995).Previous studies on aflatoxin M1 (AFM1) in the Punjab Province of Pakistan indicated that 99%, 32% and 58% of the milk samples evaluated in the respective studies exceeded the European Union and Codex Alimentarius limit (0.05 µg/L), which indicates the constant exposure of dairy products' consumers to aflatoxins (Codex Alimentarius Commission 2001;Hussain and Anwar 2008;Iqbal and Asi 2013;Sadia et al. 2012).As feed is the central source of AFM1 in cow milk, the level of aflatoxins in dairy cattle diets should also be monitored and kept to a minimum (Sadia et al. 2012).
Previous surveillance studies on contamination of dairy cattle feed in Pakistan focused mainly on aflatoxins (AFs), zearalenone (ZEN), ochratoxin A (OTA) and trichothecenes (types A and B) (Ashiq 2015;Aslam and Wynn 2015;Gallo et al. 2015;Gruber-Dorninger et al. 2019;Santos Pereira et al. 2019;Yunus et al. 2020;Akbar et al. 2020).The most relevant investigated mycotoxins include the strictly regulated AFB1 and other mycotoxins with GV addressed by the EU legislation like deoxynivalenol (DON), ZEN, fumonisins (FBs), OTA as well as T-2 and HT-2 toxins (EC 2002(EC , 2006;;Gallo et al. 2015;Gruber-Dorninger et al. 2019).Although hundreds of compounds have been considered mycotoxins, most of the relevant studies investigated a limited number of mycotoxins in agricultural commodities (Gallo et al. 2015;Cinar and Onbaşı 2019;Battilani et al. 2020).Toxicological interactions (addition, synergism, potentiation and antagonism) among mycotoxins and other fungal metabolites affect animal and human health and reproduction (Smith et al. 2016).This requires more research and risk assessment by more integrative approaches (Battilani et al. 2020).Multi-mycotoxin contamination has been evidenced at pre-harvest and post-harvest (Rasmussen et al. 2010;Nichea et al. 2015a, b;Panasiuk et al. 2019;Hajnal et al. 2020;Penagos-Tabares et al. 2021, 2022a).It has been evidenced that dairy cattle diets such as total mixed rations (TMRs) are generally contaminated with complex cocktails of dozens of mycotoxins and other fungal and plant metabolites (Awapak et al. 2021;Penagos-Tabares et al. 2022b).TMR is a "complete ration" feeding system, which is very popular worldwide on dairy farms with big herds.TMR is produced by mixing forages, by-products, cereal grains, concentrates, minerals, vitamins and additives.From this mix, animals get the nutrients needed to meet maintenance and production requirements (Bueno et al. 2020;Schingoethe 2017).
Sub-clinical disorders in dairy cows, such as disrupted rumen function or increased susceptibility to infections, might be related to the impact of complex mixtures of toxic fungal secondary metabolites (Santos and Fink-Gremmels 2014).The relevance of synergistic interactions and consequences of long-term exposure to such mycotoxin mixtures is recognized, and the importance of integrative and innovative approaches based on multi-mycotoxin analyses has been highlighted (Battilani et al. 2020).Therefore, this investigation planned to determine the frequency, co-occurrences and concentration of contamination with mycotoxins and other fungal metabolites (> 500) derived from species of Alternaria, Aspergillus, Fusarium, Penicillium and other fungi in the TMR samples of dairy cattle farms from Punjab, Pakistan.The analysis was accomplished by employing a validated multi-metabolite liquid chromatography/electrospray ionization-tandem mass spectrometric (LC/ESI-MS/ MS) method.The possible relationship of the main dietary ingredients to the dietary concentrations of mycotoxins and other metabolites was also explored.

Sampling and sample preparation
After obtainment of written authorization and consent of the farmers, TMR (n = 30) samples were collected from corporate dairy farms in Punjab, Pakistan (Fig. 1a).The herd size of the farms was over 200 Holstein-Frisian lactating cows.The farms were selected so that each of the nine administrative divisions of Punjab contributed at least three farms.Information regarding the TMR composition (main ingredients, proportions and estimated feed intake) was provided by the farmers (n = 29/30) via a personal (questionnaire-guided) interview.Each representative sample of TMR consisted of a minimum of 30 incremental samples, which were manually collected from the feed bunk with gloves directly after the serving (Fig. 1b).The final TMR sample amount was 1-1.5 kg, which was mixed and immediately vacuum-packed and stored in the dark at − 20 °C.Sampling was carried out during the period June-July of 2020.For the sample preparation, the frozen TMRs were thawed at room temperature for 24 h and airdried at 65 °C for 48 h.Then, the dried TMRs were milled to a final particle size of ≤ 0.5 mm, using the cutting mill (SM 300; Retsch GmbH, Haan, Germany) at 1500 rpm for approximately 1 min, and the remnants (> 0.5 mm) were processed using an ultra-centrifugal mill (ZM 200; Retsch GmbH, Haan, Germany) at 10,000 rpm for approximately 30 s, following the procedures described by Penagos- Tabares et al. (2022a, b).Finally, 5 g (± 0.01 g) of each homogenized TMR sample was weighed into 50-mL polypropylene conical tubes (Sarstedt, Nümbrecht, Germany) and stored at − 20 °C until posterior analysis targeting multiple mycotoxins and other fungal secondary metabolites.

Multi-mycotoxin analysis (LC/ESI-MS/MS)
Water purification was done using a Purelab Ultra system (ELGA LabWater, Celle, Germany).Glacial acetic acid (p.a.) and ammonium acetate (LC-MS grade) were purchased from Sigma-Aldrich (Vienna, Austria).HiPerSolv Chromanorm HPLC gradient grade acetonitrile was obtained from VWR Chemicals (Vienna, Austria), and LC-MS Chromasolv grade methanol was acquired from Honeywell (Seelze, Germany).Standards of fungal, plant and unspecific secondary metabolites were purchased from several commercial suppliers or obtained via a donation from different research institutions (Sulyok et al. 2020).For simultaneous multiple metabolite quantification, 5 g (± 0.01 g) of the TMR sample was extracted in 20 mL of the extraction solvent (acetonitrile/water/acetic acid 79:20:1, v/v/v) following the procedures reported by Sulyok et al. (2020).These volumes were put into the QTrap 5500 LC-MS/MS system (Applied Biosystems, Foster City, CA, USA) equipped with a Turbo V electrospray ionization (ESI) source coupled to a 1290 series UHPLC system (Agilent Technologies, Waldbronn, Germany, as described by Sulyok et al. 2020).Subsequently, quantification from external calibration by serial dilutions of a stock solution of analysed compounds was accomplished.In the end, the results were adjusted for apparent recoveries defined through spiking experiments according to Steiner et al. (2020).This analytical methodology has been validated (Steiner et al. 2020;Sulyok et al. 2020) and has been utilized to study the multi-mycotoxin occurrence in complex feedstuff matrices like silage, pastures, concentrate feed and TMR (Shimshoni et al. 2013;Nichea et al. 2015b;Reisinger et al. 2019;Awapak et al. 2021;Penagos-Tabares et al. 2021, 2022a, b).The method accuracy has been verified on a routine basis by proficiency testing organized by BIPEA (Gennevilliers, France).Satisfactory z-scores between − 2 and 2 have been obtained for > 95% of > 1700 results submitted so far.

Statistical analysis
Concentrations of mycotoxins and other fungal metabolites were presented on a dry matter basis in µg/kg.Descriptive statistics (i.e.occurrences, mean, median and range of the concentrations of mycotoxins and metabolites) were processed, considering only the positive values (x ≥ limit of detection (LOD) using Microsoft ® Excel ® .Values between the limit of quantification (LOQ) and LOD were calculated as LOQ/2.A two-tailed Spearman's correlation test was conducted to explore possible relationships between dietary compounds and levels of metabolites, as well as relationships among metabolites within each ingredient compound.For this, only data of metabolites with frequencies over 30% was studied.Spearman's correlation coefficients were considered significant at a p value < 0.05.Accordingly, the correlation coefficients were interpreted according to Hinkle et al. (2003) as follows: "very high" (0.90 up to 1.00), "high" (0.70 up to 0.90), "moderate" (0.50 up to 0.70), "low" (0.30 up to 0.50) and "negligible" (< 0.30).Low and negligible correlations were not considered during interpretation in the results' description.Linear regressions between fungal metabolites and the content of certain feed ingredients were performed to corroborate the promising relationships.The statistical analyses and graphs were performed using GraphPad Prism version 9.1 (GraphPad Software, San Diego, CA, USA).

Main dietary components
The selected farms fed TMR, which is a feeding method consisting of complete diet composed of the mixtures of forages and varying quantities of concentrate feed, byproducts and mineral supplements.The frequency and rate of the inclusion levels of the main TMR ingredients offered to lactating cows for all selected farms are shown in Table 1.The most common dietary components included were maize silage (100%), commercial concentrate (90%), corn (maize) grain (83%), soybean meal (83%), canola meal (79%), molasses (72%), wheat straw (52%), Rhodes grass hay (34%), rice polish (21%) and wheat bran (21%).Other feedstuffs, including lucerne hay, rapeseed cake, palm kernel cake, maize gluten, lucerne, sugar beet pulp, cotton seed cake and rice bran, were less frequently (< 20%) included (Table 1).Proportionally, maize silage was the most abundant dietary ingredient, with an average inclusion of 65.1% DM of the ration, varying from 41.6 to 77%.Maize grain was incorporated on an average proportion of 10.2% of the ration (DM basis), followed by commercial concentrate (8.7%), lucerne hay (6.1%) and soybean meal (6.1%).On average, the TMR samples contained 30.4% of concentrate feeds and 69.6% of forages.The forage-to-concentrate ratio (F:C) fluctuated between 52:48 and 84:16 (Table 1).

Relationship between concentrations and groups of mycotoxins and metabolites and the dietary ingredients
Positive moderate correlations (ρ > 0.5, p value < 0.001) were observed between corn grain, soybean meal and canola meal with FB1 and FB2 (Table 3).Also, the ergot alkaloid ergometrinine correlated positively with the content of molasses (ρ = 0.54, p value < 0.001).The proportion of commercial concentrate correlated negatively with the contamination levels of bikaverin (ρ = 0.54, p value < 0.001), FB1 (ρ = − 0.56, p value < 0.001), FB3 (ρ = − 0.50, p value < 0.001) and moniliformin (ρ = − 0.56, p value < 0.001).All the values of the correlation analysis, i.e. ρ correlation coefficients and p values, are available in Supplementary Table S1.The moderate correlations between dietary ingredients and some of the toxins and metabolites (like FB1, FB2, bikaverin and ergometrinine) were confirmed by regression analyses (p value < 0.05) (Supplementary Fig. S1).Via regression analysis, a proportion of commercial concentrate showed no significant negative relationship with the levels of moniliformin (p value = 0.137).

Discussion
This study describes the mixtures of mycotoxins and other fungal metabolites in the complete diets of lactating cows at corporate dairy farms in Punjab, Pakistan.The region is considered a central crop-producing province and a crucial livestock-keeping area in the country (Younas and Yaqoob 2005;Akbar et al. 2019).The presented results again demonstrated the ubiquitous presence of mycotoxin mixtures in the complete diets of dairy cows.The mixtures fluctuated from 11 to 20 mycotoxins per ration.The cocktails of mycotoxins in commercial dairy cow farms have been previously revealed using multi-metabolite approaches (Awapak et al. 2021;Penagos-Tabares et al. 2022b).It is vital to note that this study is not representative for Pakistan, where most milk production and commercialization are informal.In the country, the formal sector (which sampled farms belong) has a small market share of merely 5%.Most dairy production and commercialization remain informal (Godfrey et al. 2018).Thus, public health's additional investigation focused on Pakistan's informal dairy chain is highly advised.
This research confirmed (as expected) the presence of AFB1.This toxin, produced mainly by Aspergillus spp., is the most toxic and recurrent among the AFs and is the most potent natural hepatocarcinogenic agent in mammals.It is classified as a group 1 human carcinogen by the International Agency for Research on Cancer (IARC 2012; Marchese et al. 2018).AFB1 was detected in around 40% of the rations in the current study.This situation implicates the AFM1 contamination of milk in some of these farms, which is a global issue, particularly in developing countries (Groopman et al. 2008).AFB1 concentrations exceeding the maximum level of the EU were detected in 7% of the TMR.Additionally, precursors of AFs, such as averufin, STC and versicolorin C (Cary et al. 2006;Hsieh et al. 1973), were detected in frequencies < 35%.Regarding AFs and STC, it has been suggested that these mycotoxins can be produced pre-and post-harvest (Mo et al. 2015).Like AFs, STC is known to be carcinogenic with immunotoxin and immunomodulatory activity.The information available on exposure data of dairy cows and other animals to STC is limited (EFSA 2013;Gruber-Dorninger et al. 2017;Chuang et al. 2020).One sample was contaminated with AFM1 (1.18 µg/kg), and AFM1 is not found in plants, a fact that indicates the contamination with animal products, particularly milk or dairy products mixed in the TMR (Min et al. 2021).
Aflatoxicosis in cattle includes clinical signs such as poor weight gains, decreased feed conversion and milk production, lethargy, inappetence, ataxia and increment of hepatic enzymes and bilirubin, in addition to prolonged clotting times (Diekman and Green 1992).Cows fed with diets containing AFB1 at concentrations of 20 µg/kg Fig. 3 Scatter plots showing the co-contamination (number of metabolites/sample) in each metabolite group detected in the TMR samples from Punjab, Pakistan.The grey lines indicated the average numbers of detected metabolites per sample presented a depletion in the feed intake and milk yield.Three days after the source of AFB1 was removed, the clinical signs began to improve (Jones and Ewart 1979).Similarly, another field study, which assessed the effect of aflatoxin-contaminated corn on lactating dairy cattle, observed a decline in reproductive efficiency.After the inclusion of an aflatoxin-free diet, an increment of 25% of the milk yield was evidenced (Guthrie and Bedell 1979) as cited by Jouany and Diaz (2005).Several case reports of acute aflatoxicosis in cattle have been described.For example, a group of crossbred feeder steers fed with corn contaminated with aflatoxin at a concentration of 1.5 µg/ kg generated typical hepatic lesions.Mycotoxin residues were detected in the kidney tissue (Colvin et al. 1984).In the same way, a small herd of cattle having access to mouldy and unharvested sweet corn revealed via postmortem examinations oedema of all soft tissues and liver lesions consistent with aflatoxicosis.Weather conditions were favourable for the proliferation of Aspergillus flavus and Aspergillus parasiticus, and the contamination levels of the corn samples taken from the field contained 2365 ng of aflatoxin/g (Hall et al. 1989).Minimizing AFB1 contamination in dairy feeds needs good agricultural and management practices at pre-harvest stages, such as appropriate harvest time, maintenance of crop health by avoiding pest infestations and use of fungal-resistant varieties of crops.At post-harvest, the reduction of moisture for conserved feedstuffs, proper storage at low temperature and humidity and protection against pest infestation (by insecticides and fungicides) are advocated.Also, routine monitoring for aflatoxins in feeds using aflatoxin binders/ inactivators in feed and creating awareness among farmers on the health impacts of aflatoxins have been proposed to reduce risks (Patyal et al. 2021).
Mycotoxins from Fusarium (e.g.FB1, FB2, beauvericin and bikaverin), Alternaria (e.g.alternariolmethylether and tenuazonic acid) and Aspergillus (e.g.kojic acid, averufin and STC) along with Penicillium toxins (like mycophenolic acid) and other metabolites were recently reported in diets of dairy cattle in Thailand and Austria (Awapak et al. 2021;Penagos-Tabares et al. 2022, survey).Like in European dairy cattle diets (Penagos-Tabares et al. 2022b), Mycotoxins included in this analysis occurred in ≥ 30% of the samples the current study found Fusarium mycotoxins/metabolites as a dominant group of fungal metabolites.The occurrences and levels of ZEN, enniatins and ergot alkaloids were lower than those of diets of Austrian dairy cattle (Penagos- Tabares et al. 2022b), and fumonisin contamination was higher in the dairy diets of Pakistan (100% of the samples contaminated with at least one fumonisin).Our results also highlight the role of tenuazonic acid as the most abundant mycotoxin produced by Alternaria spp.; however, the information regarding the occurrence and toxic effects of these toxins in animals is still scarce and, therefore, health risks associated with Alternaria toxins in feeds have not yet been clarified (EFSA 2011).Regarding the occurrence of trichothecenes, with an occurrence of 40%, NIV was the only mycotoxin of this group detected.A study reported NIV as the most occurrent trichothecene (12.3%), which was usually detected co-occurring with other trichothecenes like DON, T-2, HT-2 and 3-acetyldeoxynivalenol in maize grain in Punjab, Pakistan.However, the study also reports samples contaminated only with NIV, which suggest the presence of NIV-dominant Fusarium chemotypes (Khatoon et al. 2012).
Various Penicillium-derived compounds have been previously detected in silages, such as mycophenolic acid, mycophenolic acid IV and andrastin A (Gallo et al. 2015;Penagos-Tabares et al. 2022a, b;Storm et al. 2014).OTA, contemplated in the European regulation, is produced by species of Penicillium and Aspergillus and presented a low occurrence and contamination levels, which suggest that this mycotoxin presents a negligible risk for dairy herds, in line with previous studies (Driehuis et al. 2008;Awapak et al. 2021;Penagos-Tabares et al. 2022b).Additionally, produced primarily by Aspergillus spp., but by Penicillium and Acetobacter fungi (Parrish et al. 1966), kojic acid has shown low toxicity for human macrophages and antibacterial and immunomodulatory properties (Morton et al. 1945;Kotani et al. 1976;Bashir et al. 2021).Additionally, further less-known metabolites are produced by other fungi detected in the dairy cows' diets.Some of them have antibacterial activity, for example the illicicolins (Hayakawa et al. 1971), cytochalasins (Aldridge et al. 1967;Jouda et al. 2016) and ascochlorin also known as antibiotic LL-Z1272γ and ilicicolin D (Molnár et al. 2010).The current results showed tenuazonic acid as the most abundant mycotoxin produced by Alternaria spp.; however, the information regarding the occurrence and toxic effects of these toxins in animals is still scarce and, therefore, health risks associated with Alternaria toxins in feeds have not been elucidated (EFSA 2011).
The critical factors facilitating the growth of aflatoxinproducing moulds in corn grains and silage include, among others, lack of good agricultural storage practices and unfavourable climatic conditions (Kebede et al. 2012;Frazzoli et al. 2016).The risk of aflatoxin contamination is generally higher in geographical regions with a tropical climate or a sub-tropical climate, but an extremely hot and droughty season may promote the growth of Aspergillus spp. in crops (Kebede et al. 2012).AFB1 has been reported in dairy feeds in Thailand, with an occurrence of 39% in concentrate (Awapak et al. 2021).European reports are rare; however, 61% of the TMR from Lithuanian dairy farms tested positive for AFB1 (mean: 2.42 µg/kg, range: 1.03-5.00µg/kg) (Vaičiulienė et al. 2021).The incidence  (Decastelli et al. 2007).A moderate positive correlation of molasses with ergot alkaloids (specifically with ergometrinine) can be explained because ergot can grow on sugarcane (Saccharum spp.) (Singh 1976).Molasses is dehydrated sugarcane juice; some ergot alkaloids could be found in high concentration in molasses due to its concentration during fabrication and the thermostability of the ergot alkaloids.Different metabolite profiles could result from the same genus and species depending on the high variability of strains, substrate and growing conditions (Daou et al. 2021).The diversity of mycotoxins and fungal secondary metabolites is due to the multi-commodity composition of the diets.Despite the risk associated with toxicological interactions of mycotoxins, there is hardly any regulation on their combined occurrence globally (Battilani et al. 2020;Singh 2022).This study evidenced the high occurrence of a broad number of mycotoxins (most of them not contemplated in the legislation) and other fungal secondary metabolites occurring in dairy TMR in Pakistan.Around 7% of the samples exceeded the GVs of the EU commission for AFB1.Moreover, a vast majority of mycotoxins and metabolites are emerging, as well as less-known and less-studied fungal metabolites.After the compounds derived from other fungal species were analysed, it was observed that Fusariumproduced metabolites and mycotoxins were the dominant fungal contaminants.Additionally, the data derived from Spearman's correlation test (Table 3 and Table S1) and lineal regressions (Fig. S1) show consistently that moderate positive relationships among the dietary contents of ingredients like corn grain, soybean meal and canola meal were related to increased contamination of some Fusarium mycotoxins (like FB1, FB2 and FB3) in the TMR from the province of Punjab, Pakistan.Considering the low sample size of this exploratory study, both statistical methods (correlations and linear regression) were used to explore the relationships of ration formulation (ingredients) with mycotoxin/metabolite concentrations.In contrast with studies in other regions like South America and Europe (Driehuis et al. 2008;Signorini et al. 2012;Reisinger et al. 2019;Penagos-Tabares et al. 2022a), our results do not reveal maize silage as one of the most influential feedstuffs to the mycotoxin/metabolite contamination.Among the typical forage sources, the content of maize silage was ubiquitous in the analysed rations.However, a previous study suggested the role of cottonseed cake as the contributor to around 80% of the AFB1 in diets of dairy cattle in periurban farms in Punjab (Yunus et al. 2020).It is also crucial to consider that more consistent association and relationship assessments would require a higher sample size.
Except for AFB1, which represents a risk for animals and human consumers due to AFM1 content in milk (Min et al. 2021), other detected mycotoxins correspond to a relatively low level of risk.However, the realistic scenario, the longterm exposure to multiple mycotoxins and other fungal secondary metabolites could have unpredictable effects on animal health, reproduction and productivity.However, the high co-occurrence of various mycotoxins/metabolites should be investigated because of their potential toxicological interactions (additive or synergistic effects) and long-term effects at low chronic exposure (Smith et al. 2016;Battilani et al. 2020).At the detected levels, no other mycotoxin than AFB1 was reported to have a considerable transfer of metabolites into milk and other animal products.The findings suggest that it is necessary to design effective strategies to verify the safety of feedstuffs utilized in ration formulation.More surveillance and further research based on multi-metabolite methodologies in the dairy industry in other geographic regions of Pakistan and the world, considering seasonal variation, are still strongly encouraged.More governmental interest and research are essential for this concern to ensure the offer of safe dairy products to the consumer and support animal health and the productive potential of dairy herds.
This exploratory study evidenced that the most relevant mycotoxin for public health, the carcinogenic AFB1, is occurring in diets of big commercial dairy farms (> 200 lactating dairy cows) in the province of Punjab, Pakistan.AFB1 was detected in concentrations seven times higher than the EU maximum limit, representing a severe risk to animal health and human milk consumers.No other mycotoxin than AFB1 exceeded the EU guidance values.Except for ergot alkaloids, all the groups of metabolites (i.e.derived from Alternaria spp., Aspergillus spp., Fusarium spp., Penicillium spp.and other fungi) occurred in 100% of the TMR samples.Although the detected contamination levels of single compounds (mycotoxins/metabolites) are moderately low, the effects on animal health, reproduction and productivity under the detected realistic scenario ("cocktails effect") are still unpredictable.Similar studies with higher sample size and approaching other regions are extremely advocated.Thus, future toxicological studies should address such interactions (additivity, potentiation, synergism and antagonism), as well as the long-term exposure effects of "mycotoxin mixtures".The presented results reconfirm that the monitoring and surveillance of aflatoxin M1 in dairy products in the South Asian region are essential and highly required.

Fig. 1
Fig.1Representative sampling of total mixed rations (TMRs) from dairy farms in Punjab, Pakistan.a Map of the province of Punjab, illustrating the localization of explored farms.b The representative sampling consisted of at least 30 incremental (handful) samples collected from the feeding table immediately after serving

Fig. 2
Fig. 2 Distribution of the concentration of categorized mycotoxins and other fungal metabolites detected in the TMR samples in Punjab, Pakistan.The total number of metabolites detected per group is shown in parentheses

Fig. 4
Fig. 4 Heat map of the most frequent mycotoxin/metabolite combinations (in %) detected in the TMR samples (n = 30) from Punjab, Pakistan.Mycotoxins included in this analysis occurred in ≥ 30% of the samples

Table 2
Occurrences and levels of mycotoxins and other fungal metabolites detected in total mixed rations of dairy farms in Punjab, Pakistan

Table 3
Spearman's correlation among the proportions of dietary ingredients incorporated and levels of mycotoxins/fungal metabolites detected in total mixed rations of dairy farms in Punjab, Pakistan a n = 29, one sampled farm declined to provide the information of the total mixed ration (TMR) composition *p value < 0.05, significant; **p value < 0.01, highly significant