Effect of cultivar and substrate on the efficacy of biopesticides to suppress Pythium on greenhouse crops

Oomycetes, such as Pythium, are soil-borne plant pathogens that can cause significant losses in greenhouse crop production due to their swimming zoospores and wide host range. Additionally, the increasing use of substrates that lack microbial diversity in greenhouse production creates a “biological vacuum” that can reduce the substrate’s capacity to resist microbial invasion by soil-borne diseases. The lack of competition by a natural microbial community and the environmental conditions of greenhouse production creates an ideal situation for the use of biopesticides. Biopesticides are commercial products that use beneficial microorganisms (biocontrol agents) to suppress disease and promote plant health. Greenhouse producers can utilize commercial biopesticides in addition to chemical treatments to protect plants from soil-borne pathogens. One barrier to use of biopesticides is the variability of their performance which can be attributed to differences in environmental conditions, such as plant species and substrate materials. Few studies have evaluated the effect of plant cultivar and current substrates on the efficacy of biopesticides to suppress disease in horticulture crops. The objectives of this research were to 1) develop a greenhouse-based assay to study biopesticide suppression of Pythium root rot of greenhouse grown crops, 2) evaluate the effect of plant cultivar on biopesticide efficacy in a tomato system, and 3) evaluate the implications of propagation substrate on biopesticide efficacy in cucumber and calibrachoa systems. For each of these systems, commercially available biopesticides were applied at the label rate twice during propagation. At transplant, plants were challenged with Pythium spp. or a water control. Root rot and root growth were evaluated at 21 days post infection. Findings reveal that the plant cultivars tested did not affect biopesticide efficacy, however a different cultivar panel with greater genetic diversity may affect biopesticide efficacy. There was a significant effect of propagation substrate on disease severity. Plants propagated in coconut coir had higher root disease than those propagated in Oasis®. These findings suggest that the chemical and physical properties of these substrates affect plant susceptibility to disease or pathogen activity, however further research is needed to evaluate this observation. Furthermore, there was a significant interaction between the effects of substrate and biopesticide on root rot severity in which biopesticide efficacy varied by substrate. This result suggests that substrate may affect biopesticide performance, but further research is needed to confirm these results and to understand the mechanisms behind this phenomenon. Finally, in all experiments, the commercial biopesticide Rootshield® WP suppressed root disease compared to the infested water control. These experiments provided initial data for determining the mechanisms driving variation in biopesticide performance and to improve on-farm performance and adoption

Model simplification of signal transduction pathway networks via a hybrid inference strategy

A full-scale mathematical model of cellular networks normally involves a large number of variables and parameters. How to effectively develop manageable and reliable models is crucial for effective computation, analysis and design of such systems. The aim of model simplification is to eliminate parts of a model that are unimportant for the properties of interest. In this work, a model reduction strategy via hybrid inference is proposed for signal pathway networks. It integrates multiple techniques including conservation analysis, local sensitivity analysis, principal component analysis and flux analysis to identify the reactions and variables that can be considered to be eliminated from the full-scale model. Using an I·B-NF-·B signalling pathway model as an example, simulation analysis demonstrates that the simplified model quantitatively predicts the dynamic behaviours of the network

Insights into the behaviour of systems biology models from dynamic sensitivity and identifiability analysis: a case study of an NF-kB signaling pathway

Mathematical modelling offers a variety of useful techniques to help in understanding the intrinsic behaviour of complex signal transduction networks. From the system engineering point of view, the dynamics of metabolic and signal transduction models can always be described by nonlinear ordinary differential equations (ODEs) following mass balance principles. Based on the state-space formulation, many methods from the area of automatic control can conveniently be applied to the modelling, analysis and design of cell networks. In the present study, dynamic sensitivity analysis is performed on a model of the IB-NF-B signal pathway system. Univariate analysis of the Euclidean-form overall sensitivities shows that only 8 out of the 64 parameters in the model have major influence on the nuclear NF-B oscillations. The sensitivity matrix is then used to address correlation analysis, identifiability assessment and measurement set selection within the framework of least squares estimation and multivariate analysis. It is shown that certain pairs of parameters are exactly or highly correlated to each other in terms of their effects on the measured variables. The experimental design strategy provides guidance on which proteins should best be considered for measurement such that the unknown parameters can be estimated with the best statistical precision. The whole analysis scheme we describe provides efficient parameter estimation techniques for complex cell networks

Unscented Kalman filter with parameter identifiability analysis for the estimation of multiple parameters in kinetic models

