Aflatoxin B1 contamination and fungal diversity in oilseed cakes and formulated feeds from livestock feed markets in Addis Ababa and surrounding areas
Adey Melesse Yalew, Alemayehu Teressa Negawo, Tesfaye Alemu, Jean Hanson, Michael Blummel, Nebiyat Nigusie Woldeyohannis, Yilikal Assefa Teferi, Yonas Asmare, Prasad Kodukula V.S.V., Chris S. Jones
TL;DR
This study found high levels of aflatoxin B1 and diverse fungi in livestock feeds in Addis Ababa, highlighting the need for better feed safety practices.
Contribution
The study provides new insights into fungal diversity and aflatoxin contamination levels in oilseed cakes and formulated feeds in Ethiopia.
Findings
92% of poultry ration and 72% of dairy ration exceeded aflatoxin B1 permissible levels.
Aspergillaceae was the dominant fungal family, with significant variation across feed types.
Multiple harmful fungal species were identified, including Aspergillus flavus in several feed samples.
Abstract
This study was conducted in Addis Ababa, surrounding peri-urban and urban feed markets to assess the prevalence of aflatoxin B1 (AFB1) and identify fungal diversity in livestock feeds. Two hundred and ninety samples of oil seed cakes and formulated feeds were collected from the randomly selected feed vendors. The samples were analysed for aflatoxin B1 using a Helica® Aflatoxin B1 (low matrix) Enzyme-linked immunoassay (ELISA) kit. To validate the ELISA results, a subset of randomly selected samples was further analysed using High-performance liquid chromatography-fluorescence detection (HPLC-FLD). The results indicated that 92% of poultry ration, 72% of dairy ration, 66% of linseed cake, 64% of noug seed cake, 61% of cotton seed cake, and 36% of soybean cake were contaminated with AFB1 exceeding the permissible level (20 µg/kg) in Ethiopia. To evaluate fungal diversity, genomic DNA was…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6- —https://doi.org/10.13039/501100015815Consortium of International Agricultural Research Centers
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsMycotoxins in Agriculture and Food · Indoor Air Quality and Microbial Exposure · Marine Toxins and Detection Methods
Introduction
Ethiopia is the second most populous country in Africa, with approximately 22.17% of its population living in urban areas (UN 2024). The demand for livestock products has been increasing due to population growth and increased urbanization (Menale 2019), leading to a corresponding rise in urban and peri-urban livestock production and a subsequent high demand for feed. Feed is a crucial commodity and a fundamental pillar in achieving the economic, social, and environmental goals of livestock production (Wendimu 2021). The quantity and quality of nutrients, as well as advances in feeding science, significantly impact the economic viability of animal husbandry (Bediye et al. 2018). For animals to achieve their genetically selected production potential, nutrients must be made available and partitioned towards productive functions.
Feed resources traded in local Ethiopian markets include cereals, agro-industry by-products, crop residues, hay, and green fodder (Yosef et al. 2022). Utilizing agro-industry by-products in livestock feed formulation is a common practice. Oilseed cakes, among these by-products, serve as crucial ingredients as protein sources for the growth and development of livestock (Melesse et al. 2023). Oilseeds cultivated in Ethiopia include noug seed, sesame, soybean, linseed and cottonseed (Alemaw and Gurmu 2023). Although Ethiopia has the largest livestock population in Africa (CSA 2020a), livestock productivity is still limited by feed shortages, poor feed quality (FAO 2019), and high feed costs (Yosef et al. 2022). Concentrate feed prices set by the market often disregard feed quality considerations (Melesse et al. 2023). The production, processing, handling, storage, and market-related factors significantly influence the quality and safety of feed (Negash 2020).
Fungal contamination is one of the critical factors affecting feed quality and safety. The diversity of fungi in crops is influenced by pre- and post-harvest conditions, geographical location and weather patterns (Orina et al. 2020). Filamentous fungi pose significant challenges to food security and safety because they produce mycotoxins (Avery et al. 2019). Xu et al. (2022) reported that mycotoxins in animal feeds and food products can cause mycotoxicosis in both humans and animals. Even at low doses, they can disrupt nutrient digestion, absorption, metabolism and physiology. Aflatoxins are considered among the most significant mycotoxins affecting both human and animal health globally (Marshall et al. 2020). These toxic compounds, known for their hepatotoxicity and immunosuppressive properties (IARC 2015), can severely reduce feed intake and animal productivity (Xu et al. 2022). Aflatoxins are produced primarily by Aspergillus species, which thrive in warm, humid environments and can contaminate a wide range of crops. When livestock consume aflatoxin B1 (AFB_1_) contaminated feed, toxic residues may accumulate in animal products such as milk or meat, posing potential public health risks through human consumption (Tolosa et al. 2021).
Nutritional quality of oil seed by-products has been evaluated previously in Ethiopia (Walelegne 2023; Tsehay 2021; and Abdelghani 2008). These analyses provide crucial insights into the composition and suitability of oilseed cakes and other by-products for livestock feed formulations. Concurrently, aflatoxin contamination in these by-products has been assessed (Rehrahie et al. 2018; Gizachew et al. 2016). Despite these investigations, a substantial gap persists in understanding the fungal diversity within livestock feed. Metagenomic studies, focusing specifically on fungal populations in oilseed cakes and formulated feeds, have not yet been conducted. Such research is crucial, as it can reveal the diverse array of fungal species present, including aflatoxin-producing fungi such as Aspergillus flavus, and elucidate their interactions with other fungal communities.
This study aims to provide information on the prevalence of AFB_1_ and aflatoxin producing fungal species, as well as their co-existence with other mycotoxin producing fungal populations, in selected oilseed cakes and formulated feeds used in Addis Ababa and surrounding districts. Understanding the diversity and co-existence of mycotoxin producing fungi in feeds supports targeted monitoring and regulatory measures to reduce fungal toxin risks in Ethiopia’s feed and food system.
Materials and methods
Study area
A cross-sectional study was conducted in Addis Ababa and surrounding districts (Ada’, Mojo, Adama, Addaberga, Wolmera, Sheno, Sebeta, and Sululta) (Fig. 1; Table 1) in 2019. These districts are known for small-scale commercial livestock farming in urban and peri-urban areas. Livestock feed markets, primarily located along main roads for easy transportation access, were also present in the districts. The feed markets were identified and randomly selected for this study.
Fig. 1. Livestock feed collection market locationsSource: Study area mapped by Climate, Geospatial & Biometrics Research Directorate, Ethiopian Institute of Agricultural Research (EIAR)
Sample collection and preparation
Information on the number of feed vendors in each selected district was obtained from the respective administrative focal persons. 50% of feed vendors in each districts markets were randomly selected based on their accessibility. The availability of feed types varied from one vendor to another; this led to unequal sample numbers across feed vendors and districts (Table 1). Therefore, all available feeds of interest were sampled from the feed vendors. As a result, the overall level of aflatoxin contamination was estimated from all feed samples collected across the districts. The final number of samples (290) was therefore informed by vendor distribution, feed availability, and practical field conditions. A kilogram of representative sample of similar feeds from each feed shop was thoroughly mixed according to the procedure (ISO 6497:2002) and put in a double layered kaki paper bag. Each bag was labelled with the type of feed sample, date of collection, vendor code and district of collection and sample number to allow traceability. The samples were transported to the International Livestock Research Institute (ILRI) molecular biology laboratory on the collection days and stored in a cold room (4 °C).
Table 1. Location coordinates and number of feed samples collected for aflatoxin B_1_ (AFB_1_) analysisFeed typesArea nameLatitudeLongitudeCotton seed cakeLinseed cakeNoug seed cakeSoybean cakeDairy rationPoultry rationTotalAda’8°45′21″N39°10′09″E145414937Addaberga9°24′53″N38°25′10″E0020305Adama8°32’30"N39°16’13"E14143127766Addis Ababa9°01’28"N38°44’14"E103561126Mojo8°35’19"N39°7’38’’E21186120773Sebeta8°54′14″N38°34′50″E02826018Sheno9°19′21″N39°17′21″E1111007433Sululta9°10’52"N38°45’20"E016016023Welmera9°3’41"N38°29’35"E0340209Total3853471310138290
Samples were ground in the ILRI Addis Ababa Feed Nutritional Analysis Laboratory, using a CHRISTY (Christy and Norris limited process engineers Chelmsford England) mill and passed through a 1 mm sieve. To avoid cross-contamination between milling, the mill was cleaned between each milling. To remove remaining feed from corners and gaps, a vacuum cleaner (compact, Moulinex 1200 Electronic) was used and then lint-free wipes moistened with 70% ethanol were used to wipe the inside surface. Each ground sample was well mixed and representative sub samples were then taken for AFB_1_ assay and DNA extraction.
Aflatoxin analysis
The Helica ^®^ Aflatoxin B_1_ (low matrix) Enzyme-linked immunoassay (ELISA) for quantitative detection of AFB_1_ was used to determine aflatoxin B_1_ concentrations in the samples. The AFB_1_ extraction and ELISA procedure specifically for animal feed was executed according to the manufacturer’s instructions (Helica Biosystems Inc, Low Matrix ELISA Kit).
AFB_1_ was extracted following the procedure described by Gizachew et al. (2016); a representative subsample (0.1g) of each ground feed sample was added to 10 ml of 80% acetonitrile (High performance liquid chromatography (HPLC) grade acetonitrile, Fisher Scientific, USA). The mixture was shaken in a sealed 15 ml falcon tube for 10 min, centrifuged (2,058 x g) for 5 min, then the supernatant was collected for analysis, an aliquot was diluted 1:10 in reconstituted wash buffer (Phosphate buffered saline (PBS) with 0.05% Tween 20), and a final dilution of 1:1000 was used for the calculation. To validate efficiency of extraction of this method, 30 randomly selected samples were analysed for total aflatoxin using High Performance Liquid Chromatography-Fluorescence (HPLC-FLD) detection.
AFB_1_ was analysed according to the manufacturer’s instructions. 100 µl of each standard and the prepared samples, were mixed with 200 µl of sample assay diluent making 300 µl in the mixing well. Then, 100 µl from each mixing well was transferred to a monoclonal antibody (low matrix specific for AFB_1_ antibody) coated microtiter plate well and incubated at room temperature for 30 min. The liquid was then discarded and the wells drained by tapping on absorbent paper, before the well was washed three times with wash buffer and drained on absorbent paper. Then 100 µl of horseradish peroxidase (HRP), aflatoxin enzyme conjugate, was added to the wells. After incubating for 30 min at room temperature, the liquid was decanted and the well was washed three times with any residual wash buffer drained. 100 µl of substrate solution was added to each well and incubated for 10 min in the dark at room temperature. Finally, 100 µl of stop solution was added to stop the reaction. The optical density of each standard and the samples was measured at 450 nm using a microplate plate reader ELx 800 (BioTek Instrumentals, Inc., USA).
The standard concentrations were plotted along the x-axis on a log scale. The corresponding %B/Bo values were plotted along the y-axis. Reader Fit software (Gen5 a Software Version 2.01.14) was used to fit the standard curve (R^2^ > 0.99). A Data Analysis Worksheet for the Low Matrix ELISA Kit was used to record the AFB_1_ concentration.
To confirm the validity of the AFB_1_ results analysed with ELISA, 30 samples were randomly selected (among which the highest detected samples were included) and sent to Biosciences Eastern and Central Africa (BecA-ILRI Hub) laboratory in Nairobi for HPLC-FLD analysis. HPLC was validated using a low-level standard solution in six replicates. The limits of detection and quantification were determined as 0.27 ppb and 0.91 ppb, respectively; relative standard deviation was 4.3%, and uncertainty 8.7%. Spike recoveries were 86.6–88.5% (mean 87.6 ± 0.98%) which is in the acceptable range (70–120%). The correlation of the results run by HPLC-FLD and ELISA methods was done with linear regression analysis; Coefficient of determination (R^2^), the slope and intercept were calculated using Excel scatter plots.
DNA extraction and sequencing
Fungal genomic DNA was extracted from 50 mg of ground feed samples using DNeasy^®^ Plant Mini Kit Cat. No. 69,104 and 69,106 (QIAGEN, Venlo, The Netherlands), following the manufacturer’s protocol. Each sample was lysed by adding 500 µL Buffer AP1 with 4 µL RNase A; then the mixture was incubated in a Clifton water bath at 65 °C for 20 min with mixing by inverting the tubes three times. Following lysis 130 µL of Buffer P3 was added, the lysate was cooled in ice for 5 min, then centrifuged (eppendorf centrifuge 5424 R) at 20,000 X g for 5 min. The supernatant was put through a QIAshredder spin column placed in a 2 ml collection tube by centrifuging for 2 min at 20,000 X g. The flowthrough was transferred into a 2 ml eppendorf tube and mixed by pipetting with 1.5 ml of AW1 buffer. The mixture was then loaded onto a DNeasy Mini spin column in 750 µL aliquots and centrifuged at ≥ 8000 rpm for 1 min. The column was washed by adding 500 µL of Buffer AW2 and centrifuging for 1 min at ≥ 8,000 rpm, followed by a second wash with 500 µL Buffer AW2 and centrifugation at 20,000 × g for 2 min. To elute the DNA, the spin column was transferred to a clean 1.5 mL microcentrifuge tube, and 100 µL of Buffer AE was added directly to the membrane. After a 5-minute incubation at room temperature, the column was centrifuged at 8,000 × g for 1 min. A second elution using another 100 µL of Buffer AE was performed to maximize DNA yield. Eluted DNA was stored at − 20 °C until further analysis.
The extracted genomic DNA was sent to Genohub, Admera Health, USA for library preparation and sequencing. First, the genomic DNA was quantified using a Qubit 2.0 DNA HS Assay (ThermoFisher, Massachusetts, USA) and the quality assessed using a Tapestation genomic DNA Assay (Agilent Technologies, California, USA). Fifty nanograms of the isolated genomic DNA was used for Polymerase Chain Reaction (PCR) amplification using proprietary primers (Admera Health, LLC, South Plainfield, NJ) covering the hypervariable regions, ITS1 and ITS2. Primer selection and design were undertaken to achieve a comprehensive taxonomic coverage, and spike-in controls were used to generate maximal data. The quantity of the final libraries was assessed using a Qubit 2.0 and quality was assessed by TapeStation D1000 ScreenTape. Illumina^®^ 8-ntdual-indices adaptors were used followed by sequencing on an Illumina^®^ Miseq (Illumina, California, USA) with a read length of 250 paired-end reads per sample (500 K in each direction). The fungal diversity of the top 50 taxonomies at family, genus, and species levels were identified using the QIIME2 platform integrated with a pre-trained Naïve Bayes taxonomic classifier. However, phylum level taxonomic labelling was done using the taxonomy identifier from NCBI using the species as a starting point. Rarefaction measurement was run on the basis of each species (species richness) versus operational taxonomic unit (OTU) abundance.
Statistical analysis
Descriptive statistics were used in Microsoft Excel 2013 (Microsoft Office Professional Plus 2013, Version 15.0) to present the mean aflatoxin quantity in the feed samples. Data were summarized as mean ± standard error (SE), range and median. To compare the significance of AFB_1_ contamination level between the different feed types, Wilcoxon Rank Sum (Mann-Whitney U), R version 4.5.0 was used. For metagenomic data, Qiime2 (version 2021.2) was used to filter the raw data. The final non-chimeric reads were summarized, and Dada2 was used to cluster the OTUs. The taxonomy was classified based on the file: unite-ver8-99-classifier-04.02.2020.qza. The abundance-based coverage estimator (Ace) and Chao1 index (Chao1) were used to provide an estimate of the species richness, while the Shannon, Simpson and InvSimpson indices were used to provide an indication of the species diversity. Beta-diversity was calculated and compared among the six DNA extracts. Finally, bar plots were developed at family, genus, and species levels using the R package “metacoder”.
Results
Aflatoxin B1 (AFB1) contamination of feed samples
AFB_1_ was detected in 92% of the poultry ration, 72% of the dairy ration, 66% of the linseed cake, 64% of the noug seed cake, 61% of the cotton seed cake and 36% of the soybean cake (Fig. 2). The detected concentrations exceeded the maximum allowable level set by European Union (EU) and Ethiopian standard authority (ES ISO14718 Aflatoxin B_1_ µg/kg Max. 20.0). Samples holding AFB_1_ below the threshold of 0.1 µg/kg are below detection limit (BDL) of ELISA assay.
Fig. 2. Percentage occurrence of AFB_1_ in feed samples
The average contamination level for each feed type along the standard error indicated the variability of the contamination in each feed type (Table 2). The results showed that there was a wide range of contamination level in each feed type. The highest median concentration (245.62 µg/kg) was observed in linseed cake showing that most of the samples had the highest AFB_1_ contamination level, followed by poultry ration and noug seed cake.
Table 2. Aflatoxin B1 (AFB_1_) concentration (µg/kg) in oil seed cakes and formulated feedsFeed typeMean ± SERange (Min-Max)MedianDairy ration129.4 ± 11.98BDL − 44286.47Poultry ration162.3 ± 14.34BDL − 293192.29Cotton seed cake108.8 ± 31.87BDL − 79333.54Linseed cake202.1 ± 24.54BDL − 556245.62Noug seed cake162.1 ± 24.79BDL − 518104.58Soybean cake18.9 ± 4.30BDL − 4914.78
AFB_1_ contamination in feed samples was not evenly distributed. However, comparison using the Wilcoxon rank sum test with continuity correction (Table S1) showed differences in AFB_1_ contamination level among samples with significance values ranging from P < 0.001 up to not significant.
To validate the findings, HPLC (gold standard) analysis was conducted on 30 randomly selected ELISA-tested samples revealing a strong correlation between the two methods (r (27) = 0.96, P < 0.01). This correlation was further supported by a linear regression equation (y = 0.7817x + 15.215, R² = 0.9223) (Fig. 3), confirming the reliability of ELISA for AFB_1_ quantification. However, two noug seed cake samples showed very high concentration of AFB_1_ (1708.44 µg/kg and 1549.09 µg/kg). The HPLC conformation test results for these samples were 1232.71 µg/kg and 4954.16 µg/kg, respectively. Due to their deviation from the mean result, they were not included in the summary statistics (mean, range, standard error) shown above.
Fig. 3. Correlation between HPLC and ELISA results of feed samples for aflatoxin B_1_ (µg/kg)
Fungal distribution in the selected feeds
The DNA sequence results showed that Ascomycota was the dominant phylum in all feeds, while the phyla Basidiomycota and Mucoromycota were found in smaller numbers in all the feeds except in feed S4 (Figure S1). The dominant fungal families in the feed samples included Aspergillaceae, Nectriaceae, Botyosphaeriaceae and Cladosporiaceae. Aspergillaceae were represented in 53%, 43%, 16%, 33% 10.6% and 28% of the samples S1-S6, respectively, (Figure S2). The second dominant family was Nectriaceae observed in 23%, 7%, 21%, 3%, 16.6% and 13.6% of samples from the S1-S6 feed types, respectively. The third most represented family was Botryosphaeriaceae which occurred in 4%, 20%, 1.5%, 7.6%, 4.6% and 5.8% of the samples from S1-S6. Cladosporiaceae was the fourth dominating family with up to 2.7%, 11%, 1.9%, 3%, 5% and 4.8% representation in the S1-S6 feeds.
The sequence results revealed that the genus Aspergillus was found in all feed types (Fig. 4; Table 3). This genus was detected in cotton seed cake (S1), noug cake (S2), linseed cake (S4), dairy ration (S6), soybean cake (S3) and poultry ration (S5) in 57%, 41%, 27%, 25%, 15%, and 8% of the samples respectively, (Table 3). Figure 6 shows the bar plot constructed from the abundance of different fungal genera.
Fig. 4. Bar plot showing the average percentage or abundance of each fungal genus in each feed samples in cotton seed cake (S1), noug seed cake (S2), soybean cake (S3), linseed cake (S4), poultry ration (S5) and dairy ration (S6)
Table 3. Fungi genera ranking one to seven in the selected feed samples
The DNA sequence analysis showed that Aspergillus glaucus and Fusarium lacertarum were the most represented fungal species in the feed samples (Figure S3). A. glaucus accounted for 38%, 12.6%, 21%, 20%, 5.7%, and 18% of the fungi species detected in the studied feeds S1-S6, respectively. F. lacertarum was the next most dominant fungal species in the feeds S1-S6 representing 38%, 8%, and 3% 2% 4% and 4%, respectively, of the fungal species detected. The representation of A. flavus was 1.8%, 3.8%, 0%, 0.4%, 1.2%, and 1.2% of fungal species in the feeds S1–S6, respectively.
Evaluation of diversity of fungi in feeds
A principal component analysis (PCA) using unifrac was calculated to analyse and measure the differences or dissimilarity of the types of fungal species present across the six feed types. Dissimilarity in fungal diversity across the six feed samples show component 1 and component 2 explain 30.4% and 26.1% of the variance respectively (Figure S4). Soybean cake (S3) showed the highest level of dissimilarity with a diversity index of 0.4, indicating its difference from the other five feed types. Cotton seed cake (S1) also indicated much more dissimilarity. Linseed cake (S4) showed a moderate level of dissimilarity, not as far away as S3. Other feeds, poultry ration (S5), and dairy cake (S4) showed a moderate level of dissimilarity, not as far away as S3. Other feeds, poultry ration (S5), and dairy ration (S6) indicated they contained fungal communities that are more similar than S3 and S4.
Rarefaction measurement
The rarefaction curve indicated that the sequencing depth was sufficient to include all fungal species, including those present in low numbers. It was run on the basis of each species (species richness) versus Operational Taxonomic Unit (OTU) abundance. Each of the samples contained a different species richness. The feed S1 had the highest species richness followed by S4. Feed Samples S3 and S6 showed medium richness. Feed types S2 and S5 had the lowest species richness (Figure S5).
Evaluation of richness and evenness of fungal dynamics in the different feeds
Different measures of alpha diversity matrix (Cao1, ACE, Shannon, Simpson and Invsimpson) were applied to the Operational Taxonomic Unit (OTU) abundance data of the feeds (Fig. 5).
Fig. 5. Alpha diversity Measures (Cao1, ACE, Shannon, Simpson and Invsimpson) for Cotton seed cake (S1), Noug seed cake (S2), Soybean cake (S3), Linseed cake (S4), poultry ration (S5) and dairy ration (S6)
Alpha diversity indices such as Choa 1 and ACE (species richness), Shannon, Simpson and InvSimpson (species evenness) quantified species richness in descending order. According to the results, Choa 1 indicated that S4, S3, and S5 represented a high species richness; when tested with the ACE S4, S3, and S5 demonstrated richness in the same higher to lower manner. The evenness indices Shannon, Simpson, and InvSimpson revealed that S6, S5, and S4 had the most even fungal dispersion in descending order.
Discussion
Aflatoxin B1 contamination of feeds and feed ingredients
Concerns about agricultural commodities due to fungal incidence and mycotoxin contamination have been widely recognized (Bhat and Reddy 2017). Among these mycotoxins, AFB_1_ has significant negative effects on both animal and human health (Anelich 2014). In Ethiopia, agro-industrial by-products are commonly used as feed ingredients for dairy and poultry production, including oilseed by-products (oilseed cakes) from both modern and traditional processing units (FAO 2019). Oilseed cakes are vital components in the formulation of dairy and poultry feeds due to their high protein content (Gizachew et al. 2016; Mengistu et al. 2017; Menale 2019; Demissie 2022). However, despite their nutritional benefits, oilseed cakes have also been shown to contain aflatoxin B_1_, posing potential risks to both animal and human health. Consumption of contaminated dairy feed can convert aflatoxin B_1_into aflatoxin M_1_ in milk (Pettersson 2012; Campagnollo et al. 2016), and poultry products can also be affected (Oliveira et al. 2000; Tatfo et al. 2022).
In this study, the majority of animal feed samples collected from feed markets were found to contain aflatoxin B_1_. In most cases, the concentration of aflatoxin B_1_ exceeded the maximum limit set by the Ethiopian Standards Authority (ESA 2019; ES ISO 14718: Aflatoxin B_1_µg/kg Max. 20.0), as well as the limits established by the European Union (EU 2002) and Food and Drug Administration (FDA 2019). For dairy animals, the Food and Drug Administration (FDA) sets a maximum AFB_1_ limit of 20 µg/kg, as AFB_1_ present in animal feed can be metabolised into aflatoxin M_1_ (AFM_1_), which can subsequently be passed into milk (Sect. 683.100 Action Levels for Aflatoxins in Animal food).
The study showed that 76% of dairy ration samples were contaminated with aflatoxin B_1_ above the permissible limit, with levels up to 442 µg/kg. This finding aligns with a similar study by Gizachew et al. (2016), who reported contamination levels of up to 419 µg/kg using a similar ELISA procedure. Tadele et al. (2023) found contamination levels as high as 306.9 µg/kg, and Mengesha et al. (2024) reported AFB_1_contamination in dairy feeds up to 370.51 µg/kg. This study also revealed that 92% of the poultry ration samples were contaminated with AFB_1_ above the permissible limit, with concentrations reaching up to 293 µg/kg. Similarly, Mengesha et al. (2024), using a HPLC method, found layer and broiler feeds contaminated with aflatoxin B_1_ at concentrations of up to 139 µg/kg and 148.68 µg/kg, respectively. The findings by Kassaw et al. (2022) reported an average AFB_1_ concentration in poultry feed of 70.11 µg/kg, with a maximum of 633.94 µg/kg. These results supported the findings of this study, indicating that there could be significant risk to poultry health and food safety.
Compared to other feed ingredients in the study, soybean cake exhibited the lowest levels of AFB_1_ contamination, with only 36% of samples showing amounts above 20 µg/kg. The average concentration was 18.9 ± 4.30 µg/kg, ranging from undetectable levels to 49 µg/kg. This result aligns with the findings by Motbaynor et al. (2021) where the mean AFB_1_ range recorded was from 5 to 53 µg/kg. Soybean cake showed the lowest risk of AFB_1_ contamination; in support of this, Mengesha et al. (2024) reported aflatoxin B_1_ in soybean cake at levels below 5 µg/kg. Additionally, aflatoxin concentrations in soybean meal from two regions in Ghana ranged from non-detectable to 3.27 µg/kg, and mostly below the Ghana Standard Agency’s limit of 20 µg/kg for feed ingredients (Nsiah et al. 2023). Gupta and Venkitasubramanian (1975) suggested that the low level of aflatoxin in soybean meal may be due to the small amount of zinc (0.01 mg/g) bound to phytic acid in soybeans. Li et al. (2023) also suggested that phytic acid could act as a plant-derived antifungal agent. In contrast, 66% of the cottonseed cake, 64% of linseed cake, and 61% of noug seed cake exhibited significantly higher contamination rates that often exceeded the permissible limits of AFB_1_. The average AFB_1_ concentration levels detected in the three ingredients were 108.8 ± 31.87, 202.1 ± 24.54 and 162.1 ± 24.79 µg/kg and ranging up to 793. 556, and 518 µg/kg of AFB_1_, respectively. Poultry ration and linseed cake had highest median value and were statistically not different (P = 0.29) and both are significantly more contaminated than dairy ration (P < 0.06). Dairy ration is more contaminated than cotton seed cake (P = 1.015).
Concerning noug seed cake, the result aligns with findings by Gizachew et al. (2016), who reported contamination levels between 290 and 397 µg/kg. Furthermore, Gizachew et al. (2019) demonstrated the impact of water activity and temperature on the growth of A. flavus and subsequent production of AFB_1_, with levels ranging from 203 to 282 µg/kg. More concerning and outlying, AFB_1_ concentration in two noug seed cake samples were also recorded with the levels 1,708.44 µg/kg and 1,549.09 µg/kg, with subsequent HPLC analysis detecting levels of 1,232.71 µg/kg and 4,954.16 µg/kg, respectively. These findings highlight severe contamination risks in the tested samples. In linseed cakes, Ting et al. (2020) reported similarly high levels of AFB_1_ in ground linseed, reaching concentrations of up to 3,644 µg/kg. Cracked, crushed, or ground oilseeds are more susceptible to fungal growth and aflatoxin production due to the removal of the protective seed coat allowing fungi access to readily available nutrients (Ting et al. 2020). While limited studies on cottonseed cake in Ethiopia exist, the wide variability in levels in this study aligns with studies conducted in Pakistan by Shar et al. (2020) and Yunus et al. (2020), who reported concentration levels of 89 µg/kg and 595.9 µg/kg, respectively.
Although samples can be traced to feed vendors, a key limitation of this study is that information was not collected on the pre- and post-harvest practices. Market-based sampling, while reflective of actual exposure risk, does not provide information on the origin or handling history of feed samples. Future studies should prioritize traceable sampling along the supply chain to identify critical control points to realize targeted interventions. The findings in this study serves as a baseline for further assessments and policy development to implement improved record-keeping and supply chain management in Ethiopian feed markets.
Evaluation of fungal diversity in feeds
The fungal DNA sequence showed that Ascomycota constituted the most abundant phylum in all samples. Basidiomycota and Mucoromycota were the next most abundant, but in significantly smaller quantities in all feed samples except in linseed cake where they appeared in a relatively higher abundance. Basidiomycota and Mucoromycota phyla can pose significant risks to both livestock and human health (Petrikkos et al. 2012); in this study, mould species such as Rhizopus arrhizus and Mucor circinelloides were detected under the phylum Mucoromycota. Although, we could not find any evidence that shows the production of mycotoxin, these fungi are known to cause a rare but serious disease, mucormycosis, which is caused by inhaling fungal spores or consuming contaminated feed, resulting in respiratory, digestive, or systemic infections (Petrikkos et al. 2012). Among the species in phylum Basidiomycota, the pathogens Wallemia sebi and Wallemia muriae, which are reported as a human health problem, were identified in this study. The species W.* sebi* and W. muriae cause the condition related with lung issues commonly referred to as “farmer’s lung” (Zajc and Gunde-Cimerman 2018); and the same authors reported fatal livestock toxicosis associated with contaminated hay in the UK. The impact of Ascomycota; however, is particularly concerning due to its association with mycotoxin-producing species, which have direct implications in food safety. The fungal families identified in this study were predominantly from the Ascomycota phylum. Within this group, important genera such as Aspergillus,* Fusarium*,* Penicillium*, and Alternaria were detected (Fig. 4). These genera include species known to produce mycotoxins, which can cause acute and chronic toxicity in both humans and animals (Navale et al. 2021).
The variation in fungal species found in the study may be attributed to different crop cultivation environments (Beccari et al. 2020) and storage conditions (Cao et al. 2022), and differences in nutrient content of the samples, which can influence fungal growth (Gizachew et al. 2019). Other pre-harvest and post-harvest conditions, and feed storage conditions (Zeng et al. 2022; Cao et al. 2022), could also contribute to the observed differences, particularly in A. flavus abundance. Although aflatoxin B_1_ was detected in soybean cake at a low level, A. flavus was not detected in the DNA sample. This could be due to the low abundance of the fungus, making it difficult to detect at the sequencing depth used despite the positive rarefaction result (Sanchez-Cid et al. 2022). Other significant mycotoxin-producing fungi detected in other studies include Fusarium species, such as F. proliferatum (Alshannaq 2017), F. oxysporum, and F. poae (Shabeer et al. 2021), as well as Alternaria alternata, which is known to produce mycotoxins harmful to animals (EFSA 2011). The saprophytic fungi Alternaria and Mucor spp., such as A. glaucus, Mucor circinelloides, and M. racemosus were identified in the feed samples. These fungi have been reported to have the ability to deplete and breakdown organic material, causing loss of nutritional value and feed deterioration (Thomma 2003; Buckle 1983).
In this study, although AFB_1_ was detected in high amounts, Aspergillus flavus was not the dominant fungus (0% − 3.8%) found in the studied feed fungal community. This mismatch can be explained by several factors; one reason could be antagonism or competition for nutrients by other species of fungi identified, for instance like Fusarium. proliferatum and A. flavus (Phan et al. 2021); and bacteria which are naturally found on livestock feeds (Palumbo et al. 2006; Watkinson 2016), although we did not assess the bacterial community. On the other hand, this big variation in abundance could be as a result of raw material handling (Zeng et al. 2022; Cao et al. 2022), nutrient content, storage conditions (Cao et al. 2022) and cultivation area of raw materials (Beccari et al. 2020). Khan et al. (2018) noted that the favourable environmental conditions, especially optimal temperature and humidity can favour the growth of fungi and of the synthesis of their toxic metabolites. As a result, even the small population of aflatoxin producing fungi can build up aflatoxin. Due to their chemical nature, aflatoxins exhibit high stability, (Tahir et al. 2018). Therefore, once formed these toxins persist due to their chemically stable structure, remaining intact even after the fungi themselves have been eliminated. In addition, in the sequencing experiment there could be draw backs of the experiment such as primer bias in targeting ITS regions and also library preparations that cause low coverage of certain species during the sequencing pipeline (Ross et al. 2013; Schirmer et al. 2015). Database annotation issues like fungal ITS databases (e.g., UNITE, GenBank) may have limited or inconsistent annotations for Aspergillus flavus as has been seen generally for Aspergillus species (Steenwyk et al. 2024).
This study demonstrated the possible risks associated with AFB_1_ contamination in oilseed cakes used in Ethiopian livestock feeds, which pose potential health risks to animal and human populations. The fungal contamination observed in the six feed samples aligns with the global trends where genera such as Aspergillus and Fusarium are common contaminants. A. flavus, a powerful pathogen, was identified in all feed samples except soybean cake. AFB_1_ was detected in all samples in the study, with the exception of soybean cake, mostly at higher than the internationally allowable level of 20 µg/kg. Considering the co-existence of different filamentous mycotoxin producing fungi, there is a high probability that the feeds were contaminated with different mycotoxins. Saprophytic fungal Botryosphaeriaceae species which cause feed quality deterioration were also identified in the samples. Therefore, to minimise contamination and toxicity, strict feed management measures at every step of the livestock feed value chain must be enforced. Continual study and observation are vital to understand the dynamics and nature of the contaminating fungi in livestock feed. The findings highlight the critical need for enhanced feed management practices, improved education and training for feed handlers, and stricter regulatory measures to mitigate aflatoxin contamination in Ethiopia.
This study provides insightful information that can serve as a foundation for future research on mycotoxins, feed/food safety studies that could benefit human and animal health. The DNA-based profiling used in the study offers a broad view of the fungal community present in the samples, and the concentration of aflatoxin B1 was also assessed. However, this approach didn’t allow the researchers to establish the direct correlation of the presence or amount of Aspergillus species (viable fungi) to the actual AFB_1_ being produced in the feed. To better understand how Aspergillus flavus behaves in relation to the other existing fungi (especially with dominant fungi species) in the same feed or environment, further research is needed to explore their population dynamics. Adding culture-based methods in future studies could help fill this gap by identifying which fungi are actively producing toxins, and could lead to discovering effective biological control options to reduce aflatoxin contamination.
Supplementary Information
Below is the link to the electronic supplementary material.
Supplementary Material 1
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1United Nations Department of Economic and Social Affairs (UN 2024 Revision) - Population Division. World Population. https://www.worldometers.info/world-population/ethiopia-population/