Buscar en
Revista Argentina de Microbiología
Toda la web
Inicio Revista Argentina de Microbiología Growth characterization and predictive behavior of Eurotium species in a feedstu...
Información de la revista
Vol. 53. Núm. 3.
Páginas 248-256 (Julio - Septiembre 2021)
Compartir
Compartir
Descargar PDF
Más opciones de artículo
Visitas
1775
Vol. 53. Núm. 3.
Páginas 248-256 (Julio - Septiembre 2021)
Original article
Open Access
Growth characterization and predictive behavior of Eurotium species in a feedstuff matrix
Caracterización del crecimiento y comportamiento predictivo de especies de Eurotium en matriz de alimento balanceado
Visitas
1775
Mariana Grecoa,b,
Autor para correspondencia
mariana.greco@gmail.com

Corresponding author.
, Graciela Posea,b, Alejandro Pardoa,b
a Laboratorio de Micología Molecular, Departamento de Ciencia y Tecnología, Universidad Nacional de Quilmes, Roque Sáenz Peña 352, B1876BXD, Buenos Aires, Argentina
b Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Argentina
Este artículo ha recibido

Under a Creative Commons license
Información del artículo
Resumen
Texto completo
Bibliografía
Descargar PDF
Estadísticas
Figuras (4)
Mostrar másMostrar menos
Tablas (2)
Table 1. Experimental and estimated cardinal values of water activity and temperature, growth rate and coefficients
Table 2. Estimated parameters by secondary logistic regression models
Mostrar másMostrar menos
Material adicional (2)
Abstract

Animal feeds are characterized by low water activity values. Nevertheless, fungal contamination with Eurotium species are quite common, causing nutritional depletion, spoilage and economic losses.

The aim of this work was to assess Eurotium amstelodami, E. chevalieri, E. repens and E. rubrum growth in a feed matrix at different conditions of water activity (0.71–0.97) and temperature (5, 15, 25, 30 and 37°C).

It was found that Eurotium species are able to grow in a wide range of water activity and temperature in a short period of time (7 days) and faster than in synthetic media. Rosso and probabilistic models were applied in order to determine the limiting and optimum growth conditions as well as growth probability at certain combinations of environmental factors. Both models provided an accurate fit to the cardinal parameters and good performance for growth/no growth cases.

This is the first report assessing the growth parameters of Eurotium species directly in animal feed. Data obtained in the present study is useful to predict and avoid Eurotium species growth in animal feed.

Keywords:
Eurotium
Feed matrix
Environmental factors
Water activity
Modeling
Resumen

Los alimentos balanceados se caracterizan por tener bajos valores de actividad de agua. Sin embargo, la contaminación por hongos con especies de Eurotium es bastante común y causa agotamiento nutricional, deterioro y pérdidas económicas.

El objetivo de este trabajo fue evaluar el crecimiento de Eurotium amstelodami, E. chevalieri, E. repens y E. rubrum en una matriz de alimento balanceado a diferentes condiciones de actividad de agua (0,71-0,97) y temperatura (5, 15, 25, 30 y 37°C).

Se determinó que las especies de Eurotium pueden crecer en un amplio rango de actividades de agua y temperatura en un corto período de tiempo (7 días), y a mayor velocidad que en medio sintético. Se utilizaron los modelos de Rosso y probabilísticos para determinar las condiciones de crecimiento limitantes y óptimas, así como la probabilidad de crecimiento en ciertas combinaciones de factores ambientales. Ambos modelos proporcionaron un ajuste preciso a los parámetros cardinales y una buena performance para los casos de crecimiento/sin crecimiento.

Este es el primer trabajo que evalúa los parámetros de crecimiento de las especies de Eurotium directamente en alimento balanceado. Los datos obtenidos en el presente estudio son útiles para predecir y evitar el crecimiento de especies de Eurotium en este tipo de alimentos.

Palabras clave:
Eurotium
Matriz alimento balanceado
Factores ambientales
Actividad de agua
Modelado
Texto completo
Introduction

Feedlot practices were introduced in Argentina in the 90s mainly for economic reasons. Nowadays, this practice still continues and at the same time feedstuff production has globally increased. A recent feed survey indicates that worldwide feed production has exceeded a billion metric tons4; 10.6% of feed production in the Latin American region is represented by Argentina34. Particularly, more than 14 million tons were produced in 2015, while the most recent estimates suggest that around 25 million tons will be produced by 202027. Animal feeds are intended for cattle, dairy cattle, porcine, poultry, aquaculture, rabbit and chinchilla, and more recently buffalo. Moreover, an important growth has been registered in the pet food industry in the last years. In addition to industrial feed production, small farmers also elaborate their own feed. However, in some circumstances, lack of quality assurance of raw material, improper technical production process or inappropriate storage conditions can lead to fungal spoilage development and mycotoxin production. In addition, sometimes due to lack of knowledge or negligence, despite the fact that visible fungal growth is observed, feeds are mixed with non-contaminated feed in order to dilute its contamination, thus economic losses are avoided. Therefore, feed quality decreases and becomes a risk for animals, and consequently for the complete food chain37.

