Original Research ARTICLE
Experimental Study and ANN Dual-Time Scale Perturbation Model of Electrokinetic Properties of Microbiota
- 1Key Laboratory for Agro-Ecological Processes in Subtropical Region, Hunan Research Center of Livestock and Poultry Sciences, South-Central Experimental Station of Animal Nutrition and Feed Science in the Ministry of Agriculture, Institute of Subtropical Agriculture, Chinese Academy of Sciences, Changsha, China
- 2RNASA-IMEDIR, Computer Science Faculty, University of A Coruna, A Coruña, Spain
- 3Instituto de Investigación Biomédica de A Coruña, Complexo Hospitalario Universitario de A Coruña, A Coruña, Spain
- 4Hunan Co-Innovation Center of Animal Production Safety, CICAPS, Changsha, China
- 5Institute of Biological Resources, Jiangxi Academy of Sciences, Jiangxi, China
- 6Department of Organic Chemistry II, University of the Basque Country UPV/EHU, Leioa, Spain
- 7IKERBASQUE, Basque Foundation for Science, Bilbao, Spain
The electrokinetic properties of the rumen microbiota are involved in cell surface adhesion and microbial metabolism. An in vitro study was carried out in batch culture to determine the effects of three levels of special surface area (SSA) of biomaterials and four levels of surface tension (ST) of culture medium on electrokinetic properties (Zeta potential, ξ; electrokinetic mobility, μe), fermentation parameters (volatile fatty acids, VFAs), and ST over fermentation processes (ST-a, γ). The obtained results were combined with previously published data (digestibility, D; pH; concentration of ammonia nitrogen, c(NH3-N)) to establish a predictive artificial neural network (ANN) model. Concepts of dual-time series analysis, perturbation theory (PT), and Box-Jenkins Operators were applied for the first time to develop an ANN model to predict the variations of the electrokinetic properties of microbiota. The best dual-time series Radial Basis Functions (RBR) model for ξ of rumen microbiota predicted ξ for >30,000 cases with a correlation coefficient >0.8. This model provided insight into the correlations between electrokinetic property (zeta potential) of rumen microbiota and the perturbations of physical factors (specific surface area and surface tension) of media, digestibility of substrate, and their metabolites (NH3-N, VFAs) in relation to environmental factors.
The rumen environment is characterized by a resident microbial population that rapidly colonizes and digests feed particles, thereby providing fermentation end products that are utilized by the host animal (Paracer and Ahmadjian, 2000). Enhancing feed digestibility by rumen microbiota is highly desirable for improving animal performance (Fox et al., 1995; Krause et al., 2003). Volatile fatty acids (VFAs) as a type of metabolites of feed fermentation in rumen providing metabolic energy for the microbial protein synthesis are directly absorbed by the host intestinal tract (Hungate, 1984). The physical-chemical properties of microbial cells and their environment play a vital role in microbial metabolic processes. Specifically, the specific surface area (SSA) of feed materials contributes to the biological processes of catalysis, adhesion, and digestion for microbes (Christensen et al., 1985; Shida et al., 2013; Yoda et al., 2014) as a vital factor of interface property (Garzón and Sánchez-Soto, 2015). Electrokinetic properties, such as zeta potential (ξ) and electrophoretic mobility (μe), are highly relevant to this biological process (de Wouters et al., 2015). The ξ is quantified through measuring the μe of microbes under an electric field. The relationship between ξ and μe can be expressed as Henry equation (Hunter, 1981; Kaszuba et al., 2010) presented as Equation (1).
Where ε0 represents the permittivity of free space, ερ is the dielectric constant and η is the dynamic viscosity of dispersion medium, μe is the electrophoretic mobility and ξ is the zeta potential. The ξ represents the potential difference between the dispersion medium and stationary layer of fluid attached to the dispersed particle. It implies that ξ is a key indicator of the stability of colloidal dispersions. It has been proven that lower ξ absolute value means low viability in organism cells, which leads to the aggregation of cells in a dispersion medium (Kłodzińska et al., 2010).
In previous works, the experimental outcomes and predictive models have been reported based on perturbation values of cell physical properties (permeability and hydrophobicity) of in vitro rumen microbiota (Liu et al., 2016, 2017a). This study is focused on the perturbation effects on parameters such as ST, SSA, pH, c(NH3-N), VFAs and digestibility (D) of neutral detergent fiber (NDF) under the simulated rumen environment. However, no model has been reported yet to predict the ξ effect involved in time-dependent perturbations on these types of parameters during in vitro ruminal fermentation. One of the extreme difficulties, in this case, is that this is a very heterogeneous system with input variables that have perturbations over time measured on two different time scales.
In this line of thinking, ideas from the Perturbation Theory (PT) can be used to model the previously mentioned problem. The PT models are useful to find an exactly related and simpler solution to an existing sophisticated problem. The main feature of PT is to break the sophisticated issue into “perturbation” or “solvable” components. Our group has developed different PT models that start with a known solution to a problem and add corrections considering the variations of various experimental conditions (cj) (Gonzalez-Diaz et al., 2013). In the context of Cheminformatics, the effect of these deviations has been quantified using Moving Average (MA) operators. MA operators have been used to quantify perturbations in complex systems (Duardo-Sanchez et al., 2014), and nanoparticles (Luan et al., 2014; Messina et al., 2015).
This study is aimed at predicting ξ by fitting a linear model that uses the MA operators of all parameters along with different time scales as input values. In some cases, the linear model hypothesis may fail due to the high complexities of the natural system. In agreement with the previously mentioned concepts, some non-linear Machine Learning methods can be combined with time-series analysis methods, such as two of the most known examples of each area of Autoregressive Integrated Moving Average (ARIMA) and non-linear Artificial Neural Network (ANN). ARIMA and ANN algorithms have been compared to each other in time-series problems (Turias et al., 2008). However, the ARIMA and non-linear ANN can be further integrated to seek hybrid ARIMA–ANN models for different time-series issues (Babu and Reddy, 2014). In addition, the MA component of ARIMA models, developed by Box and Jenkins, also named the Box-Jenkins Operator, is basically the same operator that has been used to pre-process data previous to the ANN analysis (Box and Jenkins, 1968). In some previous works, a combination of ARIMA and ANN models was proposed in order to deal with data applying MA methods (Babu and Reddy, 2014), and to increase the prediction of time-series issues by using the smoothing strategies (Barba et al., 2014).
The present work focuses on the integration of ANN and ARIMA as models with PT operators to seek a new model able to predict the values of ξ (ζ-potential), given a set of perturbations on many input parameters measured on two different time scales. To this end, firstly new experimental measurements have been performed in terms of the perturbations on the values of ξ and μe of the rumen microbiota over time. Perturbations of different experimental physicochemical parameters such as SSA, ST, pH, D, c(NH3-N), VFA, and γ have also been measured over time. The present experimental data and the previous published data have been combined herein to develop a new predictive model over ξ in two-time scales (t1 and t2). ANN regression algorithms have been used with two types of variables as input variables. The first one was the expected measure (EM) of the ξ values for a given set of cj. The second type of variables was the MA operators of the experimental physicochemical parameters (for further details please see the following part). The best ANN time-series model accurately predicts the output ξ values for a set of 31,104 perturbations in a given set of cj over a time span of 6–72 h.
Materials and Methods
Methods and In vitro Culture
The experiment was conducted to determine the effect of rumen culture medium at different ST levels and SSA of substrates on the electrokinetic properties of rumen microbes and fermented in vitro performances, the procedures of in vitro ruminal fermentation was described by our previous reference (Tang et al., 2008; Liu et al., 2016). In short, rumen fluid was equally obtained from 3 goats with rumen fistula before feeding, strained through 4 layers of cheesecloth, and mixed with an anaerobic fermented buffer in a ratio of 1:2 (v/v) under continuous flushing with CO2. The fermented buffer was prepared according to the description of Tang et al. (2006). About 500 ± 50 mg NDF extracted from rice straw was accurately weighed into a screw-cap serum bottle (145 mL) and full flushed with CO2. The fermented bottle, providing the inoculum of 50 mL mixed solution (39°C, pH = 6.9 ~ 7.0) of rumen fluid with buffer and 500 ± 50 mg DNF sealed with rubber stopper in full of CO2, was incubated and gently shaken in an incubator at the constant temperature of 39°C. As the limited of the experimental condition, we conducted the experimental in three times. For each time, we proceeded to all the treatments or combinational conditions (12 combinations = 4 ST × 3 SSA) at 6 or 4 fermented time-points (t1 = 6, 12, 24, 36, 48, and 72 h or t2 = 6, 12, 24, and 48 h after fermentation). The electrokinetic properties (μe, ξ) of the ruminal microbiota and the conventional parameters were determined, along with the molar concentrations of VFAs and ST during fermentation (ST-a, or γ) in two different time series (t1 and t2). The parameters of μe, ξ and γ were carried out on a time scale t1, while VFAs on a time scale t2. The experiment was conducted to measure these parameters under the conditions of 12 combinations with 4 ST (36, 43, 46, and 54 mN/m) × 3 SSA (3.27, 3.73, and 4.44 m2/g) with 3 replicates each. All experiments and animals used were provided by the Animal Care Committee, Institute of Subtropical Agriculture, the Chinese Academy of Sciences, China, with the document number (No. ISA-2012-018).
Construction of Combinatorial Conditions
In this study, an exogenous non-ionic surfactant, alkyl polyglucoside (APG) bought from Hunan Diyuan Co., Ltd., China, was provided to construct the different ST levels. The dosages of APG used in the present work were adjusted as described in our previous reports (Liu et al., 2013). The final ST levels were set to ST1, ST2, ST3, and ST4 with the values of 53.95, 46.09, 42.78, and 36.07 dynes/cm by supplying 0, 0.02, 0.05, and 0.12% (v/v) of APG, respectively. In addition, the percentages of APG used in the present study were similar to the ones from the current animal production (Yuan et al., 2010; Zeng et al., 2011). On the other hand, NDF extracted from rice straw was used as the substrate material for in vitro fermentation. According to rumen gastrointestinal digest particle distribution (Li and Jiang, 2001). NDF particles were ground into three different screen sizes with a grinder. The surface property of NDF particles expressed as specific surface areas (SSA1, SSA2, and SSA3) was determined by Surface Area Analyzer (Quadrasorb-SI, Quantachrome Inc. Florida, CA, USA), with the values of 3.37, 3.73, and 4.44 cm2/g, respectively.
Electrokinetic Properties Assay
In general, the fermentation microbial samples at each time-point were collected by centrifuging and washing with phosphate-buffered saline at least two times. The details for collecting microbial cells were described in our previous works (Liu et al., 2016). Briefly, fermented liquid (10 mL) from each bottle at the particular fermented time-point was centrifuged at 500 × rpm at 4°C for 10 min to remove the feed particles. The supernatant (2 mL) was then further centrifuged at 12,000 × rpm at 4°C for 10 min to collect the mixed microbial particles. The precipitate was full dissolved in 5 mL 4°C KNO3 with the concentration of 1 mM. The solution was centrifuged at 12,000 × rpm at 4°C for 10 min again, removed the supernatant, and then accurately diluted to 5 mL 4 C 1 mM KNO3 for 30 min. The solution (5 mL) was conducted to carry out the electrokinetic properties of RMMs assay. The electrophoretic mobility (μe) of microbial cells was measured according to a phase amplitude light scattering (PALS) described by Hayashi et al. (Hayashi et al., 2003) and Eboigbodin et al. (de Kerchove and Elimelech, 2005), using the Zeta Potential Analyzer (NanoBrook ZetaPALS, Brookhaven Instruments Corp., Holtsville, NY, USA) (Instruments, 2012) equipped with a He-Ne laser (658.0 nm) as a source of light in a high precision system (30 cycles) at 39°C. The ξ was calculated according to equation (1) by the Zeta potential analyzer.
Surface Tension during Fermentation (γ)
The ST of inoculum in the fermentation process was determined on a time scale t1. The ST of the ruminal culture medium at each point time was measured immediately using the surface tension analyzer (model K100 Tensiometer, KRÜSS GmbH, Hamburg, Germany) (Blake et al., 1957). This ST was set in the actual fermentation processes as ST-a/γ.
Determination of Volatile Fatty Acids
In current study, VFAs include acetic acid (C2:0), propionic acid (C3:0), butyric acid (C4:0), isobutyric acid (isoC4:0), valeric acid (C5:0), and isovaleric acid (isoC5:0). Every 2 mL of incubation fluid from each fermentation bottle on a time scale t2 were centrifuged at 10,000 × g and 4°C for 15 min. 1.5 mL of supernatant solution was collected, which was immediately mixed and homogenized with 0.15 mL 25% metaphosphoric acid. The mixture solution was centrifuged again at 10,000 × g at 4°C for 15 min, and the supernatant solution was collected to determine each VFA content by gas chromatography (HP5890, Agilent 5890; Agilent Technologies Co. Ltd, USA). The DB-FFAP column (Agilent, No.: 122-3232, 30 m in length with a 0.25 mm i.d. and 0.25 um thickness) was used. The parameters of this column were set as the attenuation in a nitrogen split ratio of 1:50, hydrogen flow 30 mL/min, airflow 365 mL/min, injector temperature 250°C, column temperature 150°C and detector temperature 220°C with N2 as carrier gas at a flow rate of 0.8 mL/min. The relative response factor, represented as the peak of each VFA, was calculated against a standard VFA mixture analyzed following every 10 measurements.
Database Construction and Modeling
Original Data Resources
In this work, the ξ and μe values of ruminal mixed microbes (RMMs), VFA and γ were obtained with the same initial combinatorial conditions of ST × SSA. The recorded data of pH, digestibility (D) of NDF, and c(NH3-N) were also collected from our previous works (Liu et al., 2016). All data came from two different time scale series (t1 and t2) to establish a predictive model. More specifically, for the variables of ξ, μe, γ, D, pH, and c(NH3-N), the time series was t1. For the variables of VFAs (six individual VFA), the time series was t2. The details for all data resources used to develop the predictive ANN Time Series model are shown in Figure 1. The full dataset is provided online (Liu et al., 2017b) (SM01_Variables_t1.pdf, SM02_Variables_t2.pdf, and SM03_Zp_model_dataset.xlsx).
Details of Dataset Construction and Modeling
The experiment considered the initial combinatorial conditions made up of 4 ST × 3 SSA = ST1-SSA1, ST1-SSA2 …ST4-SSA3 for all the variables, such as (a) experimental variables: ξ, μe, γ, and VFA; (b) record data: pH, c(NH3-N), and D. However, the experimental variables were measured with a different time span of 6–72 h in t1 and 6–48 h in t2. More specifically, all variables combined with the time series can be expressed as ξ(t1), μe(t1), ST-a(t1), pH(t1), c(NH3-N)(t1), and D(t1) in time series of t1; VFA(t2): C2:0(t2), C3:0(t2), C4:0(t2), isoC4:0(t2), C5:0(t2), and isoC5:0(t2) in time series of t2. Therefore, in order to study the effect of the perturbations of all variables over zeta potential (ξ) of RMMs, an integrated input dataset, made up of combinations of all the variables presented in these two different time scales (t1 and t2), needs to be assembled. There were two blocks of experimental data: (n1) No.(t1) = 12 initial cj × 3 replicates × 6 point-of-time = 216 cases; (n2) No.(t2) = 12 initial cj × 3 replicates × 4 point-of-time = 144 cases. In this work, an integrated dataset made up of No. = 31 104 (= 216 No.(t1) × 144 No.(t2)) cases was constructed, combining all variables in these two blocks of data on different time scales. The ξ value was standardized using the mean for all ξ. Next, the ξ-ANN models were developed to predict the numerical simulation ξ of rumen microbes over the environmental variables or factors in rumen fermentation ecosystem. For further details on n1 and n2 stocks experimental data (please see the online files SM01_Variables_t1.pdf, SM02_Variables_t2.pdf, SM03_Zp_model_dataset.xlsx, and SM04_Zp_dataset_std_values.sta) (Liu et al., 2017b).
In this work, it was firstly assumed that the zeta potential of RMMs might associate with the environmental factors to some extent. Herein, it was assumed that a numerical model used to map the physicochemical properties of RMMs might be beneficial to reflect its mechanism or function. Therefore, the environmental factors such as pH, D, c(NH3-N), γ, and VFAs were taken into consideration. The variable symbols Vq(t1) and Vf(t2) were defined according to two different time scales, where the subscript “q” and “f” indicate the different input variables corresponding to the different time series, respectively. Vq(t1) or Vf(t2) indicates that the input variables change with the corresponding time scale. As the variability or the “small” change in different conditions cj, we introduced the deviations (perturbations) of each input variables based on the corresponding expected values (data dispersion). Therefore, the following types of terms were introduced. The first type is the variable < ξ>, which is the Expected Measurement (EM) component used to account for the expected value of the output property marked as ξ(e), where the subscript “e” is the expected measurement. This theoretical notion has been widely used in our previous works (Ran et al., 2016). The other type refers to the Box-Jenkins Operators (or perturbation values) ΔVϕ(tk), tk representing the time scale t1 and t2, while Vϕ is Vq or Vf, respectively. In addition, the MA component was used to account for the dispersion of variables (Gonzalez-Diaz et al., 2013). In the present work, the terms < Vϕ(tk)> are the MA of the variable Vϕ(tk) in one of the original combinatorial experimental conditions with the format of 4 ST × 3 SSA.
In this sense, this is a complex theoretical model which combined the methods of the PT, MA and Time Series Analysis (Fisher, 1936). The general formula of the proposed model is shown as follows:
In this general formula, symbols such as a0, bq, and cf, are the intercept, or the coefficient of the corresponding variable. In this specific case, the ANN model is a linear additive function of all the terms mentioned above. The terms of ΔVϕ(tk) were added as corrections to the EM of zeta potential ξ(e) = a1. < ξ(tk)>, where a1 refers to the weighted value of < ξ(tk)>, and < ξ(tk)> refers to expected value of ξ calculated as the standardized value of ξ based on the mean value of each cj.
Next, the expected measurements < ξ> and its corrections ΔVϕ(tk) were obtained for the entire dataset. First the different ANN algorithms (Fisher, 1936) were trained/validated in a ratio of 3:1 (training/validation) to search for the best ξ-ANN model with STATISTICA software version 6.0 (STATISTICA, 2002). Three types of ANN algorithms were chosen to search for the best regression prediction model: Linear Neural Networks (LNN), Multilayer Perceptron (MLP), and Radial Basis Functions (RBF). In this sense, it was mandatory to use all factors, including the EM variable and all the corrections ΔVϕ(tk), such as the input variables, in order to develop the predictive ANN regression models. In so doing, our group attempted to find the contribution weighted value of each variable in the zeta potential of RMMs. Thus, the back propagation with 100 epochs and Conjugate Gradient Descent with 500 epochs were used for retraining. The details of experimental and theoretical sections used to develop the predictive model for zeta potential were presented in Figure 2.
The experimental results of ζ-potential (ξ), μe (μe) and ST during fermentation (expressed as γ) in t1 time series are shown in Table 1. In order to understand the data errors in this table, we decided to include two different boxplots (for raw and standardized data) in the Figshare repository of this paper (Liu et al., 2017b). All the ζ-potential (ξ) values of microbes were ranged from −40 mV to −20 mV. In the present study, the ζ-potential (ξ) value was not affected by SSA, ST and their interaction (P > 0.05), while it changed along the incubation time (t) (P < 0.01) by the significance analysis. In addition, the effects of incubation time (t) on ζ-potential depended on its interaction with SSA of the substrate (SSA × time interaction, P < 0.05). For SSA at 3.37 m2/g, the increase of ξ from 12 to 24 h, and decrease from 24 to 48 h was greater than that of ζ-potential at 3.73 m2/g. The ST of rumen fluid could be decreased by adding a non-ionic surfactant, and a moderate decrease of ST could increase the permeability of the microbial membrane, which also was noted in the present study. The values in Table 1 are the mean of the three experiments.
The VFA values of the present study were shown in Table 2. In order to increase the interpretability of the information in this table, we decided to include two different boxplots (for raw and standardized data) in the Supplementary Materials of this paper (Figshare repository) (Liu et al., 2017b). The total VFA concentration and the molar percentages of most individual VFA were not affected by SSA, except for the molar percentage of propionate, which was increased (P < 0.05) with the increasing of SSA. By increasing the SSA of the substrate, the percent of propionate increased, while the ratio of acetate/propionate decreased. The ratio of acetate to propionate decreased (P < 0.05) with the increasing of SSA. Total VFA concentration and the molar percentage of isovalerate decreased (P < 0.01), but the molar percentage of propionate increased (P < 0.01) with increasing ST. By decreasing the ST of rumen inoculums, the percent of propionate decreased, but isovalerate increased, this phenomenon indicates that both of the SSA of the substrate and the ST of rumen fluid could change fermentation patterns of the fiber. The values in Table 2 are the mean of the three experiments. Thus, four plots (Raw_data_Errors_Table 1.png, Raw_data_Errors_Table 2.png, Standardized_data_Errors_Table 1.png, Standardized_data_Errors_Table 2.png) and an Excel file with the plot data (Data_errors.xlsx) were included in the Figshare repository (Liu et al., 2017b).
Table 2. The VFAs concentrations with the designated initial combinatorial conditions (3 SSA × 4 ST) in time series t2 (h).
In the previous section, the output values of ξ (ζ-potential) for rumen microbiota were conducted to measure after the perturbations of the initial values of ST and SSA in this experiment. We separated all data into two sets of input variables, measured on two different time scales under the same conditions cj. In details, the physicochemical parameters per se and D were measured on the time scale t1 and the concentrations of VFAs on the time scale t2. Next, our group started to develop a general predictive model based on the data from the present work and from our previous studies (Liu et al., 2013). The ANN topologies tested included LNN, MLP, and RBF. The best model found was the RBF 14:14-794-1:1; which is able to predict the values of ξ for more than 30,000 cases with a correlation coefficient > 0.8, see details in Table 3. The models exported from STATISTICA can be downloaded as SM05_Zp.snn (Liu et al., 2017b).
Table 4 depicts the sensitivity ratio value (Var Ratio) for each input variable of the corresponding ANN model. Evidently, the value of the Var Ratio provides extremely important information for a predictive model; a greater value of Var Ratio provides a more important role in the predictive model. In addition, when the value is greater than 1, the corresponding variable plays an important role in the model, on the other hand, less than 1 indicates no or a little role for the model.
Thus, for the case of the RBF model, all variables have sensitivity ratios >1. It shows that the perturbations on the variables measured on both time scales (t1 and t2) have a significant effect on the final value of ξ (ζ-potential) to some extent. In the context of physiological properties of RMMs, our best RBF model obtained further proof of the vital role of combinatorial conditions with the different levels of SSA and ST on the electrokinetic properties of RMMs. In addition, this new model showed that the fermentation time and digestibility play the most important role in zeta potential property of RMMs. However, the pH value and ST-a (γ) in the time series of t1 showed the lowest influence role in zeta potential of RMMs. The concentration of metabolites, such as ammonia nitrogen and VFAs, also plays an important role in this electrokinetic property for rumen microbiota. In particular, the role of the propionic acid is higher than that of other metabolites (NH3-N and other components of VFAs) for the ζ-potential of rumen microbiota in this case.
The cellulose degradation is highly correlated to the activity and quantity of fibrolytic enzyme. The previous work had proved that the inclusion of a non-ionic surfactant in the culture media resulted in a significant increase in the production of extracellular celluloses and β-xylosidase in vitro production (Long and Knapp, 1991). Despite the fact that rumen microorganisms can produce some bio-surfactants to decrease the surface tension of rumen liquor, but the decreased degree is limited (Liu et al., 2013). The results obtained from the in vitro study implied that the ST could be regulated to an optimum level for rumen microorganisms to grow better. On the other hand, the material surface roughness (in particular the specific surface area property), surface stiffness and softness are the other key factors influencing cell adhesion and behavior, such as the cellular morphology, proliferation, and phenotype expression (Engler et al., 2006; Chang and Wang, 2011).
The adhesion of ruminal microbiota to the substrate is the prerequisite for bacterial colonization and proliferation (McAllister et al., 1994). The electrokinetic properties (ζ-potential and μe) are considered as the key factors in the adhesion of microorganisms on biomaterial surface (Ascencio et al., 1995). ζ-potential indicates the degree of repulsion between adjacent and similarly charged particles in the dispersion, and maintains the stability of the suspension or dispersion (Van der Biest and Vandeperre, 1999). The Zeta potential can be used as a parameter for characterizing the physicochemical properties of the bacterial cell envelope which are often used to characterize the adhesiveness of bacteria and biofilm formation (Cieśla et al., 2011). In the present study, it was found that rumen microbiota imparted negative ζ-potential values to keep the electrosteric stabilization from incipient instability (ζ-potential range from ±10 mV to ±30 mV) to moderate stability (±10 mV to ±30 mV) in suspensions at all SSA and ST levels (Hanaor et al., 2012). The cell surfaces of the rumen microorganism were negatively charged, which was consistent with the previous report (Delgado et al., 2005). It implies that the Gram-negative bacteria were the main microorganism in the medium (Zhang, 2008). In the present study, the ζ-potential of rumen microorganism unaffected by ST, which was inconsistent with the previous studies, while the ζ-potential changed with increasing bio-surfactant and non-ionic surfactant concentration for single cell or carbon black powers (Hua et al., 2003; Zhang et al., 2009; Cheng, 2014).
In the present study, the results also showed that the ζ-potential of rumen microorganisms increased with increasing fermentation time up to 24 h, and the increasing charge density on the microorganism surface decreased the resistance of microbe adhering to substrates (Li et al., 2014). The negative charge density of cell or bacteria were more adhesive to the positive-charged material (e.g., amine group-grafted polyethylene) compared to neutral- or negative-charged materials in the case of Chinese hamster ovary (CHO) cells, leukocyte adhesion, phagocyte migration, and osteoblast differentiation, etc. (Keselowsky et al., 2003; Chang and Wang, 2011). In addition, the incorporation of negative charges for cell/bacterial organisms may facilitate the adsorption of proteins which promote cell adhesion and responses (Thevenot et al., 2008) and promote alkaline phosphatase enzymatic activity (Keselowsky et al., 2005). In this study, the results may indicate that the main adhesion of microbiota to the substrates occurred within 24 h. In addition, the decreasing of the microbial ζ-potential after 24 h may be due to the decreasing of bacterial activity, or the agglomeration of bacteria as the population increasing. Soni et al. (2008) reported that the number and the shape of peaks obtained on electropherograms, which reflected the number of small aggregates originating from the individual or clustered bacterial cells, were negatively correlated with the zeta potential. Furthermore, in the late fermentation phase, dead bacterial cells may have also contributed to a lower zeta potential (Soni et al., 2008). In addition, the surface charge of bacteria was also significantly affected by the growth medium, bacterial age and bacterial surface structure (Delgado et al., 2005).
On the other hand, the VFAs are the secondary metabolites of rumen carbohydrate degradation. The iso-acids are required absolutely to stimulate the growth of several rumen bacterial species, particularly for the fibrolytic organisms (Allison, 1963; Bryant, 1973). Our results implied that the surface tension property of media and the SSA characteristic of the substrate could directly influence the composition and distribution of VFAs metabolites of rumen microbiota. This result is consistent with some other previous reports, such as the optimal surface tension value might increase the VFA production and microbial metabolic ability in the anaerobic fermentation of waste activated sludge by addition of bio-surfactants (surfactin, rhamnolipid and saponin) (Huang et al., 2015, 2016). In the fermentation of rumen microorganisms, the total VFA concentration was significantly higher with 2% Tween 80 than with no or 1% Tween 80 treatments (Kim et al., 2005). Our results showed that the molar percentage of propionate significantly increased with the ST of media, similar with the previous report by supplementation of Tween 80 (Kim et al., 2005).
In general, the PT method sets with an initial of a known exact solution for an issue or problem by adding the corresponding corrections due to the variations of different experimental conditions (Liu et al., 2015). This idea was widely used to solve many practical problems, such as the fatty acid distribution (Liu et al., 2015), goat growth yield via mRNA expression of ghrelin receptor and growth hormone receptor (Ran et al., 2016), carbon nanotubes (González-Durruthy et al., 2017), and drug-lymphocyte interactome networks (Tenorio-Borroto et al., 2016), etc.
In our previous studies, we reported that the ideas of PT combining with the ANN modeling were successfully used to predict the physic-chemical properties (permeability and hydrophobicity) of rumen microbes based on the combinatorial conditions of the ST and SSA in the processes of in vitro fermentation (Liu et al., 2016, 2017a). However, for another key factor of physic-chemical properties of rumen microbiota, the predictive model focusing on the electrokinetic property (ζ-potential, ξ) has not yet been reported. Thus, the present work focused on developing for the first time ANN models, which evaluate the variations of the ξ (ζ-potential) of the cells, present in the rumen microbiota due to perturbations of many different parameters on two different time scales. The best ANN model (RBF) for ζ-potential (ξ) of rumen microbiota was established with a correlation coefficient >0.8 for > 30,000 cases based on the input perturbations of the physicochemical parameters such as SSA, ST, pH value, digestibility of NDF, c(NH3-N) and VFAs for the dual-time scale. From a physiological view, the initial combinatorial conditions of the SSA of material and surface tension of media, as well as the fermentation time, digestibility of NDF and concentrations of metabolites (NH3-N, and VFAs) in fermentation process are highly related to the electrokinetic property (ζ-potential, ξ) of the rumen microbiota.
YL, CZ, and ZT designed the research; YL, TR, YY, and ST conducted the research; YL, CM, CF, AP, and HG analyzed the experimental original data and constructed the predictive model; YL, TR, ST, CM, HG, and ZT wrote the full manuscript. All authors approved the final manuscript.
The present study was supported by grants from the National Natural Science Foundation of China (Grant No. 31172234 and No. 31260556), the Planned Science and Technology Project of Hu-nan Province (2015NK3041), the Technology Specialty Fund for Cooperation between Jilin Province and the Chinese Academy of Sciences (2016SYHZ0022), and the Hunan Provincial Creation Development Project (Grant No. 2013TF3006). This work was also supported by the General Directorate of Culture, the Education and University Management of Xunta de Galicia (Ref. GRC2014/049), and the “Galician Network for Colorectal Cancer Research” (Ref. R2014/039). CF acknowledges a Juan de la Cierva postdoctoral fellowship grant from the Spanish Ministry of Economy, Industry and Competitiveness (Ref. FJCI-2015-26071).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors are extremely grateful to the Key Laboratory of Subtropical Agro-ecological Engineering, Institute of Subtropical Agriculture, from the Chinese Academy of Sciences (CAS) for providing the experimental materials and apparatuses to complete the present experiment.
Ascencio, F., Johansson, G., and Wadström, T. (1995). Cell-surface charge and cell-surface hydrophobicity of collagen-binding Aeromonas and Vibrio strains. Arch. Microbiol. 164, 223–230. doi: 10.1007/BF02529975
Barba, L., Rodríguez, N., and Montt, C. (2014). Smoothing strategies combined with ARIMA and neural networks to improve the forecasting of traffic accidents. Sci. World J. 2014:12. doi: 10.1155/2014/152375
Chang, H.-I., and Wang, Y. (2011). “Cell responses to surface and architecture of tissue engineering scaffolds,” in Regenerative Medicine and Tissue Engineering - Cells and Biomaterials, ed D. Eberli (Rijeka: InTech), Ch. 27.
Christensen, G. D., Simpson, W. A., Younger, J. J., Baddour, L. M., Barrett, F. F., Melton, D. M., et al. (1985). Adherence of coagulase-negative staphylococci to plastic tissue culture plates: a quantitative model for the adherence of staphylococci to medical devices. J. Clin. Microbiol. 22, 996–1006.
Cieśla, J., Bieganowski, A., Janczarek, M., and Urbanik–Sypniewska, T. (2011). Determination of the electrokinetic potential of Rhizobium leguminosarum bv trifolii Rt24.2 using Laser Doppler Velocimetry—a methodological study. J. Microbiol. Methods 85, 199–205. doi: 10.1016/j.mimet.2011.03.004
de Kerchove, A. J., and Elimelech, M. (2005). Relevance of electrokinetic theory for “Soft” particles to bacterial cells: Implications for bacterial adhesion. Langmuir 21, 6462–6472. doi: 10.1021/la047049t
Delgado, A. V., González-Caballero, F., Hunter, R. J., Koopal, L. K., and Lyklema, J. (2005). Measurement and interpretation of electrokinetic phenomena. Pure Appl. Chem. 77, 1753–1805. doi: 10.1351/pac200577101753
de Wouters, T., Jans, C., Niederberger, T., Fischer, P., and Rühs, P. A. (2015). Adhesion potential of intestinal microbes predicted by physico-chemical characterization methods. PLoS ONE 10:e0136437. doi: 10.1371/journal.pone.0136437
Duardo-Sanchez, A., Munteanu, C. R., Riera-Fernandez, P., Lopez-Diaz, A., Pazos, A., and Gonzalez-Diaz, H. (2014). Modeling complex metabolic reactions, ecological systems, and financial and legal networks with MIANN models based on Markov-Wiener node descriptors. J. Chem. Inf. Model. 54, 16–29. doi: 10.1021/ci400280n
Fox, D. G., Barry, M. C., Pitt, R. E., Roseler, D. K., and Stone, W. C. (1995). Application of the Cornell Net Carbohydrate and Protein model for cattle consuming forages. J. Anim. Sci. 73, 267–277. doi: 10.2527/1995.731267x
Garzón, E., and Sánchez-Soto, P. J. (2015). An improved method for determining the external specific surface area and the plasticity index of clayey samples based on a simplified method for non-swelling fine-grained soils. Appl. Clay Sci. 115, 97–107. doi: 10.1016/j.clay.2015.07.015
Gonzalez-Diaz, H., Arrasate, S., Gomez-SanJuan, A., Sotomayor, N., Lete, E., Besada-Porto, L., et al. (2013). General theory for multiple input-output perturbations in complex molecular systems. 1. Linear QSPR electronegativity models in physical, organic, and medicinal chemistry. Curr. Top. Med. Chem. 13, 1713–1741. doi: 10.2174/1568026611313140011
González-Durruthy, M., Castro, M., Nunes, S. M., Ventura-Lima, J., Alberici, L. C., Naal, Z., et al. (2017). QSPR/QSAR-based Perturbation Theory approach and mechanistic electrochemical assays on carbon nanotubes with optimal properties against mitochondrial Fenton reaction experimentally induced by Fe2+-overload. Carbon N.Y. 115, 312–330. doi: 10.1016/j.carbon.2017.01.002
Hanaor, D., Michelazzi, M., Leonelli, C., and Sorrell, C. C. (2012). The effects of carboxylic acids on the aqueous dispersion and electrophoretic deposition of ZrO2. J. Eur. Ceram. Soc. 32, 235–244. doi: 10.1016/j.jeurceramsoc.2011.08.015
Hayashi, H., Seiki, H., Tsuneda, S., Hirata, A., and Sasaki, H. (2003). Influence of growth phase on bacterial cell electrokinetic characteristics examined by soft particle electrophoresis theory. J. Colloid Interface Sci. 264, 565–568. doi: 10.1016/S0021-9797(03)00418-1
Hua, Z. Z., Chen, J., Lun, S. Y., and Wang, X. R. (2003). Influence of biosurfactants produced by Candida antarctica on surface properties of microorganism and biodegradation of n-alkanes. Water Res. 37, 4143–4150. doi: 10.1016/S0043-1354(03)00380-4
Huang, X., Mu, T., Shen, C., Lu, L., and Liu, J. (2016). Effects of bio-surfactants combined with alkaline conditions on volatile fatty acid production and microbial community in the anaerobic fermentation of waste activated sludge. Int. Biodeterior. Biodegradation 114, 24–30. doi: 10.1016/j.ibiod.2016.05.014
Huang, X., Shen, C., Liu, J., and Lu, L. (2015). Improved volatile fatty acid production during waste activated sludge anaerobic fermentation by different bio-surfactants. Chem. Eng. J. 264, 280–290. doi: 10.1016/j.cej.2014.11.078
Hunter, R. J. (1981). “Chapter 3 - The calculation of zeta potential,” in Zeta Potential in Colloid Science: Principles and Applications, eds R. H. Ottewill and R. L. Rowell (London; San Diego; New York; Berkeley; Boston; Sydney; Tokyo; Toronto: Academic Press), 59–124.
Instruments, M. (2012). Zeta Potential an Introduction in 30 Minutes. Available online at: www.malvern.com
Kaszuba, M., Corbett, J., Watson, F. M., and Jones, A. (2010). High-concentration zeta potential measurements using light-scattering techniques. Philos. Tran. A Math. Phys. Eng. Sci. 368, 4439–4451. doi: 10.1098/rsta.2010.0175
Keselowsky, B. G., Collard, D. M., and Garcia, A. J. (2003). Surface chemistry modulates fibronectin conformation and directs integrin binding and specificity to control cell adhesion. J. Biomed. Mater. Res. A 66, 247–259. doi: 10.1002/jbm.a.10537
Keselowsky, B. G., Collard, D. M., and García, A. J. (2005). Integrin binding specificity regulates biomaterial surface chemistry effects on cell differentiation. Proc. Natl. Acad. Sci. U.S.A. 102, 5953–5957. doi: 10.1073/pnas.0407356102
Kim, W., Gamo, Y., Yahaya, M. S., Yimiti, W., and Goto, M. (2005). Pretreatment of surfactant Tween 80 and fibrolytic enzyme influencing volatile fatty acid and methane production of a total mixed ration by mixed rumen microorganisms at in vitro. Grassl. Sci. 51, 145–152. doi: 10.1111/j.1744-697X.2005.00020.x
Kłodzińska, E., Szumski, M., Dziubakiewicz, E., Hrynkiewicz, K., Skwarek, E., Janusz, W., et al. (2010). Effect of zeta potential value on bacterial behavior during electrophoretic separation. Electrophoresis 31, 1590–1596. doi: 10.1002/elps.200900559
Krause, D. O., Denman, S. E., Mackie, R. I., Morrison, M., Rae, A. L., Attwood, G. T., et al. (2003). Opportunities to improve fiber degradation in the rumen: microbiology, ecology, and genomics. FEMS Microbiol. Rev. 27, 663–693. doi: 10.1016/S0168-6445(03)00072-X
Li, F., Weir, M. D., Chen, J., and Xu, H. H. K. (2014). Effect of charge density of bonding agent containing a new quaternary ammonium methacrylate on antibacterial and bonding properties. Dental Mater. 30, 433–441. doi: 10.1016/j.dental.2014.01.002
Li, J.-S., and Jiang, Z.-W. (2001). Distribution of food particles with different size in the digestive tract in Mongolian gazelle. Acta Zool. Sin. 47, 488–494. doi: 10.3969/j.issn.1674-5507.2001.05.002
Liu, Y., Buendia-Rodriguez, G., Penuelas-Rivas, C. G., Tan, Z., Rivas-Guevara, M., Tenorio-Borroto, E., et al. (2015). Experimental and computational studies of fatty acid distribution networks. Mol. Biosyst. 11, 2964–2977. doi: 10.1039/C5MB00325C
Liu, Y., Ran, T., Tan, Z.-L., Tang, S.-X., and Wang, P.-P. (2013). Effects of surface tension and specific surface areas on in vitro fermentation of fiber. Acta Vet. Zootech. Sin. 44, 901–910. doi: 10.11843/j.issn.0366-6964.2013.06.011
Liu, Y., Ran, T., Tenorio-Borroto, E., Tang, S., Pazos, A., Tan, Z., et al. (2016). Experimental and chemometric studies of cell membrane permeability. Chemometr. Intell. Lab. Syst. 154, 1–6. doi: 10.1016/j.chemolab.2016.03.010
Liu, Y., Tang, S., Fernandez-Lozano, C., Munteanu, C. R., Pazos, A., Yu, Y.-Z., et al. (2017a). Experimental study and Random Forest prediction model of microbiome cell surface hydrophobicity. Expert Syst. Appl. 72, 306–316. doi: 10.1016/j.eswa.2016.10.058
Liu, Y., Tang, S., and Munteanu, C. R (2017b). Full Dataset and Figures for Experimental Study and ANN Dual-Time Scale Perturbation Model of Electrokinetic Properties of Microbiota. Available online at: https://doi.org/10.6084/m9.figshare.4558243.v5 [Accessed Feb. 01 2017].
Long, K., and Knapp, J. S. (1991). The effect of junlon PW110 and tween 80 on the production of cellulolytic enzymes by Coprinus cinereus. Mycol. Res. 95, 1077–1081. doi: 10.1016/S0953-7562(09)80550-9
Luan, F., Kleandrova, V. V., Gonzalez-Diaz, H., Ruso, J. M., Melo, A., Speck-Planche, A., et al. (2014). Computer-aided nanotoxicology: assessing cytotoxicity of nanoparticles under diverse experimental conditions by using a novel QSTR-perturbation approach. Nanoscale 6, 10623–10630. doi: 10.1039/C4NR01285B
Messina, P. V., Besada-Porto, J. M., Gonzalez-Diaz, H., and Ruso, J. M. (2015). Self-assembled binary nanoscale systems: multioutput model with LFER-covariance perturbation theory and an experimental-computational study of NaGDC-DDAB micelles. Langmuir 31, 12009–12018. doi: 10.1021/acs.langmuir.5b03074
Ran, T., Liu, Y., Li, H., Tang, S., He, Z., Munteanu, C. R., et al. (2016). Gastrointestinal spatiotemporal mRNA expression of ghrelin vs. growth hormone receptor and new growth yield machine learning model based on perturbation theory. Sci. Rep. 6:30174. doi: 10.1038/srep30174
Shida, T., Koseki, H., Yoda, I., Horiuchi, H., Sakoda, H., and Osaki, M. (2013). Adherence ability of Staphylococcus epidermidis on prosthetic biomaterials: an in vitro study. Int. J. Nanomedicine 8, 3955–3961. doi: 10.2147/ijn.s51994
Soni, K. A., Balasubramanian, A. K., Beskok, A., and Pillai, S. D. (2008). Zeta potential of selected bacteria in drinking water when dead, starved, or exposed to minimal and rich culture media. Curr. Microbiol. 56, 93–97. doi: 10.1007/s00284-007-9046-z
STATISTICA (2002). Data Analysis Software System. version 6.0 ed. Availabe online at: www.statsoft.com
Tang, S., Tan, Z., Zhou, C., Jiang, H., Jiang, Y., and Sheng, L. (2006). A comparison of in vitro fermentation characteristics of different botanical fractions of mature maize stover. J. Anim. Feed Sci. 15, 505–515. doi: 10.22358/jafs/66920/2006
Tang, S. X., Tayo, G. O., Tan, Z. L., Sun, Z. H., Shen, L. X., Zhou, C. S., et al. (2008). Effects of yeast culture and fibrolytic enzyme supplementation on in vitro fermentation characteristics of low-quality cereal straws. J. Anim. Sci. 86, 1164–1172. doi: 10.2527/jas.2007-0438
Tenorio-Borroto, E., Ramirez, F. R., Vazquez Chagoyan, J. C., de Oca Jimenez, R. M., Garcia-Mera, X., and Gonzalez-Diaz, H. (2016). Experimental-theoretic approach to drug-lymphocyte interactome networks with flow cytometry and spectral moments perturbation theory. Curr. Pharm. Des. 22, 5114–5119. doi: 10.2174/1381612822666160805164308
Turias, I., González, F., Martin, M. L., and Galindo, P. (2008). Prediction models of CO, SPM and SO2 concentrations in the Campo de Gibraltar Region, Spain: a multiple comparison strategy. Environ. Monit. Assess. 143, 131–146. doi: 10.1007/s10661-007-9963-0
Yoda, I., Koseki, H., Tomita, M., Shida, T., Horiuchi, H., Sakoda, H., et al. (2014). Effect of surface roughness of biomaterials on Staphylococcus epidermidis adhesion. BMC Microbiol. 14:234. doi: 10.1186/s12866-014-0234-2
Yuan, Z. Q., Tang, S. X., Zeng, B., Wang, M., Tan, Z. L., Sun, Z. H., et al. (2010). Effects of dietary supplementation with alkyl polyglycoside, a nonionic surfactant, on nutrient digestion and ruminal fermentation in goats1. J. Anim. Sci. 88, 3984–3991. doi: 10.2527/jas.2009-2397
Zeng, B., Tan, Z., Tang, S., Han, X., Tan, C., Zhong, R., et al. (2011). Effects of alkyl polyglycoside, a nonionic surfactant, and forage-to-concentrate ratio on rumen fermentation, amino acid composition of rumen content, bacteria and plasma in goats. Arch. Anim. Nutr. 65, 229–241. doi: 10.1080/1745039X.2011.559729
Zhang, J. H. (2008). Determination and Analysis on Content of Rumen Internal Environment, Physiological and Biochemical Index of Blood in Grazing and Confinedness Cashmere Goat. Master, Inner Mongolia Agricultural University.
Keywords: electrokinetic properties, zeta potential, artificial neural networks, perturbation theory, predictive model, ruminal microbiome
Citation: Liu Y, Munteanu CR, Fernandez-Lozano C, Pazos A, Ran T, Tan Z, Yu Y, Zhou C, Tang S and González-Díaz H (2017) Experimental Study and ANN Dual-Time Scale Perturbation Model of Electrokinetic Properties of Microbiota. Front. Microbiol. 8:1216. doi: 10.3389/fmicb.2017.01216
Received: 05 April 2017; Accepted: 14 June 2017;
Published: 30 June 2017.
Edited by:Claudio Fabricio Gonzalez, University of Florida, United States
Reviewed by:Catrin Ffion Williams, Cardiff University, United Kingdom
Armen Trchounian, Yerevan State University, Armenia
Copyright © 2017 Liu, Munteanu, Fernandez-Lozano, Pazos, Ran, Tan, Yu, Zhou, Tang and González-Díaz. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shaoxun Tang, firstname.lastname@example.org