SciELO - Scientific Electronic Library Online

 
vol.56 número2Selección de un medio de cultivo para el crecimiento de Pichia guilliermondii LEVICA- 27 como aditivo activador de la fermentación ruminalEvaluación de variedades de Cenchrus purpureus tolerantes a la sequía en la región occidental de Cuba índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

  • No hay articulos citadosCitado por SciELO

Links relacionados

  • No hay articulos similaresSimilares en SciELO

Compartir


Cuban Journal of Agricultural Science

versión impresa ISSN 0864-0408versión On-line ISSN 2079-3480

Cuban J. Agric. Sci. vol.56 no.2 Mayabeque abr.-jun. 2022  Epub 15-Mar-2022

 

Biomathematics

ProGas v1.1: Program for the preprocessing and data analysis of in vitro gas production of ruminant foods.

Yaneilys García Avila*  1 
http://orcid.org/0000-0003-0126-6233

Mildrey Torres Martínez1 
http://orcid.org/0000-0001-7942-0195

R. Rodríguez Hernández1 
http://orcid.org/0000-0003-1894-4328

1Instituto de Ciencia Animal, Carretera Central km 47 ½, Apartado Postal 24, San José de las Lajas, Mayabeque, Cuba

ABSTRACT

ProGas v1.1 was developed to support researchers in the preprocessing and data analysis of in vitro gas production of ruminant foods. For the implementation of this program, Microsoft Office Excel 2007 was used. The conversion of data to volume (mL g-1 incOM) was carried out with an established equation. The duration of the gas production phases was estimated by means of the average variation rates and the critical points. A total of six sheets were designed in the program: general, enter pressures, accumulated volume, critical points, graphs and help. It is concluded that the proposed program allowed preparing the experimental data for its subsequent modeling. In addition, the methodology that was applied to estimate the critical points did not require a mathematical model to calculate the duration of the in vitro gas production phases, which facilitated the description of the degradation kinetics of the evaluated foods.

Key words: processing; speed; critical points; phases; fermentation

There are several techniques to study the nutritional value of foods in the rumen. To fulfill the hypotheses of the experimental design, in vivo techniques need sufficient numbers of animals and, therefore, large volumes of plant material (Rodríguez et al. 2020). However, in vitro techniques are less time consuming, require less food as substrate and allow greater control of experimental conditions (Arce et al. 2020). In this regard, Theodorou et al. (1994) proposed a method to study the fermentation kinetics of substrates from the in vitro gas production (IVGP) generated by their fermentation, which is measured in bottles of constant volume and at different times.

There are some programs that specialize in modeling ruminal degradation kinetics. These need a mathematical function to estimate the parameters that describe the IVGP kinetics. This is the case of the NEWAY EXCEL program, version 5.0 (Chen 1997) and the ROMENAL procedure (Correa 2004). Both use the Excel platform and are commonly applied in in situ experiments. One of their disadvantages is that they have a reduced number of mathematical models to describe the IVGP kinetics.

Regardless of the used program, the most common parameters during modeling are: asymptotic IVGP, IVGP rate, Lag phase and inflection point, among others. In relation to the inflection point, it is known that it allows identifying the moment in which the food reached the maximum degradation rate. However, the mathematical models that are most used to describe the IVGP kinetics, having a single inflection point, do not allow knowing the time in which the food stopped degradation. If it is added to this the low availability of computer programs related to IVPG, it is necessary to implement a program specialized in the preprocessing and data analysis of IVGP ruminant foods.

The Microsoft Office Excel 2007 development environment was used, which has a friendly platform, common among researchers, with numerous options, formulas and logical operators that make the implementation process easy. The conversion of the IVGP values in Pascal (Pa) to volume data in mL g-1 incOM was carried out using the equation proposed by Rodríguez et al. (2013):

Volume (mL g1 incOM) = (Pressure (Pa*103) + 4.95)/ (incOM)

where:

incOM

is the incubated organic matter.