Animal feeds are composed of the necessary nutrients that animals need at every stage of their growth to reach a proper development and maximal growth in a profitable period of time for producers. Several parameters are involved and considered to define the quality of feed pellets such as formulation, durability, starch gelatinization and protein content, particle size, exposure to high pressure steam and proper cooling, among others17. In addition, the final product stability is regulated by aw. Therefore, aw is the most critical parameter to be taken into consideration since variations tend to increase fungal proliferation11,29,46.

Fungal contamination in raw materials and animal feed has been widely reported worldwide. In fact, most of them emphasize mycotoxigenic genera, mainly Alternaria, Aspergillus, Fusarium and Penicillium5,16,26,37,42. In addition, based on the findings reported by Hassan et al.26, mixed cereals used as animal feeds are contaminated with high levels of mycotoxigenic fungi and mycotoxins. Moreover, several cases of human and animal intoxication and death caused by food and feedstuff contaminated with mycotoxins were also informed in the recent years18,21,33,38.

The presence and development of fungi cause changes in food and feed properties involving nutritional and organoleptic quality aspects. Due to their xerophilic physiology Eurotium species are able to grow at aw values as low as 0.6512,28. Furthermore, these species are also able to produce a large variety of secondary metabolites10,22,32,43. Therefore, it is of relevance to study the behavior of Eurotium species in order to prevent their growth in animal feed, and their potential production of secondary metabolites.

In recent years mathematical modeling has gained importance since they are considered helpful tools applied to avoid fungal growth and increase food shelf-life.

Previously, we have reported a high incidence of Eurotium species in raw materials and animal feedstuff and an initial study of the effect of aw and temperature on the growth of Eurotium species was performed on synthetic media23. However, considering the high incidence and magnitude of damages that this genus can cause, it is more accurate to study the effect of environmental parameters in a more realistic context thus, in this case in a feed matrix.

Therefore, the aims of this work were (i) to study the effect of aw and temperature on the growth of Eurotium species in a feed matrix, and (ii) to obtain predictive mathematical models for each species: the cardinal model of Rosso for secondary kinetic growth, and the probabilistic model for growth/no-growth interface.

Materials and methodsFungal isolates

Four Eurotium species previously isolated from rabbit and chinchilla feed were used in the present study. Eurotium amstelodami, E. chevalieri, E. repens and E. rubrum strains were previously identified by scanning electron microscopy and DNA sequencing and characterized by secondary metabolite profiling, and growth on synthetic media22,23.

Experimental design and data generation

A full factorial design was used to study the growth of four Eurotium species in a feed matrix. The effect of nine aw and five temperatures were evaluated. The aw assayed were in the range 0.71–0.97 (0.71, 0.73, 0.75, 0.79, 0.83, 0.87, 0.91, 0.95, 0.97), and the incubation temperatures were 5, 15, 25, 30 and 37°C.

Additional aw values were evaluated for some of the species: E. amstelodami at 0.77 and 0.81 aw at 15, 25, 30 and 37°C; E. repens at 0.85 and 0.89 aw and 5, 15, 25 and 30°C.

The maximum aw that the feed could reach without showing excess of water was 0.97 aw corresponding to a moisture content of 60%. Above this value there was surrounding water.

Feedstuff media and water activity adjustment

Commercial a package of rabbit and chinchilla feed was purchased from a local store. The total content was subdivided into 1kg samples, each of which was placed in plastic bags and the total number of bags were placed into the original package made of a double layer of wood paper, and then closed again. Immediately, the bag sample was sent to the Semi-Industrial Irradiation Plant (PISI) in the Ezeiza Atomic Center (Buenos Aires province, Argentina) where it was irradiated with 35kGy of gamma irradiation in order to achieve sterilization. Initial moisture content and aw were 11.06±0.43% (humid basis) and 0.687±0.014, respectively. Then, an adsorption isotherm was built. First, small samples (50-100g) were adjusted to the required aw, and then were scaled to 1kg sample size. Each kilogram of sample was put into a 3 liter glass jar previously sterilized in autoclave at 121°C for 15min and dried in a drying oven at 105°C for 3h. Later, appropriate volumes of sterile distilled water were added in aseptic conditions, mixed and stored at 4°C for 7 days with periodic mixing. Moisture content was determined according to the AOAC Official Method6 and aw was measured with LabSwift-aw Novasina. Each determination was done by triplicate, jars were vigorously shaken and samples were taken from different points. The feedstuff media were aseptically transferred into the Petri plates, a thin layer of feed was placed into them (10–15g), and previous to inoculation the plates were placed at each final incubation temperature for 24h in order to reach the incubation temperature40.

Inoculum, inoculation and incubation