In systems biology, experimentally measured parameters are not always available, necessitating the use of computationally based parameter estimation. In order to rely on estimated parameters, it is critical to first determine which parameters can be estimated for a given model and measurement set. This is done with parameter identifiability analysis. A kinetic model of the sucrose accumulation in the sugar cane culm tissue developed by Rohwer et al. was taken as a test case model. What differentiates this approach is the integration of an orthogonal-based local identifiability method into the unscented Kalman filter (UKF), rather than using the more common observability-based method which has inherent limitations. It also introduces a variable step size based on the system uncertainty of the UKF during the sensitivity calculation. This method identified 10 out of 12 parameters as identifiable. These ten parameters were estimated using the UKF, which was run 97 times. Throughout the repetitions the UKF proved to be more consistent than the estimation algorithms used for comparison

The Attributed Pi Calculus with Priorities

International audienceWe present the attributed $\pi$-calculus for modeling concurrent systems with interaction constraints depending on the values of attributes of processes. The $\pi$-calculus serves as a constraint language underlying the $\pi$-calculus. Interaction constraints subsume priorities, by which to express global aspects of populations. We present a nondeterministic and a stochastic semantics for the attributed $\pi$-calculus. We show how to encode the $\pi$-calculus with priorities and polyadic synchronization $\pi$@ and thus dynamic compartments, as well as the stochastic $\pi$-calculus with concurrent objects spico. We illustrate the usefulness of the attributed $\pi$-calculus for modeling biological systems at two particular examples: Euglena’s spatial movement in phototaxis, and cooperative protein binding in gene regulation of bacteriophage lambda. Furthermore, population-based model is supported beside individual-based modeling. A stochastic simulation algorithm for the attributed $\pi$-calculus is derived from its stochastic semantics. We have implemented a simulator and present experimental results, that confirm the practical relevance of our approach

A combined model reduction algorithm for controlled biochemical systems

Background: Systems Biology continues to produce increasingly large models of complex biochemical reaction networks. In applications requiring, for example, parameter estimation, the use of agent-based modelling approaches, or real-time simulation, this growing model complexity can present a significant hurdle. Often, however, not all portions of a model are of equal interest in a given setting. In such situations methods of model reduction offer one possible approach for addressing the issue of complexity by seeking to eliminate those portions of a pathway that can be shown to have the least effect upon the properties of interest. Methods: In this paper a model reduction algorithm bringing together the complementary aspects of proper lumping and empirical balanced truncation is presented. Additional contributions include the development of a criterion for the selection of state-variable elimination via conservation analysis and use of an ‘averaged’ lumping inverse. This combined algorithm is highly automatable and of particular applicability in the context of ‘controlled’ biochemical networks. Results: The algorithm is demonstrated here via application to two examples; an 11 dimensional model of bacterial chemotaxis in Escherichia coli and a 99 dimensional model of extracellular regulatory kinase activation (ERK) mediated via the epidermal growth factor (EGF) and nerve growth factor (NGF) receptor pathways. In the case of the chemotaxis model the algorithm was able to reduce the model to 2 state-variables producing a maximal relative error between the dynamics of the original and reduced models of only 2.8% whilst yielding a 26 fold speed up in simulation time. For the ERK activation model the algorithm was able to reduce the system to 7 state-variables, incurring a maximal relative error of 4.8%, and producing an approximately 10 fold speed up in the rate of simulation. Indices of controllability and observability are additionally developed and demonstrated throughout the paper. These provide insight into the relative importance of individual reactants in mediating a biochemical system’s input-output response even for highly complex networks. Conclusions: Through application, this paper demonstrates that combined model reduction methods can produce a significant simplification of complex Systems Biology models whilst retaining a high degree of predictive accuracy. In particular, it is shown that by combining the methods of proper lumping and empirical balanced truncation it is often possible to produce more accurate reductions than can be obtained by the use of either method in isolation

Inferring robust gene networks from expression data by a sensitivity-based incremental evolution method

<p>Abstract</p> <p>Background</p> <p>Reconstructing gene regulatory networks (GRNs) from expression data is one of the most important challenges in systems biology research. Many computational models and methods have been proposed to automate the process of network reconstruction. Inferring robust networks with desired behaviours remains challenging, however. This problem is related to network dynamics but has yet to be investigated using network modeling.</p> <p>Results</p> <p>We propose an incremental evolution approach for inferring GRNs that takes network robustness into consideration and can deal with a large number of network parameters. Our approach includes a sensitivity analysis procedure to iteratively select the most influential network parameters, and it uses a swarm intelligence procedure to perform parameter optimization. We have conducted a series of experiments to evaluate the external behaviors and internal robustness of the networks inferred by the proposed approach. The results and analyses have verified the effectiveness of our approach.</p> <p>Conclusions</p> <p>Sensitivity analysis is crucial to identifying the most sensitive parameters that govern the network dynamics. It can further be used to derive constraints for network parameters in the network reconstruction process. The experimental results show that the proposed approach can successfully infer robust GRNs with desired system behaviors.</p