The estimation of the speed and acceleration of the IVGP was based on the calculation of the mean variation rates (figure 1) of the i-th observations. This is the slope of the secant line to the function that passes through the abscissa points ti and ti+hi. Particularly, in the IVGP experiments, the h value corresponds to the difference between the consecutive sampling times. It is important to highlight that, when the measurements are not taken evenly spaced and the time between IVGP samples is lengthened, h increases and the approximations are less exacts. However, with these concepts it is possible to have an approximate image of the real speed of IVGP, especially when h is small. Similarly, acceleration can be calculated by performing the same procedure at the approximate IVGP point velocities.

Figure 1 Performance of the IVGP of ruminant foods: speed and acceleration calculated by means of the average variation rates 

To determine the time in which the maximum and minimum speed of IVGP occurred, the procedures for calculating the critical and inflection points were used, as reported by Pico and Alava (2018) and Jones (2019), who state that there is a local maximum of the function f at point a, if: f´(a) = 0 and f´´(a) < 0 and there is a local minimum at point a, if: f´(a) = 0 and f´´(a) > 0. On the other hand, if f´ is derivationable at point a and f´´(a) = 0, there is an inflection point at a.

To locate the critical points in the IVGP speed curve, estimated from the experimental data, the time intervals between ti and ti + hi were analyzed. When the estimated acceleration at ti was higher than zero, and the estimated acceleration at ti + hi was lower than zero, the velocity was assumed to have a local maximum in the interval between ti and ti + hi. On the other hand, when the estimated acceleration at ti was lower than zero and the estimated acceleration at ti + hi was higher than zero, it was estimated that the velocity had a local minimum in the interval between ti and ti + hi. As shown in figure 1, the velocity was maximum when the acceleration intercepted the abscissa axis, and when the acceleration slope was negative. At this time, there was also an inflection point for the IVGP. This methodology, although it constitutes an approximation, can be applied to other types of experimental observations, and does not require a mathematical model to determine the critical and inflection points of the data set.

To access ProGas v1.1, users must have a computer running Microsoft Windows XP or newer versions of the operating system. In addition, they must have the Microsoft Office 2007 work package or higher versions. Then, they must copy and open the ProGas v1.1 file, in Excel format (.xls). When starting the program, the general sheet appears, from which you can see the names of the rest: enter pressures, accumulated volume, critical points, graphs and help. The only ones that can be edited are the general one and the one for enter pressures. Sheets must be used sequentially.

The general sheet is where the user enters the characteristics of the experiment. That is: name of the experiment, number of treatments that comprise it, number of repetitions, replications and times. In addition, the names of each treatment must be entered, with their respective fresh (g) and incubated (mL) matter.

The functionality (enter pressures) allows the researcher to fill in the pressure values measured with the manometer in Pa at the different sampling times. In addition, the user must fill in the gas production of the controls, and enter the hours in which each measurement was made (figure 2a). Once all these values are filled in, the application calculates the accumulated gas production and applies the conversion equation proposed by Rodríguez et al. (2013). For this, ProGas v1.1 takes into account the IVGP of the corresponding controls. The result is shown in the accumulated table of gas production in mL g-1 incOM, belonging to the accumulated volume sheet (figure 2b). This sheet also shows another table with the average value of IVGP (mL g-1 incOM).

Figure 2 a) Sheet to enter pressures of the ProGas v1.1 program. b) Sheet of accumulative volume of ProGas v1.1 program. Outlet table 

The critical points sheet is the special functionality of ProGas v1.1, since it allows a better understanding of the IVGP performance over time (figure 3a). Here the average rates of change are calculated as approximate indicators of those derived from the IVGP. The time at which the critical points were reached is shown, and it is specified if it is a local maximum or minimum. In addition, it provides the two critical points that were most repeated in each treatment. Having an approximation of the time in which the minimum speed of IVGP occurred allows knowing when the food finished the degradation. This can contribute to not extending the sampling time in the experiments, since these estimates are made with the experimental data, and without the need to fit a mathematical model. In addition, they allow better characterization of the food, by providing the duration of each IVGP phase.