Heavy sporulating cultures of each strain were obtained after seven days of growth at 25°C on Czapek yeast extract agar with 20% of sucrose (CY20S). Conidia were collected from the plate surface and suspended in 5ml of a water/glycerol solution adjusted to the required aw level and containing 0.05% v/v of Tween 80. Conidia concentration was determined with a Neubauer chamber and adjusted to a final concentration of 105conidia/ml. Petri plates containing the feed substrate were centrally inoculated with 1μl of the conidia suspension using a calibrated loop. Petri plates with the same aw were enclosed in polyethylene bags with plastic glasses containing water/glycerol solutions of the same aw in order to keep constant the relative humidity inside the bags and plates; water/glycerol solutions were renewed weekly. Incubation temperatures were 5, 15, 25, 30 and 37°C and the experiment had a duration of 40 days and a maximum of 90 days for the most unfavorable conditions. For each combination of aw*T*species ten plates were used. Additional non-inoculated plates with adjusted aw were placed into the corresponding bags to asses any aw deviation during the study.

Growth assessment

Plates were observed twice a day with a binocular magnifier until growth was detected. Then, two perpendicular diameters of the colony were measured daily without opening the plates, until the colony reached the edge of the plate or until the end of the experiment.

Growth data treatment and mathematical analysisData treatment

Colony diameters were plotted against time, then linear regression was applied to obtain the growth rate (mm/day) as the slope of the line. All obtained data for each pair of (aw, T) conditions and replicates were transformed to obtain homogeneity of variance. All data were square root-transformed.

Secondary model

Secondary models are used to describe the influence of environmental factors on fungal growth39. These models are of great importance and fundamentally very useful because they can be used to predict the microbiological shelf life of a food product35. In particular, in the Rosso cardinal model all of the studied parameters have a physiological meaning19. The data set obtained from this study was tested with other secondary models (data not shown), however the best fit was obtained with the Rosso cardinal model.

Then, the influence of aw and temperature on fungal growth was studied. The Rosso cardinal model was applied to describe the effect of both environmental parameters on the radial growth rate. Each model performance was assessed by the root mean square error (RMSE) and coefficient of determination (R2).

The equation describing the model is:

The equation was fitted using the nls function (nonlinear least squares) from the R Stats Package, in R (version 3.4.0).

The parameter aw max was fixed at 0.97 since it was the maximum value that the feed can reach by addition of water without having free water/surrounding water. It is significant to mention that above 0.95 aw values there are notorious changes in feed appearance (pellet disintegration). Moreover, according to Sautour et al.41awmax should be set to one only for non-xerophilic fungi.

Growth boundary model

Probabilistic models are used when the aim is to determine if growth can happen or not under specific conditions19. The information provided by these models is very useful because it allows to predict and avoid fungal (or bacterial) growth or even toxin production.19 Then, this model instead of determining growth rate, determines whether or not growth occurs.

The growth/no growth boundaries at different aw, temperature and time were built according to the linear logistic regression analysis. A value of 1 was assigned when growth was registered and 0 value when no growth occurred. A total of ten replicates were used.

The equation for the logistic regression model is:

where P is the probability of growth (range 0–1), bi are the coefficients to be estimated, t (days) is the incubation time, aw is the water activity of the feed medium, and T is the incubation temperature (°C).

The equation was fitted using the glm function (generalized linear models) from the R Stats Package, in R (version 3.4.0).

The predicted growth/no growth boundaries interfaces for p=0.5 were calculated and plotted with Microsoft Excel 2007 Solver and Octave (GNU Project, Version 4.2.1).

The predictive power/capacity and the goodness of fit of the logistic models were assessed by pseudo R2 indices (McFadden and Nagelkerke), percentage of concordance and accuracy3,14.

Performance of probability models at probability level of 0.50 was additionally evaluated by comparison with independent data obtained from studies carried out in analog media1,2,24,30,44.

ResultsGrowth description in feed

Growth of E. amstelodami, E. chevalieri, E. repens and E. rubrum was studied at different aw and temperatures.

E. amstelodami grew in the range 0.77–0.97aw at 25 and 30°C. At 15°C the range was 0.81–0.97aw; and at 37°C the range was shorter, from 0.91 to 0.95 aw. The minimum growth rate was 0.19mm/day at 0.77aw and 25°C, while the maximum growth rate was 18.3mm/day at 0.91aw and 30°C. High growth rates were observed in the range 0.87–0.95 aw in the temperature interval of 15–37°C. Nevertheless, optimal aw and temperature were 0.91 and 30°C, respectively.

E. chevalieri growth was observed in the range 0.79–0.97aw from 15 to 37°C and at 0.75aw at 30°C. The lowest growth rate was 0.13mm/day at 0.75aw and 30°C, while the highest growth rate was 23.6mm/day at 0.91aw and 30°C. High growth rates ranging from 10.5 to 15.8mm/day were also observed in the intervals of 0.83–0.91aw at 25–37°C temperature range.

