(16 µm). Then transferred to a 250 ml flask with culture medium without nitrate, and left in
|
|
- James Rose
- 5 years ago
- Views:
Transcription
1 Oikos OIK Barreiro, A., Roy, S. and Vasconcelos, V. M Allelopathy prevents competitive exclusion and promotes phytoplankton biodiversity. Oikos doi: /oik Appendix 1 Parameterization of uptake, growth and allelopathic effect Uptake experiments were performed in batch cultures. The phytoplankton species were pre-cultured during 6 days in 50 ml flasks (one flask per species) with low nitrate concentration (16 µm). Then transferred to a 250 ml flask with culture medium without nitrate, and left in these conditions during 24 h. Then, phosphate was resupplied (20 µm) and a pulse of 10 µm of nitrate was added to the flasks. Then, nitrate concentration (three 6 ml replicates) and cell abundances were monitored every 2 h during 18 h. All these experiments were performed in a culture room where light conditions were 40 µmols photons m -2 s -1, during 24 h per day, and temperature 20 ºC. The data obtained were fitted to a Michaelis-Menten equation (see further details below). Experiments for the estimation of the growth parameters were performed in continuous cultures. Flask volume was 400 ml and dilution rate was 0.3 day -1. Initial nitrate concentration was 320 µm and phosphate 200 µm. Experiments were performed in culture room where light conditions were 40 µmols photons m -2 s -1, during 24 h per day and temperature 20 ºC. Nitrate concentrations (in three 6 ml replicates) and cell abundances were monitored every 24 h. Experiments lasted days. Data were fitted to the Monod model (see further details below). Allelopathic effect of each species was estimated by testing the effect of increasing
2 concentrations (0, 10%, 25%, 50% and 75% of final volume) of cell-free filtrate of the donor species on the growth rate of a target species during 24 h (but 72 h when Oscillatoria sp. was used as target species). Cell-free filtrate was obtained by filtering through 0.22 µm filters (Millipore Express membrane filters) medium from continuous cultures of the donor species in steady-state, grown in the same conditions as above. These filtrates were added to 3 replicated 5 ml-vials containing cell suspensions of the target species, with culture medium and saturating concentrations of nitrate (3200 µm) and phosphate (200 µm). Experiments were performed in culture room where light conditions were 40 µmols photons m -2 s -1, during 24 h per day and temperature 20 ºC. The initial abundances of each target species were set at approximately equal biovolumes. Cell abundances were estimated with Neubauer or Sedgewyck-Rafter counting chambers (depending on cell abundance) on each of the 3 replicates immediately at the initial time and after 24 or 72 h (see above). Nitrate analysis Nitrate analysis were performed in samples of medium filtered through 0.22 µm filters (Millipore Express membrane filters). Analysis was performed in a Skalar Sanplus nutrient auto-analyzer, using the method Skalar M (EPA 353.2). Model formulations In order to fit the results from the experiments shown in Figure 1, we considered a resource-competition model of two phytoplankton competing for nitrate, and included the allelopathic effects from Oscillatoria sp. by a non-liner mortality term as function of the densities of both species. The choice of the non-linear mortality is crucial to correctly represent
3 the allelopathic effects. We tested several possible options based on the literature, and also introduced new functional forms to investigate the most suitable one that can describe our experimental results. The most basic formulations for a model of this kind could be as follows (Model 1): dn dt = D N & N F ) N 1 η ) F, N 1 η, dp ) dt = F ) N γp, P ) DP ) dp, dt = F, N DP, Where N is the nitrate concentration in the culture, N 0 is the inflowing nitrate concentration, D is the dilution rate of the system, P i (being i = 1,2) are the population abundances of Ankistrodesmus falcatus and Oscillatoria sp. respectively. η i are the yield coefficients of species i. as a ratio of mass of cells to mass of nutrient, γ is the parameter denoting the allelopathic effect. The functions F i (N) are the growth functions: F / N = µμ /P / N K / + N Where µ i are the maximum growth rates of species i, and K i the half saturation constants for nitrate. A model with this formulation does not support coexistence (known from theoretical analysis), so this formulation was discarded. A second form of the model (Model 2) includes a positive relationship between the growth rate of the allelopathic species (P 2 ) and the abundance of the non-allelopathic species (P 1 ). This positive relationship could be due to nutritional exudates or alkaline phosphatases by P 1, but also to the release of nutrients from cells of P 1 that died due to allelopathy. So, the only difference in Model 2 relative to Model 1 is the growth equation of the allelopathic species.
4 The equation for the non allelopathic growth is given by: F ) N = µμ )P ) N K ) + N with the same parameters as Model 1. The growth equation for the allelopathic species is given by: F, N = µμ,p, N K, + N where: µμ, P ) = µμ,& 1 + εp ) K 4 + P ) µ 20 is the basal maximum growth rate of this species, which is affected by the abundance of the non-allelopathic species, ε is a rate associated to this process and K p the half saturation constant. This form of the model was able to fit very accurately the two extreme cases of exclusion as well, and was also able to fit coexistence. For certain parameter values and conditions, this model can fit oscillatory coexistence (S. Roy, in prep.), similar to what is shown in the Figure 1 D, E, and F (main text). However, the values of the parameters and the state variables needed to show this oscillatory coexistence differed substantially from those needed to fit the data from this manuscript. In order to be able to fit the oscillatory coexistence using a reasonable range of parameter values and initial conditions for our experiments, we needed to incorporate an additional mechanism to those previous basic model formulations. We considered the inclusion of a more biologically feasible mechanism, which is a density dependent function determining allelochemical production. Density dependent mechanisms are inherent to most of the biological rate processes. It is very reasonable as well to consider that this rate would have lower and upper activation limits. So, we incorporated a mechanism with these
5 features through two alternative functional forms. The general model formulation in both cases is as follows: dn dt = D N & N F ) N 1 η ) F, N 1 η, dp ) dt = F ) N φ ) (P, )φ, (P, )γp ) DP ) dp, dt = F, N DP, All the terms are equivalent to Model 1 except Φ 1 (P 2 ) and Φ 2 (P 2 ), which are new. In the first alternative functional form (Model 3) Φ 1 (P 2 ) is the function that determines the negative density dependent allelopathic effect with lower threshold: φ ) P, = e 9:; < + b Where a is the rate of the process and b the lower threshold. In the same Model 3, Φ 2 (P 2 ) determines the upper limit of allelopathy effect (the minimum abundance of the allelopathic species in order for the process to occur): φ, P, = max P, n & P,, 0 Where n 0 is this minimum abundance of the allelopathic species needed for the allelopatic process. This model was able to fit qualitatively and quantitatively all the outcomes we observed in our experiments. To establish that the dynamics predicted by the model is not very specific to one non-linear functional form, we further considered alternative formulations, in order to select among them, the one providing the best fits. The first of those alternatives consists in the same model structure with a different functional form for the density dependent allelopathic killing (a similar
6 function considered in Roy et al in a different context) where: Φ 1 (P 2 ) is: φ ) P, = 1 e 9D E; < e 9D <; < Where C 1 determines the upper threshold of allelochemical production, whereas C 2 affects the slope of the process. Φ 2 (P 2 ) is the same as in Model 3. We termed this new alternative model formulation Model 4. Finally, we also considered two additional model formulations, which consisted in the same basic formulations as Models 3 and 4, but including the more explicit function for the growth of the allelopathic species, as in Model 2. The first of these two models was Model 5, where, similarly as in Model 2: F ) N = µμ )P ) N K ) + N and: F, N = µμ,p, N K, + N where: µμ, P ) = µμ,& 1 + εp ) K 4 + P ) And, as in Model 3: φ ) P, φ, P, = e 9:; < + b = max P, n & P,, 0 And the second alternative formulation is Model 6, which only changes with respect to Model 5 in: φ ) P, = 1 e 9D E; < e 9D <; < Which are the same functional forms used in Model 4 for these processes.
7 There are several mechanisms that could generate oscillatory coexistence in allelopathy models. One of them consists in incorporating a fully mechanistic structure for allelopathy (including rates of allelochemical production, uptake, and degradation) and adopting a sigmoidal function for the uptake of allelochemicals (Hsu and Waltman 2004). However, this was not successful in our case, probably because the referenced model was developed for a completely different context (bacterial allelopathy via plasmids). Other features in the model formulations such as delays (Mukhopadhyay et al. 1998) and adaptation (Mougi 2013) also showed to generate oscillatory coexistence, but in Lotka-Volterra models. We discarded incorporating a nonmechanistic delay in our models such as Mukhopadhyay et al. 1998, because we were more interested in explaining the observed dynamics as much as possible with biological mechanisms. Regarding the adaptation, we considered that this is a complex process that would need more evidences to support its inclusion in our models. Model optimization and selection In order to fit a Michealis Menten equation for nitrate uptake and a Monod model of growth (data shown in Table 1, main text) we used a nonlinear least squares method implemented in the nls function from stats package in R software. The Michaelis Menten function was: V = V G:HIJ K NO M H IJK + NO M Where V maxno3 is the maximum uptake rate and H NO3 is a half saturation constant. And the Monod function: µμ G:HNO M µμ = K IJK + NO M
8 Where µ max is the maximum growth rate and K NO2 a half saturation constant. The initial parameter estimates for our tested models come from the experiments whose results are shown in Table 1, main text. From these estimates, we could set reasonable initial estimates for µ i, K i, η i, γ. Other parameters were set experimentally (D = 0.3 day -1, N 0 = 320 µm). The rest of the parameters were unknown, and their initial values were set tentatively for estimation. Parameter optimization was performed in two steps, first by applying a global optimization method, the simulated annealing, and then, with the estimates obtained from it, a local optimization method based in the Nelder Mead algorithm (Nelder and Mead 1965). The simulated annealing was implemented with the GenSA function from GenSA R package. The Nelder Mead algorithm was implemented with the optim function from the stats R package. In both global and local optimization, the objective functions were set according to a version of the Levenberg Mardquart minimization criterion (Levenberg 1944): E = YZ[ YZ) P ) 4QRS P ) TUV P ) TUV 1 CV ;E + P, 4QRS P, TUV P, TUV 1 CV ;< The parameter set providing the minimum value of E is considered the best estimate. In this formula, P ipred (being i = 1 for A. falcatus, 2 for Oscillatoria sp.) correspond to the abundances of each species predicted by the model for times t = 1 to t = n, and P iobs, the same abundances from the real data. CV Pi is the coefficient of variation of the observed population abundances. For the values of the observed population abundances, with used a smoothed series (Ellner et al. 2002). After optimization of models 3, 4, 5 and 6, we used an ad hoc version of the Akaike Information Criterion, in which the values of E obtained from the best model fits were penalized by the number of parameters in the model, according to this formula:
9 A = log10 SZ` SZ) E + 2 log10 k Where d = 1,... 9 are the experimental runs. E is the minimized value from the criterion detailed above, and k is the number of model parameters. In order to compare A between models, A models were calculated as follows: ΔA model i = A model i lowest A Another goodness of fit criterion, the relative mean error (RME) was calculated as follows: RME = YZ[ YZ) P )4QRS P )TUV P )TUV + P,4QRS P,TUV P,TUV In this formula, P ipred (being i = 1 for A. falcatus, 2 for Oscillatoria sp.) correspond to the abundances of each species predicted by the model for times t = 1 to t = n, and P iobs, the same abundances from the real data. The results from model optimization process and the values of A are shown in Table A1 considering all the experiments as a single data set and in Table A2 by using each independent experiment as a single data set. In Table A1 we can see that, for Models 3 to 6, parameters take similar values, always in the same orders of magnitude, as expected from similarly structured models, with just minor differences in some functional forms. Models 1 and 2 were not able to fit the experiments showing species coexistence (these models do not predict coexistence under any circumstances, see above). Our results showing coexistence are the main findings of this work. For this reason, they were not considered in the comparison of A and RME values with the other models. Also, these two models showed larger amplitude in the range of parameter values than models 3 to 6. Particularly for the allelopathy parameter (γ) which spans one order of magnitude. Such a large range in the value of one of the most important model parameters indicates poor performance of Models 1 and 2. This suggests that the model structure needs to include some
10 additional feature that makes it possible to modulate the effect of allelopathy without the need to vary a lot the value of this parameter. This is what models 3 to 6 do with their nonlinear functional forms. Overall, the selection criterion benefits Model 3, the simplest model (considering 3 to 6 only) together with Model 4. In our context, model selection is not a very relevant issue, since the fitting to our data could always be improved by adding more complex functional forms or additional, less relevant mechanism. Our aim here was to keep the minimum relevant mechanisms, using the simplest possible functional forms in order to generate a model that was able to fit qualitatively and quantitatively the experimental results in an acceptable way. From this point of view, we could say that models 3 or 4 are better than 5 and 6, since the latter two incorporate an additional function for growth that makes the models more complex. But, this increase in complexity does not improve the fitting substantially. However, with respect to models 3 and 4, there are no strong reasons to decide whether Model 3 is better than Model 4, despite the slightly lower value of A for Model 3. Table A2 allows us to perform the same comparison as previously for each independent experiment. This makes it easier to see the effect of each specific feature of a model for each of the outcomes of competition. Table A2 shows that the best performance (in general) of Models 3 and 4 relative to Models 5 and 6 was due to a better fit to those experiments showing exclusion of either species. If we look at the experiments showing coexistence, the performance of all the models was about the same, being 5 or 6 even slightly better than 3 or 4. The reason why, globally, Model 3 performs better (see Table A1) was because it was the most consistent model across all of the experiments.
11 Model simulation In figure 1, we simulated the models with the optimized values of the parameters and the actual initial conditions of the state variables for each individual experiment. The dilution rate of the system and the inflowing nitrate concentration were the only parameters that were not optimized, and their values were 0.3 day -1 and 320 µm respectively. References Roy, S The coevolution of two phytoplankton species on a single resource. Allelopathy as a pseudo mixotrophy. - Theor Popul Biol. 75: Roy, S Importance of allelopathy as pseudo-mixotrophy for the dynamics and diversity of phytoplankton. Biodiversity in Ecosystems - Linking Structure and Function. Eds Lo et al, InTech, UK, pp Roy, S. et al Competing effects of toxin-producing phytoplankton on overall plankton populations in the Bay of Bengal. - B Math Biol 68: Hsu, S. B. and Waltman, P A survey of mathematical models of competition with an inhibitor. - Math Biosci. 187: Mukhopadyhay, A. et al A delay differential equations model of plankton allelopathy. - Math Biosci. 149: Mougi, A Allelopathic adaptation can cause competitive coexistence. - Theor Ecol. 6: Nelder, J. A. and Mead, R A simplex method for function minimization. - Comput J. 7: Levenberg, K A Method for the Solution of Certain Non-Linear Problems in Least
12 Squares. Quart Appl Math. 2: Ellner, S.P. et al Fitting population dynamics models to time series data by gradient matching. Ecology 83:
13 Species abundance LN(cell ml 1 ) A D B E C F Time (days)
14 Scaled population abundances A(a) B(a) C(a) Time (days) Scaled A.falcatus abundance A(b) B(b) C(b) Scaled Oscillatoria sp. abundance
15
16
17 Realized peaks of P 1 within 100 days Realized peaks of P 2 within 100 days Ini$al abundance of Allelopathic species Ini$al abundance of Allelopathic species Ini$al abundance of Allelopathic species
18 Figure A1. Plots of Oscillatoria and Ankistrodesmus abundances in long term experiments shown in log-scale, which is more appropriate to see in detail when the species were excluded. A. B, C correspond to those experiments where Oscillatoria was excluded (Figure 1 A, B, C in the main text) and D, E, F correspond to those experiments where Ankistrodesmus was excluded (Figure 1 G, H, I in the main text). Figure A2. Details on the oscillatory patterns of the three experiments with intermediate Oscillatoria sp. to A. falcatus ratios: 0.11 (A(a), A(b)), 0.20 (B(a), B(b)) and 0.36 (C(a)), C(b)). Data are represented by the smoothed curves used for model fitting (see Model optimization and selection sections in this Appendix, Figure 1 in the main text). In A(a), B(a) and C(a) green color represents A. falcatus and blue Oscillatoria sp. Major period and amplitude oscillations seem to overlap with other smaller oscillations, as the model fits in Figure 1 in the main text show. Population densities in C(a) are nearly stable. Data begin in days and end in day 90. Numbers labeling the lines in plots A(b), B(b) and C(b) correspond to days. Figure A3. Results of the same model simulations from Figure 4 in the main text showing the % of occurrence of exclusion of Ankistrodesmus falcatus, P 1 (left column), coexistence (central column) and exclusion of Oscillatoria sp., P 2. The colour scale indicates increased % from blue to red. The upper row panel represents the simulations of 'Model 3' describing density-dependent allelopathic killing, and the lower row panel 'Model 4' (see Table 2). The simulations were performed with the same conditions and methods (see Methods) for Figure 4 in the main text. Spline interpolation was applied to smooth the contours of the risk of extinction and coexistence regions. Figure A4. Results of oscillation amplitudes from species population dynamics obtained from model simulations performed under the same conditions as in Figure A3 and Figure 4 in the main
19 text. Figure A5. Same as figure A3 for realized peaks shown in the dynamics of each species.
20 Table A1. Results of model optimization and criterions for model selection considering all the experimental runs together. The values of each parameter are the average of the optimized values from each individual experiment, and within brackets, the range of those values. Last three rows show the values of the different criterions for model selection (RME: relative mean error, A, and the differential with the lowest A computed). Models 1 and 2 were not considered in the model selection process due to their inability to predict either coexistence or oscillatory coexistence (see text in Appendix 1, Model optimization and selection). Parameters Model 1 Model 2 Model 3 Model 4 Model 5 Model 6 η [ ] 9500 [ ] 8314 [ ] 7663 [ ] 7915 [ ] 9266 [ ] η [ ] 3000 [ ] 2380 [ ] 2365 [ ] 1618 [ ] 2791 [ ] µ [ ] 1.4 [ ] 0.9 [ ] 1.1[ ] 1.1 [ ] 1 [ ] µ [ ] [ ] 0.44 [ ] - - µ [ ] [ ] 0.45 [ ] K [ ] 31.2 [ ] 15.4 [ ] 15.4 [ ] 15 [9-16] 15 [15-16] K [ ] [ ] 36.2 [ ] 45.1 [20-69] 28 [22-40] 41 [40-42] ε [ ] [ ] 0.04 [ ] K p [ ] [ ] 1995 [ ] γ [ ] [ ] 549 [ ] 448 [ ] 583 [ ] 487 [ ] n [ ] [ ] [ ] [ ] a [ ] [ ] - b [ ] [ ] - C [ ] [ ] C [ ] [ ] RME A A
21 Table A2. Same model selection criterions as in Table A1, but reported for each single experiment separately. The experiments are identified with the same indexes as in Figure 1 from the main manuscript. Experiment Model RME A A Model Model A Model Model Model Model Model Model B Model Model Model Model Model Model C Model Model Model Model Model Model D Model Model Model Model Model Model E Model Model Model Model Model Model F Model Model Model Model Model Model G Model Model Model Model Model Model H Model Model Model Model Model Model I Model Model Model Model
COMPETITION OF FAST AND SLOW MOVERS FOR RENEWABLE AND DIFFUSIVE RESOURCE
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 2, Number, Spring 22 COMPETITION OF FAST AND SLOW MOVERS FOR RENEWABLE AND DIFFUSIVE RESOURCE SILOGINI THANARAJAH AND HAO WANG ABSTRACT. In many studies of
More informationFeedback-mediated oscillatory coexistence in the chemostat
Feedback-mediated oscillatory coexistence in the chemostat Patrick De Leenheer and Sergei S. Pilyugin Department of Mathematics, University of Florida deleenhe,pilyugin@math.ufl.edu 1 Introduction We study
More information3.5 Competition Models: Principle of Competitive Exclusion
94 3. Models for Interacting Populations different dimensional parameter changes. For example, doubling the carrying capacity K is exactly equivalent to halving the predator response parameter D. The dimensionless
More informationModeling Microbial Populations in the Chemostat
Modeling Microbial Populations in the Chemostat Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Modeling Microbial Populations in the Chemostat MBI, June 3, 204 / 34 Outline Why
More informationInterspecific Competition
Interspecific Competition Intraspecific competition Classic logistic model Interspecific extension of densitydependence Individuals of other species may also have an effect on per capita birth & death
More informationDYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD
Journal of Sustainability Science Management Volume 10 Number 2, December 2015: 16-23 ISSN: 1823-8556 Penerbit UMT DYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD SARKER MD SOHEL
More informationImportance of Allelopathy as Peudo-Mixotrophy for the Dynamics and Diversity of Phytoplankton
Chapter Importance of Allelopathy as Peudo-Mixotrophy for the Dynamics and Diversity of Phytoplankton Shovonlal Roy Additional information is available at the end of the chapter http://dx.doi.org/0.577/59055.
More informationOrdinary Differential Equations
Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations
More informationBIOS 6150: Ecology Dr. Stephen Malcolm, Department of Biological Sciences
BIOS 6150: Ecology Dr. Stephen Malcolm, Department of Biological Sciences Week 7: Dynamics of Predation. Lecture summary: Categories of predation. Linked prey-predator cycles. Lotka-Volterra model. Density-dependence.
More informationGause, Luckinbill, Veilleux, and What to Do. Christopher X Jon Jensen Stony Brook University
Gause, Luckinbill, Veilleux, and What to Do Christopher X Jon Jensen Stony Brook University Alternative Models of Predation: Functional Responses: f (N) Prey Dependent f (N/P) Ratio Dependent Possible
More informationCompetitive Exclusion in a Discrete-time, Size-structured Chemostat Model
Competitive Exclusion in a Discrete-time, Size-structured Chemostat Model Hal L. Smith Department of Mathematics Arizona State University Tempe, AZ 85287 1804, USA E-mail: halsmith@asu.edu Xiao-Qiang Zhao
More informationBIOS 5970: Plant-Herbivore Interactions Dr. Stephen Malcolm, Department of Biological Sciences
BIOS 5970: Plant-Herbivore Interactions Dr. Stephen Malcolm, Department of Biological Sciences D. POPULATION & COMMUNITY DYNAMICS Week 10. Population models 1: Lecture summary: Distribution and abundance
More informationBIO 114 Spring Protozoan Population Ecology and Interactions
Protozoan Population Ecology and Interactions Laboratory Learning Outcomes Conceptual 1. Describe the effect of birth rate and death rate on population growth. 2. Apply the concepts in BMEs 24.1, 24.2
More informationOikos. Appendix A1. o20830
1 Oikos o20830 Valdovinos, F. S., Moisset de Espanés, P., Flores, J. D. and Ramos-Jiliberto, R. 2013. Adaptive foraging allows the maintenance of biodiversity of pollination networks. Oikos 122: 907 917.
More informationGrowth models for cells in the chemostat
Growth models for cells in the chemostat V. Lemesle, J-L. Gouzé COMORE Project, INRIA Sophia Antipolis BP93, 06902 Sophia Antipolis, FRANCE Valerie.Lemesle, Jean-Luc.Gouze@sophia.inria.fr Abstract A chemostat
More informationMathematical Modeling of Competition for Light and Nutrients Between Phytoplankton Species in a Poorly Mixed Water Column
University of Wisconsin Milwaukee UWM Digital Commons Theses and Dissertations May 2014 Mathematical Modeling of Competition for Light and Nutrients Between Phytoplankton Species in a Poorly Mixed Water
More informationField experiments on competition. Field experiments on competition. Field experiments on competition
INTERACTIONS BETWEEN SPECIES Type of interaction species 1 species 2 competition consumer-resource (pred, herb, para) mutualism detritivore-detritus (food is dead) Field experiments on competition Example
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationMicrobial Grazers Lab
Microbial Grazers Lab Objective: Measure the rate at which bacteria are consumed by predators. Overview Size based food webs Microbial loop concepts Bacterial predators Methods to assess microbial grazing
More informationdv dt Predator-Prey Models
Predator-Prey Models This is a diverse area that includes general models of consumption: Granivores eating seeds Parasitoids Parasite-host interactions Lotka-Voterra model prey and predator: V = victim
More informationCompetition. Not until we reach the extreme confines of life, in the arctic regions or on the borders of an utter desert, will competition cease
Competition Not until we reach the extreme confines of life, in the arctic regions or on the borders of an utter desert, will competition cease Darwin 1859 Origin of Species Competition A mutually negative
More informationThe influence of resource limitation on the allelopathic effect of Chlamydomonas reinhardtii on other unicellular freshwater planktonic organisms
JPR Advance Access published August 21, 2013 Journal of Plankton Research plankt.oxfordjournals.org J. Plankton Res. (2013) 0(0): 1 6. doi:10.1093/plankt/fbt080 SHORT COMMUNICATION The influence of resource
More informationBehaviour of simple population models under ecological processes
J. Biosci., Vol. 19, Number 2, June 1994, pp 247 254. Printed in India. Behaviour of simple population models under ecological processes SOMDATTA SINHA* and S PARTHASARATHY Centre for Cellular and Molecular
More informationEcology 302: Lecture VII. Species Interactions.
Ecology 302: Lecture VII. Species Interactions. (Gotelli, Chapters 6; Ricklefs, Chapter 14-15) MacArthur s warblers. Variation in feeding behavior allows morphologically similar species of the genus Dendroica
More informationCommunity Structure. Community An assemblage of all the populations interacting in an area
Community Structure Community An assemblage of all the populations interacting in an area Community Ecology The ecological community is the set of plant and animal species that occupy an area Questions
More informationShifting N:P ratios do not affect the competitive traits of dinoflagellates under present and future climate scenarios
Shifting N:P ratios do not affect the competitive traits of dinoflagellates under present and future climate scenarios Maarten De Rijcke, Gabriel Orellana, Julie Vanden Bussche, Lynn Vanhaecke, Karel A.C.
More informationDETERMINATION OF DRUG RELEASE DURING DISSOLUTION OF NICORANDIL IN TABLET DOSAGE FORM BY USING REVERSE PHASE HIGH PERFORMANCE LIQUID CHROMATOGRAPHY
CHAPTER 9 DETERMINATION OF DRUG RELEASE DURING DISSOLUTION OF NICORANDIL IN TABLET DOSAGE FORM BY USING REVERSE PHASE HIGH PERFORMANCE LIQUID CHROMATOGRAPHY CHAPTER 9 Determination of drug release during
More informationLinear Models for Regression CS534
Linear Models for Regression CS534 Prediction Problems Predict housing price based on House size, lot size, Location, # of rooms Predict stock price based on Price history of the past month Predict the
More informationCompetition-induced starvation drives large-scale population cycles in Antarctic krill
In the format provided by the authors and unedited. SUPPLEMENTARY INFORMATION VOLUME: 1 ARTICLE NUMBER: 0177 Competition-induced starvation drives large-scale population cycles in Antarctic krill Alexey
More informationSpring /30/2013
MA 138 - Calculus 2 for the Life Sciences FINAL EXAM Spring 2013 4/30/2013 Name: Sect. #: Answer all of the following questions. Use the backs of the question papers for scratch paper. No books or notes
More informationModeling and Simulation in Medicine and the Life Sciences
Frank C. Hoppensteadt Charles S. Peskin Modeling and Simulation in Medicine and the Life Sciences Second Edition With 93 Illustrations fyj Springer Contents Series Preface Preface vii ix Introduction 1
More informationMicrobial Grazers Lab
Microbial Grazers Lab Objective: Measure the rate at which bacteria are consumed by predators. Overview Size based food webs Microbial loop concepts acterial predators Methods to assess microbial grazing
More informationCommonly Used Designs
Commonly Used Designs Partial Additive Replacement Series Additive Complete Additive The Complete Additive Design replacement series design partial additive design Provides information across a range of
More informationSIMPLE MODEL Direct Binding Analysis
Neurochemistry, 56:120:575 Dr. Patrick J. McIlroy Supplementary Notes SIMPLE MODEL Direct Binding Analysis The interaction of a (radio)ligand, L, with its receptor, R, to form a non-covalent complex, RL,
More informationD. Correct! Allelopathy is a form of interference competition in plants. Therefore this answer is correct.
Ecology Problem Drill 18: Competition in Ecology Question No. 1 of 10 Question 1. The concept of allelopathy focuses on which of the following: (A) Carrying capacity (B) Limiting resource (C) Law of the
More informationLecture 20/Lab 21: Systems of Nonlinear ODEs
Lecture 20/Lab 21: Systems of Nonlinear ODEs MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Coupled ODEs: Species
More informationMetabolic trade-offs promote diversity in a model ecosystem
Metabolic trade-offs promote diversity in a model ecosystem Anna Posfai, Thibaud Taillefumier, Ben Weiner, Ned Wingreen Princeton University q-bio Rutgers University, July 25 2017 How can we explain species
More informationAllelopathic interactions in phytoplankton population ecology
International Allelopathy Society I.S.S.N: 2312-8666 Journal of Allelochemical Interactions July 2016 (1) July 2016. JAI 2 (2): 25-34 Allelopathic interactions in phytoplankton population ecology Aldo
More informationFurther analysis. Objective Based on the understanding of the LV model from last report, further study the behaviours of the LV equations parameters.
Further analysis Objective Based on the understanding of the LV model from last report, further study the behaviours of the LV equations parameters. Changes have been made Method of finding the mean is
More informationInterspecific Patterns. Interference vs. exploitative
Types of Competition Interference vs. exploitative Intraspecific vs. Interspeific Asymmetric vs. Symmetric Interspecific Patterns When two similar species coexist, there are three outcomes: Competitive
More informationMeasurement of Enzyme Activity - ALP Activity (ALP: Alkaline phosphatase)
Measurement of Enzyme Activity - ALP Activity (ALP: Alkaline phosphatase) Measurement and analysis of enzyme activity is often used in the field of life science such as medicines and foods to investigate
More informationBIOL 410 Population and Community Ecology. Predation
BIOL 410 Population and Community Ecology Predation Intraguild Predation Occurs when one species not only competes with its heterospecific guild member, but also occasionally preys upon it Species 1 Competitor
More informationThe effect of emigration and immigration on the dynamics of a discrete-generation population
J. Biosci., Vol. 20. Number 3, June 1995, pp 397 407. Printed in India. The effect of emigration and immigration on the dynamics of a discrete-generation population G D RUXTON Biomathematics and Statistics
More informationSystems Biology in Photosynthesis INTRODUCTION
1 / 26 Systems Biology in Photosynthesis INTRODUCTION Rainer Machné Institute for Theoretical Chemistry, University of Vienna, Austria PSI - Photon System Instruments, Czech Republic Brno, April, 2011
More informationSupporting Information. Methods. Equations for four regimes
Supporting Information A Methods All analytical expressions were obtained starting from quation 3, the tqssa approximation of the cycle, the derivation of which is discussed in Appendix C. The full mass
More informationBioinformatics 3. V18 Kinetic Motifs. Fri, Jan 8, 2016
Bioinformatics 3 V18 Kinetic Motifs Fri, Jan 8, 2016 Modelling of Signalling Pathways Curr. Op. Cell Biol. 15 (2003) 221 1) How do the magnitudes of signal output and signal duration depend on the kinetic
More informationBioinformatics 3! V20 Kinetic Motifs" Mon, Jan 13, 2014"
Bioinformatics 3! V20 Kinetic Motifs" Mon, Jan 13, 2014" Modelling of Signalling Pathways" Curr. Op. Cell Biol. 15 (2003) 221" 1) How do the magnitudes of signal output and signal duration depend on the
More informationThe implications of neutral evolution for neutral ecology. Daniel Lawson Bioinformatics and Statistics Scotland Macaulay Institute, Aberdeen
The implications of neutral evolution for neutral ecology Daniel Lawson Bioinformatics and Statistics Scotland Macaulay Institute, Aberdeen How is How is diversity Diversity maintained? maintained? Talk
More informationPrinciples of Synthetic Biology: Midterm Exam
Principles of Synthetic Biology: Midterm Exam October 28, 2010 1 Conceptual Simple Circuits 1.1 Consider the plots in figure 1. Identify all critical points with an x. Put a circle around the x for each
More informationEffect of Species 2 on Species 1 Competition - - Predator-Prey + - Parasite-Host + -
Community Ecology Community - a group of organisms, of different species, living in the same area Community ecology is the study of the interactions between species The presence of one species may affect
More informationCommunity Ecology. Classification of types of interspecific interactions: Effect of Species 1 on Species 2
Community Ecology Community - a group of organisms, of different species, living in the same area Community ecology is the study of the interactions between species The presence of one species may affect
More informationAdaptive dynamics in an individual-based, multi-resource chemostat model
Adaptive dynamics in an individual-based, multi-resource chemostat model Nicolas Champagnat (INRIA Nancy) Pierre-Emmanuel Jabin (Univ. Maryland) Sylvie Méléard (Ecole Polytechnique) 6th European Congress
More informationOptimization of an Adapta Kinase Assay for CAMK1
Overview Optimization of an Adapta Kinase Assay for CAMK1 This protocol describes how to perform an Adapta assay with the kinase CAMK1. To maximize the ability of the assay to detect ATP-competitive inhibitors,
More informationDynamical Systems and Chaos Part II: Biology Applications. Lecture 6: Population dynamics. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part II: Biology Applications Lecture 6: Population dynamics Ilya Potapov Mathematics Department, TUT Room TD325 Living things are dynamical systems Dynamical systems theory
More informationREVISION: POPULATION ECOLOGY 18 SEPTEMBER 2013
REVISION: POPULATION ECOLOGY 18 SEPTEMBER 2013 Lesson Description In this lesson we: Revise population ecology by working through some exam questions. Key Concepts Definition of Population A population
More informationCompetition for resources: complicated dynamics in the simple Tilman model
DOI 86/s464-5-246-6 ESEACH Open Access Competition for resources: complicated dynamics in the simple Tilman model Joost H J van Opheusden *, Lia Hemerik, Mieke van Opheusden 2 and Wopke van der Werf 2
More informationOptimal foraging and predator prey dynamics III
Theoretical Population Biology 63 (003) 69 79 http://www.elsevier.com/locate/ytpbi Optimal foraging and predator prey dynamics III Vlastimil Krˇ ivan and Jan Eisner Department of Theoretical Biology, Institute
More informationBIO S380T Page 1 Summer 2005: Exam 2
BIO S380T Page 1 Part I: Definitions. [5 points for each term] For each term, provide a brief definition that also indicates why the term is important in ecology or evolutionary biology. Where I ve provided
More informationModeling Marine Microbes: Past, Present and Prospects Mick Follows MIT
Modeling Marine Microbes: Past, Present and Prospects Mick Follows MIT What do we mean by models? Concepts Statistical relationships Mathematical descriptions Numerical models Why theory and numerical
More informationParameter Sensitivity In A Lattice Ecosystem With Intraguild Predation
Parameter Sensitivity In A Lattice Ecosystem With Intraguild Predation N. Nakagiri a, K. Tainaka a, T. Togashi b, T. Miyazaki b and J. Yoshimura a a Department of Systems Engineering, Shizuoka University,
More informationSimplicity is Complexity in Masquerade. Michael A. Savageau The University of California, Davis July 2004
Simplicity is Complexity in Masquerade Michael A. Savageau The University of California, Davis July 2004 Complexity is Not Simplicity in Masquerade -- E. Yates Simplicity is Complexity in Masquerade One
More informationIdentification of Non-linear Dynamical Systems
Identification of Non-linear Dynamical Systems Linköping University Sweden Prologue Prologue The PI, the Customer and the Data Set C: I have this data set. I have collected it from a cell metabolism experiment.
More informationPredicting the outcome of competition when fitness inequality is variable
rsos.royalsocietypublishing.org Predicting the outcome of competition when fitness inequality is variable Research Cite this article: Pedruski MT, Fussmann GF, Gonzalez A. 25 Predicting the outcome of
More informationInterspecific Competition
Use Interspecific Competition 0.8 0.6 0.4 0. 0 0 0.5.5.5 3 Resource The niche and interspecific competition Species A Use Use 0.8 0.6 0.4 0. Species B 0 0 0.5.5.5 3 0.8 0.6 0.4 0. 0 0 0.5.5.5 3 Resource
More information8 Ecosystem stability
8 Ecosystem stability References: May [47], Strogatz [48]. In these lectures we consider models of populations, with an emphasis on the conditions for stability and instability. 8.1 Dynamics of a single
More informationJOURNAL OF MATH SCIENCES -JMS- Url: Jl. Pemuda, No. 339 Kolaka Southeast Sulawesi, Indonesia
JOURNAL OF MATH SCIENCES -JMS- Url: http://usnsj.com/index.php/jms Jl. Pemuda, No. 339 Kolaka Southeast Sulawesi, Indonesia STABILITY ANALYSIS OF ONE PREY TWO PREDATOR MODEL WITH HOLLING TYPE III FUNCTIONAL
More informationGary G. Mittelbach Michigan State University
Community Ecology Gary G. Mittelbach Michigan State University Sinauer Associates, Inc. Publishers Sunderland, Massachusetts U.S.A. Brief Table of Contents 1 Community Ecology s Roots 1 PART I The Big
More informationEmergence of diversity in a biological evolution model
Journal of Physics: Conference Series PAPER OPE ACCESS Emergence of diversity in a biological evolution model To cite this article: R Wang and C Pujos 2015 J. Phys.: Conf. Ser. 604 012019 Related content
More informationx 2 F 1 = 0 K 2 v 2 E 1 E 2 F 2 = 0 v 1 K 1 x 1
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 20, Number 4, Fall 1990 ON THE STABILITY OF ONE-PREDATOR TWO-PREY SYSTEMS M. FARKAS 1. Introduction. The MacArthur-Rosenzweig graphical criterion" of stability
More informationROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 6, Number 1, Winter 1997 ROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM J. CHATTOPADHYAY, E. BERETTA AND F. SOLIMANO ABSTRACT. The present paper deals with
More informationRank-abundance. Geometric series: found in very communities such as the
Rank-abundance Geometric series: found in very communities such as the Log series: group of species that occur _ time are the most frequent. Useful for calculating a diversity metric (Fisher s alpha) Most
More informationBistability and a differential equation model of the lac operon
Bistability and a differential equation model of the lac operon Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Fall 2016 M. Macauley
More informationFeedback control for a chemostat with two organisms
Feedback control for a chemostat with two organisms Patrick De Leenheer and Hal Smith Arizona State University Department of Mathematics and Statistics Tempe, AZ 85287 email: leenheer@math.la.asu.edu,
More informationIt has become increasingly evident that nonlinear phenomena
Increased competition may promote species coexistence J. Vandermeer*, M. A. Evans*, P. Foster*, T. Höök, M. Reiskind*, and M. Wund* *Department of Ecology and Evolutionary Biology, and School of Natural
More informationName Student ID. Good luck and impress us with your toolkit of ecological knowledge and concepts!
Page 1 BIOLOGY 150 Final Exam Winter Quarter 2000 Before starting be sure to put your name and student number on the top of each page. MINUS 3 POINTS IF YOU DO NOT WRITE YOUR NAME ON EACH PAGE! You have
More informationGause s exclusion principle revisited: artificial modified species and competition
arxiv:nlin/0007038v1 [nlin.ao] 26 Jul 2000 Gause s exclusion principle revisited: artificial modified species and competition J. C. Flores and R. Beltran Universidad de Tarapacá, Departamento de Física,
More informationEcology 203, Exam III. November 16, Print name:
Ecology 203, Exam III. November 16, 2005. Print name: Read carefully. Work accurately and efficiently. The exam is worth 100 points (plus 6 extra credit points). Choose four of ten concept-exploring questions
More informationNumerical Methods of Approximation
Contents 31 Numerical Methods of Approximation 31.1 Polynomial Approximations 2 31.2 Numerical Integration 28 31.3 Numerical Differentiation 58 31.4 Nonlinear Equations 67 Learning outcomes In this Workbook
More informationNatal versus breeding dispersal: Evolution in a model system
Evolutionary Ecology Research, 1999, 1: 911 921 Natal versus breeding dispersal: Evolution in a model system Karin Johst 1 * and Roland Brandl 2 1 Centre for Environmental Research Leipzig-Halle Ltd, Department
More informationESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp DOI: /m2an:
Mathematical Modelling and Numerical Analysis ESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp. 339 344 DOI: 10.1051/m2an:2003029 PERSISTENCE AND BIFURCATION ANALYSIS
More informationStochastic models in biology and their deterministic analogues
Stochastic models in biology and their deterministic analogues Alan McKane Theory Group, School of Physics and Astronomy, University of Manchester Newton Institute, May 2, 2006 Stochastic models in biology
More informationWhat is competition? Competition among individuals. Competition: Neutral Theory vs. the Niche
Competition: Neutral Theory vs. the Niche Reading assignment: Ch. 10, GSF (especially p. 237-249) Optional: Clark 2009 9/21/09 1 What is competition? A reduction in fitness due to shared use of a limited
More informationBifurcation and Stability Analysis of a Prey-predator System with a Reserved Area
ISSN 746-733, England, UK World Journal of Modelling and Simulation Vol. 8 ( No. 4, pp. 85-9 Bifurcation and Stability Analysis of a Prey-predator System with a Reserved Area Debasis Mukherjee Department
More informationOptimal Designs for 2 k Experiments with Binary Response
1 / 57 Optimal Designs for 2 k Experiments with Binary Response Dibyen Majumdar Mathematics, Statistics, and Computer Science College of Liberal Arts and Sciences University of Illinois at Chicago Joint
More informationFrom simple rules to cycling in community assembly
OIKOS 105: 349/358, 2004 From simple rules to cycling in community assembly Sebastian J. Schreiber and Seth Rittenhouse Schreiber, S. J. and Rittenhouse, S. 2004. From simple rules to cycling in community
More information1 Overview. 2 Learning from Experts. 2.1 Defining a meaningful benchmark. AM 221: Advanced Optimization Spring 2016
AM 1: Advanced Optimization Spring 016 Prof. Yaron Singer Lecture 11 March 3rd 1 Overview In this lecture we will introduce the notion of online convex optimization. This is an extremely useful framework
More informationSUPPLEMENTARY INFORMATION
Social diversity promotes the emergence of cooperation in public goods games Francisco C. Santos 1, Marta D. Santos & Jorge M. Pacheco 1 IRIDIA, Computer and Decision Engineering Department, Université
More informationThe major interactions between two organisms in a mixed culture are : Competition Neutralism Mutualism Commensalism Amensalism Prey-predator
1 Introduction Major classes of interaction in mixed cultures Simple models describing mixed-cultures interactions Mixed culture in nature Industrial utilization of mixed cultures 2 1 The dynamic of mixed
More informationYield-Density Equations
Yield-Density Equations A General Model of Intraspecific Density Effects Yield-Density Equations Y Y wn N max ( 1+ an ) b Total yield of the population per unit area Yield-Density Equations Y w Y wn N
More informationInfluence of Noise on Stability of the Ecosystem
Commun. Theor. Phys. 60 (2013) 510 514 Vol. 60, No. 4, October 15, 2013 Influence of Noise on Stability of the Ecosystem LIU Tong-Bo ( ) and TANG Jun (» ) College of Science, China University of Mining
More informationAfter lectures by. disappearance of reactants or appearance of. measure a reaction rate we monitor the. Reaction Rates (reaction velocities): To
Revised 3/21/2017 After lectures by Dr. Loren Williams (GeorgiaTech) Protein Folding: 1 st order reaction DNA annealing: 2 nd order reaction Reaction Rates (reaction velocities): To measure a reaction
More informationMinimization of Static! Cost Functions!
Minimization of Static Cost Functions Robert Stengel Optimal Control and Estimation, MAE 546, Princeton University, 2017 J = Static cost function with constant control parameter vector, u Conditions for
More informationAggregations on larger scales. Metapopulation. Definition: A group of interconnected subpopulations Sources and Sinks
Aggregations on larger scales. Metapopulation Definition: A group of interconnected subpopulations Sources and Sinks Metapopulation - interconnected group of subpopulations sink source McKillup and McKillup
More informationExpanding brackets and factorising
CHAPTER 8 Epanding brackets and factorising 8 CHAPTER Epanding brackets and factorising 8.1 Epanding brackets There are three rows. Each row has n students. The number of students is 3 n 3n. Two students
More informationSymbol Definition Value Range (Experimental) Time of the end of residual growth 10 h 10
Table 1. Definitions and experimental values of parameters Symbol Definition Value ange (Experimental) T I Time of the end of residual growth 10 h 10 T D D D ate s s G max V max c-fed c-fed Delay in Death
More informationAnalysis of Plankton Populations
Analysis of Plankton Populations Sean Li September 21, 25 Abstract I use numerical and analytical techniques to study the dynamics of various models of plankton food webs with and without resource fluctuations.
More informationPRACTICE TEST ANSWER KEY & SCORING GUIDELINES GRADE 8 MATHEMATICS
Ohio s State Tests PRACTICE TEST ANSWER KEY & SCORING GUIDELINES GRADE 8 MATHEMATICS Table of Contents Questions 1 25: Content Summary and Answer Key... iii Question 1: Question and Scoring Guidelines...
More informationName Class Date. Section: How Organisms Interact in Communities. In the space provided, explain how the terms in each pair differ in meaning.
Section: How Organisms Interact in Communities In the space provided, explain how the terms in each pair differ in meaning 1 coevolution, secondary compounds 2 predation, parasitism Complete each statement
More informationAge (x) nx lx. Population dynamics Population size through time should be predictable N t+1 = N t + B + I - D - E
Population dynamics Population size through time should be predictable N t+1 = N t + B + I - D - E Time 1 N = 100 20 births 25 deaths 10 immigrants 15 emmigrants Time 2 100 + 20 +10 25 15 = 90 Life History
More informationVEGETATION PROCESSES IN THE PELAGIC: A MODEL FOR ECOSYSTEM THEORY
Colin S. Reynolds VEGETATION PROCESSES IN THE PELAGIC: A MODEL FOR ECOSYSTEM THEORY Introduction (Otto Kinne) Colin S. Reynolds: A Laudatio (William D. Williams) Publisher: Ecology Institute Nordbunte
More information