A comparison of approximation techniques for variance-based sensitivity analysis of biochemical reaction systems

<p>Abstract</p> <p>Background</p> <p>Sensitivity analysis is an indispensable tool for the analysis of complex systems. In a recent paper, we have introduced a thermodynamically consistent variance-based sensitivity analysis approach for studying the robustness and fragility properties of biochemical reaction systems under uncertainty in the standard chemical potentials of the activated complexes of the reactions and the standard chemical potentials of the molecular species. In that approach, key sensitivity indices were estimated by Monte Carlo sampling, which is computationally very demanding and impractical for large biochemical reaction systems. Computationally efficient algorithms are needed to make variance-based sensitivity analysis applicable to realistic cellular networks, modeled by biochemical reaction systems that consist of a large number of reactions and molecular species.</p> <p>Results</p> <p>We present four techniques, derivative approximation (DA), polynomial approximation (PA), Gauss-Hermite integration (GHI), and orthonormal Hermite approximation (OHA), for <it>analytically </it>approximating the variance-based sensitivity indices associated with a biochemical reaction system. By using a well-known model of the mitogen-activated protein kinase signaling cascade as a case study, we numerically compare the approximation quality of these techniques against traditional Monte Carlo sampling. Our results indicate that, although DA is computationally the most attractive technique, special care should be exercised when using it for sensitivity analysis, since it may only be accurate at low levels of uncertainty. On the other hand, PA, GHI, and OHA are computationally more demanding than DA but can work well at high levels of uncertainty. GHI results in a slightly better accuracy than PA, but it is more difficult to implement. OHA produces the most accurate approximation results and can be implemented in a straightforward manner. It turns out that the computational cost of the four approximation techniques considered in this paper is orders of magnitude smaller than traditional Monte Carlo estimation. Software, coded in MATLAB^®, which implements all sensitivity analysis techniques discussed in this paper, is available free of charge.</p> <p>Conclusions</p> <p>Estimating variance-based sensitivity indices of a large biochemical reaction system is a computationally challenging task that can only be addressed via approximations. Among the methods presented in this paper, a technique based on orthonormal Hermite polynomials seems to be an acceptable candidate for the job, producing very good approximation results for a wide range of uncertainty levels in a fraction of the time required by traditional Monte Carlo sampling.</p

The integrated analysis of metabolic and protein interaction networks reveals novel molecular organizing principles

Background: The study of biological interaction networks is a central theme of systems biology. Here, we investigate the relationships between two distinct types of interaction networks: the metabolic pathway map and the protein-protein interaction network (PIN). It has long been established that successive enzymatic steps are often catalyzed by physically interacting proteins forming permanent or transient multi-enzymes complexes. Inspecting high-throughput PIN data, it was shown recently that, indeed, enzymes involved in successive reactions are generally more likely to interact than other protein pairs. In our study, we expanded this line of research to include comparisons of the underlying respective network topologies as well as to investigate whether the spatial organization of enzyme interactions correlates with metabolic efficiency. Results: Analyzing yeast data, we detected long-range correlations between shortest paths between proteins in both network types suggesting a mutual correspondence of both network architectures. We discovered that the organizing principles of physical interactions between metabolic enzymes differ from the general PIN of all proteins. While physical interactions between proteins are generally dissortative, enzyme interactions were observed to be assortative. Thus, enzymes frequently interact with other enzymes of similar rather than different degree. Enzymes carrying high flux loads are more likely to physically interact than enzymes with lower metabolic throughput. In particular, enzymes associated with catabolic pathways as well as enzymes involved in the biosynthesis of complex molecules were found to exhibit high degrees of physical clustering. Single proteins were identified that connect major components of the cellular metabolism and may thus be essential for the structural integrity of several biosynthetic systems. Conclusion: Our results reveal topological equivalences between the protein interaction network and the metabolic pathway network. Evolved protein interactions may contribute significantly towards increasing the efficiency of metabolic processes by permitting higher metabolic fluxes. Thus, our results shed further light on the unifying principles shaping the evolution of both the functional (metabolic) as well as the physical interaction network