For E. repens, growth was registered from 0.91 to 0.97 aw at 5, 15, 25 and 30°C. Minimum growth rate was 0.29mm/day at 0.97aw and 5°C, whereas maximum growth occurred at 0.91aw and 25°C (20.16mm/day). At 15 and 30°C, high growth rates also occurred at both 0.91 and 0.95aw (10.2mm/day at 15°C and 16.7mm/day at 30°C; 9.3mm/day at 15°C and 6.3mm/day at 30°C, respectively).

With regard to E. rubrum, growth was observed in the interval of 0.75–0.97aw. At 0.75aw only grew at 25°C, and growth at 5°C only occurred at 0.91aw. From 0.79 to 0.97aw growth took place at 15, 25 and 30°C. The optimal aw values ranged from 0.83 to 0.95 at 25 and 30°C with growth rates above 10mm/day, except at 0.95aw and 25°C. The maximum growth rate was 15.9mm/day (0.91aw at 25°C) and the minimum growth rate was 0.16mm/day (0.75aw and 25°C). No aw upper limit of growth was observed.

It has been observed that the rate of growth in all species studied decreased at aw above 0.91 at all the temperatures tested.

Secondary model

Table 1 shows the cardinal values of environmental factors obtained using the secondary cardinal model. The RMSE was calculated to assess the performance of the predictive models; in the calculation no-growth conditions next to the growth ones were included31. In general terms, the models showed a good fit with R2 ranging from 0.885 to 0.950. Nevertheless, the RMSE values were above 1, indicating a poor fit to experimental data. Despite the fact that all the RMSE values were higher than 1, there was high agreement between the parameters estimated by the model and the experimental data. For aw min, aw opt and aw max, the model provided values that fell within the experimental domain.

Table 1.

Experimental and estimated cardinal values of water activity and temperature, growth rate and coefficients

Species  Parameter  Estimated value±SE  Experimental data  R2  RMSE 
E. amstelodamiμopt (mm/day)  17.76±0.912  18.3  0.908  1.614 
awmax  0.974±0.001  0.97     
awopt  0.913±0.003  0.91     
awmin  0.757±0.002  0.77     
Tmax (°C)  38.93±0.483  37     
Topt (°C)  30.41±0.597  30     
Tmin (°C)  5.09±0.436  15     
E. chevalieriμopt (mm/day)  24.88±6.450  23.6  0.950  2.021 
awmax  0.978±0.001  0.97     
awopt  0.909±0.002  0.91     
awmin  0.744±0.002  0.75     
Tmax (°C)  37.20±0.402  37     
Topt (°C)  34.99±2.127  30     
Tmin (°C)  5.39±0.403  15     
E. repensμopt (mm/day)  10.21±1.275  20.2  0.885  2.139 
awmax  0.990±0.006  0.97     
awopt  0.928±0.010  0.91     
awmin  0.897±0.002  0.91     
Tmax (°C)  36.96±0.249  30     
Topt (°C)  24.44±0.603  25     
Tmin (°C)  2.78±1.810     
E. rubrumμopt (mm/day)  14.50±0.605  15.9  0.896  1.687 
awmax  0.972±0.0006  0.97     
awopt  0.911±0.003  0.91     
awmin  0.727±0.002  0.75     
Tmax (°C)  37.04±0.029  37     
Topt (°C)  26.17±0.497  25     
Tmin (°C)  5.38±0.339     

For three out of four species there was agreement between the experimental aw opt and aw max and those estimated by the model (0.91 and 0.97, respectively). In addition, for Topt and Tmax, the model also showed high concordance with experimental data.

With regard to Tmin for E. repens and E. rubrum, both species were able to grow at 5°C. The estimated value was below 5°C (2.78°C) for E. repens, whereas for E. rubrum, Tmin was estimated slightly above 5°C.

With respect to optimal growth rate (μopt), there was high agreement for three out of four species. Nevertheless, there was a significant difference for E. repens, where the estimation of the model turned out to be twice as that of the experimental value.

Growth/no-growth interface model

Probabilistic models were built for each of the Eurotium species at 7, 15 and 40 days. Table 2 shows the coefficients obtained for each model.

Table 2.

Estimated parameters by secondary logistic regression models

Species  Variable  Estimated value  Standard Error 
E. amstelodamiIntercept  −481.10000  14.66 
t  0.25680  0.01026 
aw  1010.00000  31.45 
T  2.70700  0.10007 
aw2  −550.10000  17.43 
T2  −0.03314  0.001079 
aw*T  −1.16900  0.07698 
t2  −0.00296  0.0001321 
E. chevalieriIntercept  −454.000  17.180 
t  0.435  0.018 
aw  953.100  37.210 
T  1.594  0.054 
aw2  −524.000  20.840 
T2  −0.028  0.001 
t2  −0.004  0.000 
E. repensIntercept  −3220  172 
t  0.211  0.015 
aw  6840  365 
T  −1.63  0.312 
aw2  −3640  195 
T2  −0.0355  0.00166 
aw*T  3.10  0.344 
t2  −0.0017  0.00015 
E. rubrumIntercept  −51.7  2.208 
t  0.1992  0.01027 
aw  39.45  2.091 
T  2.114  0.0874 
T2  −0.03568  0.00110 
aw*T  −0.613703  0.06746 
t2  −0.002345  0.00013 

