Transcript Document
Consequences of Heterogeneous Survival Rates of an Entomopathogenic Nematode. Chris Dugaw Department of Mathematics Humboldt State University Outline Biological background Understanding Nematode Survival Experimental Setup Survival Analysis Results Discussion Entomopathogenic nematodes Insect predators, in soil or litter Can move >2 cm/day following volatiles Kills prey with symbiotic bacteria injected into host One nematode in 800K+ emerge Images courtesy of Ed Lewis, Virginia Tech Nematode life cycle http://www.bath.ac.uk/bio-sci/clarke.htm Current Uses of Nematodes as Biocontrol Agents Commodity Insect Pests Artichokes Artichoke plume moth Berries Root weevils Citrus Root weevils Cranberries Root Weevils Cranberry girdler Mushrooms Fungus gnats Ornamentals Root Weevils Wood borers Fungus gnats Turf Scarabs Mole crickets Billbugs Armyworm, Cutworm, Webworm Source: http://www.oardc.ohio-state.edu/nematodes/biologyecology.htm Study Site: the Bodega Marine Reserve The predatory nematode Heterorhabditis marelatus neudorff.de/nuetzlinge/img/hmne.jpg A natural host: the ghost moth Hepialus californicus Host larvae Adult Host Larvae infected by nematodes Ghost moth caterpillars feed on the roots of bush lupine (Lupinus arboreus) Lupine killed by ghost moth caterpillars Large-scale ghost moth outbreaks occur, killing 10,000+ mature lupines 2001 2002 The trophic cascade: predators indirectly affect producers by suppressing herbivores Entomopathogenic nematode Heterorhabditis marelatus Ghost moth Hepialus californicus Bush lupine Lupinus arboreus + Strong 1997, Strong et al. 1999 Seasonal Dynamics Wet Winter •Nematodes search for hosts •Nematode reproduction occurs •Hosts are in larval stage Dry Summer • Nematodes are inactive • Nematodes must survive • Host become adults and disperse • Host eggs are deposited on bush Seasonal Dynamics Wet Winter •Nematodes search for hosts •Nematode reproduction occurs •Hosts are in larval stage Dry Summer • Nematodes are inactive • Nematodes must survive • Host become adults and disperse • Host eggs are deposited on bush Outline Biological background Understanding Nematode Survival Experimental Setup Survival Analysis Results Discussion Experimental design 2 treatments = Lupine, Grasslands 4 replicates/treatment = 8 total replicates 50 tubes/replicate = 400 total tubes Each sampling date, removed 10 tubes/replicate = 80 total tubes/sampling date Assessed nematodes using ‘bait’ insects Sampled 3 times over a Summer Each tube - 30 g past. soil - 1100 IJ nematodes - Fine mesh covers Survival Analysis Survival Analysis Homogeneous Death Rates Survival Analysis Homogeneous Death Rates Exponential Distribution Survival Analysis Homogeneous Death Rates Heterogeneous Death Rates Exponential Distribution Survival Analysis Homogeneous Death Rates Heterogeneous Death Rates Exponential Distribution Mixed Exponential Distribution First Step: Exponential Fit Mixed Exponential Distributions Individuals have different mortality rates, k. Mixed Exponential Distributions Individuals have different mortality rates, k. Risk of death for each individual is constant over time. Mixed Exponential Distributions Individuals have different mortality rates, k. Risk of death for each individual is constant over time. The conditional distribution for individual lifespan, T, given k is exponential. Pareto Distribution: The mixed exponential you get when you assume k is gamma distributed. Pareto Distribution: The mixed exponential you get when you assume k is gamma distributed. A simple function form can be derived by integrating: fT (t ) fT |K (t | k ) f K (k )dk 0 Pareto Distribution: The mixed exponential you get when you assume k is gamma distributed. A simple function form can be derived by integrating: fT (t ) fT |K (t | k ) f K (k )dk 0 kt 1 k ke k e dk ( ) 0 Pareto Distribution: The mixed exponential you get when you assume k is gamma distributed. A simple function form can be derived by integrating: f T (t ) fT |K (t | k ) f K ( k )dk 0 kt 1 k ke k e dk ( ) 0 , (t ) 1 The distribution of survival rates shifts over time leading to a decrease in mean mortality rate. McNolty, Doyle and Hansen, Technometrics, 1980 Improvement: Pareto Fit Improvement: Pareto Fit = 0.29 = 2.77 = 0.73 = 2.77 Why is it an improvement? 1. Akaike says so: AICc = 3.46 Why is it an improvement? 1. Akaike says so: AICc = 3.46 2. Provides a greater understanding of the biological system. Why is it an improvement? 1. Akaike says so: AICc = 3.46 2. Provides a greater understanding of the biological system. 3. Allows us to quantify heterogeneity using the scale parameter, . Results: Survival in soil is heterogeneous. Results: Survival in soil is heterogeneous. Mean mortality is higher in the grasslands. (log ratio test: 2 = 0.449, df=1, p = 0.050) Results: Survival in soil is heterogeneous. Mean mortality is higher in the grasslands. (log ratio test: 2 = 0.449, df=1, p = 0.050) Heterogeneity same in the two treatments. (log ratio test: 2 = 0.279, df=1, p = 0.98) Outline Biological background Understanding Nematode Survival Experimental Setup Survival Analysis Results Discussion Feedback loop in trophic cascade Entomopathogenic nematode Heterorhabditis marelatus + non-trophic trophic Ghost moth dry season Hepialus californicus wet season + + Bush lupine Lupinus arboreus Preisser, Dugaw, et al., In Review Alternative Explanations for Observations Decreasing individual hazards Alternative Explanations for Observations Decreasing individual hazards Density Dependant Survival Future work: Apply this analysis to new experiments to assess survival and heterogeneity in different soil types. Future work: Apply this analysis to new experiments to assess survival and heterogeneity in different soil types. Compare fitted shape parameter to physical soil properties. Future work: Apply this analysis to new experiments to assess survival and heterogeneity in different soil types. Compare fitted shape parameter to physical soil properties. Incorporate heterogeneous survival into a stochastic model that includes nematode reproduction. Thanks to: Evan Preisser Mike Eng Don Strong Brian Dennis Support of: NSF UC Davis Dissertation Year Fellowship UCD Faculty Fellow