1
|
Alsahafi S, Woodcock S. Mutual inhibition in presence of a virus in continuous culture. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2021; 18:3258-3273. [PMID: 34198384 DOI: 10.3934/mbe.2021162] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
In this paper, we consider two species competing for a limiting substrate such that each species impedes the growth of the other one (Mutual inhibition) in presence of a virus inhibiting one bacterial species. A system of ordinary differential equations is proposed as a mathematical model for this competition. A detailed local qualitative analysis of the system is carried out. We proved that for a general nonlinear growth rates, the Competitive Exclusion Principle still valid, that at least one species goes extinct. For some cases where we have two locally stable equilibrium points, initial species concentrations are important in determining which is the winning species. Obtained results were confirmed by some numerical simulations using Matlab software.
Collapse
Affiliation(s)
- Salah Alsahafi
- University of Technology Sydney, School of Mathematical and Physical Sciences, Australia
| | - Stephen Woodcock
- University of Technology Sydney, School of Mathematical and Physical Sciences, Australia
| |
Collapse
|
2
|
Mtar T, Fekih-Salem R, Sari T. Interspecific density-dependent model of predator–prey relationship in the chemostat. INT J BIOMATH 2020. [DOI: 10.1142/s1793524520500862] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The objective of this study is to analyze a model of competition for one resource in the chemostat with general interspecific density-dependent growth rates, taking into account the predator–prey relationship. This relationship is characterized by the fact that the prey species promotes the growth of the predator species which in turn inhibits the growth of the first species. The model is a three-dimensional system of ordinary differential equations. With the same dilution rates, the model can be reduced to a planar system where the two models have the same local and even global behavior. The existence and stability conditions of all steady states of the reduced model in the plane are determined according to the operating parameters. Using the nullcline method, we present a geometric characterization of the existence and stability of all equilibria showing the multiplicity of coexistence steady states. The bifurcation diagrams illustrate that the steady states can appear or disappear only through saddle-node or transcritical bifurcations. Moreover, the operating diagrams describe the asymptotic behavior of this system by varying the control parameters and show the effect of the inhibition of predation on the emergence of the bistability region and the reduction until the disappearance of the coexistence region by increasing this inhibition parameter.
Collapse
Affiliation(s)
- Tahani Mtar
- University of Tunis El Manar, National Engineering School of Tunis, LAMSIN 1002, Tunis, Tunisia
| | - Radhouane Fekih-Salem
- University of Tunis El Manar, National Engineering School of Tunis, LAMSIN 1002, Tunis, Tunisia
- University of Monastir, Higher Institute of Computer Science of Mahdia, 5111 Mahdia, Tunisia
| | - Tewfik Sari
- ITAP, Univ Montpellier, INRAE, Institut Agro, Montpellier, France
| |
Collapse
|
3
|
Oriano M, Zorzetto L, Guagliano G, Bertoglio F, van Uden S, Visai L, Petrini P. The Open Challenge of in vitro Modeling Complex and Multi-Microbial Communities in Three-Dimensional Niches. Front Bioeng Biotechnol 2020; 8:539319. [PMID: 33195112 PMCID: PMC7606986 DOI: 10.3389/fbioe.2020.539319] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/29/2020] [Accepted: 09/28/2020] [Indexed: 12/03/2022] Open
Abstract
The comprehension of the underlying mechanisms of the interactions within microbial communities represents a major challenge to be faced to control their outcome. Joint efforts of in vitro, in vivo and ecological models are crucial to controlling human health, including chronic infections. In a broader perspective, considering that polymicrobial communities are ubiquitous in nature, the understanding of these mechanisms is the groundwork to control and modulate bacterial response to any environmental condition. The reduction of the complex nature of communities of microorganisms to a single bacterial strain could not suffice to recapitulate the in vivo situation observed in mammals. Furthermore, some bacteria can adapt to various physiological or arduous environments embedding themselves in three-dimensional matrices, secluding from the external environment. Considering the increasing awareness that dynamic complex and dynamic population of microorganisms (microbiota), inhabiting different apparatuses, regulate different health states and protect against pathogen infections in a fragile and dynamic equilibrium, we underline the need to produce models to mimic the three-dimensional niches in which bacteria, and microorganisms in general, self-organize within a microbial consortium, strive and compete. This review mainly focuses, as a case study, to lung pathology-related dysbiosis and life-threatening diseases such as cystic fibrosis and bronchiectasis, where the co-presence of different bacteria and the altered 3D-environment, can be considered as worst-cases for chronic polymicrobial infections. We illustrate the state-of-art strategies used to study biofilms and bacterial niches in chronic infections, and multispecies ecological competition. Although far from the rendering of the 3D-environments and the polymicrobial nature of the infections, they represent the starting point to face their complexity. The increase of knowledge respect to the above aspects could positively affect the actual healthcare scenario. Indeed, infections are becoming a serious threat, due to the increasing bacterial resistance and the slow release of novel antibiotics on the market.
Collapse
Affiliation(s)
- Martina Oriano
- Molecular Medicine Department (DMM), Center for Health Technologies (CHT), UdR INSTM, University of Pavia, Pavia, Italy
- Department of Pathophysiology and Transplantation, University of Milan, Milan, Italy
- Internal Medicine Department, Respiratory Unit and Adult Cystic Fibrosis Center, Fondazione IRCCS Ca’ Granda Ospedale Maggiore Policlinico, Milan, Italy
| | - Laura Zorzetto
- Department of Biomaterials, Max Planck Institute of Colloids and Interfaces, Potsdam, Germany
| | - Giuseppe Guagliano
- Department of Chemistry, Materials and Chemical Engineering “Giulio Natta” and UdR INSTM Politecnico di Milano, Milan, Italy
| | - Federico Bertoglio
- Molecular Medicine Department (DMM), Center for Health Technologies (CHT), UdR INSTM, University of Pavia, Pavia, Italy
- Technische Universität Braunschweig, Institute of Biochemistry, Biotechnology and Bioinformatic, Department of Biotechnology, Braunschweig, Germany
| | - Sebastião van Uden
- Department of Chemistry, Materials and Chemical Engineering “Giulio Natta” and UdR INSTM Politecnico di Milano, Milan, Italy
| | - Livia Visai
- Molecular Medicine Department (DMM), Center for Health Technologies (CHT), UdR INSTM, University of Pavia, Pavia, Italy
- Department of Occupational Medicine, Toxicology and Environmental Risks, Istituti Clinici Scientifici (ICS) Maugeri, IRCCS, Pavia, Italy
| | - Paola Petrini
- Department of Chemistry, Materials and Chemical Engineering “Giulio Natta” and UdR INSTM Politecnico di Milano, Milan, Italy
| |
Collapse
|
4
|
El Hajji M. How can inter-specific interferences explain coexistence or confirm the competitive exclusion principle in a chemostat? INT J BIOMATH 2019. [DOI: 10.1142/s1793524518501115] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
In this paper, I consider two species feeding on limiting substrate in a chemostat taking into account some possible effects of each species on the other one. System of differential equations is proposed as model of these effects with general inter-specific density-dependent growth rates. Three cases were considered. The first one for a mutual inhibitory relationship where it is proved that at most one species can survive which confirms the competitive exclusion principle. Initial concentrations of species have great importance in determination of which species is the winner. The second one for a food web relationship where it is proved that under general assumptions on the dilution rate, both species persist for any initial conditions. Finally, a third case dealing with an obligate mutualistic relationship was discussed. It is proved that initial condition has a great importance in determination of persistence or extinction of both species.
Collapse
Affiliation(s)
- Miled El Hajji
- ENIT-LAMSIN, BP. 37, 1002 Tunis-Belvédère, Tunis El Manar University, Tunis, Tunisia
- General Studies Department, College of Telecom and Electronics, Technical and Vocational Training Corporation, Jeddah 2146, Saudi Arabia
| |
Collapse
|
5
|
Vet S, de Buyl S, Faust K, Danckaert J, Gonze D, Gelens L. Bistability in a system of two species interacting through mutualism as well as competition: Chemostat vs. Lotka-Volterra equations. PLoS One 2018; 13:e0197462. [PMID: 29874266 PMCID: PMC5991418 DOI: 10.1371/journal.pone.0197462] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2018] [Accepted: 05/02/2018] [Indexed: 12/17/2022] Open
Abstract
We theoretically study the dynamics of two interacting microbial species in the chemostat. These species are competitors for a common resource, as well as mutualists due to cross-feeding. In line with previous studies (Assaneo, et al., 2013; Holland, et al., 2010; Iwata, et al., 2011), we demonstrate that this system has a rich repertoire of dynamical behavior, including bistability. Standard Lotka-Volterra equations are not capable to describe this particular system, as these account for only one type of interaction (mutualistic or competitive). We show here that the different steady state solutions can be well captured by an extended Lotka-Volterra model, which better describe the density-dependent interaction (mutualism at low density and competition at high density). This two-variable model provides a more intuitive description of the dynamical behavior than the chemostat equations.
Collapse
Affiliation(s)
- Stefan Vet
- Interuniversity Institute of Bioinformatics in Brussels (IB2), VUB-ULB, Brussels, Belgium
- Applied Physics Research Group, Vrije Universiteit Brussel (VUB), Brussels, Belgium
- Unité de Chronobiologie théorique, Université Libre de Bruxelles (ULB), Brussels, Belgium
| | - Sophie de Buyl
- Interuniversity Institute of Bioinformatics in Brussels (IB2), VUB-ULB, Brussels, Belgium
- Applied Physics Research Group, Vrije Universiteit Brussel (VUB), Brussels, Belgium
| | - Karoline Faust
- Interuniversity Institute of Bioinformatics in Brussels (IB2), VUB-ULB, Brussels, Belgium
- Laboratory of Molecular Bacteriology, KU Leuven, Leuven, Belgium
| | - Jan Danckaert
- Interuniversity Institute of Bioinformatics in Brussels (IB2), VUB-ULB, Brussels, Belgium
- Applied Physics Research Group, Vrije Universiteit Brussel (VUB), Brussels, Belgium
| | - Didier Gonze
- Interuniversity Institute of Bioinformatics in Brussels (IB2), VUB-ULB, Brussels, Belgium
- Unité de Chronobiologie théorique, Université Libre de Bruxelles (ULB), Brussels, Belgium
| | - Lendert Gelens
- Applied Physics Research Group, Vrije Universiteit Brussel (VUB), Brussels, Belgium
- Laboratory of Dynamics in Biological Systems, KU Leuven, Leuven, Belgium
| |
Collapse
|
6
|
Sari T, Wade MJ. Generalised approach to modelling a three-tiered microbial food-web. Math Biosci 2017; 291:21-37. [PMID: 28709972 PMCID: PMC5552901 DOI: 10.1016/j.mbs.2017.07.005] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2016] [Revised: 05/26/2017] [Accepted: 07/10/2017] [Indexed: 11/18/2022]
Abstract
The complexity of the anaerobic digestion process has motivated the development of complex models, such as the widely used Anaerobic Digestion Model No. 1. However, this complexity makes it intractable to identify the stability profile coupled to the asymptotic behaviour of existing steady-states as a function of conventional chemostat operating parameters (substrate inflow concentration and dilution rate). In a previous study this model was simplified and reduced to its very backbone to describe a three-tiered chlorophenol mineralising food-web, with its stability analysed numerically using consensus values for the various biological parameters of the Monod growth functions. Steady-states where all organisms exist were always stable and non-oscillatory. Here we investigate a generalised form of this three-tiered food-web, whose kinetics do not rely on the specific kinetics of Monod form. The results are valid for a large class of growth kinetics as long as they keep the signs of their derivatives. We examine the existence and stability of the identified steady-states and find that, without a maintenance term, the stability of the system may be characterised analytically. These findings permit a better understanding of the operating region of the bifurcation diagram where all organisms exist, and its dependence on the biological parameters of the model. For the previously studied Monod kinetics, we identify four interesting cases that show this dependence of the operating diagram with respect to the biological parameters. When maintenance is included, it is necessary to perform numerical analysis. In both cases we verify the discovery of two important phenomena; i) the washout steady-state is always stable, and ii) a switch in dominance between two organisms competing for hydrogen results in the system becoming unstable and a loss in viability. We show that our approach results in the discovery of an unstable operating region in its positive steady-state, where all three organisms exist, a fact that has not been reported in a previous numerical study. This type of analysis can be used to determine critical behaviour in microbial communities in response to changing operating conditions.
Collapse
Affiliation(s)
- T Sari
- Irstea, UMR Itap, Montpellier, France & Université de Haute Alsace, Laboratoire de Mathématiques, Mulhouse, France
| | - M J Wade
- School of Civil Engineering and Geosciences, Newcastle University, Newcastle-upon-Tyne NE1 7RU, United Kingdom.
| |
Collapse
|
7
|
Fekih-Salem R, Lobry C, Sari T. A density-dependent model of competition for one resource in the chemostat. Math Biosci 2017; 286:104-122. [PMID: 28212840 DOI: 10.1016/j.mbs.2017.02.007] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2016] [Revised: 02/08/2017] [Accepted: 02/13/2017] [Indexed: 11/28/2022]
Abstract
This paper deals with a two-microbial species model in competition for a single-resource in the chemostat including general intra- and interspecific density-dependent growth rates with distinct removal rates for each species. In order to understand the effects of intra- and interspecific interference, this general model is first studied by determining the conditions of existence and local stability of steady states. With the same removal rate, the model can be reduced to a planar system and then the global stability results for each steady state are derived. The bifurcations of steady states according to interspecific interference parameters are analyzed in a particular case of density-dependent growth rates which are usually used in the literature. The operating diagrams show how the model behaves by varying the operating parameters and illustrate the effect of the intra- and interspecific interference on the disappearance of coexistence region and the occurrence of bi-stability region. Concerning the small enough interspecific interference terms, we would shed light on the global convergence towards the coexistence steady state for any positive initial condition. When the interspecific interference pressure is large enough this system exhibits bi-stability where the issue of the competition depends on the initial condition.
Collapse
Affiliation(s)
- Radhouane Fekih-Salem
- Université de Tunis El Manar, ENIT, LAMSIN, BP 37, Le Belvédère, 1002 Tunis, Tunisie; Université de Monastir, ISIMa, BP 49, Av Habib Bourguiba, 5111 Mahdia, Tunisie.
| | - Claude Lobry
- Université de Nice et MODEMIC, Le Gd Palais Bt6, 2 BD de Cimiez, 06000 Nice, France.
| | - Tewfik Sari
- IRSTEA, UMR Itap, 361 rue Jean-François Breton, 34196 Montpellier, France; Université de Haute Alsace, LMIA, 4 rue des frères Lumière, 68093 Mulhouse, France.
| |
Collapse
|
8
|
A Computational Study of Amensalistic Control of Listeria monocytogenes by Lactococcus lactis under Nutrient Rich Conditions in a Chemostat Setting. Foods 2016; 5:foods5030061. [PMID: 28231156 PMCID: PMC5302389 DOI: 10.3390/foods5030061] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2016] [Revised: 08/24/2016] [Accepted: 08/29/2016] [Indexed: 11/22/2022] Open
Abstract
We study a previously introduced mathematical model of amensalistic control of the foodborne pathogen Listeria monocytogenes by the generally regarded as safe lactic acid bacteria Lactococcus lactis in a chemostat setting under nutrient rich growth conditions. The control agent produces lactic acids and thus affects pH in the environment such that it becomes detrimental to the pathogen while it is much more tolerant to these self-inflicted environmental changes itself. The mathematical model consists of five nonlinear ordinary differential equations for both bacterial species, the concentration of lactic acids, the pH and malate. The model is algebraically too involved to allow a comprehensive, rigorous qualitative analysis. Therefore, we conduct a computational study. Our results imply that depending on the growth characteristics of the medium in which the bacteria are cultured, the pathogen can survive in an intermediate flow regime but will be eradicated for slower flow rates and washed out for higher flow rates.
Collapse
|
9
|
Sari T, Harmand J. A model of a syntrophic relationship between two microbial species in a chemostat including maintenance. Math Biosci 2016; 275:1-9. [PMID: 26926240 DOI: 10.1016/j.mbs.2016.02.008] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2015] [Revised: 02/11/2016] [Accepted: 02/17/2016] [Indexed: 10/22/2022]
Abstract
Many microbial ecosystems can be seen as microbial 'food chains' where the different reaction steps can be seen as such: the waste products of the organisms at a given reaction step are consumed by organisms at the next reaction step. In the present paper we study a model of a two-step biological reaction with feedback inhibition, which was recently presented as a reduced and simplified version of the anaerobic digestion model ADM1 of the International Water Association (IWA). It is known that in the absence of maintenance (or decay) the microbial 'food chain' is stable. In a previous study, using a purely numerical approach and ADM1 consensus parameter values, it was shown that the model remains stable when decay terms are added. However, the authors could not prove in full generality that it remains true for other parameter values. In this paper we prove that introducing decay in the model preserves stability whatever its parameters values are and for a wide range of kinetics.
Collapse
Affiliation(s)
- Tewfik Sari
- Irstea, UMR ITAP, 361 rue Jean-François Breton, 34196 Montpellier, France; Université de Haute Alsace, LMIA, 4 rue des Frères Lumière, 68093 Mulhourse, France.
| | - Jérôme Harmand
- INRA, UR050, Laboratoire de Biotechnologie de l'Environnement, 102 Avenue des Etangs, Narbonne F-11100, France.
| |
Collapse
|
10
|
Wade M, Harmand J, Benyahia B, Bouchez T, Chaillou S, Cloez B, Godon JJ, Moussa Boudjemaa B, Rapaport A, Sari T, Arditi R, Lobry C. Perspectives in mathematical modelling for microbial ecology. Ecol Modell 2016. [DOI: 10.1016/j.ecolmodel.2015.11.002] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
11
|
Sari T, El Hajji M, Harmand J. The mathematical analysis of a syntrophic relationship between two microbial species in a chemostat. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2012; 9:627-645. [PMID: 22881029 DOI: 10.3934/mbe.2012.9.627] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
A mathematical model involving a syntrophic relationship between two populations of bacteria in a continuous culture is proposed. A detailed qualitative analysis is carried out as well as the analysis of the local and global stability of the equilibria. We demonstrate, under general assumptions of monotonicity which are relevant from an applied point of view, the asymptotic stability of the positive equilibrium point which corresponds to the coexistence of the two bacteria. A syntrophic relationship in the anaerobic digestion process is proposed as a real candidate for this model.
Collapse
Affiliation(s)
- Tewfik Sari
- Irstea, UMR ITAP, 361 rue Jean-Francois Breton 34196 Montpellier, France.
| | | | | |
Collapse
|