p-value<0.001.

Developed models include linear, quadratic and interaction terms. There were two exceptions since terms were not significant (p-value>0.001): aw*T for E. chevalieri (p-value=0.0079), and aw2 for E. rubrum (p-value=0.452). Despite the fact that coefficients to assess predictive power and goodness of fit were accurate in both cases, the terms aw*T and aw2 were suppressed in each respective case in order to improve model fit and avoid overfitting (backward elimination of terms).

Figures 1–4 show plots of probability of growth (p=0.50) for temperature and aw at 7, 15 and 40 days of incubation for each species. As it can be seen from the figures, growth occurred at almost all suitable conditions during the first 7 days. In general, a high percentage of logistic model agreement with the experimental data was found. The models showed a concordance percentage above 95 and accuracy above 0.929 (Table S1, supplementary material).

Figure 1.

Growth/no growth interface of E. amstelodami after 7, 15 and 40 days of incubation in a feed matrix.

P(0.50): short dashed line: at 7 days; long dashed line: at 15 days; solid line: at 40 days; white square: growth at 7 days; gray square: growth at 15 days; black square: growth at 40 days; cross: no growth.

(0,07MB).
Figure 2.

Growth/no growth interface of E. chevalieri after 7, 15 and 40 days of incubation in a feed matrix.