The graphic sheet shows the performance of the average accumulated volumes (mL g-1 incOM) of each treatment (figure 3b). Lastly, the help sheet offers an explanation of the application's functionalities to facilitate working with the ProGas v1.1 program.

Figure 3 a) Critical points sheet of the ProGas v1.1 program. Outlet table example b) Graphic sheet of the ProGas v1.1 program. Outlet table example 

The proposed program made it possible to prepare the experimental data for subsequent modeling, especially when they were obtained using the IVGP technique, described by Theodorou et al. (1994). In addition, the methodology that was applied to estimate the critical points did not require a mathematical model to calculate the duration of the IVgp phases, which facilitated the description of the degradation kinetics of the evaluated foods.

References

Arce, W., Rojas, A. & Campos, C.M. 2020. "Determinación del contenido energético de materiales forrajeros a través de la relación entre la técnica de producción de gas in vitro y la ecuación mecanicista del NRC (2001)". Nutrición Animal Tropical 14(1): 13-35, ISSN: 2115-3527. [ Links ]

Chen, B. 1997. NEWAY EXCEL Version 5.0. A utility for processing data of feed degradability and in vitro gas production. Rowett Research Institute. Aberdeen. Scotland. 29 pp. [ Links ]

Correa, H.J. 2004. "RUMENAL: procedimiento para estimar los parámetros de cinética ruminal mediante la función Solver de Microsoft Excel®". Revista Colombiana de Ciencias Pecuarias 17(3): 250-254, ISSN: 0120-0690. [ Links ]

Jones, S.R. 2019. "Students’ Application of Concavity and Inflection Points to Real-World Contexts". International Journal of Science and Mathematics Education 17: 523-544, ISSN: 1571-0068. https://doi.org/10.1007/s10763-017-9876-5. [ Links ]

Pico, U.C.D. & Alava, L.A.R. 2018. "Asistencia de softwares matemáticos en el estudio de cálculo diferencial en estudiantes de ingeniería". Revista Bases de la Ciencia3(2): 61-76, ISSN: 2588-0764. [ Links ]

Rodríguez, R., Lores, J., Gutiérrez, D., Ramírez, A., Gómez, S., Elías, A., Aldana, A.I., Moreira, O., Sarduy, L. & Jay, O. 2013. "Inclusion of the microbial additive Vitafert in the in vitro ruminal fermentation of a goat diet". Cuban Journal of Agricultural Science, 47(2): 171-178, ISSN: 2079-3480. [ Links ]

Rodríguez, R., Ontivero, Y., García, Y., Sosa, D. & Gómez, S. 2020. "Empleo del tubérculo de boniato (Ipomoea batatas L.) y la cepa Lactobacillus pentosus LB-31 como aditivos a ensilajes mixtos para rumiantes". Livestock Research for Rural Development, 32(117), ISSN: 0121-3784. Available: http://www.lrrd.org/lrrd32/7/rodri32117.html. [ Links ]

Theodorou, M.K., Williams, B.A., Dhanoa, M.S., McAllan, A.B. & France, J. 1994. "A simple gas production method using a pressure transducer to determine the fermentation kinetics of ruminant feeds". Animal Feed Science and Technology, 48: 185-197, ISSN: 0377-8401. https://doi.org/10.1016/0377-8401(94)90171-6. [ Links ]

Received: December 02, 2021; Accepted: March 15, 2022

*Email: feygarciaa@ica.co.cu

Conflict of interest: There is no conflict of interest between the authors.

Authors contribution: Yaneilys García Avila: Conceptualization, Investigation, Data curation, Formal analysis and Writing - original draft. Mildrey Torres Martínez: Formal analysis and Writing - original draft. R. Rodríguez Hernández: Data curation, Formal analysis and Writing - original draft.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License