P (0.50 short dashed line: at 7 days; long dashed line: at 15 days; solid line: at 40 days; white square: growth at 7 days; gray square: growth at 15 days; black square: growth at 40 days; cross: no growth.

(0,07MB).
Figure 3.

Growth/no growth interface of E. repens after 7, 15 and 40 days of incubation in a feed matrix.

P (0.50): short dashed line: at 7 days; long dashed line: at 15 days; solid line: at 40 days; white square: growth at 7 days; gray square: growth at 15 days; black square: growth at 40 days; cross: no growth.

(0,07MB).
Figure 4.

Growth/no growth interface of E. rubrum after 7, 15 and 40 days of incubation in a feed matrix.

P (0.50): short dashed line: at 7 days; long dashed line: at 15 days; solid line at 40 days; white square: growth at 7 days; gray square: growth at 15 days; black square: growth at 40 days; cross: no growth.

(0,07MB).

Figures S1–S4 (supplementary material) show plot probability of growth vs. temperature for different values of aw and incubation times of 7, 15 and 40 days. These figures were also constructed by using the equations obtained by the model, and are complementary to Figures 1–4.

It is important to notice that the probability curves are almost centered at the experimental optimal temperature of growth for each species (except for E. repens, whose probability curve is centered at a lower temperature), and show how the probability of growth increases with time. It is worth noting that for E. chevalieri and E. rubrum this probability increment is notorious when comparing figures for 7 and 40 days respectively (Figs. S2 and S4).

In addition, the validation of the models was also performed with external data from other authors1,2,24,30,44. This comparison is shown in Table S2 (supplementary material), where the probability of growth of the models is shown on the “P model” column. Despite the fact that some disagreements were observed (8 false-positives and 5 false-negatives out of a total of 70 samples), the models predicted acceptably considering that the validation points were located at the domain-model boundaries.

Discussion

Growth and probabilistic models were developed for Eurotium species isolated from animal feed. It was found that in most of the evaluated conditions of aw and temperature, growth happened in a shorter period of time than in previously assayed synthetic media (MEA)23. Synthetic media are richer than natural feed and food, the latter provide a more real condition while synthetic media might induce an overestimated growth7. However, some discrepancies were found in the present study since in the feed matrix the temperatures required for growth at 0.75aw were 25 and 30°C (for E. rubrum and E. chevalieri, respectively), and no growth was detected at 15°C, while in synthetic media the lowest aw value at which growth occurred was 0.75 at 15 and 25°C. In addition, the temperature range of growth in feed was wider than in synthetic media, while the aw range of growth was more limited. Thus, we consider that the most appropriate choice is to perform the analysis in a culture media as close as possible to the original source from which molds were isolated. With regard to this topic, it was determined that the temperature and aw at which fungal growth occurred in animal feed were more in accordance with several authors1,13,15,47. However, differences might be due to the nutritional composition of the substrates and geographical origin7,20,36.

In relation to the cardinal models, they provided accurate estimations for aw max and aw opt, μopt (except for E. repens), Tmax and Topt for all the four species studied (except Topt for E. chevalieri). On the other hand, the parameters aw min and Tmin were underestimated in all cases with the exception of aw min for E. chevalieri which was closer to the experimental observation.

Rosso and Robinson39 estimated E. amstelodami and E. chevalieri radial growth rate, aw min and aw opt from external data, which was obtained in a laboratory medium at 25°C and glycerol as substrate to control aw (among other solutes). These authors determined a radial growth rate of 2.52mm/day, 0.647 aw min and 0.974aw opt for E. amstelodami; while for E. chevalieri the radial growth rate was 2.664mm/day, 0.705aw min and 0.955aw opt. Once again, our findings differ widely since the radial growth rate in feed was about ten times higher than any informed value, aw min was higher than those and aw opt was lower than the values reported by the authors. However, Guynot et al.25 reported a mean of maximum growth rates estimated by linear regression (at 0.90aw) comparable with the informed growth rates in the present study. In addition, Suhr and Nielsen44 reported that E. repens was able to grow in rye bread medium at 0.80aw (controlled by glycerol) and exhibited visible growth before 30 days. Moreover, these authors reported that the shortest time for spoilage was observed at the highest aw (0.95) within a week44. In our study, E. repens was not able to grow below 0.91aw in animal feed and all the four species were able to grow at the highest aw (0.97) in less than a week's time.

The RMSE and R2 were the parameters used to evaluate the goodness of fit and the fitting capacity of the models, respectively. The first parameter measures the deviation between observed and predicted values, while the second parameter indicates how useful the explanatory variables are in predicting the response variable9. A good fit is indicated by RMSE values close to zero, while a proper fitting capacity (R2) is indicated by values close to one.

In the present work, despite the fact of having inoculated molds directly on the feed pellets, the obtained R2 values were very high (>0.885), comparable to values usually obtained from well-defined synthetic media45. Otherwise, RMSE values were too high to indicate a proper fit (2.139>RMSE>1.614). Similarly, Basak and Guha8 informed an R2 range from 0.980 to 0.991 while RMSE ranged from 3.5 to 5.3 in a Penicillium expansum spore germination study.

Probabilistic models developed for each Eurotium spp. showed agreement between observed and predicted data. This agreement was supported by the high values of the obtained coefficients. The pseudo R2 indices (McFadden and Nagelkerke) were accurate enough to indicate a good fit for each model. In addition, accuracy and concordance percentage for all developed models were also high, ranging from 95.58 to 97.81% and from 0.929 to 0.945, respectively. This agreement was also supported graphically, although probability models failed in some cases. Nevertheless, this kind of information can be usefully applied to prevent Eurotium spp. growth in animal feed during storage. As many of the conditions which favored growth occurred during the first 7 days, a safe storage from E. amstelodami growth will be achieved below 0.77aw regardless of the temperature (p of growth<0.10). To prevent E. chevalieri growth, aw above 0.79 (p of growth<0.10) should be avoided. In the case of E. repens aw should be kept below 0.91, and below 0.73aw for E. rubrum.

In conclusion, the predictive models applied in the present study were able to satisfactorily describe and predict the effect of aw and temperature on growth and the probability of growth of Eurotium species in a feed matrix. Since these models were constructed in a real feed matrix they can be used as predictive tools to improve food safety, optimize food storage and prevent fungal spoilage. Thereby, important economic losses can be avoided and feed quality can be improved over time.

To our knowledge this is the first report assessing the growth parameters of Eurotium species directly in animal feed.

Conflict of interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was financially supported by Universidad Nacional de Quilmes and Agencia Nacional de Promoción Científica y Tecnológica.

The authors also thank Gonzalo Pacheco for programming assistance.

Appendix A
Supplementary data

The following are the supplementary data to this article:

References
[1]
M. Abellana, J. Benedí, V. Sanchis, A.J. Ramos.
Water activity and temperature effects on germination and growth of Eurotium amstelodami, E. chevalieri and E. herbariorum isolates from bakery products.
J Appl Microbiol, 87 (1999), pp. 371-380
[2]
M. Abellana, X. Magrí, V. Sanchis, A.J. Ramos.
Water activity and temperature effects on growth of Eurotium amstelodami, E. chevalieri and E. herbariorum on a sponge cake analogue.
Int J Food Microbiol, 52 (1999), pp. 97-103
[3]
L. Aldars-García, M. Berman, J. Ortiz, A.J. Ramos, S. Marín.
Probability models for growth and aflatoxin B1 production as affected by intraspecies variability in Aspergillus flavus.
Food Microbiol, (2018), pp. 166-175
[4]
Alltech [Internet].
7th Annual Alltech Global Feed Survey.
(2018),
[5]
I.F.M. Almeida, H.M.L. Martins, S.M.O. Santos, M.S. Freitas, J.M.G.N. Costa, F.M.A. Bernardo.
Mycobiota and aflatoxin B1 in feed for farmed sea bass (Dicentrarchus labrax).
Toxins, 3 (2011), pp. 163-171
[6]
AOAC.
AOAC official method 934.01.1999. Moisture in animal feed.
Official methods of analysis of AOAC International, 16th ed., AOAC International, (1999), pp. 4-1
[7]
A. Astoreca, G. Vaamonde, A. Dalcero, A.J. Ramos, S. Marín.
Modelling the effect of temperature and water activity of Aspergillus flavus isolates from corn.
Int J Food Microbiol, 156 (2012), pp. 60-67
[8]
S. Basak, P. Guha.
Modelling the effect of essential oil of betel leaf (Piper betle L.) on germination, growth, and apparent lag time of Penicillium expansum on semi-synthetic media.
Int J Food Microbiol, 215 (2015), pp. 171-178
[9]
V. Bewick, L. Cheek, J. Ball.
Statistics Review 14: logistic regression.
Crit Care, 9 (2005), pp. 112-118
[10]
L. Butinar, P. Zalar, J.C. Frisvad, N. Gunde-Cimerman.
The genus Eurotium – members of indigenous fungal community in hypersaline waters of salterns.
FEMS Microb Ecol, 51 (2005), pp. 155-166
[11]
B. Carter, A. Fontana.
Water activity: the key to pet food.
Application Note, Decagon Devices, (2008),
[12]
J.E.L. Corry.
Relationships of water activity to fungal growth.
Food and beverage mycology, pp. 51-88
[13]
S. Dagnas, B. Onno, J.M. Membré.
Modeling growth of three bakery product spoilage molds as a function of water activity, temperature and pH.
Int J Food Microbiol, 186 (2014), pp. 95-104
[14]
N. Deschuyffeleer, A. Vermeulen, J. Daelman, E. Castelein, M. Eeckhout, F. Devlieghere.
Modelling of the growth/no growth interface of Wallemia sebi and Eurotium herbariorum as a function of pH, aw and ethanol concentration.
Int J Food Microbiol, 192 (2015), pp. 77-85
[15]
K.H. Domsch, W. Gams, T.H. Anderson.
Compendium of soil fungi 2.
Academic Press Publisher, (1980),
[16]
C.N. Ezekiel, W.A. Abia, I.M. Ogara, M. Sulyok, B. Warth, R. Krska.
Fate of mycotoxins in two popular traditional cereal-based beverages (kunu-zaki and pito) from rural Nigeria.
LWT Food Sci Technol, 60 (2015), pp. 137-141
[17]
Farahat MH [Internet].
Good quality feed pellets: Does it make sense? All About Feed.
(2015),
[18]
M.C. Ferreras, J. Benavides, C. García-Pariente, L. Delgado, M. Fuertes, M. Muñoz, J.F. García-Marín, V. Pérez.
Acute and chronic disease associated with naturally occurring T-2 mycotoxicosis in sheep.
J Comp Pathol, 148 (2013), pp. 236-242
[19]
D. García, A.J. Ramos, V. Sanchis, S. Marín.
Predicting mycotoxins in foods: a review.
Food Microbiol, 26 (2009), pp. 757-769
[20]
D. García, A.J. Ramos, V. Sanchis, S. Marín.
Is intraspecific variability of growth and mycotoxin production dependent on environmental conditions?. A study with Aspergillus carbonarius isolates.
Int J Food Microbiol, 144 (2011), pp. 432-439
[21]
F. Giannitti, E. Odriozola, C.A. Margineda, E. Fernández, L. Cámpora, N. Weber, G. Clemente, J.P. García, C.M. Campero.
Leucoencefalomalacia equina por pastoreo de maíz contaminado con fumonisinas en Argentina.
Rev Vet Argent, 28 (2011), pp. 1-6
[22]
M. Greco, M. Kemppainen, G. Pose, A. Pardo.
Taxonomic characterization and secondary metabolite profiling of Aspergillus section Aspergillus contaminating feeds and feedstuffs.
Toxins, 7 (2015), pp. 3512-3537
[23]
M. Greco, A. Pardo, G. Pose, A. Patriarca.
Effect of water activity and temperature on the growth of Eurotium species isolated from animal feeds.
Rev Iberoam Micol, 35 (2018), pp. 39-48
[24]
M.E. Guynot, A.J. Ramos, D. Sala, V. Sanchis, S. Marín.
Combined effects of weak acid preservatives, pH and water activity on growth of Eurotium species on a sponge cake.
Int J Food Microbiol, 76 (2002), pp. 39-46
[25]
M.E. Guynot, S. Marín, V. Sanchis, A.J. Ramos.
An attempt to optimize potassium sorbate use to preserve low pH (4.5-5.5) intermediate moisture bakery products by modelling Eurotium spp., Aspergillus spp. and Penicillium corylophilum growth.
Int J Food Microbiol, 101 (2005), pp. 169-177
[26]
Z.U. Hassan, R.F. Al-Thani, Q. Migheli, S. Jaoua.
Detection of toxigenic mycobiota and mycotoxins in cereal feed market.
Food Control, 84 (2018), pp. 389-394
[27]
InfoCampo [Internet].
(2014),
[28]
J.M. Jay.
Intrinsic parameters of foods that affect microbial growth.
Modern food microbiology, Chapman and Hall, (1992), pp. 38-62
[29]
J.A. Lowe, S.J. Kershaw.
Water activity-moisture content relationship as a predictive indicator for control of spoilage in commercial pet diet components.
Anim Feed Sci Technol, 56 (1995), pp. 187-194
[30]
S. Marín, M. Abellana, M. Rubinat, V. Sanchis, A.J. Ramos.
Efficacy of sorbates on the control of the growth of Eurotium species in bakery products with near neutral pH.
Int J Food Microbiol, 87 (2003), pp. 251-258
[31]
S. Marín, A.J. Ramos, V. Sanchis.
Modelling Aspergillus flavus growth and alfatoxins production in pistachio nuts.
Food Microbiol, 32 (2012), pp. 378-388
[32]
B. Mouhamadou, L. Sage, S. Pérgion, V. Séguin, V. Bouchart, P. Legendre, M. Caillat, H. Yamouni, D. Garon.
Molecular screening of xerophilic Aspergillus strains producing mycophenolic acid.
Fungal Biol, 121 (2017), pp. 103-111
[33]
M.R. Moyano, A.M. Molina, A.J. Lora, J. Mendez, A. Rueda.
Tremorgenic mycotoxicosis caused by Paspalum paspaloides (Michx.) Scribner infected by Claviceps paspali: a case report.
Vet Med, 55 (2010), pp. 336-338
[34]
NutriNews [Internet].
Informe sobre la producción de Alimentación Animal en LATAM.
(2016),
[35]
C. O’Mahony, D.L. Seman.
Modeling the microbiological shelf life of foods and beverages.
The stability and shelf life of food, 2nd ed., pp. 253-289
[36]
E.Z. Panagou, S. Chelonas, I. Chatzipavlidis, G.J.E. Nychas.
Modelling the effect of temperature and water activity on the growth rate and growth/no growth interface of Byssochlamys fulva and Byssochlamys nivea.
Food Microbiol, 27 (2010), pp. 618-627
[37]
R.C. Radziejewska, K. Stuper, T. Szablewski.
Micoflora and mycotoxin contamination in poultry feed mixtures from western Poland.
Ann Agric Environ Med, 20 (2013), pp. 30-35
[38]
R. Rivero, S. Collazo, E. Ugarte, F. Nan, Y. Mazzollini, C. Matto.
Leucoencefalomalacia en equinos en el litoral oeste de Uruguay.
Veterinaria, 49 (2013), pp. 31-39
[39]
L. Rosso, T.P. Robinson.
A cardinal model to describe the effect of water activity on the growth of moulds.
Int J Food Microbiol, 63 (2001), pp. 265-273
[40]
S. Samapundo, F. Devlieghere, B. De Meulenaer, J. Debevere.
Growth kinetics of cultures from single spores of Aspergillus flavus and Fusarium verticillioides on yellow dent corn meal.
Food Microbiol, 24 (2007), pp. 336-345
[41]
A. Sautour, P. Dantigny, C. Divies, M. Bensoussan.
A temperature-type model for describing the relationship between fungal growth and water activity.
Int J Food Microbiol, 67 (2001), pp. 63-69
[42]
P.M. Scott, W. Zhao, S. Feng, B.P.Y. Lau.
Alternaria toxins alternariol and alternariol monomethyl ether in grain foods in Canada.
Mycotoxin Res, 28 (2012), pp. 261-266
[43]
G.J. Slack, E. Puniani, J.C. Frisvad, R.A. Samson.
Secondary metabolites from Eurotium species. Aspergillus calidoustus and A. insuetus common in Canadian homes with a review of their chemistry and biological activities.
Mycol Res, 113 (2009), pp. 480-490
[44]
K.I. Suhr, P.V. Nielsen.
Effect of weak acid preservatives on growth of bakery product spoilage fungi at different water activities and pH values.
Int J Food Microbiol, 95 (2004), pp. 67-78
[45]
C.C. Tassou, P.I. Natskoulis, N. Magan, E.Z. Panagou.
Effect of temperature and water activity on growth and ochratoxin A production boundaries of two Aspergillus carbonarius isolates on a simulated grape juice medium.
J Appl Microbiol, 107 (2009), pp. 257-268
[46]
Timmons RA [Internet].
Water activity as a tool for predicting and controlling the stability of pet foods.
Engormix, (2007),
[47]
K.A. Wheeler, A.D. Hocking.
Water relations of Paecilomyces variotii, Eurotium amstelodami, Aspergillus candidus and Aspergillus sydowii, xerophilic fungi isolated from Indonesian dried fish.
Int J Food Microbiol, 7 (1988), pp. 73-78
Copyright © 2020. Asociación Argentina de Microbiología
Opciones de artículo
Herramientas
Material suplementario
es en pt

¿Es usted profesional sanitario apto para prescribir o dispensar medicamentos?

Are you a health professional able to prescribe or dispense drugs?

Você é um profissional de saúde habilitado a prescrever ou dispensar medicamentos